Simulations of intracellular calcium release dynamics in response to a high-intensity, ultrashort electric pulse

Numerical simulations for electrically induced, intracellular calcium release from the endoplasmic reticulum are reported. A two-step model is used for self-consistency. Distributed electrical circuit representation coupled with the Smoluchowski equation yields the ER membrane nanoporation for calcium outflow based on a numerical simulation. This is combined with the continuum Li-Rinzel model and drift diffusion for calcium dynamics. Our results are shown to be in agreement with reported calcium release data. A modest increase (rough doubling) of the cellular calcium is predicted in the absence of extra-cellular calcium. In particular, the applied field of 15 kV/cm with 60 ns pulse duration makes for a strong comparison. No oscillations are predicted and the net recovery period of about 5 min are both in agreement with published experimental results. A quantitative explanation for the lack of such oscillatory behavior, based on the density dependent calcium fluxes, is also provided.