Localization of sodium channels in intercalated disks modulates cardiac conduction

Acute Adenoviral Infection Elicits an Arrhythmogenic Substrate Prior to Myocarditis

BACKGROUND

Viral cardiac infection represents a significant clinical challenge encompassing several etiological agents, disease stages, complex presentation, and a resulting lack of mechanistic understanding. Myocarditis is a major cause of sudden cardiac death in young adults, where current knowledge in the field is dominated by later disease phases and pathological immune responses. However, little is known regarding how infection can acutely induce an arrhythmogenic substrate before significant immune responses. Adenovirus is a leading cause of myocarditis, but due to species specificity, models of infection are lacking, and it is not understood how adenoviral infection may underlie sudden cardiac arrest. Mouse adenovirus type-3 was previously reported as cardiotropic, yet it has not been utilized to understand the mechanisms of cardiac infection and pathology.

METHODS

We have developed mouse adenovirus type-3 infection as a model to investigate acute cardiac infection and molecular alterations to the infected heart before an appreciable immune response or gross cardiomyopathy.

RESULTS

Optical mapping of infected hearts exposes decreases in conduction velocity concomitant with increased Cx43Ser368 phosphorylation, a residue known to regulate gap junction function. Hearts from animals harboring a phospho-null mutation at Cx43Ser368 are protected against mouse adenovirus type-3-induced conduction velocity slowing. Additional to gap junction alterations, patch clamping of mouse adenovirus type-3-infected adult mouse ventricular cardiomyocytes reveals prolonged action potential duration as a result of decreased IK1 and IKs current density. Turning to human systems, we find human adenovirus type-5 increases phosphorylation of Cx43Ser368 and disrupts synchrony in human induced pluripotent stem cell-derived cardiomyocytes, indicating common mechanisms with our mouse whole heart and adult cardiomyocyte data.

CONCLUSIONS

Together, these findings demonstrate that adenoviral infection creates an arrhythmogenic substrate through direct targeting of gap junction and ion channel function in the heart. Such alterations are known to precipitate arrhythmias and likely contribute to sudden cardiac death in acutely infected patients.

Determinants of electrical propagation and propagation block in Arrhythmogenic Cardiomyopathy

Gap junction and ion channel remodeling occur early in Arrhythmogenic Cardiomyopathy (ACM), but their pathogenic consequences have not been elucidated. Here, we identified the arrhythmogenic substrate, consisting of propagation slowing and conduction block, in ACM models expressing two different desmosomal gene variants. Neonatal rat ventricular myocytes were transduced to express variants in genes encoding desmosomal proteins plakoglobin or plakophilin-2. Studies were performed in engineered cells and anisotropic tissues to quantify changes in conduction velocity, formation of unidirectional propagation, cell-cell electrical coupling, and ion currents. Conduction velocity decreased by 71% and 63% in the two ACM models. SB216763, an inhibitor of glycogen synthase kinase-3 beta, restored conduction velocity to near normal levels. Compared to control, both ACM models showed greater propensity for unidirectional conduction block, which increased further at greater stimulation frequencies. Cell-cell electrical conductance measured in cell pairs was reduced by 86% and 87% in the two ACM models. Computer modeling showed close correspondence between simulated and experimentally determined changes in conduction velocity. The simulation identified that reduced cell-cell electrical coupling was the dominant factor leading to slow conduction, while the combination of reduced cell-cell electrical coupling, reduced sodium current and inward rectifier potassium current explained the development of unidirectional block. Expression of two different ACM variants markedly reduced cell-cell electrical coupling and conduction velocity, and greatly increased the likelihood of developing unidirectional block - both key features of arrhythmogenesis. This study provides the first quantitative analysis of cellular electrophysiological changes leading to the substrate of reentrant arrhythmias in early stage ACM.

Unraveling Impacts of Chamber-Specific Differences in Intercalated Disc Ultrastructure and Molecular Organization on Cardiac Conduction

BACKGROUND

Propagation of action potentials through the heart coordinates the heartbeat. Thus, intercalated discs, specialized cell-cell contact sites that provide electrical and mechanical coupling between cardiomyocytes, are an important target for study. Impaired propagation leads to arrhythmias in many pathologies, where intercalated disc remodeling is a common finding, hence the importance and urgency of understanding propagation dependence on intercalated disc structure. Conventional modeling approaches cannot predict changes in propagation elicited by perturbations that alter intercalated disc ultrastructure or molecular organization, because of lack of quantitative structural data at subcellular through nano scales.

OBJECTIVES

This study sought to quantify intercalated disc structure at these spatial scales in the healthy adult mouse heart and relate them to chamber-specific properties of propagation as a precursor to understanding the effects of pathological intercalated disc remodeling.

METHODS

Using super-resolution light microscopy, electron microscopy, and computational image analysis, we provide here the first ever systematic, multiscale quantification of intercalated disc ultrastructure and molecular organization.

RESULTS

By incorporating these data into a rule-based model of cardiac tissue with realistic intercalated disc structure, and comparing model predictions of electrical propagation with experimental measures of conduction velocity, we reveal that atrial intercalated discs can support faster conduction than their ventricular counterparts, which is normally masked by interchamber differences in myocyte geometry. Further, we identify key ultrastructural and molecular organization features underpinning the ability of atrial intercalated discs to support faster conduction.

CONCLUSIONS

These data provide the first stepping stone to elucidating chamber-specific effects of pathological intercalated disc remodeling, as occurs in many arrhythmic diseases.

Extracellular Perinexal Separation Is a Principal Determinant of Cardiac Conduction

BACKGROUND

Cardiac conduction is understood to occur through gap junctions. Recent evidence supports ephaptic coupling as another mechanism of electrical communication in the heart. Conduction via gap junctions predicts a direct relationship between conduction velocity (CV) and bulk extracellular resistance. By contrast, ephaptic theory is premised on the existence of a biphasic relationship between CV and the volume of specialized extracellular clefts within intercalated discs such as the perinexus. Our objective was to determine the relationship between ventricular CV and structural changes to micro- and nanoscale extracellular spaces.

METHODS

Conduction and Cx43 (connexin43) protein expression were quantified from optically mapped guinea pig whole-heart preparations perfused with the osmotic agents albumin, mannitol, dextran 70 kDa, or dextran 2 MDa. Peak sodium current was quantified in isolated guinea pig ventricular myocytes. Extracellular resistance was quantified by impedance spectroscopy. Intercellular communication was assessed in a heterologous expression system with fluorescence recovery after photobleaching. Perinexal width was quantified from transmission electron micrographs.

RESULTS

CV primarily in the transverse direction of propagation was significantly reduced by mannitol and increased by albumin and both dextrans. The combination of albumin and dextran 70 kDa decreased CV relative to albumin alone. Extracellular resistance was reduced by mannitol, unchanged by albumin, and increased by both dextrans. Cx43 expression and conductance and peak sodium currents were not significantly altered by the osmotic agents. In response to osmotic agents, perinexal width, in order of narrowest to widest, was albumin with dextran 70 kDa; albumin or dextran 2 MDa; dextran 70 kDa or no osmotic agent, and mannitol. When compared in the same order, CV was biphasically related to perinexal width.

CONCLUSIONS

Cardiac conduction does not correlate with extracellular resistance but is biphasically related to perinexal separation, providing evidence that the relationship between CV and extracellular volume is determined by ephaptic mechanisms under conditions of normal gap junctional coupling.

Fibroblasts in heart scar tissue directly regulate cardiac excitability and arrhythmogenesis

After heart injury, dead heart muscle is replaced by scar tissue. Fibroblasts can electrically couple with myocytes, and changes in fibroblast membrane potential can lead to myocyte excitability, which suggests that fibroblast-myocyte coupling in scar tissue may be responsible for arrhythmogenesis. However, the physiologic relevance of electrical coupling of myocytes and fibroblasts and its impact on cardiac excitability in vivo have never been demonstrated. We genetically engineered a mouse that expresses the optogenetic cationic channel ChR2 (H134R) exclusively in cardiac fibroblasts. After myocardial infarction, optical stimulation of scar tissue elicited organ-wide cardiac excitation and induced arrhythmias in these animals. Complementing computational modeling with experimental approaches, we showed that gap junctional and ephaptic coupling, in a synergistic yet functionally redundant manner, excited myocytes coupled to fibroblasts.

Nonmyocytes as electrophysiological contributors to cardiac excitation and conduction

Although cardiac action potential (AP) generation and propagation have traditionally been attributed exclusively to cardiomyocytes (CM), other cell types in the heart are also capable of forming electrically conducting junctions. Interactions between CM and nonmyocytes (NM) enable and modulate each other's activity. This review provides an overview of the current understanding of heterocellular electrical communication in the heart. Although cardiac fibroblasts were initially thought to be electrical insulators, recent studies have demonstrated that they form functional electrical connections with CM in situ. Other NM, such as macrophages, have also been recognized as contributing to cardiac electrophysiology and arrhythmogenesis. Novel experimental tools have enabled the investigation of cell-specific activity patterns in native cardiac tissue, which is expected to yield exciting new insights into the development of novel or improved diagnostic and therapeutic strategies.

Sodium channel subpopulations with distinct biophysical properties and subcellular localization enhance cardiac conduction

Sodium (Na+) current is responsible for the rapid depolarization of cardiac myocytes that triggers the cardiac action potential upstroke. Recent studies have illustrated the presence of multiple pools of Na+ channels with distinct biophysical properties and subcellular localization, including clustering of channels at the intercalated disk and along the lateral membrane. Computational studies predict that Na+ channel clusters at the intercalated disk can regulate cardiac conduction via modulation of the narrow intercellular cleft between electrically coupled myocytes. However, these studies have primarily focused on the redistribution of Na+ channels between intercalated disk and lateral membranes and have not considered the distinct biophysical properties of the Na+ channel subpopulations. In this study, we use computational modeling to simulate computational models of single cardiac cells and one-dimensional cardiac tissues and predict the function of distinct Na+ channel subpopulations. Single-cell simulations predict that a subpopulation of Na+ channels with shifted steady-state activation and inactivation voltage dependency promotes an earlier action potential upstroke. In cardiac tissues that account for distinct subcellular spatial localization, simulations predict that shifted Na+ channels contribute to faster and more robust conduction in response to changes in tissue structure (i.e., cleft width), gap junctional coupling, and rapid pacing rates. Simulations predict that the intercalated disk-localized shifted Na+ channels contribute proportionally more to total Na+ charge than lateral membrane-localized Na+ channels. Importantly, our work supports the hypothesis that Na+ channel redistribution may be a critical mechanism by which cells can respond to perturbations to support fast and robust conduction.

Diversity of cells and signals in the cardiovascular system

This white paper is the outcome of the seventh UC Davis Cardiovascular Research Symposium on Systems Approach to Understanding Cardiovascular Disease and Arrhythmia. This biannual meeting aims to bring together leading experts in subfields of cardiovascular biomedicine to focus on topics of importance to the field. The theme of the 2022 Symposium was 'Cell Diversity in the Cardiovascular System, cell-autonomous and cell-cell signalling'. Experts in the field contributed their experimental and mathematical modelling perspectives and discussed emerging questions, controversies, and challenges in examining cell and signal diversity, co-ordination and interrelationships involved in cardiovascular function. This paper originates from the topics of formal presentations and informal discussions from the Symposium, which aimed to develop a holistic view of how the multiple cell types in the cardiovascular system integrate to influence cardiovascular function, disease progression and therapeutic strategies. The first section describes the major cell types (e.g. cardiomyocytes, vascular smooth muscle and endothelial cells, fibroblasts, neurons, immune cells, etc.) and the signals involved in cardiovascular function. The second section emphasizes the complexity at the subcellular, cellular and system levels in the context of cardiovascular development, ageing and disease. Finally, the third section surveys the technological innovations that allow the interrogation of this diversity and advancing our understanding of the integrated cardiovascular function and dysfunction.

Efficient, cell-based simulations of cardiac electrophysiology; The Kirchhoff Network Model (KNM)

Mathematical models based on homogenized representation of cardiac tissue have greatly improved our understanding of cardiac electrophysiology. However, these models are too coarse to investigate the dynamics at the level of the myocytes since the cells are not present in homogenized models. Recently, fine scale models have been proposed to allow for cell-level resolution of the dynamics, but these models are too computationally expensive to be used in applications like whole heart simulations of large animals. To address this issue, we propose a model that balances computational demands and physiological accuracy. The model is founded on Kirchhoff's current law, and represents every myocyte in the tissue. This allows specific properties to be assigned to individual cardiomyocytes, and other cell types like fibroblasts can be added to the model in an accurate manner while keeping the computing efforts reasonable.

Ephaptic Coupling as a Resolution to the Paradox of Action Potential Wave Speed and Discordant Alternans Spatial Scales in the Heart

Previous computer simulations have suggested that existing models of action potential wave propagation in the heart are not consistent with observed wave propagation behavior. Specifically, computer models cannot simultaneously reproduce the rapid wave speeds and small spatial scales of discordant alternans patterns measured experimentally in the same simulation. The discrepancy is important, because discordant alternans can be a key precursor to the development of abnormal and dangerous rapid rhythms in the heart. In this Letter, we show that this paradox can be resolved by allowing so-called ephaptic coupling to play a primary role in wave front propagation in place of conventional gap-junction coupling. With this modification, physiological wave speeds and small discordant alternans spatial scales both occur with gap-junction resistance values that are more in line with those observed in experiments. Our theory thus also provides support to the hypothesis that ephaptic coupling plays an important role in normal wave propagation.

Role of ephaptic coupling in discordant alternans domain sizes and action potential propagation in the heart

Discordant alternans, the spatially out-of-phase alternation of the durations of propagating action potentials in the heart, has been linked to the onset of fibrillation, a major cardiac rhythm disorder. The sizes of the regions, or domains, within which these alternations are synchronized are critical in this link. However, computer models employing standard gap junction-based coupling between cells have been unable to reproduce simultaneously the small domain sizes and rapid action potential propagation speeds seen in experiments. Here we use computational methods to show that rapid wave speeds and small domain sizes are possible when a more detailed model of intercellular coupling that accounts for so-called ephaptic effects is used. We provide evidence that the smaller domain sizes are possible, because different coupling strengths can exist on the wavefronts, for which both ephaptic and gap-junction coupling are involved, in contrast to the wavebacks, where only gap-junction coupling plays an active role. The differences in coupling strength are due to the high density of fast-inward (sodium) channels known to localize on the ends of cardiac cells, which are only active (and thus engage ephaptic coupling) during wavefront propagation. Thus, our results suggest that this distribution of fast-inward channels, as well as other factors responsible for the critical involvement of ephaptic coupling in wave propagation, including intercellular cleft spacing, play important roles in increasing the vulnerability of the heart to life-threatening tachyarrhythmias. Our results, combined with the absence of short-wavelength discordant alternans domains in standard gap-junction-dominated coupling models, also provide evidence that both gap-junction and ephaptic coupling are critical in wavefront propagation and waveback dynamics.

Unraveling Chamber-specific Differences in Intercalated Disc Ultrastructure and Molecular Organization and Their Impact on Cardiac Conduction

During each heartbeat, the propagation of action potentials through the heart coordinates the contraction of billions of individual cardiomyocytes and is thus, a critical life process. Unsurprisingly, intercalated discs, which are cell-cell contact sites specialized to provide electrical and mechanical coupling between adjacent cardiomyocytes, have been the focus of much investigation. Slowed or disrupted propagation leads to potentially life-threatening arrhythmias in a wide range of pathologies, where intercalated disc remodeling is a common finding. Hence, the importance and urgency of understanding intercalated disc structure and its influence on action potential propagation. Surprisingly, however, conventional modeling approaches cannot predict changes in propagation elicited by perturbations that alter intercalated disc ultrastructure or molecular organization, owing to lack of quantitative structural data at subcellular through nano scales. In order to address this critical gap in knowledge, we sought to quantify intercalated disc structure at these finer spatial scales in the healthy adult mouse heart and relate them to function in a chamber-specific manner as a precursor to understanding the impacts of pathological intercalated disc remodeling. Using super-resolution light microscopy, electron microscopy, and computational image analysis, we provide here the first ever systematic, multiscale quantification of intercalated disc ultrastructure and molecular organization. By incorporating these data into a rule-based model of cardiac tissue with realistic intercalated disc structure, and comparing model predictions of electrical propagation with experimental measures of conduction velocity, we reveal that atrial intercalated discs can support faster conduction than their ventricular counterparts, which is normally masked by inter-chamber differences in myocyte geometry. Further, we identify key ultrastructural and molecular organization features underpinning the ability of atrial intercalated discs to support faster conduction. These data provide the first stepping stone to elucidating chamber-specific impacts of pathological intercalated disc remodeling, as occurs in many arrhythmic diseases.

Subcellular diversity of Nav1.5 in cardiomyocytes: distinct functions, mechanisms and targets

In cardiomyocytes, the rapid depolarisation of the membrane potential is mediated by the α-subunit of the cardiac voltage-gated Na⁺ channel (Na_V 1.5), encoded by the gene SCN5A. This ion channel allows positively charged Na⁺ ions to enter the cardiomyocyte, resulting in the fast upstroke of the action potential and is therefore crucial for cardiac excitability and electrical propagation. This essential role is underscored by the fact that dysfunctional Na_V 1.5 is associated with high risk for arrhythmias and sudden cardiac death. However, development of therapeutic interventions regulating Na_V 1.5 has been limited due to the complexity of Na_V 1.5 structure and function and its diverse roles within the cardiomyocyte. In particular, research from the last decade has provided us with increased knowledge on the subcellular distribution of Na_V 1.5 as well as the proteins which it interacts with in distinct cardiomyocyte microdomains. We here review these insights, detailing the potential role of Na_V 1.5 within subcellular domains as well as its dysfunction in the setting of arrhythmia disorders. We furthermore provide an overview of current knowledge on the pathways involved in (microdomain-specific) trafficking of Na_V 1.5, and their potential as novel targets. Unravelling the complexity of Na_V 1.5 (dys)function may ultimately facilitate the development of therapeutic strategies aimed at preventing lethal arrhythmias. This is not only of importance for pathophysiological conditions where sodium current is specifically decreased within certain subcellular regions, such as in arrhythmogenic cardiomyopathy and Duchenne muscular dystrophy, but also for other acquired and inherited disorders associated with Na_V 1.5.

Role of Connexin 43 phosphorylation on Serine-368 by PKC in cardiac function and disease

Intercellular communication mediated by gap junction channels and hemichannels composed of Connexin 43 (Cx43) is vital for the propagation of electrical impulses through cardiomyocytes. The carboxyl terminal tail of Cx43 undergoes various post-translational modifications including phosphorylation of its Serine-368 (S368) residue. Protein Kinase C isozymes directly phosphorylate S368 to alter Cx43 function and stability through inducing conformational changes affecting channel permeability or promoting internalization and degradation to reduce intercellular communication between cardiomyocytes. Recent studies have implicated this PKC/Cx43-pS368 circuit in several cardiac-associated diseases. In this review, we describe the molecular and cellular basis of PKC-mediated Cx43 phosphorylation and discuss the implications of Cx43 S368 phosphorylation in the context of various cardiac diseases, such as cardiomyopathy, as well as the therapeutic potential of targeting this pathway.

Ephaptic Coupling Is a Mechanism of Conduction Reserve During Reduced Gap Junction Coupling

Many cardiac pathologies are associated with reduced gap junction (GJ) coupling, an important modulator of cardiac conduction velocity (CV). However, the relationship between phenotype and functional expression of the connexin GJ family of proteins is controversial. For example, a 50% reduction of GJ coupling has been shown to have little impact on myocardial CV due to a concept known as conduction reserve. This can be explained by the ephaptic coupling (EpC) theory whereby conduction is maintained by a combination of low GJ coupling and increased electrical fields generated in the sodium channel rich clefts between neighboring myocytes. At the same time, low GJ coupling may also increase intracellular charge accumulation within myocytes, resulting in a faster transmembrane potential rate of change during depolarization (dV/dt_max) that maintains macroscopic conduction. To provide insight into the prevalence of these two phenomena during pathological conditions, we investigated the relationship between EpC and charge accumulation within the setting of GJ remodeling using multicellular simulations and companion perfused mouse heart experiments. Conduction along a fiber of myocardial cells was simulated for a range of GJ conditions. The model incorporated intercellular variations, including GJ coupling conductance and distribution, cell-to-cell separation in the intercalated disc (perinexal width—W_P), and variations in sodium channel distribution. Perfused heart studies having conditions analogous to those of the simulations were performed using wild type mice and mice heterozygous null for the connexin gene Gja1. With insight from simulations, the relative contributions of EpC and charge accumulation on action potential parameters and conduction velocities were analyzed. Both simulation and experimental results support a common conclusion that low GJ coupling decreases and narrowing W_P increases the rate of the AP upstroke when sodium channels are densely expressed at the ends of myocytes, indicating that conduction reserve is more dependent on EpC than charge accumulation during GJ uncoupling.

Mechanisms of arrhythmia termination during acute myocardial ischemia: Role of ephaptic coupling and complex geometry of border zone

Myocardial ischemia occurs when blood flow to the heart is reduced, preventing the heart muscle from receiving enough oxygen required for survival. Several anatomical and electrophysiological changes occur at the ischemic core (IC) and border zone (BZ) during myocardial ischemia, for example, gap junctional remodeling, changes in ionic channel kinetics and electrophysiologic changes in cell excitability, which promote the development of cardiac arrhythmia. Ephaptic coupling (EpC), which is an electrical field effect developed in the shared cleft space between adjacent cells, has been suggested to rescue the conduction when gap junctions are impaired, such as myocardial ischemia. In this manuscript, we explored the impact of EpC, electrophysiological and anatomical components of myocardial ischemia on reentry termination during non-ischemic and ischemic condition. Our results indicated that EpC and BZ with complex geometry have opposite effects on the reentry termination. In particular, the presence of homogeneous EpC terminates reentry, whereas BZ with complex geometry alone facilitates reentry by producing wave break-up and alternating conduction block. The reentry is terminated in the presence of homogeneous or heterogeneous EpC despite the presence of complex geometry of the BZ, independent of the location of BZ. The inhibition of reentry can be attributed to a current-to-load mismatch. Our results points to an antiarrhythmic role of EpC and a pro-arrhythmic role of BZ with complex geometry.

Simulation of atrial fibrillation in a non-ohmic propagation model with dynamic gap junctions

Gap junctions exhibit nonlinear electrical properties that have been hypothesized to be relevant to arrhythmogenicity in a structurally remodeled tissue. Large-scale implementation of gap junction dynamics in 3D propagation models remains challenging. We aim to quantify the impact of nonlinear diffusion during episodes of arrhythmias simulated in a left atrial model. Homogenization of conduction properties in the presence of nonlinear gap junctions was performed by generalizing a previously developed mathematical framework. A monodomain model was solved in which conductivities were time-varying and depended on transjunctional potentials. Gap junction conductances were derived from a simplified Vogel-Weingart model with first-order gating and adjustable time constant. A bilayer interconnected cable model of the left atrium with 100  μm resolution was used. The diffusion matrix was recomputed at each time step according to the state of the gap junctions. Sinus rhythm and atrial fibrillation episodes were simulated in remodeled tissue substrates. Slow conduction was induced by reduced coupling and by diffuse or stringy fibrosis. Simulations starting from the same initial conditions were repeated with linear and nonlinear gap junctions. The discrepancy in activation times between the linear and nonlinear diffusion models was quantified. The results largely validated the linear approximation for conduction velocities >20 cm/s. In very slow conduction substrates, the discrepancy accumulated over time during atrial fibrillation, eventually leading to qualitative differences in propagation patterns, while keeping the descriptive statistics, such as cycle lengths, unchanged. The discrepancy growth rate was increased by impaired conduction, fibrosis, conduction heterogeneity, lateral uncoupling, fast gap junction time constant, and steeper action potential duration restitution.

Automaticity in ventricular myocyte cell pairs with ephaptic and gap junction coupling

Spontaneous electrical activity, or automaticity, in the heart is required for normal physiological function. However, irregular automaticity, in particular, originating from the ventricles, can trigger life-threatening cardiac arrhythmias. Thus, understanding mechanisms of automaticity and synchronization is critical. Recent work has proposed that excitable cells coupled via a shared narrow extracellular cleft can mediate coupling, i.e., ephaptic coupling, that promotes automaticity in cell pairs. However, the dynamics of these coupled cells incorporating both ephaptic and gap junction coupling has not been explored. Here, we show that automaticity and synchronization robustly emerges via a Hopf bifurcation from either (i) increasing the fraction of inward rectifying potassium channels (carrying the I_K1 current) at the junctional membrane or (ii) by decreasing the cleft volume. Furthermore, we explore how heterogeneity in the fraction of potassium channels between coupled cells can produce automaticity of both cells or neither cell, or more rarely in only one cell (i.e., automaticity without synchronization). Interestingly, gap junction coupling generally has minor effects, with only slight changes in regions of parameter space of automaticity. This work provides insight into potentially new mechanisms that promote spontaneous activity and, thus, triggers for arrhythmias in ventricular tissue.

Computational modeling of aberrant electrical activity following remuscularization with intramyocardially injected pluripotent stem cell-derived cardiomyocytes

Acute engraftment arrhythmias (EAs) remain a serious complication of remuscularization therapy. Preliminary evidence suggests that a focal source underlies these EAs stemming from the automaticity of immature pluripotent stem cell-derived cardiomyocytes (PSC-CMs) in nascent myocardial grafts. How these EAs arise though during early engraftment remains unclear. In a series of in silico experiments, we probed the origin of EAs-exploring aspects of altered impulse formation and altered impulse propagation within nascent PSC-CM grafts and at the host-graft interface. To account for poor gap junctional coupling during early PSC-CM engraftment, the voltage dependence of gap junctions and the possibility of ephaptic coupling were incorporated. Inspired by cardiac development, we also studied the contributions of another feature of immature PSC-CMs, circumferential sodium channel (NaCh) distribution in PSC-CMs. Ectopic propagations emerged from nascent grafts of immature PSC-CMs at a rate of <96 bpm. Source-sink effects dictated this rate and contributed to intermittent capture between host and graft. Moreover, ectopic beats emerged from dynamically changing sites along the host-graft interface. The latter arose in part because circumferential NaCh distribution in PSC-CMs contributed to preferential conduction slowing and block of electrical impulses from host to graft myocardium. We conclude that additional mechanisms, in addition to focal ones, contribute to EAs and recognize that their relative contributions are dynamic across the engraftment process.

Protein and mRNA Quantification in Small Samples of Human-Induced Pluripotent Stem Cell-Derived Cardiomyocytes in 96-Well Microplates

We describe a method for protein quantification and for mRNA quantification in small sample quantities of human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs). Demonstrated here is how the capillary-based protein detection system Wes™ by ProteinSimple and the Power SYBR™ Green Cells-to-CT™ Kit by Invitrogen can be applied to individual samples in a 96-well microplate format and thus made compatible with high-throughput (HT) cardiomyocyte assays. As an example of the usage, we illustrate that Cx43 protein and GJA1 mRNA levels in hiPSC-CMs are enhanced when the optogenetic actuator, channelrodopsin-2 (ChR2), is genetically expressed in them. Instructions are presented for cell culture and lysate preparations from hiPSC-CMs, along with optimized parameter settings and experimental protocol steps. Strategies to optimize primary antibody concentrations as well as ways for signal normalization are discussed, i.e., antibody multiplexing and total protein assay. The sensitivity of both the Wes and Cells-to-CT kit enables protein and mRNA quantification in a HT format, which is important when dealing with precious small samples. In addition to being able to handle small cardiomyocyte samples, these streamlined and semi-automated processes enable quick mechanistic analysis.

Hypernatremia and intercalated disc edema synergistically exacerbate long-QT syndrome type 3 phenotype

Cardiac voltage-gated sodium channel gain-of-function prolongs repolarization in the long-QT syndrome type 3 (LQT3). Previous studies suggest that narrowing the perinexus within the intercalated disc, leading to rapid sodium depletion, attenuates LQT3-associated action potential duration (APD) prolongation. However, it remains unknown whether extracellular sodium concentration modulates APD prolongation during sodium channel gain-of-function. We hypothesized that elevated extracellular sodium concentration and widened perinexus synergistically prolong APD in LQT3. LQT3 was induced with sea anemone toxin (ATXII) in Langendorff-perfused guinea pig hearts (n = 34). Sodium concentration was increased from 145 to 160 mM. Perinexal expansion was induced with mannitol or the sodium channel β1-subunit adhesion domain antagonist (βadp1). Epicardial ventricular action potentials were optically mapped. Individual and combined effects of varying clefts and sodium concentrations were simulated in a computational model. With ATXII, both mannitol and βadp1 significantly widened the perinexus and prolonged APD, respectively. The elevated sodium concentration alone significantly prolonged APD as well. Importantly, the combination of elevated sodium concentration and perinexal widening synergistically prolonged APD. Computational modeling results were consistent with animal experiments. Concurrently elevating extracellular sodium and increasing intercalated disc edema prolongs repolarization more than the individual interventions alone in LQT3. This synergistic effect suggests an important clinical implication that hypernatremia in the presence of cardiac edema can markedly increase LQT3-associated APD prolongation. Therefore, to our knowledge, this is the first study to provide evidence of a tractable and effective strategy to mitigate LQT3 phenotype by means of managing sodium levels and preventing cardiac edema in patients.NEW & NOTEWORTHY This is the first study to demonstrate that the long-QT syndrome type 3 (LQT3) phenotype can be exacerbated or concealed by regulating extracellular sodium concentrations and/or the intercalated disc separation. The animal experiments and computational modeling in the current study reveal a critically important clinical implication: sodium dysregulation in the presence of edema within the intercalated disc can markedly increase the risk of arrhythmia in LQT3. These findings strongly suggest that maintaining extracellular sodium within normal physiological limits may be an effective and inexpensive therapeutic option for patients with congenital or acquired sodium channel gain-of-function diseases.

Cellular Size, Gap Junctions, and Sodium Channel Properties Govern Developmental Changes in Cardiac Conduction

Electrical conduction in cardiac ventricular tissue is regulated via sodium (Na⁺) channels and gap junctions (GJs). We and others have recently shown that Na⁺channels preferentially localize at the site of cell-cell junctions, the intercalated disc (ID), in adult cardiac tissue, facilitating coupling via the formation of intercellular Na⁺nanodomains, also termed ephaptic coupling (EpC). Several properties governing EpC vary with age, including Na⁺channel and GJ expression and distribution and cell size. Prior work has shown that neonatal cardiomyocytes have immature IDs with Na⁺channels and GJs diffusively distributed throughout the sarcolemma, while adult cells have mature IDs with preferentially localized Na⁺channels and GJs. In this study, we perform an in silico investigation of key age-dependent properties to determine developmental regulation of cardiac conduction. Simulations predict that conduction velocity (CV) biphasically depends on cell size, depending on the strength of GJ coupling. Total cell Na⁺channel conductance is predictive of CV in cardiac tissue with high GJ coupling, but not correlated with CV for low GJ coupling. We find that ephaptic effects are greatest for larger cells with low GJ coupling typically associated with intermediate developmental stages. Finally, simulations illustrate how variability in cellular properties during different developmental stages can result in a range of possible CV values, with a narrow range for both neonatal and adult myocardium but a much wider range for an intermediate developmental stage. Thus, we find that developmental changes predict associated changes in cardiac conduction.

Initiation and entrainment of multicellular automaticity via diffusion limited extracellular domains

Electrically excitable cells often spontaneously and synchronously depolarize in vitro and in vivo preparations. It remains unclear how cells entrain and autorhythmically activate above the intrinsic mean activation frequency of isolated cells with or without pacemaking mechanisms. Recent studies suggest that cyclic ion accumulation and depletion in diffusion-limited extracellular volumes modulate electrophysiology by ephaptic mechanisms (nongap junction or synaptic coupling). This report explores how potassium accumulation and depletion in a restricted extracellular domain induces spontaneous action potentials in two different computational models of excitable cells without gap junctional coupling: Hodgkin-Huxley and Luo-Rudy. Importantly, neither model will spontaneously activate on its own without external stimuli. Simulations demonstrate that cells sharing a diffusion-limited extracellular compartment can become autorhythmic and entrained despite intercellular electrical heterogeneity. Autorhythmic frequency is modulated by the cleft volume and potassium fluxes through the cleft. Additionally, inexcitable cells can suppress or induce autorhythmic activity in an excitable cell via a shared cleft. Diffusion-limited shared clefts can also entrain repolarization. Critically, this model predicts a mechanism by which diffusion-limited shared clefts can initiate, entrain, and modulate multicellular automaticity in the absence of gap junctions.

Cardiac radiotherapy induces electrical conduction reprogramming in the absence of transmural fibrosis

Cardiac radiotherapy (RT) may be effective in treating heart failure (HF) patients with refractory ventricular tachycardia (VT). The previously proposed mechanism of radiation-induced fibrosis does not explain the rapidity and magnitude with which VT reduction occurs clinically. Here, we demonstrate in hearts from RT patients that radiation does not achieve transmural fibrosis within the timeframe of VT reduction. Electrophysiologic assessment of irradiated murine hearts reveals a persistent supraphysiologic electrical phenotype, mediated by increases in NaV1.5 and Cx43. By sequencing and transgenic approaches, we identify Notch signaling as a mechanistic contributor to NaV1.5 upregulation after RT. Clinically, RT was associated with increased NaV1.5 expression in 1 of 1 explanted heart. On electrocardiogram (ECG), post-RT QRS durations were shortened in 13 of 19 patients and lengthened in 5 patients. Collectively, this study provides evidence for radiation-induced reprogramming of cardiac conduction as a potential treatment strategy for arrhythmia management in VT patients.

Localization of Na+ channel clusters in narrowed perinexi of gap junctions enhances cardiac impulse transmission via ephaptic coupling: a model study

Abstract

It has been proposed that when gap junctional coupling is reduced in cardiac tissue, action potential propagation can be supported via ephaptic coupling, a mechanism mediated by negative electric potentials occurring in narrow intercellular clefts of intercalated discs (IDs). Recent studies showed that sodium (Na⁺) channels form clusters near gap junction plaques in nanodomains called perinexi, where the ID cleft is even narrower. To examine the electrophysiological relevance of Na⁺ channel clusters being located in perinexi, we developed a 3D finite element model of two longitudinally abutting cardiomyocytes, with a central Na⁺ channel cluster on the ID membranes. When this cluster was located in the perinexus of a closely positioned gap junction plaque, varying perinexal width greatly modulated impulse transmission from one cell to the other, with narrow perinexi potentiating ephaptic coupling. This modulation occurred via the interplay of Na⁺ currents, extracellular potentials in the cleft and patterns of current flow within the cleft. In contrast, when the Na⁺ channel cluster was located remotely from the gap junction plaque, this modulation by perinexus width largely disappeared. Interestingly, the Na⁺ current in the ID membrane of the pre‐junctional cell switched from inward to outward during excitation, thus contributing ions to the activating channels on the post‐junctional ID membrane. In conclusion, these results indicate that the localization of Na⁺ channel clusters in the perinexi of gap junction plaques is crucial for ephaptic coupling, which is furthermore greatly modulated by perinexal width. These findings are relevant for a comprehensive understanding of cardiac excitation.

Key points

Ephaptic coupling is a cardiac conduction mechanism involving nanoscale‐level interactions between the sodium (Na⁺) current and the extracellular potential in narrow intercalated disc clefts.

When gap junctional coupling is reduced, ephaptic coupling acts in conjunction with the classical cardiac conduction mechanism based on gap junctional current flow.

In intercalated discs, Na⁺ channels form clusters that are preferentially located in the periphery of gap junction plaques, in nanodomains known as perinexi, but the electrophysiological role of these perinexi has never been examined.

In our new 3D finite element model of two cardiac cells abutting each other with their intercalated discs, a Na⁺ channel cluster located inside a narrowed perinexus facilitated impulse transmission via ephaptic coupling.

Our simulations demonstrate the role of narrowed perinexi as privileged sites for ephaptic coupling in pathological situations when gap junctional coupling is decreased.

Targeting the Microtubule EB1-CLASP2 Complex Modulates NaV1.5 at Intercalated Discs

[Figure: see text].

Intercalated disk nanoscale structure regulates cardiac conduction

The intercalated disk (ID) is a specialized subcellular region that provides electrical and mechanical connections between myocytes in the heart. The ID has a clearly defined passive role in cardiac tissue, transmitting mechanical forces and electrical currents between cells. Recent studies have shown that Na+ channels, the primary current responsible for cardiac excitation, are preferentially localized at the ID, particularly within nanodomains such as the gap junction-adjacent perinexus and mechanical junction-associated adhesion-excitability nodes, and that perturbations of ID structure alter cardiac conduction. This suggests that the ID may play an important, active role in regulating conduction. However, the structures of the ID and intercellular cleft are not well characterized and, to date, no models have incorporated the influence of ID structure on conduction in cardiac tissue. In this study, we developed an approach to generate realistic finite element model (FEM) meshes replicating nanoscale of the ID structure, based on experimental measurements from transmission electron microscopy images. We then integrated measurements of the intercellular cleft electrical conductivity, derived from the FEM meshes, into a novel cardiac tissue model formulation. FEM-based calculations predict that the distribution of cleft conductances is sensitive to regional changes in ID structure, specifically the intermembrane separation and gap junction distribution. Tissue-scale simulations predict that ID structural heterogeneity leads to significant spatial variation in electrical polarization within the intercellular cleft. Importantly, we found that this heterogeneous cleft polarization regulates conduction by desynchronizing the activation of postjunctional Na+ currents. Additionally, these heterogeneities lead to a weaker dependence of conduction velocity on gap junctional coupling, compared with prior modeling formulations that neglect or simplify ID structure. Further, we found that disruption of local ID nanodomains can either slow or enhance conduction, depending on gap junctional coupling strength. Our study therefore suggests that ID nanoscale structure can play a significant role in regulating cardiac conduction.

Coupling between cardiac cells-An important determinant of electrical impulse propagation and arrhythmogenesis

Cardiac arrhythmias are an important cause of sudden cardiac death-a devastating manifestation of many underlying causes, such as heart failure and ischemic heart disease leading to ventricular tachyarrhythmias and ventricular fibrillation, and atrial fibrillation causing cerebral embolism. Cardiac electrical propagation is a main factor in the initiation and maintenance of cardiac arrhythmias. In the heart, gap junctions are the basic unit at the cellular level that host intercellular low-resistance channels for the diffusion of ions and small regulatory molecules. The dual voltage clamp technique enabled the direct measurement of electrical conductance between cells and recording of single gap junction channel openings. The rapid turnover of gap junction channels at the intercalated disk implicates a highly dynamic process of trafficking and internalization of gap junction connexons. Recently, non-canonical roles of gap junction proteins have been discovered in mitochondria function, cytoskeletal organization, trafficking, and cardiac rescue. At the tissue level, we explain the concepts of linear propagation and safety factor based on the model of linear cellular structure. Working myocardium is adequately represented as a discontinuous cellular network characterized by cellular anisotropy and connective tissue heterogeneity. Electrical propagation in discontinuous cellular networks reflects an interplay of three main factors: cell-to-cell electrical coupling, flow of electrical charge through the ion channels, and the microscopic tissue structure. This review provides a state-of-the-art update of the cardiac gap junction channels and their role in cardiac electrical impulse propagation and highlights a combined approach of genetics, cell biology, and physics in modern cardiac electrophysiology.

Connexins in the Heart: Regulation, Function and Involvement in Cardiac Disease

Connexins are a family of transmembrane proteins that play a key role in cardiac physiology. Gap junctional channels put into contact the cytoplasms of connected cardiomyocytes, allowing the existence of electrical coupling. However, in addition to this fundamental role, connexins are also involved in cardiomyocyte death and survival. Thus, chemical coupling through gap junctions plays a key role in the spreading of injury between connected cells. Moreover, in addition to their involvement in cell-to-cell communication, mounting evidence indicates that connexins have additional gap junction-independent functions. Opening of unopposed hemichannels, located at the lateral surface of cardiomyocytes, may compromise cell homeostasis and may be involved in ischemia/reperfusion injury. In addition, connexins located at non-canonical cell structures, including mitochondria and the nucleus, have been demonstrated to be involved in cardioprotection and in regulation of cell growth and differentiation. In this review, we will provide, first, an overview on connexin biology, including their synthesis and degradation, their regulation and their interactions. Then, we will conduct an in-depth examination of the role of connexins in cardiac pathophysiology, including new findings regarding their involvement in myocardial ischemia/reperfusion injury, cardiac fibrosis, gene transcription or signaling regulation.

Mechanisms underlying age-associated manifestation of cardiac sodium channel gain-of-function

Cardiac action potentials are initiated by sodium ion (Na+) influx through voltage-gated Na+ channels. Na+ channel gain-of-function (GOF) can arise in inherited conditions due to mutations in the gene encoding the cardiac Na+ channel, such as Long QT syndrome type 3 (LQT3). LQT3 can be a "concealed" disease, as patients with LQT3-associated mutations can remain asymptomatic until later in life; however, arrhythmias can also arise early in life in LQT3 patients, demonstrating a complex age-associated manifestation. We and others recently demonstrated that cardiac Na+ channels preferentially localize at the intercalated disc (ID) in adult cardiac tissue, which facilitates ephaptic coupling and formation of intercellular Na+ nanodomains that regulate pro-arrhythmic early afterdepolarization (EAD) formation in tissue with Na+ channel GOF. Several properties related to ephaptic coupling vary with age, such as cell size and Na+ channel and gap junction (GJ) expression and distribution: neonatal cells have immature IDs, with Na+ channels and GJs primarily diffusively distributed, while adult myocytes have mature IDs with preferentially localized Na+ channels and GJs. Here, we perform an in silico study varying critical age-dependent parameters to investigate mechanisms underlying age-associated manifestation of Na+ channel GOF in a model of guinea pig cardiac tissue. Simulations predict that total Na+ current conductance is a critical factor in action potential duration (APD) prolongation. We find a complex cell size/ Na+ channel expression relationship: increases in cell size (without concurrent increases in Na+ channel expression) suppress EAD formation, while increases in Na+ channel expression (without concurrent increases in cell size) promotes EAD formation. Finally, simulations with neonatal and early age-associated parameters predict normal APD with minimal dependence on intercellular cleft width; however, variability in cellular properties can lead to EADs presenting in early developmental stages. In contrast, for adult-associated parameters, EAD formation is highly dependent on cleft width, consistent with a mechanism underlying the age-associated manifestation of the Na+ channel GOF.

The conduction velocity-potassium relationship in the heart is modulated by sodium and calcium

The relationship between cardiac conduction velocity (CV) and extracellular potassium (K⁺) is biphasic, with modest hyperkalemia increasing CV and severe hyperkalemia slowing CV. Recent studies from our group suggest that elevating extracellular sodium (Na⁺) and calcium (Ca²⁺) can enhance CV by an extracellular pathway parallel to gap junctional coupling (GJC) called ephaptic coupling that can occur in the gap junction adjacent perinexus. However, it remains unknown whether these same interventions modulate CV as a function of K⁺. We hypothesize that Na⁺, Ca²⁺, and GJC can attenuate conduction slowing consequent to severe hyperkalemia. Elevating Ca²⁺ from 1.25 to 2.00 mM significantly narrowed perinexal width measured by transmission electron microscopy. Optically mapped, Langendorff-perfused guinea pig hearts perfused with increasing K⁺ revealed the expected biphasic CV-K⁺ relationship during perfusion with different Na⁺ and Ca²⁺ concentrations. Neither elevating Na⁺ nor Ca²⁺ alone consistently modulated the positive slope of CV-K⁺ or conduction slowing at 10-mM K⁺; however, combined Na⁺ and Ca²⁺ elevation significantly mitigated conduction slowing at 10-mM K⁺. Pharmacologic GJC inhibition with 30-μM carbenoxolone slowed CV without changing the shape of CV-K⁺ curves. A computational model of CV predicted that elevating Na⁺ and narrowing clefts between myocytes, as occur with perinexal narrowing, reduces the positive and negative slopes of the CV-K⁺ relationship but do not support a primary role of GJC or sodium channel conductance. These data demonstrate that combinatorial effects of Na⁺ and Ca²⁺ differentially modulate conduction during hyperkalemia, and enhancing determinants of ephaptic coupling may attenuate conduction changes in a variety of physiologic conditions.

Identifying Drug Response by Combining Measurements of the Membrane Potential, the Cytosolic Calcium Concentration, and the Extracellular Potential in Microphysiological Systems

Cardiomyocytes derived from human induced pluripotent stem cells (hiPSC-CMs) offer a new means to study and understand the human cardiac action potential, and can give key insight into how compounds may interact with important molecular pathways to destabilize the electrical function of the heart. Important features of the action potential can be readily measured using standard experimental techniques, such as the use of voltage sensitive dyes and fluorescent genetic reporters to estimate transmembrane potentials and cytosolic calcium concentrations. Using previously introduced computational procedures, such measurements can be used to estimate the current density of major ion channels present in hiPSC-CMs, and how compounds may alter their behavior. However, due to the limitations of optical recordings, resolving the sodium current remains difficult from these data. Here we show that if these optical measurements are complemented with observations of the extracellular potential using multi electrode arrays (MEAs), we can accurately estimate the current density of the sodium channels. This inversion of the sodium current relies on observation of the conduction velocity which turns out to be straightforwardly computed using measurements of extracellular waves across the electrodes. The combined data including the membrane potential, the cytosolic calcium concentration and the extracellular potential further opens up for the possibility of accurately estimating the effect of novel drugs applied to hiPSC-CMs.

Anisotropic Cardiac Conduction

Anisotropy is the property of directional dependence. In cardiac tissue, conduction velocity is anisotropic and its orientation is determined by myocyte direction. Cell shape and size, excitability, myocardial fibrosis, gap junction distribution and function are all considered to contribute to anisotropic conduction. In disease states, anisotropic conduction may be enhanced, and is implicated, in the genesis of pathological arrhythmias. The principal mechanism responsible for enhanced anisotropy in disease remains uncertain. Possible contributors include changes in cellular excitability, changes in gap junction distribution or function and cellular uncoupling through interstitial fibrosis. It has recently been demonstrated that myocyte orientation may be identified using diffusion tensor magnetic resonance imaging in explanted hearts, and multisite pacing protocols have been proposed to estimate myocyte orientation and anisotropic conduction in vivo. These tools have the potential to contribute to the understanding of the role of myocyte disarray and anisotropic conduction in arrhythmic states.

Vascular endothelial growth factor promotes atrial arrhythmias by inducing acute intercalated disk remodeling

Atrial fibrillation (AF) is the most common arrhythmia and is associated with inflammation. AF patients have elevated levels of inflammatory cytokines known to promote vascular leak, such as vascular endothelial growth factor A (VEGF). However, the contribution of vascular leak and consequent cardiac edema to the genesis of atrial arrhythmias remains unknown. Previous work suggests that interstitial edema in the heart can acutely promote ventricular arrhythmias by disrupting ventricular myocyte intercalated disk (ID) nanodomains rich in cardiac sodium channels (Na_V1.5) and slowing cardiac conduction. Interestingly, similar disruption of ID nanodomains has been identified in atrial samples from AF patients. Therefore, we tested the hypothesis that VEGF-induced vascular leak can acutely increase atrial arrhythmia susceptibility by disrupting ID nanodomains and slowing atrial conduction. Treatment of murine hearts with VEGF (30–60 min, at clinically relevant levels) prolonged the electrocardiographic P wave and increased susceptibility to burst pacing-induced atrial arrhythmias. Optical voltage mapping revealed slower atrial conduction following VEGF treatment (10 ± 0.4 cm/s vs. 21 ± 1 cm/s at baseline, p < 0.05). Transmission electron microscopy revealed increased intermembrane spacing at ID sites adjacent to gap junctions (GJs; 64 ± 9 nm versus 17 ± 1 nm in controls, p < 0.05), as well as sites next to mechanical junctions (MJs; 63 ± 4 nm versus 27 ± 2 nm in controls, p < 0.05) in VEGF–treated hearts relative to controls. Importantly, super-resolution microscopy and quantitative image analysis revealed reorganization of Na_V1.5 away from dense clusters localized near GJs and MJs to a more diffuse distribution throughout the ID. Taken together, these data suggest that VEGF can acutely predispose otherwise normal hearts to atrial arrhythmias by dynamically disrupting Na_V1.5-rich ID nanodomains and slowing atrial conduction. These data highlight inflammation-induced vascular leak as a potential factor in the development and progression of AF.

Specific decreasing of Na+ channel expression on the lateral membrane of cardiomyocytes causes fatal arrhythmias in Brugada syndrome

Reduced cardiac sodium (Na+) channel current (INa) resulting from the loss-of-function of Na+ channel is a major cause of lethal arrhythmias in Brugada syndrome (BrS). Inspired by previous experimental studies which showed that in heart diseases INa was reduced along with expression changes in Na+ channel within myocytes, we hypothesized that the local decrease in INa caused by the alteration in Na+ channel expression in myocytes leads to the occurrence of phase-2 reentry, the major triggering mechanism of lethal arrhythmias in BrS. We constructed in silico human ventricular myocardial strand and ring models, and examined whether the Na+ channel expression changes in each myocyte cause the phase-2 reentry in BrS. Reducing Na+ channel expression in the lateral membrane of each myocyte caused not only the notch-and-dome but also loss-of-dome type action potentials and slowed conduction, both of which are typically observed in BrS patients. Furthermore, the selective reduction in Na+ channels on the lateral membrane of each myocyte together with spatial tissue heterogeneity of Na+ channel expression caused the phase-2 reentry and phase-2 reentry-mediated reentrant arrhythmias. Our data suggest that the BrS phenotype is strongly influenced by expression abnormalities as well as genetic abnormalities of Na+ channels.

Modeling the Interactions Between Sodium Channels Provides Insight Into the Negative Dominance of Certain Channel Mutations

BACKGROUND

Nav1.5 cardiac Na+ channel mutations can cause arrhythmogenic syndromes. Some of these mutations exert a dominant negative effect on wild-type channels. Recent studies showed that Na+ channels can dimerize, allowing coupled gating. This leads to the hypothesis that allosteric interactions between Na+ channels modulate their function and that these interactions may contribute to the negative dominance of certain mutations.

METHODS

To investigate how allosteric interactions affect microscopic and macroscopic channel function, we developed a modeling paradigm in which Markovian models of two channels are combined. Allosteric interactions are incorporated by modifying the free energies of the composite states and/or barriers between states.

RESULTS

Simulations using two generic 2-state models (C-O, closed-open) revealed that increasing the free energy of the composite states CO/OC leads to coupled gating. Simulations using two 3-state models (closed-open-inactivated) revealed that coupled closings must also involve interactions between further composite states. Using two 6-state cardiac Na+ channel models, we replicated previous experimental results mainly by increasing the energies of the CO/OC states and lowering the energy barriers between the CO/OC and the CO/OO states. The channel model was then modified to simulate a negative dominant mutation (Nav1.5 p.L325R). Simulations of homodimers and heterodimers in the presence and absence of interactions showed that the interactions with the variant channel impair the opening of the wild-type channel and thus contribute to negative dominance.

CONCLUSION

Our new modeling framework recapitulates qualitatively previous experimental observations and helps identifying possible interaction mechanisms between ion channels.

Disease Phenotypes and Mechanisms of iPSC-Derived Cardiomyocytes From Brugada Syndrome Patients With a Loss-of-Function SCN5A Mutation

Brugada syndrome (BrS) is one of the major causes of sudden cardiac death in young people, while the underlying mechanisms are not completely understood. Here, we investigated the pathophysiological phenotypes and mechanisms using induced pluripotent stem cell (iPSC)-derived cardiomyocytes (CMs) from two BrS patients (BrS-CMs) carrying a heterozygous SCN5A mutation p.S1812X. Compared to CMs derived from healthy controls (Ctrl-CMs), BrS-CMs displayed a 50% reduction of I_Na density, a 69.5% reduction of Na_V1.5 expression, and the impaired localization of Na_V1.5 and connexin 43 (Cx43) at the cell surface. BrS-CMs exhibited reduced action potential (AP) upstroke velocity and conduction slowing. The I_to in BrS-CMs was significantly augmented, and the I_CaL window current probability was increased. Our data indicate that the electrophysiological mechanisms underlying arrhythmia in BrS-CMs may involve both depolarization and repolarization disorders. Cilostazol and milrinone showed dramatic inhibitions of I_to in BrS-CMs and alleviated the arrhythmic activity, suggesting their therapeutic potential for BrS patients.

Syncytium cell growth increases Kir2.1 contribution in human iPSC-cardiomyocytes

Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) enable cardiotoxicity testing and personalized medicine. However, their maturity is of concern, including relatively depolarized resting membrane potential and more spontaneous activity compared with adult cardiomyocytes, implicating low or lacking inward rectifier potassium current (I_k1). Here, protein quantification confirms Kir2.1 expression in hiPSC-CM syncytia, albeit several times lower than in adult heart tissue. We find that hiPSC-CM culture density influences Kir2.1 expression at the mRNA level (potassium inwardly rectifying channel subfamily J member 2) and at the protein level and its associated electrophysiology phenotype. Namely, all-optical cardiac electrophysiology and pharmacological treatments reveal reduction of spontaneous and irregular activity and increase in action potential upstroke in denser cultures. Blocking I_k1-like currents with BaCl₂ increased spontaneous frequency and blunted action potential upstrokes during pacing in a dose-dependent manner only in the highest-density cultures, in line with I_k1’s role in regulating the resting membrane potential. Our results emphasize the importance of syncytial growth of hiPSC-CMs for more physiologically relevant phenotype and the power of all-optical electrophysiology to study cardiomyocytes in their multicellular setting.

NEW & NOTEWORTHY We identify cell culture density and cell-cell contact as an important factor in determining the expression of a key ion channel at the transcriptional and the protein levels, KCNJ2/Kir2.1, and its contribution to the electrophysiology of human induced pluripotent stem cell-derived cardiomyocytes. Our results indicate that studies on isolated cells, out of tissue context, may underestimate the cellular ion channel properties being characterized.

Single-Molecule Localization of the Cardiac Voltage-Gated Sodium Channel Reveals Different Modes of Reorganization at Cardiomyocyte Membrane Domains

BACKGROUND

Mutations in the gene encoding the cardiac voltage-gated sodium channel Na_v1.5 cause various cardiac arrhythmias. This variety may arise from different determinants of Na_v1.5 expression between cardiomyocyte domains. At the lateral membrane and T-tubules, Na_v1.5 localization and function remain insufficiently characterized.

METHODS

We used novel single-molecule localization microscopy and computational modeling to define nanoscale features of Na_v1.5 localization and distribution at the lateral membrane, the lateral membrane groove, and T-tubules in cardiomyocytes from wild-type (N=3), dystrophin-deficient (mdx; N=3) mice, and mice expressing C-terminally truncated Na_v1.5 (ΔSIV; N=3). We moreover assessed T-tubules sodium current by recording whole-cell sodium currents in control (N=5) and detubulated (N=5) wild-type cardiomyocytes.

RESULTS

We show that Na_v1.5 organizes as distinct clusters in the groove and T-tubules which density, distribution, and organization partially depend on SIV and dystrophin. We found that overall reduction in Na_v1.5 expression in mdx and ΔSIV cells results in a nonuniform redistribution with Na_v1.5 being specifically reduced at the groove of ΔSIV and increased in T-tubules of mdx cardiomyocytes. A T-tubules sodium current could, however, not be demonstrated.

CONCLUSIONS

Na_v1.5 mutations may site-specifically affect Na_v1.5 localization and distribution at the lateral membrane and T-tubules, depending on site-specific interacting proteins. Future research efforts should elucidate the functional consequences of this redistribution.

Identification of potential candidate genes and pathways in atrioventricular nodal reentry tachycardia by whole-exome sequencing

Background

Atrioventricular nodal reentry tachycardia (AVNRT) is the most common manifestation of paroxysmal supraventricular tachycardia (PSVT). Increasing data have indicated familial clustering and participation of genetic factors in AVNRT, and no pathogenic genes related to AVNRT have been reported.

Methods

Whole‐exome sequencing (WES) was performed in 82 patients with AVNRT and 100 controls. Reference genes, genome‐wide association analysis, gene‐based collapsing, and pathway enrichment analysis were performed. A protein‐protein interaction (PPI) network was then established; WES database in the UK Biobank and one only genetic study of AVNRT in Denmark were used for external data validation.

Results

Among 95 reference genes, 126 rare variants in 48 genes were identified in the cases (minor allele frequency < 0.001). Gene‐based collapsing analysis and pathway enrichment analysis revealed six functional pathways related to AVNRT as with neuronal system/neurotransmitter release cycles and ion channel/cardiac conduction among the top 30 enriched pathways, and then 36 candidate pathogenic genes were selected. By combining with PPI analysis, 10 candidate genes were identified, including RYR2, NOS1, SCN1A, CFTR, EPHB4, ROBO1, PRKAG2, MMP2, ASPH, and ABCC8. From the UK Biobank database, 18 genes from candidate genes including SCN1A, PRKAG2, NOS1, and CFTR had rare variants in arrhythmias, and the rare variants in PIK3CB, GAD2, and HIP1R were in patients with PSVT. Moreover, one rare variant of RYR2 (c.4652A > G, p.Asn1551Ser) in our study was also detected in the Danish study. Considering the gene functional roles and external data validation, the most likely candidate genes were SCN1A, PRKAG2, RYR2, CFTR, NOS1, PIK3CB, GAD2, and HIP1R.

Conclusion

The preliminary results first revealed potential candidate genes such as SCN1A, PRKAG2, RYR2, CFTR, NOS1, PIK3CB, GAD2, and HIP1R, and the pathways mediated by these genes, including neuronal system/neurotransmitter release cycles or ion channels/cardiac conduction, might be involved in AVNRT.

Intercellular Sodium Regulates Repolarization in Cardiac Tissue with Sodium Channel Gain of Function

In cardiac myocytes, action potentials are initiated by an influx of sodium (Na⁺) ions via voltage-gated Na⁺ channels. Na⁺ channel gain of function (GOF), arising in both inherited conditions associated with mutation in the gene encoding the Na⁺ channel and acquired conditions associated with heart failure, ischemia, and atrial fibrillation, enhance Na⁺ influx, generating a late Na⁺ current that prolongs action potential duration (APD) and triggering proarrhythmic early afterdepolarizations (EADs). Recent studies have shown that Na⁺ channels are highly clustered at the myocyte intercalated disk, facilitating formation of Na⁺ nanodomains in the intercellular cleft between cells. Simulations from our group have recently predicted that narrowing the width of the intercellular cleft can suppress APD prolongation and EADs in the presence of Na⁺ channel mutations because of increased intercellular cleft Na⁺ ion depletion. In this study, we investigate the effects of modulating multiple extracellular spaces, specifically the intercellular cleft and bulk interstitial space, in a novel computational model and experimentally via osmotic agents albumin, dextran 70, and mannitol. We perform optical mapping and transmission electron microscopy in a drug-induced (sea anemone toxin, ATXII) Na⁺ channel GOF isolated heart model and modulate extracellular spaces via osmotic agents. Single-cell patch-clamp experiments confirmed that the osmotic agents individually do not enhance late Na⁺ current. Both experiments and simulations are consistent with the conclusion that intercellular cleft narrowing or expansion regulates APD prolongation; in contrast, modulating the bulk interstitial space has negligible effects on repolarization. Thus, we predict that intercellular cleft Na⁺ nanodomain formation and collapse critically regulates cardiac repolarization in the setting of Na⁺ channel GOF.

Real-Time In Vivo Imaging of Mouse Left Ventricle Reveals Fluctuating Movements of the Intercalated Discs

Myocardial contraction is initiated by action potential propagation through the conduction system of the heart. It has been thought that connexin 43 in the gap junctions (GJ) within the intercalated disc (ID) provides direct electric connectivity between cardiomyocytes (electronic conduction). However, recent studies challenge this view by providing evidence that the mechanosensitive cardiac sodium channels Na_v1.5 localized in perinexii at the GJ edge play an important role in spreading action potentials between neighboring cells (ephaptic conduction). In the present study, we performed real-time confocal imaging of the CellMask-stained ID in the living mouse heart in vivo. We found that the ID structure was not rigid. Instead, we observed marked flexing of the ID during propagation of contraction from cell to cell. The variation in ID length was between ~30 and ~42 μm (i.e., magnitude of change, ~30%). In contrast, tracking of α-actinin-AcGFP revealed a comparatively small change in the lateral dimension of the transitional junction near the ID (i.e., magnitude of change, ~20%). The present findings suggest that, when the heart is at work, mechanostress across the perinexii may activate Na_v1.5 by promoting ephaptic conduction in coordination with electronic conduction, and, thereby, efficiently transmitting excitation-contraction coupling between cardiomyocytes.

Interplay between ephaptic coupling and complex geometry of border zone during acute myocardial ischemia: Effect on arrhythmogeneity

Acute myocardial ischemia is an imbalance between myocardial blood supply and demand, which is caused by the cessation of blood flow within the heart resulting from an obstruction in one of the major coronary arteries. A severe blockage may result in a region of nonperfused tissue known as ischemic core (IC). As a result, a border zone (BZ) between perfused and nonperfused regions is created due to differences in blood and oxygen supplies. Recent experimental findings reveal a complex "finger-like" geometry in BZ; however, its effect on arrhythmogenicity is not clear. Ephaptic coupling, which relies on the intercalated disk between cell ends, has been suggested to play an active role in mediating intercellular electrical communication when gap junctions are impaired. In this paper, we explored the interplay between ephaptic coupling and the geometry of BZ on action potential propagation across the ischemic region. Our study shows that ephaptic coupling can greatly suppress the occurrence of a conduction block, which points to its beneficial effect. The beneficial effect of ephaptic coupling is more evident in BZ with the "finger-like" geometry. In addition, the complex geometry of BZ, i.e., more frequent, deeper, and wider "fingers," promotes the conduction through the ischemic region. In contrast, the larger size of IC impedes the cardiac conduction across the ischemic region. Our results also show that ephaptic coupling promotes the impact of the complex geometry of BZ on signal propagation; however, it inhibits the impact of IC size.

Non-ohmic tissue conduction in cardiac electrophysiology: Upscaling the non-linear voltage-dependent conductance of gap junctions

Gap junctions are key mediators of intercellular communication in cardiac tissue, and their function is vital to sustaining normal cardiac electrical activity. Conduction through gap junctions strongly depends on the hemichannel arrangement and transjunctional voltage, rendering the intercellular conductance highly non-Ohmic, particularly under steady-state regimes of conduction. Despite this marked non-linear behavior, current tissue-level models of cardiac conduction are rooted in the assumption that gap-junctions conductance is constant (Ohmic), which results in inaccurate predictions of electrical propagation, particularly in the low junctional-coupling regime observed under pathological conditions. In this work, we present a novel non-Ohmic homogenization model (NOHM) of cardiac conduction that is suitable to tissue-scale simulations. Using non-linear homogenization theory, we develop a conductivity model that seamlessly upscales the voltage-dependent conductance of gap junctions, without the need of explicitly modeling gap junctions. The NOHM model allows for the simulation of electrical propagation in tissue-level cardiac domains that accurately resemble that of cell-based microscopic models for a wide range of junctional coupling scenarios, recovering key conduction features at a fraction of the computational complexity. A unique feature of the NOHM model is the possibility of upscaling the response of non-symmetric gap-junction conductance distributions, which result in conduction velocities that strongly depend on the direction of propagation, thus allowing to model the normal and retrograde conduction observed in certain regions of the heart. We envision that the NOHM model will enable organ-level simulations that are informed by sub- and inter-cellular mechanisms, delivering an accurate and predictive in-silico tool for understanding the heart function. Codes are available for download at https://github.com/dehurtado/NonOhmicConduction.

The heart relies on the propagation of electrical impulses that are mediated gap junctions, whose conduction properties vary depending on the transjunctional voltage. Despite this non-linear feature, current mathematical models assume that cardiac tissue behaves like an Ohmic (linear) material, thus delivering inaccurate results when simulated in a computer. Here we present a novel mathematical multiscale model that explicitly includes the non-Ohmic response of gap junctions in its predictions. Our results show that the proposed model recovers important conduction features modulated by gap junctions at a fraction of the computational complexity. This contribution represents an important step towards constructing computer models of a whole heart that can predict organ-level behavior in reasonable computing times.

Rat caval vein myocardium undergoes changes in conduction characteristics during postnatal ontogenesis

The electrophysiological properties of the superior vena cava (SVC) myocardium, which is considered a minor source of atrial arrhythmias, were studied in this study during postnatal development. Conduction properties were investigated in spontaneously active and electrically paced SVC preparations obtained from 7-60-day-old male Wistar rats using optical mapping and microelectrode techniques. The presence of high-conductance connexin 43 (Cx43) was evaluated in SVC cross-sections using immunofluorescence. It was found that SVC myocardium is excitable, electrically coupled with the atrial tissue, and conducts excitation waves at all stages of postnatal development. However, the conduction velocity (CV) of excitation and action potential (AP) upstroke velocity in SVC were significantly lower in neonatal than in adult animals and increased with postnatal maturation. Connexins Cx43 were identified in both neonatal and adult rat SVC myocardium; however, the abundance of Cx43 was significantly less in neonates. The gap junction uncoupler octanol affected conduction more profound in the neonatal than in adult SVC. We demonstrated for the first time that the conduction characteristics of SVC myocardium change from a slow-conduction (nodal) to a high-conduction (working) phenotype during postnatal ontogenesis. An age-related CV increase may occur due to changes of AP characteristics, electrical coupling, and Cx43 presence in SVC cardiomyocyte membranes. Observed changes may contribute to the low proarrhythmicity of adult caval vein cardiac tissue, while pre- or postnatal developmental abnormalities that delay the establishment of the working conduction phenotype may facilitate SVC proarrhythmia.

Properties of cardiac conduction in a cell-based computational model

The conduction of electrical signals through cardiac tissue is essential for maintaining the function of the heart, and conduction abnormalities are known to potentially lead to life-threatening arrhythmias. The properties of cardiac conduction have therefore been the topic of intense study for decades, but a number of questions related to the mechanisms of conduction still remain unresolved. In this paper, we demonstrate how the so-called EMI model may be used to study some of these open questions. In the EMI model, the extracellular space, the cell membrane, the intracellular space and the cell connections are all represented as separate parts of the computational domain, and the model therefore allows for study of local properties that are hard to represent in the classical homogenized bidomain or monodomain models commonly used to study cardiac conduction. We conclude that a non-uniform sodium channel distribution increases the conduction velocity and decreases the time delays over gap junctions of reduced coupling in the EMI model simulations. We also present a theoretical optimal cell length with respect to conduction velocity and consider the possibility of ephaptic coupling (i.e. cell-to-cell coupling through the extracellular potential) acting as an alternative or supporting mechanism to gap junction coupling. We conclude that for a non-uniform distribution of sodium channels and a sufficiently small intercellular distance, ephaptic coupling can influence the dynamics of the sodium channels and potentially provide cell-to-cell coupling when the gap junction connection is absent.

The electrochemical wave traversing the heart during every beat is essential for cardiac pumping function and supply of blood to the body. Understanding the stability of this wave is crucial to understanding how lethal arrhythmias are generated. Despite this importance, our knowledge of the physical determinants of wave propagation are still evolving. One particular challenge has been the lack of accurate mathematical models of conduction at the cellular level. Because cardiac muscle is an electrical syncytium, in which direct charge transfer between cells drives wave propagation, classical bidomain and monodomain tissue models employ a homogenized approximation of this process. This approximation is not valid at the length scale of single cells, and prevents any analysis of how cellular structures impact cardiac conduction. Instead, so-called microdomain models must be used for these questions. Here we utilize a recently developed modelling framework that is well suited to represent small collections of cells. By applying this framework, we show that concentration of sodium channels at the longitudinal borders of myocytes accelerates cardiac conduction. We also demonstrate that when juxtaposed cells are sufficiently close, this non-uniform distribution induces large ephaptic currents, which contribute to intercellular coupling.

Modulating cardiac conduction during metabolic ischemia with perfusate sodium and calcium in guinea pig hearts

We previously demonstrated that altering extracellular sodium (Na_o) and calcium (Ca_o) can modulate a form of electrical communication between cardiomyocytes termed "ephaptic coupling" (EpC), especially during loss of gap junction coupling. We hypothesized that altering Na_o and Ca_o modulates conduction velocity (CV) and arrhythmic burden during ischemia. Electrophysiology was quantified by optically mapping Langendorff-perfused guinea pig ventricles with modified Na_o (147 or 155 mM) and Ca_o (1.25 or 2.0 mM) during 30 min of simulated metabolic ischemia (pH 6.5, anoxia, aglycemia). Gap junction-adjacent perinexal width ( W_P), a candidate cardiac ephapse, and connexin (Cx)43 protein expression and Cx43 phosphorylation at S368 were quantified by transmission electron microscopy and Western immunoblot analysis, respectively. Metabolic ischemia slowed CV in hearts perfused with 147 mM Na_o and 2.0 mM Ca_o; however, theoretically increasing EpC with 155 mM Na_o was arrhythmogenic, and CV could not be measured. Reducing Ca_o to 1.25 mM expanded W_P, as expected during ischemia, consistent with reduced EpC, but attenuated CV slowing while delaying arrhythmia onset. These results were further supported by osmotically reducing W_P with albumin, which exacerbated CV slowing and increased early arrhythmias during ischemia, whereas mannitol expanded W_P, permitted conduction, and delayed the onset of arrhythmias. Cx43 expression patterns during the various interventions insufficiently correlated with observed CV changes and arrhythmic burden. In conclusion, decreasing perfusate calcium during metabolic ischemia enhances perinexal expansion, attenuates conduction slowing, and delays arrhythmias. Thus, perinexal expansion may be cardioprotective during metabolic ischemia. NEW & NOTEWORTHY This study demonstrates, for the first time, that modulating perfusate ion composition can alter cardiac electrophysiology during simulated metabolic ischemia.

Conduction in cardiac tissue. Historical reflections

Two hypotheses have been proposed to explain propagation of the action potential in heart. According to the gap junction hypothesis local short-circuit currents pass from the proximal depolarized cell to the distal inactive cell via gap junctions and are responsible for the depolarization of the distal cell. In the ephapse hypothesis the depolarization of the proximal cell generates an electrical field in the narrow cleft between cells resulting in depolarization beyond threshold of the distal cell. Measurements of length constant, free diffusion of 42 K, local currents between cells, existence of high-conductance gap junctions led to the conclusion that heart muscle is a functional syncytium. Propagation of the action potential, however, is not uniform but anisotropic and discontinuous; it can be also unidirectional. These findings are strong arguments in favor of the gap junction thesis. They do not exclude, as predicted by theoretical calculations, that in conditions of an abnormal fall in gap junction conductance ephaptic conduction takes over. In this last case, definitive experimental confirmation is still required. See also: https://doi.org/10.14814/phy2.13861 & https://doi.org/10.14814/phy2.13862.

The role of the gap junction perinexus in cardiac conduction: Potential as a novel anti-arrhythmic drug target

Cardiovascular disease remains the single largest cause of natural death in the United States, with a significant cause of mortality associated with cardiac arrhythmias. Presently, options for treating and preventing myocardial electrical dysfunction, including sudden cardiac death, are limited. Recent studies have indicated that conduction of electrical activation in the heart may have an ephaptic component, wherein intercellular coupling occurs via electrochemical signaling across narrow extracellular clefts between cardiomyocytes. The perinexus is a 100-200 nm-wide stretch of closely apposed membrane directly adjacent to connexin 43 gap junctions. Electron and super-resolution microscopy studies, as well as biochemical analyses, have provided evidence that perinexal nanodomains may be candidate structures for facilitating ephaptic coupling. This work has included characterization of the perinexus as a region of close inter-membrane contact between cardiomyocytes (<30 nm) containing dense clusters of voltage-gated sodium channels. Here, we review what is known about perinexal structure and function and the potential that the perinexus may have novel and pivotal roles in disorders of cardiac conduction. Of particular interest is the prospect that cell adhesion mediated by the cardiac sodium channel β subunit (Scn1b) may be a novel anti-arrhythmic target.