Simulating cardiac Ca2+ release units: effects of RyR cluster size and Ca2+ buffers on diastolic Ca2+ leak

Leak of Ca2+ out of the cardiac sarcoplasmic reticulum (SR) via ryanodine receptors (RyRs) during diastole is vital to regulate SR Ca2+ levels. This leak can become deleterious when large spontaneous RyR-mediated Ca2+ release events evoke proarrhythmic Ca2+ waves that can lead to delayed after-depolarizations. Here, we model diastolic SR Ca2+ leak at individual SR Ca2+ 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 Ca2+ leak, with bigger arrays producing larger and more frequent Ca2+ release events. Moreover, big RyR arrays are more susceptible to small changes in the levels of dyadic Ca2+ buffers. Such changes in buffering shift Ca2+ leak from small Ca2+ 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 Ca2+ loads, we find further evidence for the hypothesis that SR Ca2+ leak by RyR arrays can undergo a sudden phase transition.