The velocity of calcium waves is expected to depend non-monotoneously on the density of the calcium release units.

In this paper we develop a reaction-diffusion system describing the calcium dynamics in an agarose gel system with resuspended vesicles from the sarcoplasmic reticulum (SR vesicles). We focus on a simple model: compared with living cells (e.g. cardiac myocytes) an important property of the agarose gel system is the absence of the sarcolemma and the spatial separation of the calcium release units (CRUs). Our model includes the kinetics of ryanodine sensitive receptors (RyRs), the activity of the SERCA pumps and the diffusion of free calcium. We describe numerical simulations which show a biphasic relationship between the density of the CRUs and the propagation velocity of spreading waves. The non-monotony can be explained by changes in the amplitude of the local calcium concentration. We formulate implications for the in vitro system which could be verified in future experiments.