Period doubling of calcium spike firing in a model of a Purkinje cell dendrite

Recordings from cerebellar Purkinje cell dendrites have revealed that in response to sustained current injection, the cell firing pattern can move from tonic firing of Ca²⁺ spikes to doublet firing and even to quadruplet firing or more complex firing. These firing patterns are not modified substantially if Na⁺ currents are blocked. We show that the experimental results can be viewed as a slow transition of the neuronal dynamics through a period-doubling bifurcation. To further support this conclusion and to understand the underlying mechanism that leads to doublet firing, we develop and study a simple, one-compartment model of Purkinje cell dendrite. The neuron can also exhibit quadruplet and chaotic firing patterns that are similar to the firing patterns that some of the Purkinje cells exhibit experimentally. The effects of parameters such as temperature, applied current, and potassium reversal potential in the model resemble their effects in experiments. The model dynamics involve three time scales. Ca²⁺- dependent K⁺ currents, with intermediate time scales, are responsible for the appearance of doublet firing, whereas a very slow hyperpolarizing current transfers the neuron from tonic to doublet firing. We use the fast-slow analysis to separate the effects of the three time scales. Fast-slow analysis of the neuronal dynamics, with the activation variable of the very slow, hyperpolarizing current considered as a parameter, reveals that the transitions occurs via a cascade of period-doubling bifurcations of the fast and intermediate subsystem as this slow variable increases. We carry out another analysis, with the Ca²⁺ concentration considered as a parameter, to investigate the conditions for the generation of doublet firing in systems with one effective variable with intermediate time scale, in which the rest state of the fast subsystem is terminated by a saddle-node bifurcation. We find that the scenario of period doubling in these systems can occur only if (1) the time scale of the intermediate variable (here, the decay rate of the calcium concentration) is slow enough in comparison with the interspike interval of the tonic firing at the transition but is not too slow and (2) there is a bistability of the fast subsystem of the spike-generating variables.