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

Phosphorylation of RyR2 simultaneously expands the dyad and rearranges the tetramers

We have previously demonstrated that type II ryanodine receptors (RyR2) tetramers can be rapidly rearranged in response to a phosphorylation cocktail. The cocktail modified downstream targets indiscriminately, making it impossible to determine whether phosphorylation of RyR2 was an essential element of the response. Here, we used the β-agonist isoproterenol and mice homozygous for one of the following clinically relevant mutations: S2030A, S2808A, S2814A, or S2814D. We measured the length of the dyad using transmission electron microscopy (TEM) and directly visualized RyR2 distribution using dual-tilt electron tomography. We found that the S2814D mutation, by itself, significantly expanded the dyad and reorganized the tetramers, suggesting a direct link between the phosphorylation state of the tetramer and its microarchitecture. S2808A and S2814A mutant mice, as well as wild types, had significant expansions of their dyads in response to isoproterenol, while S2030A mutants did not. In agreement with functional data from these mutants, S2030 and S2808 were necessary for a complete β-adrenergic response, unlike S2814 mutants. Additionally, all mutants had unique effects on the organization of their tetramer arrays. Lastly, the correlation of structural with functional changes suggests that tetramer-tetramer contacts play an important functional role. We thus conclude that both the size of the dyad and the arrangement of the tetramers are linked to the state of the channel tetramer and can be dynamically altered by a β-adrenergic receptor agonist.

PHOSPHORYLATION OF RyR2 SIMULTANEOUSLY EXPANDS THE DYAD AND REARRANGES THE TETRAMERS

We have previously demonstrated that type II ryanodine receptors (RyR2) tetramers can be rapidly rearranged in response to a phosphorylation cocktail. The cocktail modified downstream targets indiscriminately making it impossible to determine whether phosphorylation of RyR2 was an essential element of the response. We therefore used the β-agonist isoproterenol and mice with one of the homozygous mutations, S2030A +/+ , S2808A +/+ , S2814A +/+ , or S2814D +/+ , to address this question and to elucidate the role of these clinically relevant mutations. We measured the length of the dyad using transmission electron microscopy (TEM) and directly visualized RyR2 distribution using dual-tilt electron tomography. We found that: 1) The S2814D mutation, by itself, significantly expanded the dyad and reorganized the tetramers suggesting a direct link between the phosphorylation state of the tetramer and the microarchitecture. 2) All of the wild-type, as well as the S2808A and S2814A mice, had significant expansions of their dyads in response to ISO, while S2030A did not. 3) In agreement with functional data from the same mutants, S2030 and S2808 were necessary for a complete β-adrenergic response, whereas S2814 was not. 4) All the mutated residues had unique effects on the organization of their tetramer arrays. 5) The correlation of structure with function suggests that tetramer-tetramer contacts play an important functional role. We conclude that both the size of the dyad and the arrangement of the tetramers are linked to the state of the channel tetramer and can be dynamically altered by a β-adrenergic receptor agonist.

Summary

Analysis of RyR2 mutants suggests a direct link between the phosphorylation state of the channel tetramer and the microarchitecture of the dyad. All phosphorylation site mutations produced significant and unique effects on the structure of the dyad and its response to isoproterenol.

Calcium oscillation on homogeneous and heterogeneous networks of ryanodine receptor

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

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.

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

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