The Subcellular Distribution of Ryanodine Receptors and L-Type Ca²⁺ Channels Modulates Ca²⁺-Transient Properties and Spontaneous Ca²⁺-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.