The RyR2 channel is a tetramer composed of four identical subunits that gate a shared central pore. Structural studies show that these subunits are tightly packed and physically interact at their interfaces (Cabra et al., 2016; Woll and Van Petegem, 2022). Motivated by these observations and by our earlier modeling work (Greene et al., 2023), we describe a simplified RyR2 model in which each subunit can exist in one of two conformational states. A central assumption of the model is that adjacent subunits interact energetically, such that conformational mismatches are penalized. In particular, if a subunit transitions from closed to open while its neighbors remain closed, the resulting mismatch produces an energetic cost. This is implemented by modifying the transition rate so that the opening rate of a closed subunit is reduced by a multiplicative factor when neighbors are closed and increased when neighbors are open. An additional assumption is that a channel is considered functionally open when three or more subunits are in the open state. This requirement was established in our previous study (Greene et al., 2023), where we demonstrated that RyR2 channels must exhibit cooperative gating to remain reliably shut at diastolic Ca concentrations. Without this constraint, channels cannot maintain stable closure during rest.

To apply this model, we represent the RyR2 channel as four subunits labeled s i , with i = 1 , 2 , 3 , 4 , and periodic boundary conditions. Each subunit can be in one of two states: closed s i = − 1 or open s i = + 1 . This configuration ensures that every subunit interacts with its two immediate neighbors, consistent with the geometry of the RyR2 channel. Transitions between subunit states are governed by rates that depend on the conformational states of neighboring subunits. The forward (opening) and backward (closing) transition rates for subunit i are given by k i − 1 → + 1 = k f · γ i ,   k i + 1 → − 1 = k b γ i − 1 , (1) where k f and k b are the intrinsic forward and backward rates in the absence of coupling. The factor γ i captures the influence of neighboring subunits and is defined as: γ i = exp δ · n i ,   n i = s i − 1 + s i + 1   . (2)

Here, δ = β J is a dimensionless parameter representing the strength of coupling between channels, and where J is the energy penalty. Also, β = 1 / k B T where k B is Boltzman’s constant, and T is the temperature. These local interaction rules form the basis of a thermodynamically consistent model that accounts for cooperative interactions within individual RyR2 channels.

Recent Cryo-EM imaging studies due to Cabra et al. (2016) reveals that RyR2 channels physically interact at their interfaces through specific molecular contacts between adjacent channels. Their high-resolution structural analysis identifies two distinct types of inter-channel arrangements: an “adjoining” configuration (Figure 1A) and an “oblique” configuration (Figure 1B). In the adjoining arrangement, channels are aligned in regular rows and columns, with horizontal neighbors coupling through their edge-to-edge contacts. In the oblique arrangement (Figure 1B), channels adopt a staggered geometry with interactions between subunits 1–3 and 2–4 only. Cabra et al. observed that native RyR2 clusters typically form disordered structures that represent combinations of these two geometric arrangements rather than uniform arrays of a single type. To model these experimentally observed arrangements, we will consider cluster models based on both interaction types.

To model heterogeneous arrays, we denote s i c as the state of subunit i in channel c . Similarly we will denote s j c ′ and s k c ″ to be the state of subunits j and k in neighboring channels c ′ and c ″ respectively. The coupling factor for subunit i in channel c is given by γ i = exp δ s i − 1 c + s i + 1 c + σ ′ s j c ′ + σ ″ s k c ″ , (3) where σ ′ and σ ″ represents the subunit interaction strength between subunit i in channel c and subunits j and k in neighboring channels c ′ and c ″ . For simplicity, we will assume σ ′ = σ ″ = σ in this study. However, it is likely that in realistic scenarios these parameters may differ due to variations in channel properties. With this formulation the transition rates have the same form as in Equation 1 and remain consistent with the principle of detailed balance.

To model inter-RyR2 interactions within a cluster, we arrange RyR2 channels either with the oblique or adjoining configurations shown in Figure 1. Each channel occupies a grid position labeled by coordinates n x , n y , and consists of four subunits labeled s i . Each subunit is uniquely identified by the triplet n x , n y , s i , which denotes the lattice position and the state of the tetramer. The model incorporates two distinct types of coupling interactions. First, intra-channel coupling occurs between adjacent subunits within the same RyR2 tetramer, as described in Equation 2. Second, inter-channel coupling links specific subunits from neighbouring channels according to Equation 3. In the adjoining arrangement, subunits form connections with two nearest-neighbour channels, reflecting the side-by-side contact geometry. In the oblique arrangement, each subunit couples to only one nearest neighbour, corresponding to the staggered interaction pattern. To capture more realistic cluster arrangements, we will also consider disordered configurations that reflect the structural heterogeneity observed in cardiac myocytes. This computational framework will allow us to systematically investigate how the loss of structural integrity affects Ca release dynamics.

Super-resolution imaging studies reveal that RyR2 cluster sizes span a wide range. The most recent work by Hou et al. (2015) reports many small clusters containing roughly 20–50 channels (spanning approximately 150–250 nm in diameter), while earlier studies by Galice et al. (2018) observed larger aggregates that could extend to several hundred RyR2s, with typical averages closer to 50–100 channels (spanning approximately 250–350 nm in diameter). Although the absolute estimates vary across techniques, the consistent finding is that ventricular myocytes contain a heterogeneous population of clusters, ranging from solitary channels and small groups to large assemblies. The functional behaviour of these clusters is strongly shaped by Ca diffusion in the narrow dyadic cleft. Ca diffuses rapidly in the cytosol, with an effective diffusion coefficient in the range D ≈ 100 − 500 μ m 2 / s (Bers and Peskoff, 1991; Swietach et al., 2010). This wide range is due to the presence of multiple buffers and the structural complexity of the intracellular space. Given an RyR2 open time of roughly τ ∼ 1 m s , the expected diffusion length ( d = 4 D τ ) is in the range 600 − 1000 n m . This length scale exceeds the size of even the largest RyR2 clusters ≤ 400 n m , implying that Ca released from a single open channel rapidly equilibrates across the entire cluster. As a result, RyR2s within a cluster are functionally coupled through the shared local Ca signal, independently of direct molecular contact.

Since diffusion is much faster than channel gating, we invoke the rapid diffusion approximation, which assumes Ca is uniform within the dyad. Under this approximation, the local Ca concentration is given by C a = c o + g · n o , (4) where c 0 is the diastolic Ca concentration in the bulk cytosol, n o represents the number of open RyR2 channels within the cluster, and g represents the Ca concentration increase in the dyadic cleft due to a single open RyR2 channel (Bers and Peskoff, 1991; Langer and Peskoff, 1996; Asfaw et al., 2013). Several computational studies of CICR have argued that [Ca] must reach relatively high concentrations in the dyadic junction. Cannell et al. (2013) argued that minimum dyadic [Ca] of ∼ 50 μ M was required to maintain CICR in rat ventricular myocytes, while Shen et al. (2022) used half-maximal RyR2 activation concentrations ( K d ) ranging from 25 − 45 μ M in their models of RyR cluster dispersion during β-adrenergic stimulation. In this study we will set g = 25 μ M , so that 1-3 open channels produce dyadic [Ca] in the range of 25–75 μM, consistent with the concentrations required to support CICR in these studies.

To determine the transition rates in our model, we assume first-order kinetics with respect to cytosolic Ca. Specifically, the forward rate for subunit opening is given by k f = k f o C a , where k f o is a binding rate constant and C a is the local Ca concentration. We set k f o = 0.005   μ M · m s − 1 to ensure that RyR2 channels remain predominantly closed at a diastolic concentratonof C a ∼ 0.1 μ M . This choice is consistent with experimental observations, see for example Shen et al. (2022), showing that spontaneous Ca sparks are extremely rare at baseline ( ∼ 0.1 s p a r k s / 100 μ m / s ). Although this forward rate is substantially lower than binding rates measured in bilayer recordings (Mukherjee et al., 2012; Xu and Meissner, 2004), the difference reflects the highly regulated environment of intact cells. In particular, physiological M g 2 + occupies the Ca-activation sites at resting [Ca] and must be displaced by Ca to allow channel opening, effectively reducing the Ca sensitivity and stabilizing the closed state compared to isolated channels in bilayers (Laver et al., 1997; Walweel and Laver, 2015; Zahradníková et al., 2025). The backward rate for subunit closing is Ca-independent and fixed at k b = 0.3   m s − 1 (Xu and Meissner, 2004), which corresponds to an average open time of approximately 3 m s . Subunit cooperativity is captured by a dimensionless mismatch penalty δ = 0.5 , which stabilizes the closed conformation at rest without imposing such a large energetic cost that the channel becomes unresponsive at elevated Ca. A pore conducts only when at least three of the four subunits are open, a majority rule that preserves observed closed-time statistics and prevents spurious openings (Greene et al., 2023). Collectively, these parameters produce physiologically plausible single-channel kinetics and cluster-level behavior across the physiological Ca range.

The Gillespie algorithm is implemented as a stochastic simulation method to model the exact temporal evolution of RyR2 channel gating within clusters (Gillespie, 2007). At each time step, the algorithm calculates the total reaction rate across all subunits in the system based on their current states and local coupling environments. The time to the next reaction is drawn from an exponential distribution with parameter equal to the total rate, ensuring proper stochastic timing. A specific reaction is then selected using weighted random selection proportional to individual subunit rates, after which the chosen subunit transitions between closed and open states. This process continues iteratively, with the system state analyzed after each reaction to count the number of open channels.

To study the stochastic dynamics of RyR2 clusters, we first simulated the time evolution of a 10 × 10 array of RyR2 channels using our stochastic simulation algorithm. To monitor the cluster, we keep track of the number of open RyR2 channels, denoted as n o , which is defined as the number of channels where 3 or more subunits are in the open state. Figure 2 shows n o as function of time for three simulation runs where the cytosolic Ca concentration is fixed at c o = 6.0 μ M . Each independent simulation run is shown as a red, black, and blue line. In these runs the cluster stays mostly closed and occasionally one or two channels transition to the open state. After a waiting period a large fluctuation pushes n o above a critical threshold, leading to a rapid cooperative avtivation of the entire cluster. Once this transition occurs, the number of open channels grows explosively from a few n o ∼ 0 − 1 to nearly all the channels in the cluster n o ∼ 100 . This behaviour arises because the local Ca concentration is dependent on n o according to Equation 4. Thus, the opening rate of channels increases in a strongly nonlinear fashion due to subunit cooperativity within and between RyR2 channels.

The rapid increase in n o corresponds to the initiation of a Ca spark. Since the excitation occurred because of stochastic fluctuations, we will refer to these events as spontaneous Ca sparks. After activation, RyR2 channels proceed to inactivate and return to the closed state. However, this inactivation process occurs on a much shorter timescale than the waiting time to a spontaneous Ca spark and is not modelled in this study. To quantify the spontaneous Ca spark initiation process, we define the activation time T , shown in Figure 2 for the simulation run corresponding to the black line, as the time required for the cluster to undergo the initial transition from the closed state to a fully activated cluster. Specifically, T is measured from t = 0 , when all RyR2 channels in the cluster are initially closed (all subunits are in the state s i = − 1 ), to the moment when n o exceeds a critical threshold n c . In general, n c should be picked to be larger than a typical fluctuation, and in this study we will use n c = 5 . However, we note here that our results are independent of n c providing n c is larger than the baseline fluctuations. For example, in the simulation shown in Figure 2, any n c ≥ 5 gives the same waiting time since once n c = 5 is reached then excitation proceeds with high likelihood. This approach provides a consistent measure of the waiting time before spontaneous Ca spark initiation, which determines the average frequency of spontaneous Ca sparks in the cell.

We next investigated how the mean waiting time T to a spontaneous Ca spark depends on the dyadic Ca concentration c o . This concentration is typically c o ∼ 0.1 μ M during rest, but will increase substantially when there is a local opening of a nearby LCC channel. Figure 3A shows the relationship between T and c o for the oblique configuration of the cluster. Here, we vary c o from 5.0 μ M to 120 μ M , which represents the typical Ca concentrations at a dyadic junction that occurs during LCC activation. Three values of the inter-subunit coupling parameter σ are examined: σ = 0 (no coupling), σ = 0.5 (moderate coupling), and σ = 1 (strong coupling). The results show that T changes exponentially with the concentration c o . At high diastolic Ca ( c o ∼ 100 μ M ) all three cases have comparable firing times in the range 0.5 m s − 2 m s . In contrast, at lower Ca concentrations the waiting time increases exponentially with decreasing c o . In particular, when subunit interactions between channels is engaged, with σ = 1.0 , the waiting time increases from ∼ 2 m s at c o ∼ 120 μ M , to ∼ 10 5 m s at c o ∼ 5.0 μ M . Figure 3B shows analogous results for the adjoining configuration. The qualitative features are the same as for the oblique case, but the dependence on σ is more pronounced, since each subunit has two contacts rather than just one. In Figures 3C,D we compute T , for both the oblique and adjoining configurations respectively, for cluster sizes of 5 × 5 , 10 × 10 , and 20 × 20 . These results shows that the dependence of T on c o is similar across the wide range of cluster sizes that are expected to be found in cardiac myocytes.

To further explore the role of coupling, we examined the dependence of T on σ at two fixed Ca concentrations. Figure 4A shows results at a resting concentration of c o = 5.0 μ M . In this regime, the mean waiting time measured for both the oblique and adjoining configurations grows almost three orders of magnitude as σ is increased. In contrast, at a diastolic Ca concentration of c o = 100 μ M , shown in Figure 4B, T increases only modestly over a few milliseconds. Together, these findings identify inter-channel coupling as a powerful mechanism for tuning cluster stability at low Ca concentrations. At high Ca, clusters activate rapidly and reliably, with little dependence on coupling strength. At low Ca, however, even modest increases in σ dramatically stabilize the closed state by exponentially prolonging the waiting time T . This dual behaviour ensures that RyR2 clusters remain quiescent at rest, while still capable of fast activation when triggered during excitation–contraction coupling.

Experimental studies have revealed that RyR2 clusters undergo fractionation and lose their structural integrity in diseased states. To investigate the functional consequences of cluster fractionation, we developed a computational model that simulates the progressive disruption of RyR2 cluster connectivity. Our approach implements stochastic bond breaking between adjacent RyR2 channels within the cluster. Specifically, we assign a probability p that a connection between two nearest-neighbor subunits remains intact. The parameter p serves as a measure of cluster integrity, where p = 1 corresponds to a fully connected cluster with all inter-channel connections preserved, while p = 0 represents a completely dissociated cluster with all connections broken. For each value of p , we generate an ensemble of connectivity matrices that define the network topology and compute the average T over the ensemble. In Figure 5, we show the relationship between cluster connectivity and spontaneous activation by plotting T versus the connection probability p at a Ca concentration of c o = 2.0 μ M . The results show that the waiting time T grows exponentially with increasing connection probability p . In particular, for the case of the adjoining configuration, we find that an uncoupled cluster of 10 × 10 RyR2s has a mean firing time of T ∼ 10 3 m s . However, as p is increased to just p ∼ 0.5 , the waiting time rises to T ∼ 10 6 m s . This, result shows that the spontaneous Ca spark frequency is exponentially sensitive to the probability that neighboring RyR2s are connected.

Experimental studies reveal that RyR2 clusters within the junctional sarcoplasmic reticulum exhibit highly heterogeneous spatial arrangements. To investigate how this structural diversity affects functional properties, we have also developed a preferential attachment growth model that generates realistic cluster morphologies. The preferential attachment algorithm sequentially places RyR2 channels on a two-dimensional lattice. For each new channel placement, the probability of selecting an empty site i , j is proportional to 1 + n e α , where n e represents the number of occupied nearest-neighbor sites and the exponent α is referred to as the clustering parameter. This formulation ensures that sites adjacent to existing channels are preferentially selected, with the strength of this preference controlled by α . When α = 0 , channel placement is uniform and random, while increasing values promote increasingly compact cluster formation.

Figure 6A demonstrates the range of cluster structures generated by varying the clustering parameter α . At α = 0 , RyR2 channels are scattered randomly on the grid. At intermediate values α = 5 more cohesive structures start to from, while retaining some fragmentation. A high clustering parameter α = 8 generates compact arrangements with a few dispersed channels. To examine how cluster structural heterogeneity affects spontaneous Ca spark dynamics, we applied the preferential attachment algorithm by placing 50 RyR2 channels on a 20 × 20 lattice domain. Given that individual RyR2 channels are approximately 30 n m × 30 n m in size (Cabra et al., 2016), this domain corresponds to approximately 600 n m in diameter, roughly the spatial extent over which Ca is strongly diffusively coupled. For each value of the clustering parameter α , we generated 50 independent cluster configurations. Once the configuration is generated the channels are coupled according to the oblique architecture. We then average the mean first passage time T over all configurations. Figure 6B shows T vs. α at a fixed c o = 2.5 μ M . The results show that T increases substantially as α is varied in the range 0–10. At low values ( α ≈ 0 ), corresponding to dispersed clusters with weak inter-channel connectivity, the mean waiting time is approximately 5 × 10 3 m s . As α increases, producing tightly organized clusters, T increases to approximately 10 5 m s for the oblique configuration, and 10 6 for the adjoining configuration. This effect saturates above α ≈ 7 since all clusters generated for larger α are close to fully compact and yield the same T as the fully connected square lattice. This exponential relationship, in the range 0–7, demonstrates that cluster compactness is a critical determinant of spontaneous Ca release dynamics, with tightly organized structures providing exponentially greater stability during diastole.

In this study, we developed a computational model of an RyR2 cluster in which each channel is composed of four interacting subunits, with channels arranged in a configuration consistent with structures observed in cryo-EM and super-resolution imaging. Our goal was to understand how inter-subunit cooperativity influences Ca leak in cardiac cells. The first result of the model is that the cluster exhibits a sharp contrast in behavior across different Ca concentration regimes. At Ca concentrations larger than 100 μ M RyR2 clusters respond rapidly to Ca signals on sub-millisecond timescales. This fast response is largely insensitive to the strength of channel-to-channel coupling. However, at lower Ca concentrations (below ∼ 5 μ M ), the system behaves very differently. In this case, the mean waiting time to a spontaneous Ca spark becomes exponentially sensitive to the strength of subunit coupling. Specifically, at a Ca concentration of 5.0 μ M , increasing the coupling strength σ between RyR2 subunits from 0 to 0.5 (for the adjoining arrangement), leads to an almost 4 orders of magnitude increase in the mean waiting time. This behavior highlights the functional significance of inter-subunit coupling: it preserves fast responsiveness to LCC openings while ensuring cluster stability at lower Ca concentrations. An important consequence of this finding is that diastolic Ca leak, which is determined by how often spontaneous Ca sparks occur, will be exponentially sensitive to the coupling between RyR2 within clusters. Therefore, even small changes in coupling strength can produce orders-of-magnitude differences in leak rate. These results highlights the critical role of RyR2 coupling in regulating Ca cycling homeostasis.

Experimental evidence from cryo-EM and super-resolution microscopy shows that RyR2 clusters in cardiac cells are highly irregular, forming heterogeneous structures rather than perfect lattices. To examine how this heterogeneity shapes Ca signaling, we extended our modeling framework to include disordered cluster geometries and systematically varied their structural integrity. In one approach, we introduced a coupling parameter p , defined as the probability that a subunit remains linked to its nearest neighbor. This parameter represents the progressive loss of inter-subunit connectivity. As p decreased from 1 (fully intact) to 0 (fully fragmented), the mean waiting time to a spontaneous Ca spark shortened by approximately 3 orders of magnitude. These results shows that Ca spark frequency is exponentially sensitive to even small disruptions in cluster connectivity. In parallel, we used a preferential attachment algorithm to generate clusters with heterogeneous shapes resembling those observed experimentally. Here, the clustering parameter α controlled the likelihood that new channels attached to existing groups, thereby tuning the degree of structural cohesion. Increasing α produced a substantial increase in spark waiting time, demonstrating that more tightly connected clusters are exponentially more stable at lower Ca concentrations. Together, these results reveal that heterogeneity in RyR2 cluster organization plays a central role in regulating the timing of spontaneous Ca sparks.

Recent high-resolution imaging studies have revealed that RyR2 clusters lose their structural integrity in heart failure, with important functional consequences for Ca handling. Sheard et al. (2022) used enhanced expansion microscopy to demonstrate that RyR2 clusters in failing myocytes frequently exhibit a frayed appearance, where small groups of channels detach from the main cluster, particularly in regions depleted of the structural protein junctophilin-2. Similarly, Kolstad et al. (2018) showed that post-infarction heart failure is characterized by RyR2 cluster dispersion and fragmentation, resulting in smaller, more numerous cluster fragments with reduced inter-channel connectivity. Our computational modelling provides a mechanistic explanation for why these structural changes change the Ca leak rate. The key insight from our model is that spontaneous Ca spark frequency exhibits exponential sensitivity to cluster structural integrity. Mechanistically, intact inter-subunit coupling within the cluster provides stabilizing interactions that maintain channels in the closed state during diastole through cooperative inhibition. When coupling is weakened by fragmentation, this stabilizing effect is lost, making the cluster far more susceptible to stochastic activation. The key new insight in this study is that this relationship is exponential, and explains why even modest cluster remodeling, such as partial fragmentation or fraying, can produce the orders of magnitude increases in diastolic Ca leak observed in heart failure. Thus, cluster architectural integrity emerges as a critical control mechanism that drives Ca leak.

Ca diffusion plays a crucial role in determining the functional coupling between RyR2 channels within the dyadic cleft. With a diffusion coefficient of approximately 100 − 500   μ m 2 / s and an RyR2 open time of roughly 1 m s , Ca diffuses roughly 600 n m before channel closure. Since this diffusion length exceeds the diameter of even the largest RyR2 clusters ( ≤ 400 n m ), Ca released from any open channel rapidly equilibrates across the entire dyadic space, creating strong functional coupling through the shared Ca microdomain. Importantly, when clusters become frayed or fragmented into separate pieces, the disjoint subclusters often remain within the same Ca microdomain and are therefore still functionally coupled through diffusion. However, these separated channels lose the stabilizing nearest-neighbor interactions that maintain closure during diastole. A small detached subcluster can thus spontaneously open without the energetic penalty of intact inter-channel coupling yet still elevate local Ca sufficiently to recruit the larger cluster. This asymmetry between lost stabilization and preserved diffusive coupling explains why frayed clusters exhibit exponentially faster spontaneous Ca spark frequencies.

The exponential sensitivity of Ca leak to RyR2 coupling revealed in our study carries significant therapeutic implications. Because spark frequency rises exponentially as coupling weakens, even subtle structural disruptions can translate into exponentially large increases in spontaneous Ca release frequency. This nonlinearity also works in the opposite direction, so that modest improvements in cluster integrity could produce dramatic reductions in Ca leak. These results shows that structural stability itself is a valuable therapeutic target. Strategies that preserve or restore RyR2 subunit coupling, such as enhancing the stabilizing action of accessory proteins like FKBP12.6 (Wehrens et al., 2003), offer a direct way to reinforce this stability. Other proteins that regulate RyR2 positioning or inter-channel contacts could likewise be harnessed to maintain cluster integrity (Reynol et al., 2016). Such approaches differ fundamentally from traditional interventions aimed at modifying channel gating kinetics or expression levels, which may have limited efficacy if the underlying cluster architecture remains compromised. By identifying inter-subunit coupling as a central determinant of spontaneous Ca release, our study highlights a specific and tunable structural feature of the RyR2 complex that could be targeted to suppress pathological Ca leak.

Our model has several important limitations. First, we do not incorporate the complex regulatory mechanisms known to modulate RyR2 gating in cardiac myocytes, including inhibition by physiological Mg concentrations and modulation by regulatory proteins such as CaMKII, calmodulin, and FKBP12.6 (Marx and Marks, 2013; Walweel and Laver, 2015; Zahradníková et al., 2025; Rokita and Anderson, 2012). These factors will alter the spark waiting times. However, they are unlikely to eliminate the exponential dependence on cluster connectivity that emerges from cooperative gating between channels. Second, we model RyR2 clusters using idealized geometric arrangements of either pure oblique or pure adjoining configurations, whereas native clusters in cardiac cells exhibit heterogeneous combinations of both interaction types. Our results indicate that mixed geometries will yield intermediate waiting times, with the precise value depending on the relative proportion of each interaction type. However, because both pure configurations exhibit exponential sensitivity to coupling strength and connectivity, we expect this relationship to persist regardless of geometric details. Third, our framework does not include channel inactivation, SR Ca depletion, or competing Ca clearance mechanisms such as the sodium calcium exchanger and the sarcoplasmic reticulum Ca ATPase (SERCA). These processes determine spark termination and recovery. However, they occur on timescales much shorter than the diastolic waiting times we investigate and therefore do not affect the spark initiation dynamics that are the focus of this study. While future work incorporating more detailed RyR2 kinetics and realistic cluster geometries will refine quantitative predictions, the central finding that spontaneous Ca spark frequency is exponentially sensitive to cluster structural integrity follows directly from the cooperativity between RyR2 channels and should remain robust.