Direct stochastic simulation of Ca²⁺ motion in Xenopus eggs

The release of important intracellular ions has been widely modeled using two approaches, namely, (1) Fickian diffusion, in which sometimes tensorial diffusion coefficients are used to fit observed temporally varying concentrations of calcium, and (2) cellular automata, which produce a set of localized finite difference equations that result in complex global behavior. Here, we take a different approach, employing some assumed, a priori, distribution of ion-binding proteins in the cell, and some assumed biochemical capture and release characteristics to explain ionic motion, and ultimately, distribution. We study several scenarios for ion distribution, based on differences in binder action and distribution. The numbers and strengths of ion binders, spatial variation in inositol 1,4,5-triphosphate concentration, together with the escalating distribution of ionic diffusion speed, are found to be key factors leading to concavity in the Ca²⁺ wave shape. We also offer an explanation for geometrical effects on previously observed ion diffusion speeds in the cellular cortex of the Xenopus laevis egg during fertilization, based on an angle-of-view correction.