On stochastic differential equation models for ion channel noise in Hodgkin-Huxley neurons

TL;DR

This paper demonstrates that certain stochastic differential equation models can accurately approximate ion channel noise in Hodgkin-Huxley neurons, providing a simpler yet effective alternative to complex Markov chain models.

Contribution

It shows that channel-based SDE models can capture ion channel noise effects more accurately than subunit-based models, clarifying the conditions for effective approximation.

Findings

Channel-based SDE models replicate channel noise distribution.

Subunit-based SDE models fail to capture noise effects accurately.

Markov chain models can be effectively approximated by SDEs.

Abstract

The random transitions of ion channels between conducting and non-conducting states generate a source of internal fluctuations in a neuron, known as channel noise. The standard method for modeling fluctuations in the states of ion channels uses continuous-time Markov chains nonlinearly coupled to a differential equation for voltage. Beginning with the work of Fox and Lu, there have been attempts to generate simpler models that use stochastic differential equation (SDEs) to approximate the stochastic spiking activity produced by Markov chain models. Recent numerical investigations, however, have raised doubts that SDE models can preserve the stochastic dynamics of Markov chain models. We analyze three SDE models that have been proposed as approximations to the Markov chain model: one that describes the states of the ion channels and two that describe the states of the ion channel subunits. We show that the former channel-based approach can capture the distribution of channel noise and its effect on spiking in a Hodgkin-Huxley neuron model to a degree not previously demonstrated, but the latter two subunit-based approaches cannot. Our analysis provides intuitive and mathematical explanations for why this is the case: the temporal correlation in the channel noise is determined by the combinatorics of bundling subunits into channels, and the subunit-based approaches do not correctly account for this structure. Our study therefore confirms and elucidates the findings of previous numerical investigations of subunit-based SDE models. Moreover, it presents the first evidence that Markov chain models of the nonlinear, stochastic dynamics of neural membranes can be accurately approximated by SDEs. This finding opens a door to future modeling work using SDE techniques to further illuminate the effects of ion channel fluctuations on electrically active cells.