Construction of calcium release sites in cardiac myocytes

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.

Mechano-Electric Coupling and Arrhythmogenic Current Generation in a Computational Model of Coupled Myocytes

A wide range of arrhythmogenic phenotypes have been associated with heterogeneous mechanical dyskinesis. Pro-arrhythmic effects are often associated with dysregulated intra-cellular calcium handling, especially via the development of intra- and inter-cellular calcium waves. Experimental evidence suggests that mechanical strain can contribute to the generation and maintenance of these calcium waves via a variety of mechano-electric coupling mechanisms. Most model studies of mechano-electric coupling mechanisms have been focused on mechano-sensitive ion channels, even though experimental studies have shown that intra- and inter-cellular calcium waves triggered by mechanical perturbations are likely to be more prevalent pro-arrhythmic mechanisms in the diseased heart. A one-dimensional strongly coupled computational model of electromechanics in rabbit ventricular cardiomyocytes showed that specific myocyte stretch sequences can modulate the susceptibility threshold for delayed after-depolarizations. In simulations of mechanically-triggered calcium waves in cardiomyocytes coupled to fibroblasts, susceptibility to calcium wave propagation was reduced as the current through the gap junction caused current drain from the myocytes. In 1D multi-cellular arrays coupled via gap junctions, mechanically-induced waves may contribute to synchronizing arrhythmogenic calcium waves and after-depolarizations.

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.

The Subcellular Distribution of Ryanodine Receptors and L-Type Ca2+ Channels Modulates Ca2+-Transient Properties and Spontaneous Ca2+-Release Events in Atrial Cardiomyocytes

Spontaneous Ca²⁺-release events (SCaEs) from the sarcoplasmic reticulum play crucial roles in the initiation of cardiac arrhythmias by promoting triggered activity. However, the subcellular determinants of these SCaEs remain incompletely understood. Structural differences between atrial and ventricular cardiomyocytes, e.g., regarding the density of T-tubular membrane invaginations, may influence cardiomyocyte Ca²⁺-handling and the distribution of cardiac ryanodine receptors (RyR2) has recently been shown to undergo remodeling in atrial fibrillation. These data suggest that the subcellular distribution of Ca²⁺-handling proteins influences proarrhythmic Ca²⁺-handling abnormalities. Here, we employ computational modeling to provide an in-depth analysis of the impact of variations in subcellular RyR2 and L-type Ca²⁺-channel distributions on Ca²⁺-transient properties and SCaEs in a human atrial cardiomyocyte model. We incorporate experimentally observed RyR2 expression patterns and various configurations of axial tubules in a previously published model of the human atrial cardiomyocyte. We identify an increased SCaE incidence for larger heterogeneity in RyR2 expression, in which SCaEs preferentially arise from regions of high local RyR2 expression. Furthermore, we show that the propagation of Ca²⁺ waves is modulated by the distance between RyR2 bands, as well as the presence of experimentally observed RyR2 clusters between bands near the lateral membranes. We also show that incorporation of axial tubules in various amounts and locations reduces Ca²⁺-transient time to peak. Furthermore, selective hyperphosphorylation of RyR2 around axial tubules increases the number of spontaneous waves. Finally, we present a novel model of the human atrial cardiomyocyte with physiological RyR2 and L-type Ca²⁺-channel distributions that reproduces experimentally observed Ca²⁺-handling properties. Taken together, these results significantly enhance our understanding of the structure-function relationship in cardiomyocytes, identifying that RyR2 and L-type Ca²⁺-channel distributions have a major impact on systolic Ca²⁺ transients and SCaEs.

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.

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.

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

In mammalian cardiac myocytes, the elementary calcium releases triggered by step voltage stimuli manifest either as solitary or as twin spikes that vary widely in kinetics and amplitude for unknown reasons. Here we examined the variability of calcium spikes measured using line-scanning confocal microscopy in patch-clamped rat ventricular myocytes. Amplitude distributions of the single and of the first of twin spikes were broader than those of the second spikes. All could be best approximated by a sum of a few elementary Gaussian probability distribution functions. The latency distributions of the single and the first spikes were identical, much shorter and less variable than those of the second spikes. The multimodal distribution of spike amplitudes and the probability of occurrence of twin spikes were stochastically congruent with activation of only a few of the many RyR2 channels present in the release site cluster. The occurrence of twin release events was rare due to refractoriness of release, induced with a probability proportional to the number of RyR2s activated in the primary release event. We conclude that the variability of the elementary calcium release events supports a calcium signalling mechanism that arises from stochastics of RyR2 gating and from inactivation of local origin.