A Computational Model of the Temperature-dependent Changes in Firing Patterns in Aplysia Neurons

We performed experiments using Aplysia neurons to identify the mechanism underlying the changes in the firing patterns in response to temperature changes. When the temperature was gradually increased from $11^{\circ}C$ to $31^{\circ}C$ the firing patterns changed sequentially from the silent state to beating, doublets, beating-chaos, bursting-chaos, square-wave bursting, and bursting-oscillation patterns. When the temperature was decreased over the same temperature range, these sequential changes in the firing patterns reappeared in reverse order. To simulate this entire range of spiking patterns we modified nonlinear differential equations that Chay and Lee made using temperature-dependent scaling factors. To refine the equations, we also analyzed the spike pattern changes in the presence of potassium channel blockers. Based on the solutions of these equations and potassium channel blocker experiments, we found that, as temperature increases, the maximum value of the potassium channel relaxation time constant, ${\tau}_n(t)$ increases, but the maximum value of the probabilities of openings for activation of the potassium channels, n(t) decreases. Accordingly, the voltage-dependent potassium current is likely to play a leading role in the temperature-dependent changes in the firing patterns in Aplysia neurons.