A Ca(2+) spark arises when a cluster of sarcoplasmic reticulum (SR) channels (ryanodine receptors or RyRs) opens to release calcium in a locally regenerative manner. Normally triggered by Ca(2+) influx across the sarcolemmal or transverse tubule membrane neighboring the cluster, the Ca(2+) spark has been shown to be the elementary Ca(2+) signaling event of excitation-contraction coupling in heart muscle. However, the question of how the Ca(2+) spark terminates remains a central, unresolved issue. Here we present a new model, "sticky cluster," of SR Ca(2+) release that simulates Ca(2+) spark behavior and enables robust Ca(2+) spark termination. Two newly documented features of RyR behavior have been incorporated in this otherwise simple model: "coupled gating" and an opening rate that depends on SR lumenal [Ca(2+)]. Using a Monte Carlo method, local Ca(2+)-induced Ca(2+) release from clusters containing between 10 and 100 RyRs is modeled. After release is triggered, Ca(2+) flux from RyRs diffuses into the cytosol and binds to intracellular buffers and the fluorescent Ca(2+) indicator fluo-3 to produce the model Ca(2+) spark. Ca(2+) sparks generated by the sticky cluster model resemble those observed experimentally, and Ca(2+) spark duration and amplitude are largely insensitive to the number of RyRs in a cluster. As expected from heart cell investigation, the spontaneous Ca(2+) spark rate in the model increases with elevated cytosolic or SR lumenal [Ca(2+)]. Furthermore, reduction of RyR coupling leads to prolonged model Ca(2+) sparks just as treatment with FK506 lengthens Ca(2+) sparks in heart cells. This new model of Ca(2+) spark behavior provides a "proof of principle" test of a new hypothesis for Ca(2+) spark termination and reproduces critical features of Ca(2+) sparks observed experimentally.

We construct a detailed mathematical model for Ca2+ regulation in the ventricular myocyte that includes novel descriptions of subcellular mechanisms based on recent experimental findings: 1) the Keizer-Levine model for the ryanodine receptor (RyR), which displays adaptation at elevated Ca2+; 2) a model for the L-type Ca2+ channel that inactivates by mode switching; and 3) a restricted subspace into which the RyRs and L-type Ca2+ channels empty and interact via Ca2+. We add membrane currents from the Luo-Rudy Phase II ventricular cell model to our description of Ca2+ handling to formulate a new model for ventricular action potentials and Ca2+ regulation. The model can simulate Ca2+ transients during an action potential similar to those seen experimentally. The subspace [Ca2+] rises more rapidly and reaches a higher level (10-30 microM) than the bulk myoplasmic Ca2+ (peak [Ca2+]i approximately 1 microM). Termination of sarcoplasmic reticulum (SR) Ca2+ release is predominately due to emptying of the SR, but is influenced by RyR adaptation. Because force generation is roughly proportional to peak myoplasmic Ca2+, we use [Ca2+]i in the model to explore the effects of pacing rate on force generation. The model reproduces transitions seen in force generation due to changes in pacing that cannot be simulated by previous models. Simulation of such complex phenomena requires an interplay of both RyR adaptation and the degree of SR Ca2+ loading. This model, therefore, shows improved behavior over existing models that lack detailed descriptions of subcellular Ca2+ regulatory mechanisms.

Rationale Transverse tubules (TTs) couple electric surface signals to remote intracellular Ca2+ release units (CRUs). Diffraction-limited imaging studies have proposed loss of TT components as disease mechanism in heart failure (HF). Objectives Objectives were to develop quantitative super-resolution strategies for live-cell imaging of TT membranes in intact cardiomyocytes and to show that TT structures are progressively remodeled during HF development, causing early CRU dysfunction. Methods and Results Using stimulated emission depletion (STED) microscopy, we characterized individual TTs with nanometric resolution as direct readout of local membrane morphology 4 and 8 weeks after myocardial infarction (4pMI and 8pMI). Both individual and network TT properties were investigated by quantitative image analysis. The mean area of TT cross sections increased progressively from 4pMI to 8pMI. Unexpectedly, intact TT networks showed differential changes. Longitudinal and oblique TTs were significantly increased at 4pMI, whereas transversal components appeared decreased. Expression of TT-associated proteins junctophilin-2 and caveolin-3 was significantly changed, correlating with network component remodeling. Computational modeling of spatial changes in HF through heterogeneous TT reorganization and RyR2 orphaning (5000 of 20 000 CRUs) uncovered a local mechanism of delayed subcellular Ca2+ release and action potential prolongation. Conclusions This study introduces STED nanoscopy for live mapping of TT membrane structures. During early HF development, the local TT morphology and associated proteins were significantly altered, leading to differential network remodeling and Ca2+ release dyssynchrony. Our data suggest that TT remodeling during HF development involves proliferative membrane changes, early excitation-contraction uncoupling, and network fracturing.