Size Matters: Ryanodine Receptor Cluster Size Heterogeneity Potentiates Calcium Waves

Role of ryanodine receptor cooperativity in Ca2+-wave-mediated triggered activity in cardiomyocytes

Ca2+ waves are known to trigger delayed after-depolarizations that can cause malignant cardiac arrhythmias. However, modelling Ca2+ waves using physiologically realistic models has remained a major challenge. Existing models with low Ca2+ sensitivity of ryanodine receptors (RyRs) necessitate large release currents, leading to an unrealistically large Ca2+ transient amplitude incompatible with the experimental observations. Consequently, current physiologically detailed models of delayed after-depolarizations resort to unrealistic cell architectures to produce Ca2+ waves with a normal Ca2+ transient amplitude. Here, we address these challenges by incorporating RyR cooperativity into a physiologically detailed model with a realistic cell architecture. We represent RyR cooperativity phenomenologically through a Hill coefficient within the sigmoid function of RyR open probability. Simulations in permeabilized myocytes with high Ca2+ sensitivity reveal that a sufficiently large Hill coefficient is required for Ca2+ wave propagation via the fire-diffuse-fire mechanism. In intact myocytes, propagating Ca2+ waves can occur only within an intermediate Hill coefficient range. Within this range, the spark rate is neither too low, enabling Ca2+ wave propagation, nor too high, allowing for the maintenance of a high sarcoplasmic reticulum load during diastole of the action potential. Moreover, this model successfully replicates other experimentally observed manifestations of Ca2+-wave-mediated triggered activity, including phase 2 and phase 3 early after-depolarizations and high-frequency voltage-Ca2+ oscillations. These oscillations feature an elevated take-off potential with depolarization mediated by the L-type Ca2+ current. The model also sheds light on the roles of luminal gating of RyRs and the mobile buffer ATP in the genesis of these arrhythmogenic phenomena. KEY POINTS: Existing mathematical models of Ca2+ waves use an excessively large Ca2+-release current or unrealistic diffusive coupling between release units. Our physiologically realistic model, using a Hill coefficient in the ryanodine receptor (RyR) gating function to represent RyR cooperativity, addresses these limitations and generates organized Ca2+ waves at Hill coefficients ranging from ∼5 to 10, as opposed to the traditional value of 2. This range of Hill coefficients gives a spark rate neither too low, thereby enabling Ca2+ wave propagation, nor too high, allowing for the maintenance of a high sarcoplasmic reticulum load during the plateau phase of the action potential. Additionally, the model generates Ca2+-wave-mediated phase 2 and phase 3 early after-depolarizations, and coupled membrane voltage with Ca2+ oscillations mediated by the L-type Ca2+ current. This study suggests that pharmacologically targeting RyR cooperativity could be a promising strategy for treating cardiac arrhythmias linked to Ca2+-wave-mediated triggered activity.

Structure-Function Relationship of the Ryanodine Receptor Cluster Network in Sinoatrial Node Cells

The rate of spontaneous action potentials (APs) generated by sinoatrial node cells (SANC) is regulated by local Ca2+ release (LCR) from the sarcoplasmic reticulum via Ca2+ release channels (ryanodine receptors, RyRs). LCR events propagate and self-organize within the network of RyR clusters (Ca release units, CRUs) via Ca-induced-Ca-release (CICR) that depends on CRU sizes and locations: While larger CRUs generate stronger release signals, the network's topology governs signal diffusion and propagation. This study used super-resolution structured illumination microscopy to image the 3D network of CRUs in rabbit SANC. The peripheral CRUs formed a spatial mesh, reflecting the cell surface geometry. Two distinct subpopulations of CRUs were identified within each cell, with size distributions conforming to a two-component Gamma mixture model. Furthermore, neighboring CRUs exhibited repulsive behavior. Functional properties of the CRU network were further examined in a novel numerical SANC model developed using our experimental data. Model simulations revealed that heterogeneities in both CRU sizes and locations facilitate CICR and increase the AP firing rate in a cooperative manner. However, these heterogeneities reduce the effect of β-adrenergic stimulation in terms of its relative change in AP firing rate. The presence of heterogeneities in both sizes and locations allows SANC to reach higher absolute AP firing rates during β-adrenergic stimulation. Thus, the CICR facilitation by heterogeneities in CRU sizes and locations regulates and optimizes cardiac pacemaker cell operation under various physiological conditions. Dysfunction of this optimization could be a key factor in heart rate reserve decline in aging and disease.

Modulation of heteromeric glycine receptor function through high concentration clustering

Ion channels are targeted by many drugs for treating neurological, musculoskeletal, renal and other diseases. These drugs bind to and alter the function of individual channels to achieve desired therapeutic effects. However, many ion channels function in high concentration clusters in their native environment. It is unclear if and how clustering modulates ion channel function. Human heteromeric glycine receptors (GlyRs) are the major inhibitory neurotransmitter receptors in the spinal cord and are active targets for developing chronic pain medications. We show that the α2β heteromeric GlyR assembles with the master postsynaptic scaffolding gephyrin (GPHN) into micron-sized clustered at the plasma membrane after heterologous expression. The inhibitory trans- synaptic adhesion protein neuroligin-2 (NL2) further increases both the cluster sizes and GlyR concentration. The apparent glycine affinity increases monotonically as a function of GlyR concentration but not with cluster size. We also show that ligand re-binding to adjacent GlyRs alters kinetics but not chemical equilibrium. A positively charged N- terminus sequence of the GlyR β subunit was further identified essential for glycine affinity modulation through clustering. Taken together, we propose a mechanism where clustering enhances local electrostatic potential, which in turn concentrates ions and ligands, modulating the function of GlyR. This mechanism is likely universal across ion channel clusters found ubiquitously in biology and provides new perspectives in possible pharmaceutical development.

Ca2+ dynamics in interstitial cells: foundational mechanisms for the motor patterns in the gastrointestinal tract

The gastrointestinal (GI) tract displays multiple motor patterns that move nutrients and wastes through the body. Smooth muscle cells (SMCs) provide the forces necessary for GI motility, but interstitial cells, electrically coupled to SMCs, tune SMC excitability, transduce inputs from enteric motor neurons, and generate pacemaker activity that underlies major motor patterns, such as peristalsis and segmentation. The interstitial cells regulating SMCs are interstitial cells of Cajal (ICC) and PDGF receptor (PDGFR)α+ cells. Together these cells form the SIP syncytium. ICC and PDGFRα+ cells express signature Ca2+-dependent conductances: ICC express Ca2+-activated Cl- channels, encoded by Ano1, that generate inward current, and PDGFRα+ cells express Ca2+-activated K+ channels, encoded by Kcnn3, that generate outward current. The open probabilities of interstitial cell conductances are controlled by Ca2+ release from the endoplasmic reticulum. The resulting Ca2+ transients occur spontaneously in a stochastic manner. Ca2+ transients in ICC induce spontaneous transient inward currents and spontaneous transient depolarizations (STDs). Neurotransmission increases or decreases Ca2+ transients, and the resulting depolarizing or hyperpolarizing responses conduct to other cells in the SIP syncytium. In pacemaker ICC, STDs activate voltage-dependent Ca2+ influx, which initiates a cluster of Ca2+ transients and sustains activation of ANO1 channels and depolarization during slow waves. Regulation of GI motility has traditionally been described as neurogenic and myogenic. Recent advances in understanding Ca2+ handling mechanisms in interstitial cells and how these mechanisms influence motor patterns of the GI tract suggest that the term "myogenic" should be replaced by the term "SIPgenic," as this review discusses.

Clustering properties of the cardiac ryanodine receptor in health and heart failure

The cardiac ryanodine receptor (RyR2) is an intracellular Ca2+ release channel vital for the function of the heart. Physiologically, RyR2 is triggered to release Ca2+ from the sarcoplasmic reticulum (SR) which enables cardiac contraction; however, spontaneous Ca2+ leak from RyR2 has been implicated in the pathophysiology of heart failure (HF). RyR2 channels have been well documented to assemble into clusters within the SR membrane, with the organisation of RyR2 clusters recently gaining interest as a mechanism by which the occurrence of pathological Ca2+ leak is regulated, including in HF. In this review, we explain the terminology relating to key nanoscale RyR2 clustering properties as both single clusters and functionally grouped Ca2+ release units, with a focus on the advancements in super-resolution imaging approaches which have enabled the detailed study of cluster organisation. Further, we discuss proposed mechanisms for modulating RyR2 channel organisation and the debate regarding the potential impact of cluster organisation on Ca2+ leak activity. Finally, recent experimental evidence investigating the nanoscale remodelling and functional alterations of RyR2 clusters in HF is discussed with consideration of the clinical implications.

Elementary intracellular Ca signals approximated as a transition of release channel system from a metastable state

Cardiac muscle contraction is initiated by an elementary Ca signal (called Ca spark) which is achieved by collective action of Ca release channels in a cluster. The mechanism of this synchronization remains uncertain. We approached Ca spark activation as an emergent phenomenon of an interactive system of release channels. We constructed a weakly lumped Markov chain that applies an Ising model formalism to such release channel clusters and probable open channel configurations and demonstrated that spark activation is described as a system transition from a metastable to an absorbing state, analogous to the pressure required to overcome surface tension in bubble formation. This yielded quantitative estimates of the spark generation probability as a function of various system parameters. We performed numerical simulations to find spark probabilities as a function of sarcoplasmic reticulum Ca concentration, obtaining similar values for spark activation threshold as our analytic model, as well as those reported in experimental studies. Our parametric sensitivity analyses also showed that the spark activation threshold decreased as Ca sensitivity of RyR activation and RyR cluster size increased.

Ca2+ spark latency and control of intrinsic Ca2+ release dyssynchrony in rat cardiac ventricular muscle cells

Cardiac excitation-contraction coupling (ECC) depends on Ca2+ release from intracellular stores via ryanodine receptors (RyRs) triggered by L-type Ca2+ channels (LCCs). Uncertain numbers of RyRs and LCCs form 'couplons' whose activation produces Ca2+ sparks, which summate to form a cell-wide Ca2+ transient that switches on contraction. Voltage (Vm) changes during the action potential (AP) and stochasticity in channel gating should create variability in Ca2+ spark timing, but Ca2+ transient wavefronts have remarkable uniformity. To examine how this is achieved, we measured the Vm-dependence of evoked Ca2+ spark probability (Pspark) and latency over a wide voltage range in rat ventricular cells. With depolarising steps, Ca2+ spark latency showed a U-shaped Vm-dependence, while repolarising steps from 50 mV produced Ca2+ spark latencies that increased monotonically with Vm. A computer model based on reported channel gating and geometry reproduced our experimental data and revealed a likely RyR:LCC stoichiometry of ∼ 5:1 for the Ca2+ spark initiating complex (IC). Using the experimental AP waveform, the model revealed a high coupling fidelity (Pcpl ∼ 0.5) between each LCC opening and IC activation. The presence of ∼ 4 ICs per couplon reduced Ca2+ spark latency and increased Pspark to match experimental data. Variability in AP release timing is less than that seen with voltage steps because the AP overshoot and later repolarization decrease Pspark due to effects on LCC flux and LCC deactivation respectively. This work provides a framework for explaining the Vm- and time-dependence of Pspark, and indicates how ion channel dispersion in disease can contribute to dyssynchrony in Ca2+ release.

The ryanodine receptor microdomain in cardiomyocytes

The ryanodine receptor type 2 (RyR) is a key player in Ca2+ handling during excitation-contraction coupling. During each heartbeat, RyR channels are responsible for linking the action potential with the contractile machinery of the cardiomyocyte by releasing Ca2+ from the sarcoplasmic reticulum. RyR function is fine-tuned by associated signalling molecules, arrangement in clusters and subcellular localization. These parameters together define RyR function within microdomains and are subject to disease remodelling. This review describes the latest findings on RyR microdomain organization, the alterations with disease which result in increased subcellular heterogeneity and emergence of microdomains with enhanced arrhythmogenic potential, and presents novel technologies that guide future research to study and target RyR channels within specific microdomains.

Unveiling the intricacies of intracellular Ca2+ regulation in the heart

Recent studies have provided valuable insight into the key mechanisms contributing to the spatiotemporal regulation of intracellular Ca2+ release and Ca2+ signaling in the heart. In this research highlight, we focus on the latest findings published in Biophysical Journal examining the structural organization of Ca2+ handling proteins and assessing the functional aspects of intracellular Ca2+ regulation in health and the detrimental consequences of Ca2+ dysregulation in disease. These important studies pave the way for future mechanistic investigations and multiscale understanding of Ca2+ signaling in the heart.

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.

Calcium oscillation on homogeneous and heterogeneous networks of ryanodine receptor

Calcium oscillation is an important calcium homeostasis, imbalance of which is the key mechanism of initiation and progression of many major diseases. The formation and maintenance of calcium homeostasis are closely related to the spatial distribution of calcium channels on endoplasmic reticulum, whose complex structure was unveiled by recent observations with superresolution imaging techniques. In the current paper, a theoretical framework is established by abstracting the spatial distribution of the calcium channels as a nonlinear biological complex network with calcium channels as nodes and Ca^{2+} as edges. A dynamical model for a ryanodine receptor (RyR) is adopted to investigate the effect of spatial distribution on calcium oscillation. The mean-field model can be well reproduced from the complete graph and dense Erdös-Rényi network. The synchronization of RyRs is found important to generate a global calcium oscillation. Below a critical density of the Erdös-Rényi or BaraBási-Albert network, the amplitude and interspike interval decrease rapidly with the end of disappearance of oscillation due to the desynchronization. The clique graph with a cluster structure cannot produce a global oscillation due to the failure of synchronization between clusters. A more realistic geometric network is constructed in a two-dimensional plane based on the experimental information about the RyR arrangement of clusters and the frequency distribution of cluster sizes. Different from the clique graph, the global oscillation can be generated with reasonable parameters on the geometric network. The simulation also suggests that existence of small clusters and rogue RyRs plays an important role in the maintenance of global calcium oscillation through keeping synchronization between large clusters. Such results support the heterogeneous distribution of RyRs with different-size clusters, which is helpful to understand recent observations with superresolution nanoscale imaging techniques. The current theoretical framework can also be extent to investigate other phenomena in calcium signal transduction.

The role of subunit cooperativity on ryanodine receptor 2 calcium signaling

The ryanodine receptor type 2 (RyR2) is composed of four subunits that control calcium (Ca) release in cardiac cells. RyR2 serves primarily as a Ca sensor and can respond to rapid sub-millisecond pulses of Ca while remaining shut at resting concentrations. However, it is not known how the four subunits interact for the RyR2 to function as an effective Ca sensor. To address this question, and to understand the role of subunit cooperativity in Ca-mediated signal transduction, we have developed a computational model of the RyR2 composed of four interacting subunits. We first analyze the statistical properties of a single RyR2 tetramer, where each subunit can exist in a closed or open conformation. Our findings indicate that the number of subunits in the open state is a crucial parameter that dictates RyR2 kinetics. We find that three or four open subunits are required for the RyR2 to harness cooperative interactions to respond to sub-millisecond changes in Ca, while at the same time remaining shut at the resting Ca levels in the cardiac cell. If the required number of open subunits is lowered to one or two, the RyR2 cannot serve as a robust Ca sensor, as the large cooperativity required to stabilize the closed state prevents channel activation. Using this four-subunit model, we analyze the kinetics of Ca release from a RyR2 cluster. We show that the closure of a cluster of RyR2 channels is highly sensitive to the balance of cooperative interactions between closed and open subunits. Based on this result, we analyze how specific interactions between RyR2 subunits can induce persistent Ca leak from the sarcoplasmic reticulum (SR), which is believed to be arrhythmogenic. Thus, these results provide a framework to analyze how a pharmacologic or genetic modification of RyR2 subunit cooperativity can induce abnormal Ca cycling that can potentially lead to life-threatening arrhythmias.

Enhanced calcium release at specialised surface sites compensates for reduced t-tubule density in neonatal sheep atrial myocytes

Cardiac myocytes rely on transverse (t)-tubules to facilitate a rapid rise in calcium throughout the cell. However, despite their importance in triggering synchronous Ca²⁺ release, t-tubules are highly labile structures. They develop postnatally, increase in density during exercise training and are lost in diseases such as heart failure (HF). In the majority of settings, an absence of t-tubules decreases function. Here we show that despite reduced t-tubule density due to immature t-tubules, the newborn atrium is highly specialised to maintain Ca²⁺ release. To compensate for fewer t-tubules triggering a central rise in Ca²⁺, Ca²⁺ release at sites on the cell surface is enhanced in the newborn, exceeding that at all Ca²⁺ release sites in the adult. Using electron and super resolution microscopy to investigate myocyte ultrastructure, we found that newborn atrial cells had enlarged surface sarcoplasmic reticulum and larger, more closely spaced surface and central ryanodine receptor clusters. We suggest that these adaptations mediate enhanced Ca²⁺ release at the sarcolemma and aid propagation to compensate for reduced t-tubule density in the neonatal atrium.

Biological noise is a key determinant of the reproducibility and adaptability of cardiac pacemaking and EC coupling

Each heartbeat begins with the generation of an action potential in pacemaking cells in the sinoatrial node. This signal triggers contraction of cardiac muscle through a process termed excitation-contraction (EC) coupling. EC coupling is initiated in dyadic structures of cardiac myocytes, where ryanodine receptors in the junctional sarcoplasmic reticulum come into close apposition with clusters of CaV1.2 channels in invaginations of the sarcolemma. Cooperative activation of CaV1.2 channels within these clusters causes a local increase in intracellular Ca2+ that activates the juxtaposed ryanodine receptors. A salient feature of healthy cardiac function is the reliable and precise beat-to-beat pacemaking and amplitude of Ca2+ transients during EC coupling. In this review, we discuss recent discoveries suggesting that the exquisite reproducibility of this system emerges, paradoxically, from high variability at subcellular, cellular, and network levels. This variability is attributable to stochastic fluctuations in ion channel trafficking, clustering, and gating, as well as dyadic structure, which increase intracellular Ca2+ variance during EC coupling. Although the effects of these large, local fluctuations in function and organization are sometimes negligible at the macroscopic level owing to spatial-temporal summation within and across cells in the tissue, recent work suggests that the "noisiness" of these intracellular Ca2+ events may either enhance or counterintuitively reduce variability in a context-dependent manner. Indeed, these noisy events may represent distinct regulatory features in the tuning of cardiac contractility. Collectively, these observations support the importance of incorporating experimentally determined values of Ca2+ variance in all EC coupling models. The high reproducibility of cardiac contraction is a paradoxical outcome of high Ca2+ signaling variability at subcellular, cellular, and network levels caused by stochastic fluctuations in multiple processes in time and space. This underlying stochasticity, which counterintuitively manifests as reliable, consistent Ca2+ transients during EC coupling, also allows for rapid changes in cardiac rhythmicity and contractility in health and disease.

Preserved cardiac performance and adrenergic response in a rabbit model with decreased ryanodine receptor 2 expression

Ryanodine receptor 2 (RyR2) is an ion channel in the heart responsible for releasing into the cytosol most of the Ca2+ required for contraction. Proper regulation of RyR2 is critical, as highlighted by the association between channel dysfunction and cardiac arrhythmia. Lower RyR2 expression is also observed in some forms of heart disease; however, there is limited information on the impact of this change on excitation-contraction (e-c) coupling, Ca2+-dependent arrhythmias, and cardiac performance. We used a constitutive knock-out of RyR2 in rabbits (RyR2-KO) to assess the extent to which a stable decrease in RyR2 expression modulates Ca2+ handling in the heart. We found that homozygous knock-out of RyR2 in rabbits is embryonic lethal. Remarkably, heterozygotes (KO+/-) show ~50% loss of RyR2 protein without developing an overt phenotype at the intact animal and whole heart levels. Instead, we found that KO+/- myocytes show (1) remodeling of RyR2 clusters, favoring smaller groups in which channels are more densely arranged; (2) lower Ca2+ spark frequency and amplitude; (3) slower rate of Ca2+ release and mild but significant desynchronization of the Ca2+ transient; and (4) a significant decrease in the basal phosphorylation of S2031, likely due to increased association between RyR2 and PP2A. Our data show that RyR2 deficiency, although remarkable at the molecular and subcellular level, has only a modest impact on global Ca2+ release and is fully compensated at the whole-heart level. This highlights the redundancy of RyR2 protein expression and the plasticity of the e-c coupling apparatus.

Nanoscale Organization, Regulation, and Dynamic Reorganization of Cardiac Calcium Channels

The architectural specializations and targeted delivery pathways of cardiomyocytes ensure that L-type Ca²⁺ channels (CaV1.2) are concentrated on the t-tubule sarcolemma within nanometers of their intracellular partners the type 2 ryanodine receptors (RyR2) which cluster on the junctional sarcoplasmic reticulum (jSR). The organization and distribution of these two groups of cardiac calcium channel clusters critically underlies the uniform contraction of the myocardium. Ca²⁺ signaling between these two sets of adjacent clusters produces Ca²⁺ sparks that in health, cannot escalate into Ca²⁺ waves because there is sufficient separation of adjacent clusters so that the release of Ca²⁺ from one RyR2 cluster or supercluster, cannot activate and sustain the release of Ca²⁺ from neighboring clusters. Instead, thousands of these Ca²⁺ release units (CRUs) generate near simultaneous Ca²⁺ sparks across every cardiomyocyte during the action potential when calcium induced calcium release from RyR2 is stimulated by depolarization induced Ca²⁺ influx through voltage dependent CaV1.2 channel clusters. These sparks summate to generate a global Ca²⁺ transient that activates the myofilaments and thus the electrical signal of the action potential is transduced into a functional output, myocardial contraction. To generate more, or less contractile force to match the hemodynamic and metabolic demands of the body, the heart responds to β-adrenergic signaling by altering activity of calcium channels to tune excitation-contraction coupling accordingly. Recent accumulating evidence suggests that this tuning process also involves altered expression, and dynamic reorganization of CaV1.2 and RyR2 channels on their respective membranes to control the amplitude of Ca²⁺ entry, SR Ca²⁺ release and myocardial function. In heart failure and aging, altered distribution and reorganization of these key Ca²⁺ signaling proteins occurs alongside architectural remodeling and is thought to contribute to impaired contractile function. In the present review we discuss these latest developments, their implications, and future questions to be addressed.

Mechanisms and Physiological Implications of Cooperative Gating of Ion Channels Clusters

Ion channels play a central role in the regulation of nearly every cellular process. Dating back to the classic 1952 Hodgkin-Huxley model of the generation of the action potential, ion channels have always been thought of as independent agents. A myriad of recent experimental findings exploiting advances in electrophysiology, structural biology, and imaging techniques, however, have posed a serious challenge to this long-held axiom as several classes of ion channels appear to open and close in a coordinated, cooperative manner. Ion channel cooperativity ranges from variable-sized oligomeric cooperative gating in voltage-gated, dihydropyridine-sensitive Ca_v1.2 and Ca_v1.3 channels to obligatory dimeric assembly and gating of voltage-gated Na_v1.5 channels. Potassium channels, transient receptor potential channels, hyperpolarization cyclic nucleotide-activated channels, ryanodine receptors (RyRs), and inositol trisphosphate receptors (IP₃Rs) have also been shown to gate cooperatively. The implications of cooperative gating of these ion channels range from fine tuning excitation-contraction coupling in muscle cells to regulating cardiac function and vascular tone, to modulation of action potential and conduction velocity in neurons and cardiac cells, and to control of pace-making activity in the heart. In this review, we discuss the mechanisms leading to cooperative gating of ion channels, their physiological consequences and how alterations in cooperative gating of ion channels may induce a range of clinically significant pathologies.

The Physiology and Pathophysiology of T-Tubules in the Heart

In cardiomyocytes, invaginations of the sarcolemmal membrane called t-tubules are critically important for triggering contraction by excitation-contraction (EC) coupling. These structures form functional junctions with the sarcoplasmic reticulum (SR), and thereby enable close contact between L-type Ca²⁺ channels (LTCCs) and Ryanodine Receptors (RyRs). This arrangement in turn ensures efficient triggering of Ca²⁺ release, and contraction. While new data indicate that t-tubules are capable of exhibiting compensatory remodeling, they are also widely reported to be structurally and functionally compromised during disease, resulting in disrupted Ca²⁺ homeostasis, impaired systolic and/or diastolic function, and arrhythmogenesis. This review summarizes these findings, while highlighting an emerging appreciation of the distinct roles of t-tubules in the pathophysiology of heart failure with reduced and preserved ejection fraction (HFrEF and HFpEF). In this context, we review current understanding of the processes underlying t-tubule growth, maintenance, and degradation, underscoring the involvement of a variety of regulatory proteins, including junctophilin-2 (JPH2), amphiphysin-2 (BIN1), caveolin-3 (Cav3), and newer candidate proteins. Upstream regulation of t-tubule structure/function by cardiac workload and specifically ventricular wall stress is also discussed, alongside perspectives for novel strategies which may therapeutically target these mechanisms.

Increasing SERCA function promotes initiation of calcium sparks and breakup of calcium waves

KEY POINTS

Increasing sarcoplasmic/endoplasmic reticulum calcium ATPase (SERCA) pump activity enhances sarcoplasmic reticulum calcium (Ca) load, which increases both ryanodine receptor opening and driving force of Ca release flux. Both of these effects promote Ca spark formation and wave propagation. However, increasing SERCA activity also accelerates local cytosolic Ca decay as the wave front travels to the next cluster, which limits wave propagation. As a result, increasing SERCA pump activity has a biphasic effect on the propensity of arrhythmogenic Ca waves, but a monotonic effect to increase Ca spark frequency and amplitude.

ABSTRACT

Waves of sarcoplasmic reticulum (SR) calcium (Ca) release can cause arrhythmogenic afterdepolarizations in cardiac myocytes. Ca waves propagate when Ca sparks at one Ca release unit (CRU) recruit new Ca sparks in neighbouring CRUs. Under normal conditions, Ca sparks are too small to recruit neighbouring Ca sparks where Ca sensitivity is also low. However, under pathological conditions such as a Ca overload or ryanodine receptor (RyR) sensitization, Ca sparks can be larger and propagate more readily as macro-sparks or full Ca waves. Increasing SERCA pump activity promotes SR Ca load, which promotes RyR opening and increases driving force of the Ca release flux from SR to cytosol, promoting Ca waves. However, high sarcoplasmic/endoplasmic reticulum calcium ATPase (SERCA) activity can also decrease local cytosolic [Ca] as it approaches the next CRU, thereby reducing wave appearance and propagation. In this study, we use a physiologically detailed model of subcellular Ca cycling and experiments in phospholamban-knockout mice, to show how Ca waves are initiated and propagate and how different conditions contribute to the generation and propagation of Ca waves. We show that reducing diffusive coupling between Ca sparks by increasing SERCA activity prevents Ca waves by reducing [Ca] at the next CRU, as do Ca buffers, low intra-SR Ca diffusion and distance between CRUs. Increasing SR Ca uptake rate has a biphasic effect on Ca wave propagation; initially it enhances Ca spark probability and amplitude and CRU coupling, thereby promoting arrhythmogenic Ca wave propagation, but at higher levels SR Ca uptake can abort those arrhythmogenic Ca waves.

NADPH Oxidases (NOX): An Overview from Discovery, Molecular Mechanisms to Physiology and Pathology

The reactive oxygen species (ROS)-producing enzyme NADPH oxidase (NOX) was first identified in the membrane of phagocytic cells. For many years, its only known role was in immune defense, where its ROS production leads to the destruction of pathogens by the immune cells. NOX from phagocytes catalyzes, via one-electron trans-membrane transfer to molecular oxygen, the production of the superoxide anion. Over the years, six human homologs of the catalytic subunit of the phagocyte NADPH oxidase were found: NOX1, NOX3, NOX4, NOX5, DUOX1, and DUOX2. Together with the NOX2/gp91phox component present in the phagocyte NADPH oxidase assembly itself, the homologs are now referred to as the NOX family of NADPH oxidases. NOX are complex multidomain proteins with varying requirements for assembly with combinations of other proteins for activity. The recent structural insights acquired on both prokaryotic and eukaryotic NOX open new perspectives for the understanding of the molecular mechanisms inherent to NOX regulation and ROS production (superoxide or hydrogen peroxide). This new structural information will certainly inform new investigations of human disease. As specialized ROS producers, NOX enzymes participate in numerous crucial physiological processes, including host defense, the post-translational processing of proteins, cellular signaling, regulation of gene expression, and cell differentiation. These diversities of physiological context will be discussed in this review. We also discuss NOX misregulation, which can contribute to a wide range of severe pathologies, such as atherosclerosis, hypertension, diabetic nephropathy, lung fibrosis, cancer, or neurodegenerative diseases, giving this family of membrane proteins a strong therapeutic interest.

Human Atrial Fibrillation Is Not Associated With Remodeling of Ryanodine Receptor Clusters

The release of Ca²⁺ by ryanodine receptor (RyR2) channels is critical for cardiac function. However, abnormal RyR2 activity has been linked to the development of arrhythmias, including increased spontaneous Ca²⁺ release in human atrial fibrillation (AF). Clustering properties of RyR2 have been suggested to alter the activity of the channel, with remodeling of RyR2 clusters identified in pre-clinical models of AF and heart failure. Whether such remodeling occurs in human cardiac disease remains unclear. This study aimed to investigate the nanoscale organization of RyR2 clusters in AF patients – the first known study to examine this potential remodeling in diseased human cardiomyocytes. Right atrial appendage from cardiac surgery patients with paroxysmal or persistent AF, or without AF (non-AF) were examined using super-resolution (dSTORM) imaging. Significant atrial dilation and cardiomyocyte hypertrophy was observed in persistent AF patients compared to non-AF, with these two parameters significantly correlated. Interestingly, the clustering properties of RyR2 were remarkably unaltered in the AF patients. No significant differences were identified in cluster size (mean ∼18 RyR2 channels), density or channel packing within clusters between patient groups. The spatial organization of clusters throughout the cardiomyocyte was also unchanged across the groups. RyR2 clustering properties did not significantly correlate with patient characteristics. In this first study to examine nanoscale RyR2 organization in human cardiac disease, these findings indicate that RyR2 cluster remodeling is not an underlying mechanism contributing to altered channel function and subsequent arrhythmogenesis in human AF.

Simulating cardiac Ca2+ release units: effects of RyR cluster size and Ca2+ buffers on diastolic Ca2+ leak

Leak of Ca²⁺ out of the cardiac sarcoplasmic reticulum (SR) via ryanodine receptors (RyRs) during diastole is vital to regulate SR Ca²⁺ levels. This leak can become deleterious when large spontaneous RyR-mediated Ca²⁺ release events evoke proarrhythmic Ca²⁺ waves that can lead to delayed after-depolarizations. Here, we model diastolic SR Ca²⁺ leak at individual SR Ca²⁺ release sites using computer simulations of RyR arrays like those in the dyadic cleft. The results show that RyR arrays size has a significant effect on SR Ca²⁺ leak, with bigger arrays producing larger and more frequent Ca²⁺ release events. Moreover, big RyR arrays are more susceptible to small changes in the levels of dyadic Ca²⁺ buffers. Such changes in buffering shift Ca²⁺ leak from small Ca²⁺ release events (involving few open RyRs) to larger events (with many open RyRs). Moreover, by analyzing a large parameter space of possible buffering and SR Ca²⁺ loads, we find further evidence for the hypothesis that SR Ca²⁺ leak by RyR arrays can undergo a sudden phase transition.

Dual regulation by subcellular calcium heterogeneity and heart rate variability on cardiac electromechanical dynamics

Heart rate constantly varies under physiological conditions, termed heart rate variability (HRV), and in clinical studies, low HRV is associated with a greater risk of cardiac arrhythmias. Prior work has shown that HRV influences the temporal patterns of electrical activity, specifically the formation of pro-arrhythmic alternans, a beat-to-beat alternation in the action potential duration (APD), or intracellular calcium (Ca) levels. We previously showed that HRV may be anti-arrhythmic by disrupting APD and Ca alternations in a homogeneous cardiac myocyte. Here, we expand on our previous work, incorporating variation in subcellular Ca handling (also known to influence alternans) into a nonlinear map model of a cardiac myocyte composed of diffusively coupled Ca release units (CRUs). Ca-related parameters and initial conditions of each CRU are varied to mimic subcellular Ca heterogeneity, and a stochastic pacing sequence reproduces HRV. We find that subcellular Ca heterogeneity promotes the formation of spatially discordant subcellular alternans patterns, which decreases whole cell Ca and APD alternation for low and moderate HRV, while high subcellular Ca heterogeneity and HRV both promote electromechanical desynchronization. Finally, we find that for low and moderate HRV, both the specific subcellular Ca-related parameters and the pacing sequences influence measures of electromechanical dynamics, while for high HRV, these measures depend predominantly on the pacing sequence. Our results suggest that pro-arrhythmic subcellular discordant alternans tend to form for low levels of HRV, while high HRV may be anti-arrhythmic due to mitigated influence from subcellular Ca heterogeneity and desynchronization of APD from Ca instabilities.

Sacubitrilat reduces pro-arrhythmogenic sarcoplasmic reticulum Ca2+ leak in human ventricular cardiomyocytes of patients with end-stage heart failure

Aims

Inhibition of neprilysin and angiotensin II receptor by sacubitril/valsartan (Val) (LCZ696) reduces mortality in heart failure (HF) patients compared with sole inhibition of renin–angiotensin system. Beneficial effects of increased natriuretic peptide levels upon neprilysin inhibition have been proposed, whereas direct effects of sacubitrilat (Sac) (LBQ657) on myocardial Ca²⁺ cycling remain elusive.

Methods and results

Confocal microscopy (Fluo‐4 AM) was used to investigate pro‐arrhythmogenic sarcoplasmic reticulum (SR) Ca²⁺ leak in freshly isolated murine and human ventricular cardiomyocytes (CMs) upon Sac (40 μmol/L)/Val (13 μmol/L) treatment. The concentrations of Sac and Val equalled plasma concentrations of LCZ696 treatment used in PARADIGM‐HF trial. Epifluorescence microscopy measurements (Fura‐2 AM) were performed to investigate effects on systolic Ca²⁺ release, SR Ca²⁺ load, and Ca²⁺‐transient kinetics in freshly isolated murine ventricular CMs. The impact of Sac on myocardial contractility was evaluated using in toto‐isolated, isometrically twitching ventricular trabeculae from human hearts with end‐stage HF. Under basal conditions, the combination of Sac/Val did not influence diastolic Ca²⁺‐spark frequency (CaSpF) nor pro‐arrhythmogenic SR Ca² leak in isolated murine ventricular CMs (n CMs/hearts = 80/7 vs. 100/7, P = 0.91/0.99). In contrast, Sac/Val treatment reduced CaSpF by 35 ± 9% and SR Ca²⁺ leak by 45 ± 9% in CMs put under catecholaminergic stress (isoproterenol 30 nmol/L, n = 81/7 vs. 62/7, P < 0.001 each). This could be attributed to Sac, as sole Sac treatment also reduced both parameters by similar degrees (reduction of CaSpF by 57 ± 7% and SR Ca²⁺ leak by 76 ± 5%; n = 101/4 vs. 108/4, P < 0.01 each), whereas sole Val treatment did not. Systolic Ca²⁺ release, SR Ca²⁺ load, and Ca²⁺‐transient kinetics including SERCA activity (k_SERCA) were not compromised by Sac in isolated murine CMs (n = 41/6 vs. 39/6). Importantly, the combination of Sac/Val and Sac alone also reduced diastolic CaSpF and SR Ca²⁺ leak (reduction by 74 ± 7%) in human left ventricular CMs from patients with end‐stage HF (n = 71/8 vs. 78/8, P < 0.05 each). Myocardial contractility of human ventricular trabeculae was not acutely affected by Sac treatment as the developed force remained unchanged over a time course of 30 min (n trabeculae/hearts = 3/3 vs. 4/3).

Conclusion

This study demonstrates that neprilysin inhibitor Sac directly improves Ca²⁺ homeostasis in human end‐stage HF by reducing pro‐arrhythmogenic SR Ca²⁺ leak without acutely affecting systolic Ca²⁺ release and inotropy. These effects might contribute to the mortality benefits observed in the PARADIGM‐HF trial.

Cardiolipin Structure and Oxidation Are Affected by Ca2+ at the Interface of Lipid Bilayers

Ca²⁺-overload contributes to the oxidation of mitochondrial membrane lipids and associated events such as the permeability transition pore (MPTP) opening. Numerous experimental studies about the Ca²⁺/cardiolipin (CL) interaction are reported in the literature, but there are few studies in conjunction with theoretical approaches based on ab initio calculations. In the present study, the lipid fraction of the inner mitochondrial membrane was modeled as POPC/CL large unilamellar vesicles (LUVs). POPC/CL and, comparatively, POPC, and CL LUVs were challenged by singlet molecular oxygen using the anionic porphyrin TPPS4 as a photosensitizer and by free radicals produced by Fe²⁺-citrate. Calcium ion favored both types of lipid oxidation in a lipid composition-dependent manner. In membranes containing predominantly or exclusively POPC, Ca²⁺ increased the oxidation at later reaction times while the oxidation of CL membranes was exacerbated at the early times of reaction. Considering that Ca²⁺ interaction affects the lipid structure and packing, density functional theory (DFT) calculations were applied to the Ca²⁺ association with totally and partially protonated and deprotonated CL, in the presence of water. The interaction of totally and partially protonated CL head groups with Ca²⁺ decreased the intramolecular P-P distance and increased the hydrophobic volume of the acyl chains. Consistently with the theoretically predicted effect of Ca²⁺ on CL, in the absence of pro-oxidants, giant unilamellar vesicles (GUVs) challenged by Ca²⁺ formed buds and many internal vesicles. Therefore, Ca²⁺ induces changes in CL packing and increases the susceptibility of CL to the oxidation promoted by free radicals and excited species.

Cardiac ryanodine receptor distribution is dynamic and changed by auxiliary proteins and post-translational modification

The effects of the immunophilins, FKBP12 and FKBP12.6, and phosphorylation on type II ryanodine receptor (RyR2) arrangement and function were examined using correlation microscopy (line scan confocal imaging of Ca²⁺ sparks and dual-tilt electron tomography) and dSTORM imaging of permeabilized Wistar rat ventricular myocytes. Saturating concentrations (10 µmol/L) of either FKBP12 or 12.6 significantly reduced the frequency, spread, amplitude and Ca²⁺ spark mass relative to control, while the tomograms revealed both proteins shifted the tetramers into a largely side-by-side configuration. Phosphorylation of immunophilin-saturated RyR2 resulted in structural and functional changes largely comparable to phosphorylation alone. dSTORM images of myocyte surfaces demonstrated that both FKBP12 and 12.6 significantly reduced RyR2 cluster sizes, while phosphorylation, even of immunophilin-saturated RyR2, increased them. We conclude that both RyR2 cluster size and the arrangement of tetramers within clusters is dynamic and respond to changes in the cellular environment. Further, these changes affect Ca²⁺ spark formation.

Disruption of Ca2+i Homeostasis and Connexin 43 Hemichannel Function in the Right Ventricle Precedes Overt Arrhythmogenic Cardiomyopathy in Plakophilin-2-Deficient Mice

BACKGROUND

Plakophilin-2 (PKP2) is classically defined as a desmosomal protein. Mutations in PKP2 associate with most cases of gene-positive arrhythmogenic right ventricular cardiomyopathy. A better understanding of PKP2 cardiac biology can help elucidate the mechanisms underlying arrhythmic and cardiomyopathic events consequent to PKP2 deficiency. Here, we sought to capture early molecular/cellular events that can act as nascent arrhythmic/cardiomyopathic substrates.

METHODS

We used multiple imaging, biochemical and high-resolution mass spectrometry methods to study functional/structural properties of cells/tissues derived from cardiomyocyte-specific, tamoxifen-activated, PKP2 knockout mice (PKP2cKO) 14 days post-tamoxifen injection, a time point preceding overt electrical or structural phenotypes. Myocytes from right or left ventricular free wall were studied separately.

RESULTS

Most properties of PKP2cKO left ventricular myocytes were not different from control; in contrast, PKP2cKO right ventricular (RV) myocytes showed increased amplitude and duration of Ca²⁺ transients, increased Ca²⁺ in the cytoplasm and sarcoplasmic reticulum, increased frequency of spontaneous Ca²⁺ release events (sparks) even at comparable sarcoplasmic reticulum load, and dynamic Ca²⁺ accumulation in mitochondria. We also observed early- and delayed-after transients in RV myocytes and heightened susceptibility to arrhythmias in Langendorff-perfused hearts. In addition, ryanodine receptor 2 in PKP2cKO-RV cells presented enhanced Ca²⁺ sensitivity and preferential phosphorylation in a domain known to modulate Ca²⁺ gating. RNAseq at 14 days post-tamoxifen showed no relevant difference in transcript abundance between RV and left ventricle, neither in control nor in PKP2cKO cells. Instead, we found an RV-predominant increase in membrane permeability that can permit Ca²⁺ entry into the cell. Connexin 43 ablation mitigated the membrane permeability increase, accumulation of cytoplasmic Ca²⁺, increased frequency of sparks and early stages of RV dysfunction. Connexin 43 hemichannel block with GAP19 normalized [Ca²⁺]_i homeostasis. Similarly, protein kinase C inhibition normalized spark frequency at comparable sarcoplasmic reticulum load levels.

CONCLUSIONS

Loss of PKP2 creates an RV-predominant arrhythmogenic substrate (Ca²⁺ dysregulation) that precedes the cardiomyopathy; this is, at least in part, mediated by a Connexin 43-dependent membrane conduit and repressed by protein kinase C inhibitors. Given that asymmetric Ca²⁺ dysregulation precedes the cardiomyopathic stage, we speculate that abnormal Ca²⁺ handling in RV myocytes can be a trigger for gross structural changes observed at a later stage.