] in the myoplasm or SR can influence the frequency, amplitude, and kinetics of local events. In ventricular myocytes, for instance, propagating Ca 2ϩ waves emerge when spontaneous Ca 2ϩ sparks trigger additional sparks in a regenerative fashion (5). Similarly, changes in channel gating at the level of individual RyRs can immediately affect the production of local events and, over time, influence bulk myoplasmic and SR [Ca 2ϩ ]. An increase in RyR opening will cause a gradual decrease in SR [Ca 2ϩ ] (46), whereas inhibition of RyR opening, over time, will lead to elevated SR [Ca 2ϩ ] (21). This principle has been elegantly demonstrated by Lukyanenko and coworkers (35) in quiescent rat ventricular myocytes treated with caffeine (to sensitize RyRs) or tetracaine (to inhibit RyRs).In heart cells, changes in Ca 2ϩ signaling due to altered RyR activity are currently receiving considerable attention because of close links to disease (13, 48). In particular, catecholaminergic polymorphic ventricular tachycardia (CPVT), an inherited disorder associated with a dramatic increase in arrhythmia risk, results from mutations in either RyRs or calsequestrin, a SR Ca 2ϩ buffer protein that associates with and modulates RyRs (20,33). Experiments in vitro have shown that CPVT-causing mutations usually increase the open probability of the RyR, resulting in a hyperactive or "leaky" channel (12,28,31). Studies (1, 36) have also suggested that leaky RyRs are characteristic of several experimental heart failure models. Thus, a quantitative understanding of how changes in RyR gating influence local and global Ca 2ϩ responses can provide insights into disease pathophysiology and can potentially suggest novel therapies.The difference in spatial scales between local and global Ca 2ϩ signals, however, creates significant challenges for the development of mechanistic mathematical models. In particular, gating of RyRs depends on both myoplasmic and SR [Ca 2ϩ ] (30), and concentrations within clusters during local events can be dramatically different from the bulk concentrations. In addition, because of the relatively small number of RyRs responsible for Ca 2ϩ sparks (6), the stochastic gating of these channels must be considered when simulating local events. Previous studies (9,27,42,44) have used Monte Carlo simulation methods to investigate the stochastic triggering of Ca 2ϩ sparks, but these have generally treated myoplasmic and bulk SR [Ca 2ϩ ] as fixed boundary conditions. Conversely, modeling studies (4, 10, 40) focusing on cellular Ca 2ϩ transients have usually used representations of SR Ca 2ϩ release that do not account for the stochastic nature of the local events. Attempting to simulate Ca 2ϩ signaling at both spatial scales simultaneously is a daunting prospect because the stochastic behavior of thousands of local events must be considered to determine the effects on the bulk concentrations. As a result, only a few studies (16, 17, 49, 50) have attempted to capture both phenomena.