Ryanodine receptor cluster fragmentation and redistribution in persistent atrial fibrillation enhance calcium release

Super-resolution microscopy as a drug discovery tool

At the turn of the century a fundamental resolution barrier in fluorescence microscopy known as the diffraction limit was broken, giving rise to the field of super-resolution microscopy. Subsequent nanoscopic investigation with visible light revolutionised our understanding of how previously unknown molecular features give rise to the emergent behaviour of cells. It transpires that the devil is in these fine molecular details, and essential nanoscale processes were found everywhere researchers chose to look. Now, after nearly two decades, super-resolution microscopy has begun to address previously unmet challenges in the study of human disease and is poised to become a pivotal tool in drug discovery.

Restoring Atrial T-Tubules Augments Systolic Ca Upon Recovery From Heart Failure

BACKGROUND

Transverse (t)-tubules drive the rapid and synchronous Ca2+ rise in cardiac myocytes. The virtual complete atrial t-tubule loss in heart failure (HF) decreases Ca2+ release. It is unknown if or how atrial t-tubules can be restored and how this affects systolic Ca2+.

METHODS

HF was induced in sheep by rapid ventricular pacing and recovered following termination of rapid pacing. Serial block-face scanning electron microscopy and confocal imaging were used to study t-tubule ultrastructure. Function was assessed using patch clamp, Ca2+, and confocal imaging. Candidate proteins involved in atrial t-tubule recovery were identified by western blot and expressed in rat neonatal ventricular myocytes to determine if they altered t-tubule structure.

RESULTS

Atrial t-tubules were lost in HF but reappeared following recovery from HF. Recovered t-tubules were disordered, adopting distinct morphologies with increased t-tubule length and branching. T-tubule disorder was associated with mitochondrial disorder. Recovered t-tubules were functional, triggering Ca2+ release in the cell interior. Systolic Ca2+, ICa-L, sarcoplasmic reticulum Ca2+ content, and sarcoendoplasmic reticulum Ca2+ ATPase function were restored following recovery from HF. Confocal microscopy showed fragmentation of ryanodine receptor staining and movement away from the z-line in HF, which was reversed following recovery from HF. Acute detubulation, to remove recovered t-tubules, confirmed their key role in restoration of the systolic Ca2+ transient, the rate of Ca2+ removal, and the peak L-type Ca2+ current. The abundance of telethonin and myotubularin decreased during HF and increased during recovery. Transfection with these proteins altered the density and structure of tubules in neonatal myocytes. Myotubularin had a greater effect, increasing tubule length and branching, replicating that seen in the recovery atria.

CONCLUSIONS

We show that recovery from HF restores atrial t-tubules, and this promotes recovery of ICa-L, sarcoplasmic reticulum Ca2+ content, and systolic Ca2+. We demonstrate an important role for myotubularin in t-tubule restoration. Our findings reveal a new and viable therapeutic strategy.

Impaired Intracellular Calcium Buffering Contributes to the Arrhythmogenic Substrate in Atrial Myocytes From Patients With Atrial Fibrillation

BACKGROUND

Alterations in the buffering of intracellular Ca2+, for which myofilament proteins play a key role, have been shown to promote cardiac arrhythmia. It is interesting that although studies report atrial myofibrillar degradation in patients with persistent atrial fibrillation (persAF), the intracellular Ca2+ buffering profile in persAF remains obscure. Therefore, we aimed to investigate the intracellular buffering of Ca2+ and its potential arrhythmogenic role in persAF.

METHODS

Transmembrane Ca2+ fluxes (patch-clamp) and intracellular Ca2+ signaling (fluo-3-acetoxymethyl ester) were recorded simultaneously in myocytes from right atrial biopsies of sinus rhythm (Ctrl) and patients with persAF, alongside human atrial subtype induced pluripotent stem cell-derived cardiac myocytes (iPSC-CMs). Protein levels were quantified by immunoblotting of human atrial tissue and induced pluripotent stem cell-derived cardiac myocytes. Mouse whole heart and atrial electrophysiology were measured on a Langendorff system.

RESULTS

Cytosolic Ca2+ buffering was decreased in atrial myocytes of patients with persAF because of a depleted amount of Ca2+ buffers. In agreement, protein levels of selected Ca2+ binding myofilament proteins, including cTnC (cardiac troponin C), a major cytosolic Ca2+ buffer, were significantly lower in patients with persAF. Small interfering RNA (siRNA)-mediated knockdown of cTnC (si-cTNC) in atrial iPSC-CM phenocopied the reduced cytosolic Ca2+ buffering observed in persAF. Si-cTnC treated atrial iPSC-CM exhibited a higher predisposition to spontaneous Ca2+ release events and developed action potential alternans at low stimulation frequencies. Last, indirect reduction of cytosolic Ca2+ buffering using blebbistatin in an ex vivo mouse whole heart model increased vulnerability to tachypacing-induced atrial arrhythmia, validating the direct mechanistic link between impaired cytosolic Ca2+ buffering and atrial arrhythmogenesis.

CONCLUSIONS

Our findings suggest that loss of myofilament proteins, particularly reduced cTnC protein levels, causes diminished cytosolic Ca2+ buffering in persAF, thereby potentiating the occurrence of spontaneous Ca2+ release events and atrial fibrillation susceptibility. Strategies targeting intracellular buffering may represent a promising therapeutic lead in persAF management.

BIN1 knockdown rescues systolic dysfunction in aging male mouse hearts

Cardiac dysfunction is a hallmark of aging in humans and mice. Here we report that a two-week treatment to restore youthful Bridging Integrator 1 (BIN1) levels in the hearts of 24-month-old mice rejuvenates cardiac function and substantially reverses the aging phenotype. Our data indicate that age-associated overexpression of BIN1 occurs alongside dysregulated endosomal recycling and disrupted trafficking of cardiac CaV1.2 and type 2 ryanodine receptors. These deficiencies affect channel function at rest and their upregulation during acute stress. In vivo echocardiography reveals reduced systolic function in old mice. BIN1 knockdown using an adeno-associated virus serotype 9 packaged shRNA-mBIN1 restores the nanoscale distribution and clustering plasticity of ryanodine receptors and recovers Ca2+ transient amplitudes and cardiac systolic function toward youthful levels. Enhanced systolic function correlates with increased phosphorylation of the myofilament protein cardiac myosin binding protein-C. These results reveal BIN1 knockdown as a novel therapeutic strategy to rejuvenate the aging myocardium.

Cardiac multiscale bioimaging: from nano- through micro- to mesoscales

Cardiac multiscale bioimaging is an emerging field that aims to provide a comprehensive understanding of the heart and its functions at various levels, from the molecular to the entire organ. It combines both physiologically and clinically relevant dimensions: from nano- and micrometer resolution imaging based on vibrational spectroscopy and high-resolution microscopy to assess molecular processes in cardiac cells and myocardial tissue, to mesoscale structural investigations to improve the understanding of cardiac (patho)physiology. Tailored super-resolution deep microscopy with advanced proteomic methods and hands-on experience are thus strategically combined to improve the quality of cardiovascular research and support future medical decision-making by gaining additional biomolecular information for translational and diagnostic applications.

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.

Cardiomyocyte ryanodine receptor clusters expand and coalesce after application of isoproterenol

Earlier work has shown that ventricular ryanodine receptors (RyR2) within a cluster rearrange on phosphorylation as well as with a number of other stimuli. Using dSTORM, we investigated the effects of 300 nmol/liter isoproterenol on RyR2 clusters. In rat ventricular cardiomyocytes, there was a symmetrical enlargement of RyR2 cluster areas, a decrease in the edge-to-edge nearest neighbor distance, and distribution changes that suggested movement to increase the cluster areas by coalescence. The surface area covered by the phosphorylated clusters was significantly greater than in the control cells, as was the cluster density. This latter change was accompanied by a decreased cluster fragmentation, implying that new tetramers were brought into the sarcoplasmic reticulum. We propose a possible mechanism to explain these changes. We also visualized individual RyR2 tetramers and confirmed our earlier electron-tomographic finding that the tetramers are in a disorganized but non-random array occupying about half of the cluster area. Multiclusters, cluster groups defined by the maximum distance between their members, were analyzed for various distances. At 100 nm, the areas occupied by the multiclusters just exceeded those of the single clusters, and more than half of the multiclusters had only a single subcluster that could initiate a spark. Phosphorylation increased the size of the multiclusters, markedly so for distances >100 nm. There was no relationship between the number of subclusters in a group and the area covered by it. We conclude that isoproterenol induces rapid, significant, changes in the molecular architecture of excitation-contraction coupling.

Physiology of intracellular calcium buffering

Calcium signaling underlies much of physiology. Almost all the Ca2+ in the cytoplasm is bound to buffers, with typically only ∼1% being freely ionized at resting levels in most cells. Physiological Ca2+ buffers include small molecules and proteins, and experimentally Ca2+ indicators will also buffer calcium. The chemistry of interactions between Ca2+ and buffers determines the extent and speed of Ca2+ binding. The physiological effects of Ca2+ buffers are determined by the kinetics with which they bind Ca2+ and their mobility within the cell. The degree of buffering depends on factors such as the affinity for Ca2+, the Ca2+ concentration, and whether Ca2+ ions bind cooperatively. Buffering affects both the amplitude and time course of cytoplasmic Ca2+ signals as well as changes of Ca2+ concentration in organelles. It can also facilitate Ca2+ diffusion inside the cell. Ca2+ buffering affects synaptic transmission, muscle contraction, Ca2+ transport across epithelia, and the killing of bacteria. Saturation of buffers leads to synaptic facilitation and tetanic contraction in skeletal muscle and may play a role in inotropy in the heart. This review focuses on the link between buffer chemistry and function and how Ca2+ buffering affects normal physiology and the consequences of changes in disease. As well as summarizing what is known, we point out the many areas where further work is required.

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.

Mechanisms of spontaneous Ca2+ release-mediated arrhythmia in a novel 3D human atrial myocyte model: II. Ca2+ -handling protein variation

Disruption of the transverse-axial tubule system (TATS) in diseases such as heart failure and atrial fibrillation occurs in combination with changes in the expression and distribution of key Ca2+ -handling proteins. Together this ultrastructural and ionic remodelling is associated with aberrant Ca2+ cycling and electrophysiological instabilities that underlie arrhythmic activity. However, due to the concurrent changes in TATs and Ca2+ -handling protein expression and localization that occur in disease it is difficult to distinguish their individual contributions to the arrhythmogenic state. To investigate this, we applied our novel 3D human atrial myocyte model with spatially detailed Ca2+ diffusion and TATS to investigate the isolated and interactive effects of changes in expression and localization of key Ca2+ -handling proteins and variable TATS density on Ca2+ -handling abnormality driven membrane instabilities. We show that modulating the expression and distribution of the sodium-calcium exchanger, ryanodine receptors and the sarcoplasmic reticulum (SR) Ca2+ buffer calsequestrin have varying pro- and anti-arrhythmic effects depending on the balance of opposing influences on SR Ca2+ leak-load and Ca2+ -voltage relationships. Interestingly, the impact of protein remodelling on Ca2+ -driven proarrhythmic behaviour varied dramatically depending on TATS density, with intermediately tubulated cells being more severely affected compared to detubulated and densely tubulated myocytes. This work provides novel mechanistic insight into the distinct and interactive consequences of TATS and Ca2+ -handling protein remodelling that underlies dysfunctional Ca2+ cycling and electrophysiological instability in disease. KEY POINTS: In our companion paper we developed a 3D human atrial myocyte model, coupling electrophysiology and Ca2+ handling with subcellular spatial details governed by the transverse-axial tubule system (TATS). Here we utilize this model to mechanistically examine the impact of TATS loss and changes in the expression and distribution of key Ca2+ -handling proteins known to be remodelled in disease on Ca2+ homeostasis and electrophysiological stability. We demonstrate that varying the expression and localization of these proteins has variable pro- and anti-arrhythmic effects with outcomes displaying dependence on TATS density. Whereas detubulated myocytes typically appear unaffected and densely tubulated cells seem protected, the arrhythmogenic effects of Ca2+ handling protein remodelling are profound in intermediately tubulated cells. Our work shows the interaction between TATS and Ca2+ -handling protein remodelling that underlies the Ca2+ -driven proarrhythmic behaviour observed in atrial fibrillation and may help to predict the effects of antiarrhythmic strategies at varying stages of ultrastructural remodelling.

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.

Correlative super-resolution analysis of cardiac calcium sparks and their molecular origins in health and disease

Rapid release of calcium from internal stores via ryanodine receptors (RyRs) is one of the fastest types of cytoplasmic second messenger signalling in excitable cells. In the heart, rapid summation of the elementary events of calcium release, 'calcium sparks', determine the contraction of the myocardium. We adapted a correlative super-resolution microscopy protocol to correlate sub-plasmalemmal spontaneous calcium sparks in rat right ventricular myocytes with the local nanoscale RyR2 positions. This revealed a steep relationship between the integral of a calcium spark and the sum of the local RyR2s. Segmentation of recurring spark sites showed evidence of repeated and triggered saltatory activation of multiple local RyR2 clusters. In myocytes taken from failing right ventricles, RyR2 clusters themselves showed a dissipated morphology and fragmented (smaller) clusters. They also featured greater heterogeneity in both the spark properties and the relationship between the integral of the calcium spark and the local ensemble of RyR2s. While fragmented (smaller) RyR2 clusters were rarely observed directly underlying the larger sparks or the recurring spark sites, local interrogation of the channel-to-channel distances confirmed a clear link between the positions of each calcium spark and the tight, non-random clustering of the local RyR2 in both healthy and failing ventricles.

InsP3R-RyR channel crosstalk augments sarcoplasmic reticulum Ca2+ release and arrhythmogenic activity in post-MI pig cardiomyocytes

Ca²⁺ transients (CaT) underlying cardiomyocyte (CM) contraction require efficient Ca²⁺ coupling between sarcolemmal Ca²⁺ channels and sarcoplasmic reticulum (SR) ryanodine receptor Ca²⁺ channels (RyR) for their generation; reduced coupling in disease contributes to diminished CaT and arrhythmogenic Ca²⁺ events. SR Ca²⁺ release also occurs via inositol 1,4,5-trisphosphate receptors (InsP₃R) in CM. While this pathway contributes negligeably to Ca²⁺ handling in healthy CM, rodent studies support a role in altered Ca²⁺ dynamics and arrhythmogenic Ca²⁺ release involving InsP₃R crosstalk with RyRs in disease. Whether this mechanism persists in larger mammals with lower T-tubular density and coupling of RyRs is not fully resolved. We have recently shown an arrhythmogenic action of InsP₃-induced Ca²⁺ release (IICR) in end stage human heart failure, often associated with underlying ischemic heart disease (IHD). How IICR contributes to early stages of disease is however not determined but highly relevant. To access this stage, we chose a porcine model of IHD, which shows substantial remodelling of the area adjacent to the infarct. In cells from this region, IICR preferentially augmented Ca²⁺ release from non-coupled RyR clusters that otherwise showed delayed activation during the CaT. IICR in turn synchronised Ca²⁺ release during the CaT but also induced arrhythmogenic delayed afterdepolarizations and action potentials. Nanoscale imaging identified co-clustering of InsP₃Rs and RyRs, thereby allowing Ca²⁺-mediated channel crosstalk. Mathematical modelling supported and further delineated this mechanism of enhanced InsP₃R-RyRs coupling in MI. Our findings highlight the role of InsP₃R-RyR channel crosstalk in Ca²⁺ release and arrhythmia during post-MI remodelling.

Live-cell photo-activated localization microscopy correlates nanoscale ryanodine receptor configuration to calcium sparks in cardiomyocytes

Ca2+ sparks constitute the fundamental units of Ca2+ release in cardiomyocytes. Here we investigate how ryanodine receptors (RyRs) collectively generate these events by employing a transgenic mouse with a photo-activated label on RyR2. This allowed correlative imaging of RyR localization, by super-resolution Photo-Activated Localization Microscopy, and Ca2+ sparks, by high-speed imaging. Two populations of Ca2+ sparks were observed: stationary events and "travelling" events that spread between neighbouring RyR clusters. Travelling sparks exhibited up to 8 distinct releases, sourced from local or distal junctional sarcoplasmic reticulum. Quantitative analyses showed that sparks may be triggered by any number of RyRs within a cluster, and that acute β-adrenergic stimulation augments intra-cluster RyR recruitment to generate larger events. In contrast, RyR "dispersion" during heart failure facilitates the generation of travelling sparks. Thus, RyRs cooperatively generate Ca2+ sparks in a complex, malleable fashion, and channel organization regulates the propensity for local propagation of Ca2+ release.

Regional Differences in Ca2+ Signaling and Transverse-Tubules across Left Atrium from Adult Sheep

Cardiac excitation-contraction coupling can be different between regions of the heart. Little is known at the atria level, specifically in different regions of the left atrium. This is important given the role of cardiac myocytes from the pulmonary vein sleeves, which are responsible for ectopic activity during atrial fibrillation. In this study, we present a new method to isolate atrial cardiac myocytes from four different regions of the left atrium of a large animal model, sheep, highly relevant to humans. Using collagenase/protease we obtained calcium-tolerant atrial cardiac myocytes from the epicardium, endocardium, free wall and pulmonary vein regions. Calcium transients were slower (time to peak and time to decay) in free wall and pulmonary vein myocytes compared to the epicardium and endocardium. This is associated with lower t-tubule density. Overall, these results suggest regional differences in calcium transient and t-tubule density across left atria, which may play a major role in the genesis of atrial fibrillation.

IP3R activity increases propensity of RyR-mediated sparks by elevating dyadic [Ca2＋]

Calcium (Ca^2＋) plays a critical role in the excitation contraction coupling (ECC) process that mediates the contraction of cardiomyocytes during each heartbeat. While ryanodine receptors (RyRs) are the primary Ca^2＋ channels responsible for generating the cell-wide Ca^2＋ transients during ECC, Ca^2＋ release, via inositol 1,4,5-trisphosphate (IP₃) receptors (IP₃Rs) are also reported in cardiomyocytes to elicit ECC-modulating effects. Recent studies suggest that the localization of IP₃Rs at dyads grant their ability to modify the occurrence of Ca^2＋ sparks (elementary Ca^2＋ release events that constitute cell wide Ca^2＋ releases associated with ECC) which may underlie their modulatory influence on ECC. Here, we aim to uncover the mechanism by which dyad-localized IP₃Rs influence Ca^2＋ spark dynamics. To this end, we developed a mathematical model of the dyad that incorporates the behaviour of IP₃Rs, in addition to RyRs, to reveal the impact of their activity on local Ca^2＋ handling and consequent Ca^2＋ spark occurrence and its properties. Consistent with published experimental data, our model predicts that the propensity for Ca^2＋ spark formation increases in the presence of IP₃R activity. Our simulations support the hypothesis that IP₃Rs elevate Ca^2＋ in the dyad, sensitizing proximal RyRs towards activation and hence Ca^2＋ spark formation. The stochasticity of IP₃R gating is an important aspect of this mechanism. However, dyadic IP₃R activity lowers the Ca^2＋ available in the junctional sarcoplasmic reticulum (JSR) for release, thus resulting in Ca^2＋ sparks with similar durations but lower amplitudes.

Understanding Calmodulin Variants Affecting Calcium-Dependent Inactivation of L-Type Calcium Channels through Whole-Cell Simulation of the Cardiac Ventricular Myocyte

Mutations in the calcium-sensing protein calmodulin (CaM) have been linked to two cardiac arrhythmia diseases, Long QT Syndrome 14 (LQT14) and Catecholaminergic Polymorphic Ventricular Tachycardia Type 4 (CPVT4), with varying degrees of severity. Functional characterization of the CaM mutants most strongly associated with LQT14 show a clear disruption of the calcium-dependent inactivation (CDI) of the L-Type calcium channel (LCC). CPVT4 mutants on the other hand are associated with changes in their affinity to the ryanodine receptor. In clinical studies, some variants have been associated with both CPVT4 and LQT15. This study uses simulations in a model for excitation-contraction coupling in the rat ventricular myocytes to understand how LQT14 variant might give the functional phenotype similar to CPVT4. Changing the CaM-dependent transition rate by a factor of 0.75 corresponding to the D96V variant and by a factor of 0.90 corresponding to the F142L or N98S variants, in a physiologically based stochastic model of the LCC prolonger, the action potential duration changed by a small amount in a cardiac myocyte but did not disrupt CICR at 1, 2, and 4 Hz. Under beta-adrenergic simulation abnormal excitation-contraction coupling was observed above 2 Hz pacing for the mutant CaM. The same conditions applied under beta-adrenergic stimulation led to the rapid onset of arrhythmia in the mutant CaM simulations. Simulations with the LQT14 mutations under the conditions of rapid pacing with beta-adrenergic stimulation drives the cardiac myocyte toward an arrhythmic state known as Ca2+ overload. These simulations provide a mechanistic link to a disease state for LQT14-associated mutations in CaM to yield a CPVT4 phenotype. The results show that small changes to the CaM-regulated inactivation of LCC promote arrhythmia and underscore the significance of CDI in proper heart function.

Direct proteomic and high-resolution microscopy biopsy analysis identifies distinct ventricular fates in severe aortic stenosis

The incidence of aortic valve stenosis (AS), the most common reason for aortic valve replacement (AVR), increases with population ageing. While untreated AS is associated with high mortality, different hemodynamic subtypes range from normal left-ventricular function to severe heart failure. However, the molecular nature underlying four different AS subclasses, suggesting vastly different myocardial fates, is unknown. Here, we used direct proteomic analysis of small left-ventricular biopsies to identify unique protein expression profiles and subtype-specific AS mechanisms. Left-ventricular endomyocardial biopsies were harvested from patients during transcatheter AVR, and inclusion criteria were based on echocardiographic diagnosis of severe AS and guideline-defined AS-subtype classification: 1) normal ejection fraction (EF)/high-gradient; 2) low EF/high-gradient; 3) low EF/low-gradient; and 4) paradoxical low-flow/low-gradient AS. Samples from non-failing donor hearts served as control. We analyzed 25 individual left-ventricular biopsies by data-independent acquisition mass spectrometry (DIA-MS), and 26 biopsies by histomorphology and cardiomyocytes by STimulated Emission Depletion (STED) superresolution microscopy. Notably, DIA-MS reliably detected 2273 proteins throughout each individual left-ventricular biopsy, of which 160 proteins showed significant abundance changes between AS-subtype and non-failing samples including the cardiac ryanodine receptor (RyR2). Hierarchical clustering segregated unique proteotypes that identified three hemodynamic AS-subtypes. Additionally, distinct proteotypes were linked with AS-subtype specific differences in cardiomyocyte hypertrophy. Furthermore, superresolution microscopy of immunolabeled biopsy sections showed subcellular RyR2-cluster fragmentation and disruption of the functionally important association with transverse tubules, which occurred specifically in patients with systolic dysfunction and may hence contribute to depressed left-ventricular function in AS.

Three-dimensional visualization of the cardiac ryanodine receptor clusters and the molecular-scale fraying of dyads

Clusters of ryanodine receptor calcium channels (RyRs) form the primary molecular machinery of intracellular calcium signalling in cardiomyocytes. While a range of optical super-resolution microscopy techniques have revealed the nanoscale structure of these clusters, the three-dimensional (3D) nanoscale topologies of the clusters have remained mostly unresolved. In this paper, we demonstrate the exploitation of molecular-scale resolution in enhanced expansion microscopy (EExM) along with various 2D and 3D visualization strategies to observe the topological complexities, geometries and molecular sub-domains within the RyR clusters. Notably, we observed sub-domains containing RyR-binding protein junctophilin-2 (JPH2) occupying the central regions of RyR clusters in the deeper interior of the myocytes (including dyads), while the poles were typically devoid of JPH2, lending to a looser RyR arrangement. By contrast, peripheral RyR clusters exhibited variable co-clustering patterns and ratios between RyR and JPH2. EExM images of dyadic RyR clusters in right ventricular (RV) myocytes isolated from rats with monocrotaline-induced RV failure revealed hallmarks of RyR cluster fragmentation accompanied by breaches in the JPH2 sub-domains. Frayed RyR patterns observed adjacent to these constitute new evidence that the destabilization of the RyR arrays inside the JPH2 sub-domains may seed the primordial foci of dyad remodelling observed in heart failure. This article is part of the theme issue 'The cardiomyocyte: new revelations on the interplay between architecture and function in growth, health, and disease'.

Prolonged β-adrenergic stimulation disperses ryanodine receptor clusters in cardiomyocytes

Ryanodine receptors (RyRs) exhibit dynamic arrangements in cardiomyocytes, and we previously showed that ‘dispersion’ of RyR clusters disrupts Ca²⁺ homeostasis during heart failure (HF) (Kolstad et al., eLife, 2018). Here, we investigated whether prolonged β-adrenergic stimulation, a hallmark of HF, promotes RyR cluster dispersion and examined the underlying mechanisms. We observed that treatment of healthy rat cardiomyocytes with isoproterenol for 1 hr triggered progressive fragmentation of RyR clusters. Pharmacological inhibition of Ca²⁺/calmodulin-dependent protein kinase II (CaMKII) reversed these effects, while cluster dispersion was reproduced by specific activation of CaMKII, and in mice with constitutively active Ser2814-RyR. A similar role of protein kinase A (PKA) in promoting RyR cluster fragmentation was established by employing PKA activation or inhibition. Progressive cluster dispersion was linked to declining Ca²⁺ spark fidelity and magnitude, and slowed release kinetics from Ca²⁺ propagation between more numerous RyR clusters. In healthy cells, this served to dampen the stimulatory actions of β-adrenergic stimulation over the longer term and protect against pro-arrhythmic Ca²⁺ waves. However, during HF, RyR dispersion was linked to impaired Ca²⁺ release. Thus, RyR localization and function are intimately linked via channel phosphorylation by both CaMKII and PKA, which, while finely tuned in healthy cardiomyocytes, underlies impaired cardiac function during pathology.

Molecular, Subcellular, and Arrhythmogenic Mechanisms in Genetic RyR2 Disease

The ryanodine receptor (RyR2) has a critical role in controlling Ca²⁺ release from the sarcoplasmic reticulum (SR) throughout the cardiac cycle. RyR2 protein has multiple functional domains with specific roles, and four of these RyR2 protomers are required to form the quaternary structure that comprises the functional channel. Numerous mutations in the gene encoding RyR2 protein have been identified and many are linked to a wide spectrum of arrhythmic heart disease. Gain of function mutations (GoF) result in a hyperactive channel that causes excessive spontaneous SR Ca²⁺ release. This is the predominant cause of the inherited syndrome catecholaminergic polymorphic ventricular tachycardia (CPVT). Recently, rare hypoactive loss of function (LoF) mutations have been identified that produce atypical effects on cardiac Ca²⁺ handling that has been termed calcium release deficiency syndrome (CRDS). Aberrant Ca²⁺ release resulting from both GoF and LoF mutations can result in arrhythmias through the Na⁺/Ca²⁺ exchange mechanism. This mini-review discusses recent findings regarding the role of RyR2 domains and endogenous regulators that influence RyR2 gating normally and with GoF/LoF mutations. The arrhythmogenic consequences of GoF/LoF mutations will then be discussed at the macromolecular and cellular level.

Nanoscale organization of ryanodine receptor distribution and phosphorylation pattern determines the dynamics of calcium sparks

Super-resolution imaging techniques have provided a better understanding of the relationship between the nanoscale organization and function of ryanodine receptors (RyRs) in cardiomyocytes. Recent data have indicated that this relationship is disrupted in heart failure (HF), as RyRs are dispersed into smaller and more numerous clusters. However, RyRs are also hyperphosphorylated in this condition, and this is reported to occur preferentially within the cluster centre. Thus, the combined impact of RyR relocalization and sensitization on Ca2+ spark generation in failing cardiomyocytes is likely complex and these observations suggest that both the nanoscale organization of RyRs and the pattern of phosphorylated RyRs within clusters could be critical determinants of Ca2+ spark dynamics. To test this hypothesis, we used computational modeling to quantify the relationships between RyR cluster geometry, phosphorylation patterns, and sarcoplasmic reticulum (SR) Ca2+ release. We found that RyR cluster disruption results in a decrease in spark fidelity and longer sparks with a lower amplitude. Phosphorylation of some RyRs within the cluster can play a compensatory role, recovering healthy spark dynamics. Interestingly, our model predicts that such compensation is critically dependent on the phosphorylation pattern, as phosphorylation localized within the cluster center resulted in longer Ca2+ sparks and higher spark fidelity compared to a uniformly distributed phosphorylation pattern. Our results strongly suggest that both the phosphorylation pattern and nanoscale RyR reorganization are critical determinants of Ca2+ dynamics in HF.

Structural Insight Into Ryanodine Receptor Channelopathies

The ryanodine receptors (RyRs) are large cation-selective ligand-gated channels that are expressed in the sarcoplasmic reticulum (SR) membrane. They mediate the controlled release of Ca²⁺ from SR and play an important role in many cellular processes. The mutations in RyRs are associated with several skeletal muscle and cardiac conditions, including malignant hyperthermia (MH), central core disease (CCD), catecholaminergic polymorphic ventricular tachycardia (CPVT), and arrhythmogenic right ventricular dysplasia (ARVD). Recent breakthroughs in structural biology including cryo-electron microscopy (EM) and X-ray crystallography allowed the determination of a number of near-atomic structures of RyRs, including wildtype and mutant structures as well as the structures in complex with different modulating molecules. This allows us to comprehend the physiological gating and regulatory mechanisms of RyRs and the underlying pathological mechanisms of the disease-causing mutations. In this review, based on the insights gained from the available high-resolution structures of RyRs, we address several questions: 1) what are the gating mechanisms of different RyR isoforms; 2) how RyRs are regulated by multiple channel modulators, including ions, small molecules, and regulatory proteins; 3) how do disease-causing mutations affect the structure and function of RyRs; 4) how can these structural information aid in the diagnosis of the related diseases and the development of pharmacological therapies.

CMYA5 establishes cardiac dyad architecture and positioning

Cardiac excitation-contraction coupling requires dyads, the nanoscopic microdomains formed adjacent to Z-lines by apposition of transverse tubules and junctional sarcoplasmic reticulum. Disruption of dyad architecture and function are common features of diseased cardiomyocytes. However, little is known about the mechanisms that modulate dyad organization during cardiac development, homeostasis, and disease. Here, we use proximity proteomics in intact, living hearts to identify proteins enriched near dyads. Among these proteins is CMYA5, an under-studied striated muscle protein that co-localizes with Z-lines, junctional sarcoplasmic reticulum proteins, and transverse tubules in mature cardiomyocytes. During cardiac development, CMYA5 positioning adjacent to Z-lines precedes junctional sarcoplasmic reticulum positioning or transverse tubule formation. CMYA5 ablation disrupts dyad architecture, dyad positioning at Z-lines, and junctional sarcoplasmic reticulum Ca2+ release, leading to cardiac dysfunction and inability to tolerate pressure overload. These data provide mechanistic insights into cardiomyopathy pathogenesis by demonstrating that CMYA5 anchors junctional sarcoplasmic reticulum to Z-lines, establishes dyad architecture, and regulates dyad Ca2+ release.

β2-adrenergic stimulation potentiates spontaneous calcium release by increasing signal mass and co-activation of ryanodine receptor clusters

AIMS

It is unknown how β-adrenergic stimulation affects calcium dynamics in individual RyR2 clusters and leads to the induction of spontaneous calcium waves. To address this, we analysed spontaneous calcium release events in green fluorescent protein (GFP)-tagged RyR2 clusters.

METHODS

Cardiomyocytes from mice with GFP-tagged RyR2 or human right atrial tissue were subjected to immunofluorescent labelling or confocal calcium imaging.

RESULTS

Spontaneous calcium release from single RyR2 clusters induced 91.4% ± 2.0% of all calcium sparks while 8.0% ± 1.6% were caused by release from two neighbouring clusters. Sparks with two RyR2 clusters had 40% bigger amplitude, were 26% wider, and lasted 35% longer at half maximum. Consequently, the spark mass was larger in two- than one-cluster sparks with a median and interquartile range for the cumulative distribution of 15.7 ± 20.1 vs 7.6 ± 5.7 a.u. (P < .01). β2-adrenergic stimulation increased RyR2 phosphorylation at s2809 and s2815, tripled the fraction of two- and three-cluster sparks, and significantly increased the spark mass. Interestingly, the amplitude and mass of the calcium released from a RyR2 cluster were proportional to the SR calcium load, but the firing rate was not. The spark mass was also higher in 33 patients with atrial fibrillation than in 36 without (22.9 ± 23.4 a.u. vs 10.7 ± 10.9; P = .015).

CONCLUSIONS

Most sparks are caused by activation of a single RyR2 cluster at baseline while β-adrenergic stimulation doubles the mass and the number of clusters per spark. This mimics the shift in the cumulative spark mass distribution observed in myocytes from patients with atrial fibrillation.

Multi-Scale Computational Modeling of Spatial Calcium Handling From Nanodomain to Whole-Heart: Overview and Perspectives

Regulation of intracellular calcium is a critical component of cardiac electrophysiology and excitation-contraction coupling. The calcium spark, the fundamental element of the intracellular calcium transient, is initiated in specialized nanodomains which co-locate the ryanodine receptors and L-type calcium channels. However, calcium homeostasis is ultimately regulated at the cellular scale, by the interaction of spatially separated but diffusively coupled nanodomains with other sub-cellular and surface-membrane calcium transport channels with strong non-linear interactions; and cardiac electrophysiology and arrhythmia mechanisms are ultimately tissue-scale phenomena, regulated by the interaction of a heterogeneous population of coupled myocytes. Recent advances in imaging modalities and image-analysis are enabling the super-resolution reconstruction of the structures responsible for regulating calcium homeostasis, including the internal structure of nanodomains themselves. Extrapolating functional and imaging data from the nanodomain to the whole-heart is non-trivial, yet essential for translational insight into disease mechanisms. Computational modeling has important roles to play in relating structural and functional data at the sub-cellular scale and translating data across the scales. This review covers recent methodological advances that enable image-based modeling of the single nanodomain and whole cardiomyocyte, as well as the development of multi-scale simulation approaches to integrate data from nanometer to whole-heart. Firstly, methods to overcome the computational challenges of simulating spatial calcium dynamics in the nanodomain are discussed, including image-based modeling at this scale. Then, recent whole-cell models, capable of capturing a range of different structures (such as the T-system and mitochondria) and cellular heterogeneity/variability are discussed at two different levels of discretization. Novel methods to integrate the models and data across the scales and simulate stochastic dynamics in tissue-scale models are then discussed, enabling elucidation of the mechanisms by which nanodomain remodeling underlies arrhythmia and contractile dysfunction. Perspectives on model differences and future directions are provided throughout.

Image-Driven Modeling of Nanoscopic Cardiac Function: Where Have We Come From, and Where Are We Going?

Complementary developments in microscopy and mathematical modeling have been critical to our understanding of cardiac excitation-contraction coupling. Historically, limitations imposed by the spatial or temporal resolution of imaging methods have been addressed through careful mathematical interrogation. Similarly, limitations imposed by computational power have been addressed by imaging macroscopic function in large subcellular domains or in whole myocytes. As both imaging resolution and computational tractability have improved, the two approaches have nearly merged in terms of the scales that they can each be used to interrogate. With this review we will provide an overview of these advances and their contribution to understanding ventricular myocyte function, including exciting developments over the last decade. We specifically focus on experimental methods that have pushed back limits of either spatial or temporal resolution of nanoscale imaging (e.g., DNA-PAINT), or have permitted high resolution imaging on large cellular volumes (e.g., serial scanning electron microscopy). We also review the progression of computational approaches used to integrate and interrogate these new experimental data sources, and comment on near-term advances that may unify understanding of the underlying biology. Finally, we comment on several outstanding questions in cardiac physiology that stand to benefit from a concerted and complementary application of these new experimental and computational methods.

Magnesium Ions Moderate Calcium-Induced Calcium Release in Cardiac Calcium Release Sites by Binding to Ryanodine Receptor Activation and Inhibition Sites

Ryanodine receptor channels at calcium release sites of cardiac myocytes operate on the principle of calcium-induced calcium release. In vitro experiments revealed competition of Ca2+ and Mg2+ in the activation of ryanodine receptors (RyRs) as well as inhibition of RyRs by Mg2+. The impact of RyR modulation by Mg2+ on calcium release is not well understood due to the technical limitations of in situ experiments. We turned instead to an in silico model of a calcium release site (CRS), based on a homotetrameric model of RyR gating with kinetic parameters determined from in vitro measurements. We inspected changes in the activity of the CRS model in response to a random opening of one of 20 realistically distributed RyRs, arising from Ca2+/Mg2+ interactions at RyR channels. Calcium release events (CREs) were simulated at a range of Mg2+-binding parameters at near-physiological Mg2+ and ATP concentrations. Facilitation of Mg2+ binding to the RyR activation site inhibited the formation of sparks and slowed down their activation. Impeding Mg-binding to the RyR activation site enhanced spark formation and speeded up their activation. Varying Mg2+ binding to the RyR inhibition site also dramatically affected calcium release events. Facilitation of Mg2+ binding to the RyR inhibition site reduced the amplitude, relative occurrence, and the time-to-end of sparks, and vice versa. The characteristics of CREs correlated dose-dependently with the effective coupling strength between RyRs, defined as a function of RyR vicinity, single-channel calcium current, and Mg-binding parameters of the RyR channels. These findings postulate the role of Mg2+ in calcium release as a negative modulator of the coupling strength among RyRs in a CRS, translating to damping of the positive feedback of the calcium-induced calcium-release mechanism.

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.

A Stochastic Spatiotemporal Model of Rat Ventricular Myocyte Calcium Dynamics Demonstrated Necessary Features for Calcium Wave Propagation

Calcium (Ca2+) plays a central role in the excitation and contraction of cardiac myocytes. Experiments have indicated that calcium release is stochastic and regulated locally suggesting the possibility of spatially heterogeneous calcium levels in the cells. This spatial heterogeneity might be important in mediating different signaling pathways. During more than 50 years of computational cell biology, the computational models have been advanced to incorporate more ionic currents, going from deterministic models to stochastic models. While periodic increases in cytoplasmic Ca2+ concentration drive cardiac contraction, aberrant Ca2+ release can underly cardiac arrhythmia. However, the study of the spatial role of calcium ions has been limited due to the computational expense of using a three-dimensional stochastic computational model. In this paper, we introduce a three-dimensional stochastic computational model for rat ventricular myocytes at the whole-cell level that incorporate detailed calcium dynamics, with (1) non-uniform release site placement, (2) non-uniform membrane ionic currents and membrane buffers, (3) stochastic calcium-leak dynamics and (4) non-junctional or rogue ryanodine receptors. The model simulates spark-induced spark activation and spark-induced Ca2+ wave initiation and propagation that occur under conditions of calcium overload at the closed-cell condition, but not when Ca2+ levels are normal. This is considered important since the presence of Ca2+ waves contribute to the activation of arrhythmogenic currents.

Nanoscale Organisation of Ryanodine Receptors and Junctophilin-2 in the Failing Human Heart

The disrupted organisation of the ryanodine receptors (RyR) and junctophilin (JPH) is thought to underpin the transverse tubule (t-tubule) remodelling in a failing heart. Here, we assessed the nanoscale organisation of these two key proteins in the failing human heart. Recently, an advanced feature of the t-tubule remodelling identified large flattened t-tubules called t-sheets, that were several microns wide. Previously, we reported that in the failing heart, the dilated t-tubules up to ~1 μm wide had increased collagen, and we hypothesised that the t-sheets would also be associated with collagen deposits. Direct stochastic optical reconstruction microscopy (dSTORM), confocal microscopy, and western blotting were used to evaluate the cellular distribution of excitation-contraction structures in the cardiac myocytes from patients with idiopathic dilated cardiomyopathy (IDCM) compared to myocytes from the non-failing (NF) human heart. The dSTORM imaging of RyR and JPH found no difference in the colocalisation between IDCM and NF myocytes, but there was a higher colocalisation at the t-tubule and sarcolemma compared to the corbular regions. Western blots revealed no change in the JPH expression but did identify a ~50% downregulation of RyR (p = 0.02). The dSTORM imaging revealed a trend for the smaller t-tubular RyR clusters (~24%) and reduced the t-tubular RyR cluster density (~35%) that resulted in a 50% reduction of t-tubular RyR tetramers in the IDCM myocytes (p < 0.01). Confocal microscopy identified the t-sheets in all the IDCM hearts examined and found that they are associated with the reticular collagen fibres within the lumen. However, the size and density of the RyR clusters were similar in the myocyte regions associated with t-sheets and t-tubules. T-tubule remodelling is associated with a reduced RyR expression that may contribute to the reduced excitation-contraction coupling in the failing human heart.

Modulated Calcium Homeostasis and Release Events Under Atrial Fibrillation and Its Risk Factors: A Meta-Analysis

Background: Atrial fibrillation (AF) is associated with calcium (Ca²⁺) handling remodeling and increased spontaneous calcium release events (SCaEs). Nevertheless, its exact mechanism remains unclear, resulting in suboptimal primary and secondary preventative strategies.

Methods: We searched the PubMed database for studies that investigated the relationship between SCaEs and AF and/or its risk factors. Meta-analysis was used to examine the Ca²⁺ mechanisms involved in the primary and secondary AF preventative groups.

Results: We included a total of 74 studies, out of the identified 446 publications from inception (1982) until March 31, 2020. Forty-five were primary and 29 were secondary prevention studies for AF. The main Ca²⁺ release events, calcium transient (standardized mean difference (SMD) = 0.49; I² = 35%; confidence interval (CI) = 0.33–0.66; p < 0.0001), and spark amplitude (SMD = 0.48; I² = 0%; CI = −0.98–1.93; p = 0.054) were enhanced in the primary diseased group, while calcium transient frequency was increased in the secondary group. Calcium spark frequency was elevated in both the primary diseased and secondary AF groups. One of the key cardiac currents, the L-type calcium current (I_CaL) was significantly downregulated in primary diseased (SMD = −1.07; I² = 88%; CI = −1.94 to −0.20; p < 0.0001) and secondary AF groups (SMD = −1.28; I² = 91%; CI = −2.04 to −0.52; p < 0.0001). Furthermore, the sodium–calcium exchanger (I_NCX) and NCX1 protein expression were significantly enhanced in the primary diseased group, while only NCX1 protein expression was shown to increase in the secondary AF studies. The phosphorylation of the ryanodine receptor at S2808 (pRyR-S2808) was significantly elevated in both the primary and secondary groups. It was increased in the primary diseased and proarrhythmic subgroups (SMD = 0.95; I² = 64%; CI = 0.12–1.79; p = 0.074) and secondary AF group (SMD = 0.66; I² = 63%; CI = 0.01–1.31; p < 0.0001). Sarco/endoplasmic reticulum Ca²⁺-ATPase (SERCA) expression was elevated in the primary diseased and proarrhythmic drug subgroups but substantially reduced in the secondary paroxysmal AF subgroup.

Conclusions: Our study identified that I_CaL is reduced in both the primary and secondary diseased groups. Furthermore, pRyR-S2808 and NCX1 protein expression are enhanced. The remodeling leads to elevated Ca²⁺ functional activities, such as increased frequencies or amplitude of Ca²⁺ spark and Ca²⁺ transient. The main difference identified between the primary and secondary diseased groups is SERCA expression, which is elevated in the primary diseased group and substantially reduced in the secondary paroxysmal AF subgroup. We believe our study will add new evidence to AF mechanisms and treatment targets.

Optimising Large Animal Models of Sustained Atrial Fibrillation: Relevance of the Critical Mass Hypothesis

BACKGROUND

Large animal models play an important role in our understanding of the pathophysiology of atrial fibrillation (AF). Our aim was to determine whether prospectively collected baseline variables could predict the development of sustained AF in sheep, thereby reducing the number of animals required in future studies. Our hypothesis was that the relationship between atrial dimensions, refractory periods and conduction velocity (otherwise known as the critical mass hypothesis) could be used for the first time to predict the development of sustained AF.

METHODS

Healthy adult Welsh mountain sheep underwent a baseline electrophysiology study followed by implantation of a neurostimulator connected via an endocardial pacing lead to the right atrial appendage. The device was programmed to deliver intermittent 50 Hz bursts of 30 s duration over an 8-week period whilst sheep were monitored for AF.

RESULTS

Eighteen sheep completed the protocol, of which 28% developed sustained AF. Logistic regression analysis showed only fibrillation number (calculated using the critical mass hypothesis as the left atrial diameter divided by the product of atrial conduction velocity and effective refractory period) was associated with an increased likelihood of developing sustained AF (Ln Odds Ratio 26.1 [95% confidence intervals 0.2-52.0] p = 0.048). A receiver-operator characteristic curve showed this could be used to predict which sheep developed sustained AF (C-statistic 0.82 [95% confidence intervals 0.59-1.04] p = 0.04).

CONCLUSION

The critical mass hypothesis can be used to predict sustained AF in a tachypaced ovine model. These findings can be used to optimise the design of future studies involving large animals.

Impaired Binding to Junctophilin-2 and Nanostructural Alteration in CPVT Mutation

[Figure: see text].

Cx43 hemichannel microdomain signaling at the intercalated disc enhances cardiac excitability

Cx43, a major cardiac connexin, forms precursor hemichannels that accrue at the intercalated disc to assemble as gap junctions. While gap junctions are crucial for electrical conduction in the heart, little is known about the potential roles of hemichannels. Recent evidence suggests that inhibiting Cx43 hemichannel opening with Gap19 has antiarrhythmic effects. Here, we used multiple electrophysiology, imaging, and super-resolution techniques to understand and define the conditions underlying Cx43 hemichannel activation in ventricular cardiomyocytes, their contribution to diastolic Ca2+ release from the sarcoplasmic reticulum, and their impact on electrical stability. We showed that Cx43 hemichannels were activated during diastolic Ca2+ release in single ventricular cardiomyocytes and cardiomyocyte cell pairs from mice and pigs. This activation involved Cx43 hemichannel Ca2+ entry and coupling to Ca2+ release microdomains at the intercalated disc, resulting in enhanced Ca2+ dynamics. Hemichannel opening furthermore contributed to delayed afterdepolarizations and triggered action potentials. In single cardiomyocytes, cardiomyocyte cell pairs, and arterially perfused tissue wedges from failing human hearts, increased hemichannel activity contributed to electrical instability compared with nonfailing rejected donor hearts. We conclude that microdomain coupling between Cx43 hemichannels and Ca2+ release is a potentially novel, targetable mechanism of cardiac arrhythmogenesis in heart failure.

In silico simulations reveal that RYR distribution affects the dynamics of calcium release in cardiac myocytes

The dyads of cardiac myocytes contain ryanodine receptors (RYRs) that generate calcium sparks upon activation. To test how geometric factors of RYR distribution contribute to the formation of calcium sparks, which cannot be addressed experimentally, we performed in silico simulations on a large set of models of calcium release sites (CRSs). Our models covered the observed range of RYR number, density, and spatial arrangement. The calcium release function of CRSs was modeled by RYR openings, with an open probability dependent on concentrations of free Ca2+ and Mg2+ ions, in a rapidly buffered system, with a constant open RYR calcium current. We found that simulations of spontaneous sparks by repeatedly opening one of the RYRs in a CRS produced three different types of calcium release events (CREs) in any of the models. Transformation of simulated CREs into fluorescence signals yielded calcium sparks with characteristics close to the observed ones. CRE occurrence varied broadly with the spatial distribution of RYRs in the CRS but did not consistently correlate with RYR number, surface density, or calcium current. However, it correlated with RYR coupling strength, defined as the weighted product of RYR vicinity and calcium current, so that CRE characteristics of all models followed the same state-response function. This finding revealed the synergy between structure and function of CRSs in shaping dyad function. Lastly, rearrangements of RYRs simulating hypothetical experiments on splitting and compaction of a dyad revealed an increased propensity to generate spontaneous sparks and an overall increase in calcium release in smaller and more compact dyads, thus underlying the importance and physiological role of RYR arrangement in cardiac myocytes.

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.

IP3-mediated Ca2+ release regulates atrial Ca2+ transients and pacemaker function by stimulation of adenylyl cyclases

Inositol trisphosphate (IP₃) is a Ca²⁺-mobilizing second messenger shown to modulate atrial muscle contraction and is thought to contribute to atrial fibrillation. Cellular pathways underlying IP₃ actions in cardiac tissue remain poorly understood, and the work presented here addresses the question whether IP₃-mediated Ca²⁺ release from the sarcoplasmic reticulum is linked to adenylyl cyclase activity including Ca²⁺-stimulated adenylyl cyclases (AC1 and AC8) that are selectively expressed in atria and sinoatrial node (SAN). Immunocytochemistry in guinea pig atrial myocytes identified colocalization of type 2 IP₃ receptors with AC8, while AC1 was located in close vicinity. Intracellular photorelease of IP₃ by UV light significantly enhanced the amplitude of the Ca²⁺ transient (CaT) evoked by electrical stimulation of atrial myocytes (31 ± 6% increase 60 s after photorelease, n = 16). The increase in CaT amplitude was abolished by inhibitors of adenylyl cyclases (MDL-12,330) or protein kinase A (H89), showing that cAMP signaling is required for this effect of photoreleased IP₃. In mouse, spontaneously beating right atrial preparations, phenylephrine, an α-adrenoceptor agonist with effects that depend on IP₃-mediated Ca²⁺ release, increased the maximum beating rate by 14.7 ± 0.5%, n = 10. This effect was substantially reduced by 2.5 µmol/L 2-aminoethyl diphenylborinate and abolished by a low dose of MDL-12,330, observations which are again consistent with a functional interaction between IP₃ and cAMP signaling involving Ca²⁺ stimulation of adenylyl cyclases in the SAN pacemaker. Understanding the interaction between IP₃ receptor pathways and Ca²⁺-stimulated adenylyl cyclases provides important insights concerning acute mechanisms for initiation of atrial arrhythmias.

NEW & NOTEWORTHY This study provides evidence supporting the proposal that IP₃ signaling in cardiac atria and sinoatrial node involves stimulation of Ca²⁺-activated adenylyl cyclases (AC1 and AC8) by IP₃-evoked Ca²⁺ release from junctional sarcoplasmic reticulum. AC8 and IP₃ receptors are shown to be located close together, while AC1 is nearby. Greater understanding of these novel aspects of the IP₃ signal transduction mechanism is important for future study in atrial physiology and pathophysiology, particularly atrial fibrillation.

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.

Enhanced expansion microscopy to measure nanoscale structural and biochemical remodeling in single cells

Resolution is a key feature in microscopy which allows the visualization of the fine structure of cells. Much of the life processes within these cells depend on the three-dimensional (3D) complexity of these structures. Optical super-resolution microscopies are currently the preferred choice of molecular and cell biologists who seek to visualize the organization of specific protein species at the nanometer scale. Traditional super-resolution microscopy techniques have often been limited by sample thickness, axial resolution, specialist optical instrumentation and computationally-demanding software for assembling the images. In this chapter we detail the protocol, "enhanced expansion microscopy" (EExM), which combines X10 expansion microscopy with Airyscan confocal microscopy. EExM enables 15nm lateral (and 35nm axial) resolution, and is a relatively cheap, accessible option allowing single protein resolution for the non-specialist optical microscopists. We illustrate how EExM has been utilized for mapping the 3D topology of intracellular protein arrays at sample depths which are not always compatible with some of the traditional super-resolution techniques. We demonstrate that antibody markers can recognize and map post-translational modifications of individual proteins in addition to their 3D positions. Finally, we discuss the current uncertainties and validations in EExM which include the isotropy in gel expansion and assessment of the expansion factor (of resolution improvement).

Two-variable nullcline analysis of ionic general equilibrium predicts calcium homeostasis in ventricular myocytes

Ventricular contraction is roughly proportional to the amount of calcium released from the Sarcoplasmic Reticulum (SR) during systole. While it is rather straightforward to measure calcium levels and contractibility under different physiological conditions, the complexity of calcium handling during systole and diastole has made the prediction of its release at steady state impossible. Here we approach the problem analyzing the evolution of intracellular and extracellular calcium fluxes during a single beat which is away from homeostatic balance. Using an in-silico subcellular model of rabbit ventricular myocyte, we show that the high dimensional nonlinear problem of finding the steady state can be reduced to a two-variable general equilibrium condition where pre-systolic calcium level in the cytosol and in the SR must fulfill simultaneously two different equalities. This renders calcium homeostasis as a problem that can be studied in terms of its equilibrium structure, leading to precise predictions of steady state from single-beat measurements. We show how changes in ion channels modify the general equilibrium, as shocks would do in general equilibrium macroeconomic models. This allows us to predict when an enhanced entrance of calcium in the cell reduces its contractibility and explain why SERCA gene therapy, a change in calcium handling to treat heart failure, might fail to improve contraction even when it successfully increases SERCA expression.

Cardiomyocytes, upon voltage excitation, release calcium, which leads to cell contraction. However, under some pathological conditions, calcium handling is impaired. Recently, SERCA gene therapy, whose aim is to improve Ca²⁺ sequestration by the Sarcoplasmic Reticulum (SR), has failed to improve the prognosis of patients with Heart Failure. This, together with recent counterintuitive results in calcium handling, has highlighted the need for a framework to understand calcium homeostasis across species and pathologies. We show here that the proper framework is a general equilibrium approach of two independent variables. The development of this framework allows us to find a possible mechanism for the failure of SERCA gene therapy even when it manages to increase Ca SERCA expression.

Conditional Up-Regulation of SERCA2a Exacerbates RyR2-Dependent Ventricular and Atrial Arrhythmias

Ryanodine receptor 2 (RyR2) and SERCA2a are two major players in myocyte calcium (Ca) cycling that are modulated physiologically, affected by disease and thus considered to be potential targets for cardiac disease therapy. However, how RyR2 and SERCA2a influence each others’ activities, as well as the primary and secondary consequences of their combined manipulations remain controversial. In this study, we examined the effect of acute upregulation of SERCA2a on arrhythmogenesis by conditionally overexpressing SERCA2a in a mouse model featuring hyperactive RyR2s due to ablation of calsequestrin 2 (CASQ2). CASQ2 knock-out (KO) mice were crossbred with doxycycline (DOX)-inducible SERCA2a transgenic mice to generate KO-TG mice. In-vivo ECG studies have shown that induction of SERCA2a (DOX+) overexpression markedly exacerbated both ventricular and atrial arrhythmias in vivo, compared with uninduced KO-TG mice (DOX-). Consistent with that, confocal microscopy in both atrial and ventricular myocytes demonstrated that conditional upregulation of SERCA2a enhanced the rate of occurrence of diastolic Ca release events. Additionally, deep RNA sequencing identified 17 downregulated genes and 5 upregulated genes in DOX+ mice, among which Ppp1r13l, Clcn1, and Agt have previously been linked to arrhythmias. Our results suggest that conditional upregulation of SERCA2a exacerbates hyperactive RyR2-mediated arrhythmias by further elevating diastolic Ca release.

Cardiomyocyte calcium handling in health and disease: Insights from in vitro and in silico studies

Calcium (Ca²⁺) plays a central role in cardiomyocyte excitation-contraction coupling. To ensure an optimal electrical impulse propagation and cardiac contraction, Ca²⁺ levels are regulated by a variety of Ca²⁺-handling proteins. In turn, Ca²⁺ modulates numerous electrophysiological processes. Accordingly, Ca²⁺-handling abnormalities can promote cardiac arrhythmias via various mechanisms, including the promotion of afterdepolarizations, ion-channel modulation and structural remodeling. In the last 30 years, significant improvements have been made in the computational modeling of cardiomyocyte Ca²⁺ handling under physiological and pathological conditions. However, numerous questions involving the Ca²⁺-dependent regulation of different macromolecular complexes, cross-talk between Ca²⁺-dependent regulatory pathways operating over a wide range of time scales, and bidirectional interactions between electrophysiology and mechanics remain to be addressed by in vitro and in silico studies. A better understanding of disease-specific Ca²⁺-dependent proarrhythmic mechanisms may facilitate the development of improved therapeutic strategies. In this review, we describe the fundamental mechanisms of cardiomyocyte Ca²⁺ handling in health and disease, and provide an overview of currently available computational models for cardiomyocyte Ca²⁺ handling. Finally, we discuss important uncertainties and open questions about cardiomyocyte Ca²⁺ handling and highlight how synergy between in vitro and in silico studies may help to answer several of these issues.

Calcium Signaling in Cardiomyocyte Function

Rhythmic increases in intracellular Ca²⁺ concentration underlie the contractile function of the heart. These heart muscle-wide changes in intracellular Ca²⁺ are induced and coordinated by electrical depolarization of the cardiomyocyte sarcolemma by the action potential. Originating at the sinoatrial node, conduction of this electrical signal throughout the heart ensures synchronization of individual myocytes into an effective cardiac pump. Ca²⁺ signaling pathways also regulate gene expression and cardiomyocyte growth during development and in pathology. These fundamental roles of Ca²⁺ in the heart are illustrated by the prevalence of altered Ca²⁺ homeostasis in cardiovascular diseases. Indeed, heart failure (an inability of the heart to support hemodynamic needs), rhythmic disturbances, and inappropriate cardiac growth all share an involvement of altered Ca²⁺ handling. The prevalence of these pathologies, contributing to a third of all deaths in the developed world as well as to substantial morbidity makes understanding the mechanisms of Ca²⁺ handling and dysregulation in cardiomyocytes of great importance.

Calcium Buffering in the Heart in Health and Disease

Changes of intracellular Ca²⁺ concentration regulate many aspects of cardiac myocyte function. About 99% of the cytoplasmic calcium in cardiac myocytes is bound to buffers, and their properties will therefore have a major influence on Ca²⁺ signaling. This article considers the fundamental properties and identities of the buffers and how to measure them. It reviews the effects of buffering on the systolic Ca²⁺ transient and how this may change physiologically, and in heart failure and both atrial and ventricular arrhythmias, as well. It is concluded that the consequences of this strong buffering may be more significant than currently appreciated, and a fuller understanding is needed for proper understanding of cardiac calcium cycling and contractility.

Remodeling Promotes Proarrhythmic Disruption of Calcium Homeostasis in Failing Atrial Myocytes

It is well known that heart failure (HF) typically coexists with atrial fibrillation (AF). However, until now, no clear mechanism has been established that relates HF to AF. In this study, we apply a multiscale computational framework to establish a mechanistic link between atrial myocyte structural remodeling in HF and AF. Using a spatially distributed model of calcium (Ca) signaling, we show that disruption of the spatial relationship between L-type Ca channels (LCCs) and ryanodine receptors results in markedly increased Ca content of the sarcoplasmic reticulum (SR). This increase in SR load is due to changes in the balance between Ca entry via LCCs and Ca extrusion due to the sodium-calcium exchanger after an altered spatial relationship between these signaling proteins. Next, we show that the increased SR load in atrial myocytes predisposes these cells to subcellular Ca waves that occur during the action potential (AP) and are triggered by LCC openings. These waves are common in atrial cells because of the absence of a well-developed t-tubule system in most of these cells. This distinct spatial architecture allows for the presence of a large pool of orphaned ryanodine receptors, which can fire and sustain Ca waves during the AP. Finally, we incorporate our atrial cell model in two-dimensional tissue simulations and demonstrate that triggered wave generation in cells leads to electrical waves in tissue that tend to fractionate to form wavelets of excitation. This fractionation is driven by the underlying stochasticity of subcellular Ca waves, which perturbs AP repolarization and consequently induces localized conduction block in tissue. We outline the mechanism for this effect and argue that it may explain the propensity for atrial arrhythmias in HF.