Complex oscillations and waves of calcium in pancreatic acinar cells

We perform a bifurcation analysis of a model of Ca²⁺ wave propagation in the basal region of pancreatic acinar cells. The model we consider was first presented in Sneyd et al. [J. Sneyd, K. Tsaneva-Atanasova, J.I.E. Bruce, S.V. Straub, D.R. Giovannucci, D.I. Yule, A model of calcium waves in pancreatic and parotid acinar cells, Biophys. J. 85 (2003) 1392–1405], where a partial bifurcation analysis was given of the model in the absence of diffusion. We obtain more complete information about bifurcations of the diffusionless model via numerical studies, then analyse the spatially extended model by numerical investigation of the travelling wave equations and direct numerical solution of the model equations. We find solitary waves in the model equations arising from homoclinic bifurcations in the travelling wave equations. The solitary waves exist and appear to be stable for a significant interval of the primary bifurcation parameter (i.e., the concentration of inositol trisphosphate) but are eventually replaced by irregular spatio-temporal behaviour. The homoclinic bifurcations are related to a number of complicated mathematical structures in the travelling wave equations, including an anomalous homoclinic-Hopf bifurcation, heteroclinic bifurcations between an equilibrium and a periodic orbit, and homoclinic bifurcations of periodic orbits.