 Finite calcium-store release model

For an original mathematical model, imagine a sustained, imposed release pathway, with no store refill during the plotted interval. Let $S$ be releasable store amount of substance, $X$ excess freely dissolved cytosolic calcium amount of substance, and $Q$ the cumulative amount removed from $X$ into a separate sink. Choose $S(0)=1\mathrm{amol}$, $X(0)=Q(0)=0$, cytosolic volume $V_c=1\mathrm{pL}$, release rate constant $k=0.4{\mathrm{s}}^{-1}$ and clearance rate constant $h=0.8{\mathrm{s}}^{-1}$. These are stipulated values, not physiological measurements or a model of IP3 receptor gating.

Assume fixed volume, uniform mixing, no binding, no reverse transfer, no outside entry and first-order loss from each donating compartment. ThenThe last equality is conservation of mass for the tracked calcium across the store, mobile excess and sink together; it does not assert conservation inside a biological cell if the sink represents export. A constant background $c_b=0.1\mu \mathrm{M}$ is held separate from this excess bookkeeping. Since $1\mathrm{amol}/1\mathrm{pL}=1\mu \mathrm{M}$, the response is $c=c_b+X/V_c$.

Solving these ordinary differential equations givesRelease $kS$ declines from the start, yet $X$ initially rises because delivery exceeds clearance. Its peak occurs when $kS=hX$, at $t_{\ast }=\log (h/k)/(h-k)=\log (2)/0.4\approx 1.73\mathrm{s}$. There $X=0.25\mathrm{amol}$ and $c=0.35\mu \mathrm{M}$. Ultimately $S$ and $X$ tend to zero, $Q$ tends to $1\mathrm{amol}$ and $c$ returns to $c_b$. The delayed peak is a consequence of material balance, even though this model has no delayed channel activation.

This distinction guides interpretation: a large released amount need not mean a large simultaneous free cytosolic concentration. Clearance and volume intervene. Real biological cells add binding, store refilling, IP3 turnover and IP3 receptor regulation; the two plotted curves isolate one finite-store effect rather than explaining every calcium response.