Calcium homeostasis in a local/global whole cell model of permeabilized ventricular myocytes with a Langevin description of stochastic calcium release

Population density approaches to modeling local control of Ca(2+)-induced Ca(2+) release in cardiac myocytes can be used to construct minimal whole cell models that accurately represent heterogeneous local Ca(2+) signals. Unfortunately, the computational complexity of such "local/global" whole cell models scales with the number of Ca(2+) release unit (CaRU) states, which is a rapidly increasing function of the number of ryanodine receptors (RyRs) per CaRU. Here we present an alternative approach based on a Langevin description of the collective gating of RyRs coupled by local Ca(2+) concentration ([Ca(2+)]). The computational efficiency of this approach no longer depends on the number of RyRs per CaRU. When the RyR model is minimal, Langevin equations may be replaced by a single Fokker-Planck equation, yielding an extremely compact and efficient local/global whole cell model that reproduces and helps interpret recent experiments that investigate Ca(2+) homeostasis in permeabilized ventricular myocytes. Our calculations show that elevated myoplasmic [Ca(2+)] promotes elevated network sarcoplasmic reticulum (SR) [Ca(2+)] via SR Ca(2+)-ATPase-mediated Ca(2+) uptake. However, elevated myoplasmic [Ca(2+)] may also activate RyRs and promote stochastic SR Ca(2+) release, which can in turn decrease SR [Ca(2+)]. Increasing myoplasmic [Ca(2+)] results in an exponential increase in spark-mediated release and a linear increase in nonspark-mediated release, consistent with recent experiments. The model exhibits two steady-state release fluxes for the same network SR [Ca(2+)] depending on whether myoplasmic [Ca(2+)] is low or high. In the later case, spontaneous release decreases SR [Ca(2+)] in a manner that maintains robust Ca(2+) sparks.