Calcium spike variability in cardiac myocytes results from activation of small cohorts of ryanodine receptor 2 channels

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.

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.

Structure and Function of the Human Ryanodine Receptors and Their Association with Myopathies-Present State, Challenges, and Perspectives

Cardiac arrhythmias are serious, life-threatening diseases associated with the dysregulation of Ca2+ influx into the cytoplasm of cardiomyocytes. This dysregulation often arises from dysfunction of ryanodine receptor 2 (RyR2), the principal Ca2+ release channel. Dysfunction of RyR1, the skeletal muscle isoform, also results in less severe, but also potentially life-threatening syndromes. The RYR2 and RYR1 genes have been found to harbor three main mutation “hot spots”, where mutations change the channel structure, its interdomain interface properties, its interactions with its binding partners, or its dynamics. In all cases, the result is a defective release of Ca2+ ions from the sarcoplasmic reticulum into the myocyte cytoplasm. Here, we provide an overview of the most frequent diseases resulting from mutations to RyR1 and RyR2, briefly review some of the recent experimental structural work on these two molecules, detail some of the computational work describing their dynamics, and summarize the known changes to the structure and function of these receptors with particular emphasis on their N-terminal, central, and channel domains.

Structural variability of dyads relates to calcium release in rat ventricular myocytes

Cardiac excitation-contraction coupling relies on dyads, the intracellular calcium synapses of cardiac myocytes, where the plasma membrane contacts sarcoplasmic reticulum and where electrical excitation triggers calcium release. The morphology of dyads and dynamics of local calcium release vary substantially. To better understand the correspondence between the structure and the functionality of dyads, we estimated incidences of structurally different dyads and of kinetically different calcium release sites and tested their responsiveness to experimental myocardial injury in left ventricular myocytes of rats. According to the structure of dyads estimated in random electron microscopic images of myocardial tissue, the dyads were sorted into 'compact' or 'loose' types. The calcium release fluxes, triggered at local calcium release sites in patch-clamped ventricular myocytes and recorded by laser scanning confocal fluorescence microscopy, were decomposed into 'early' and 'late' components. ANOVA tests revealed very high correlation between the relative amplitudes of early and late calcium release flux components and the relative occurrences of compact and loose dyads in the control and in the injured myocardium. This finding ascertained the relationship between the structure of dyads and the functionality of calcium release sites and the responsiveness of calcium release sites to physical load in cardiac myocytes.

Stochastic and deterministic approaches to modelling calcium release in cardiac myocytes at different spatial arrangements of ryanodine receptors

Calcium release sites (CRSs) play a key role in excitation-contraction coupling of cardiac myocytes. Recent studies based on electron tomography and super-resolution imaging revealed that CRSs are not completely filled with ryanodine receptors (RyRs) and that the spatial arrangement of RyRs is neither uniform nor static. In this work, we studied the effect of spatial arrangement of RyRs on RyR activation using simulations based on Monte Carlo (MC) and mean-field (MF) approaches. Both approaches showed that activation of RyRs is sensitive to the arrangement of RyRs in the CRS. However, the MF simulations did not reproduce results of MC simulations for non-compact CRSs, suggesting that the approximations used in the MF approach are not suitable for simulation studies of RyRs arrangements observed experimentally. MC simulations revealed the importance of realistic spatial arrangement of RyRs for adequate modelling of calcium release in cardiac myocytes.

Calcium Signaling and Contractility in Cardiac Myocyte of Wolframin Deficient Rats

Wolframin (Wfs1) is a membrane protein of the sarco/endoplasmic reticulum. Wfs1 mutations are responsible for the Wolfram syndrome, characterized by diabetic and neurological symptoms. Although Wfs1 is expressed in cardiac muscle, its role in this tissue is not clear. We have characterized the effect of invalidation of Wfs1 on calcium signaling-related processes in isolated ventricular myocytes of exon5-Wfs1 deficient rats (Wfs1-e5/-e5) before the onset of overt disease. Calcium transients and contraction were measured in field-stimulated isolated myocytes using confocal microscopy with calcium indicator fluo-3 AM and sarcomere length detection. Calcium currents and their calcium release-dependent inactivation were characterized in whole-cell patch-clamp experiments. At 4 months, Wfs1-e5/-e5 animals were euglycemic, and echocardiographic examination revealed fully compensated cardiac function. In field-stimulated isolated ventricular myocytes, both the amplitude and the duration of contraction of Wfs1-e5/-e5 animals were elevated relative to control Wfs1+/+ littermates. Increased contractility of myocytes resulted largely from prolonged cytosolic calcium transients. Neither the amplitude of calcium currents nor their voltage dependence of activation differed between the two groups. Calcium currents in Wfs1-e5/-e5 myocytes showed a larger extent of inactivation by short voltage prepulses applied to selectively induce calcium release-dependent inactivation of calcium current. Neither the calcium content of the sarcoplasmic reticulum, measured by application of 20 mmol/l caffeine, nor the expression of SERCA2, determined from Western blots, differed significantly in myocytes of Wfs1-e5/-e5 animals compared to control ones. These experiments point to increased duration of calcium release in ventricular myocytes of Wfs1-e5/-e5 animals. We speculate that the lack of functional wolframin might cause changes leading to upregulation of RyR2 channels resulting in prolongation of channel openings and/or a delay in termination of calcium release.

Reconstruction of membrane current by deconvolution and its application to membrane capacitance measurements in cardiac myocytes

Correct detection of membrane currents under whole-cell patch–clamp conditions is limited by the transfer function of a recording system. The low-pass output filter of a recording amplifier alters the time course of membrane current and causes errors in relevant descriptors. To solve these problems, we developed a practical procedure for reconstruction of the time course of membrane currents based on deconvolution of recorded currents in frequency domain. The procedure was tested on membrane capacitance estimates from current responses to step voltage pulses. The reconstructed current responses, in contrast to original current records, could be described exactly by an adequate impedance model of a recorded cell. The reconstruction allowed to increase the accuracy and reliability of membrane capacitance measurements in wide range of cell sizes and to suppress the cross-talk errors well below the noise. Moreover, it allowed resolving the instabilities in recording conditions arising from parasitic capacitance and seal resistance variation. Complex tests on hardware models, on simulated data sets, and on living cells confirmed robustness and reliability of the deconvolution procedure. The aptitude of the method was demonstrated in isolated rat cardiac myocytes by recording of spontaneous vesicular events, by discerning the formation of a fusion pore, and by revealing artefacts due to unstable seal resistance.

An adaptation of astronomical image processing enables characterization and functional 3D mapping of individual sites of excitation-contraction coupling in rat cardiac muscle

In beating cardiomyocytes, synchronized localized Ca²⁺ transients from thousands of active excitation-contraction coupling sites (ECC couplons) comprising plasma and sarcoplasmic reticulum membrane calcium channels are important determinants of the heart's performance. Nevertheless, our knowledge about the properties of ECC couplons is limited by the lack of appropriate experimental and analysis strategies. We designed CaCLEAN to untangle the fundamental characteristics of ECC couplons by combining the astronomer's CLEAN algorithm with known properties of calcium diffusion. CaCLEAN empowers the investigation of fundamental properties of ECC couplons in beating cardiomyocytes without pharmacological interventions. Upon examining individual ECC couplons at the nanoscopic level, we reveal their roles in the negative amplitude-frequency relationship and in β-adrenergic stimulation, including decreasing and increasing firing reliability, respectively. CaCLEAN combined with 3D confocal imaging of beating cardiomyocytes provides a functional 3D map of active ECC couplons (on average, 17,000 per myocyte). CaCLEAN will further enlighten the ECC-couplon-remodelling processes that underlie cardiac diseases.

Ryanodine receptor gating controls generation of diastolic calcium waves in cardiac myocytes

Calcium waves can form and propagate at low frequencies of spontaneous calcium sparks if the calcium dependence of spark frequency is sufficiently steep, or the number of open RyRs is sufficiently large.

The role of cardiac ryanodine receptor (RyR) gating in the initiation and propagation of calcium waves was investigated using a mathematical model comprising a stochastic description of RyR gating and a deterministic description of calcium diffusion and sequestration. We used a one-dimensional array of equidistantly spaced RyR clusters, representing the confocal scanning line, to simulate the formation of calcium sparks. Our model provided an excellent description of the calcium dependence of the frequency of diastolic calcium sparks and of the increased tendency for the production of calcium waves after a decrease in cytosolic calcium buffering. We developed a hypothesis relating changes in the propensity to form calcium waves to changes of RyR gating and tested it by simulation. With a realistic RyR gating model, increased ability of RyR to be activated by Ca²⁺ strongly increased the propensity for generation of calcium waves at low (0.05–0.1-µM) calcium concentrations but only slightly at high (0.2–0.4-µM) calcium concentrations. Changes in RyR gating altered calcium wave formation by changing the calcium sensitivity of spontaneous calcium spark activation and/or the average number of open RyRs in spontaneous calcium sparks. Gating changes that did not affect RyR activation by Ca²⁺ had only a weak effect on the propensity to form calcium waves, even if they strongly increased calcium spark frequency. Calcium waves induced by modulating the properties of the RyR activation site could be suppressed by inhibiting the spontaneous opening of the RyR. These data can explain the increased tendency for production of calcium waves under conditions when RyR gating is altered in cardiac diseases.

Membrane Potential and Calcium Dynamics in Beta Cells from Mouse Pancreas Tissue Slices: Theory, Experimentation, and Analysis

Beta cells in the pancreatic islets of Langerhans are precise biological sensors for glucose and play a central role in balancing the organism between catabolic and anabolic needs. A hallmark of the beta cell response to glucose are oscillatory changes of membrane potential that are tightly coupled with oscillatory changes in intracellular calcium concentration which, in turn, elicit oscillations of insulin secretion. Both membrane potential and calcium changes spread from one beta cell to the other in a wave-like manner. In order to assess the properties of the abovementioned responses to physiological and pathological stimuli, the main challenge remains how to effectively measure membrane potential and calcium changes at the same time with high spatial and temporal resolution, and also in as many cells as possible. To date, the most wide-spread approach has employed the electrophysiological patch-clamp method to monitor membrane potential changes. Inherently, this technique has many advantages, such as a direct contact with the cell and a high temporal resolution. However, it allows one to assess information from a single cell only. In some instances, this technique has been used in conjunction with CCD camera-based imaging, offering the opportunity to simultaneously monitor membrane potential and calcium changes, but not in the same cells and not with a reliable cellular or subcellular spatial resolution. Recently, a novel family of highly-sensitive membrane potential reporter dyes in combination with high temporal and spatial confocal calcium imaging allows for simultaneously detecting membrane potential and calcium changes in many cells at a time. Since the signals yielded from both types of reporter dyes are inherently noisy, we have developed complex methods of data denoising that permit for visualization and pixel-wise analysis of signals. Combining the experimental approach of high-resolution imaging with the advanced analysis of noisy data enables novel physiological insights and reassessment of current concepts in unprecedented detail.

Quantitative analysis of calcium spikes in noisy fluorescent background

Intracellular calcium signals are studied by laser-scanning confocal fluorescence microscopy. The required spatio-temporal resolution makes description of calcium signals difficult because of the low signal-to-noise ratio. We designed a new procedure of calcium spike analysis based on their fitting with a model. The accuracy and precision of calcium spike description were tested on synthetic datasets generated either with randomly varied spike parameters and Gaussian noise of constant amplitude, or with constant spike parameters and Gaussian noise of various amplitudes. Statistical analysis was used to evaluate the performance of spike fitting algorithms. The procedure was optimized for reliable estimation of calcium spike parameters and for dismissal of false events. A new algorithm was introduced that corrects the acquisition time of pixels in line-scan images that is in error due to sequential acquisition of individual pixels along the space coordinate. New software was developed in Matlab and provided for general use. It allows interactive dissection of temporal profiles of calcium spikes from x-t images, their fitting with predefined function(s) and acceptance of results on statistical grounds, thus allowing efficient analysis and reliable description of calcium signaling in cardiac myocytes down to the in situ function of ryanodine receptors.