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 Ca²⁺ and Mg²⁺ 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.

Temperature influences many aspects of cardiac excitation-contraction coupling, in particular, hypothermia increases the open probability (P_open) of cardiac sarcoplasmic reticulum (SR) Ca²⁺-release channels (ryanodine-sensitive RyR channels) rising the SR Ca²⁺ load in mammalian myocytes. However, to the best of our knowledge, no theoretical models are available for that effect. Traditional Markov chain models do not provide a reasonable molecular mechanistic insight on the origin of the temperature effects. Here in the paper we address a simple physically clear electron-conformational model to describe the RyR gating and argue that a synergetic effect of external thermal fluctuation forces (Gaussian–Markovian noise) and internal friction via the temperature stimulation/suppression of the open–close RyR tunneling probability can be considered as a main contributor to temperature effects on the RyR gating. Results of the computer modeling allowed us to successfully reproduce all the temperature effects observed for an isolated RyR gating in vitro under reducing the temperature: increase in P_open and mean open time without any significant effect on mean closed