- Poster presentation
- Open Access
Fast and accurate representations of stochastic ion channel fluctuations
BMC Neuroscience volume 16, Article number: P258 (2015)
Random ion channel gating is an important source of noise at the single neuron level. For mesoscale ion channel population sizes, fast, accurate representation of channel noise fluctuations remains an important challenge. We highlight recent progress in three areas.
In  we present an exact stochastic simulation algorithm, which takes into account the time dependence of the voltage sensitive transitions due to rapid voltage changes during action potentials. The exact algorithm is similar to a widely used approximate stochastic simulation algorithm, in which transition rates are held fixed during the intervals between channel state transitions. We compare the algorithms and show that they are inequivalent in a strong sense, meaning that sample paths diverge when driven with identical Poisson processes, leading to different precise firing times. But the stationary histograms produced are practically indistinguishable for modest channel numbers (circa N>100), indicating weak equivalence.
For channel numbers in the hundreds or higher, numerical stochastic differential equations (SDE) algorithms based on the system size expansion can be significantly faster than simulations based on discrete Markov chain simulations, while retaining reasonable accuracy. For large networks, however, even SDE based simulations become costly, particularly as more complex channel gating schemes are introduced. Schmandt and Galán introduced a stochastic shielding approximation as a fast, accurate way of simulating stochastic ion channel kinetics . In the SDE representation, each edge in the graph generates both a mean population flux between adjacent nodes, and a fluctuation about the mean. Only the fluctuations arising from edges connecting functionally distinct states directly affect fluctuations in the observed behavior of the cell. In  we analyze stochastic shielding both for the HH sodium and potassium gating models, and for an ensemble of random graphs. We derive a quantitative measure of edge importance related to the eigenvalue/eigenvector decomposition of the graph Laplacian matrix.
Channel noise makes a regularly spiking neuron a stochastic oscillator. The classical definition for the asymptotic phase of an oscillator breaks down when stochasticity is taken into account. Alternative definitions of the ``phase'' have been based on the mean first passage time property of a system of isochronal surfaces  and in terms of the spectral decomposition of a the adjoint Kolmogorov operator [5, 6]. In the vanishing noise limit, the spectrally derived isochrons approach those of the underlying mean field system, if the latter has a finite period limit cycle.
Together, these results expand our analytic, numerical, and conceptual tools for understanding the effects of random ion channel gating in conductance based neural models.
Anderson David F, Bard Ermentrout, Peter J Thomas: Stochastic representations of ion channel kinetics and exact stochastic simulation of neuronal dynamics. J Comput Neurosci. 2014, 38 (1): 67-82.
Schmandt Nicolaus T, Roberto F Galán: Stochastic-shielding approximation of Markov chains and its application to efficiently simulate random ion-channel gating. Phys rev lett. 2012, 109 (11): 118101-
Schmidt Deena R, Peter J Thomas: Measuring edge importance: a quantitative analysis of the stochastic shielding approximation for random processes on graphs. J Math Neurosci. 2014, 4 (1): 6-
Schwabedal Justus TC, Arkady Pikovsky: Phase description of stochastic oscillations. Phys rev lett. 2013, 110 (20): 204102-
Thomas Peter J, Benjamin Lindner: Asymptotic Phase for Stochastic Oscillators. Phys rev lett. 2014, 113 (25): 254101-
Potoyan Davit A, Peter G Wolynes: On the dephasing of genetic oscillators. Proc Natl Acad Sci. 2011, 111 (6): 2391-2396.
Simons Foundation (#259837 to PJT), Council for the International Exchange of Scholars (CIES), Mathematical Biosciences Institute (NSF-DMS 0931642), NSF EF-1038677, DMS-1413770, DMS-1318832, DMS-1219754, Bundesministerium für Bildung und Forschung (FKZ: 01GQ1001A), Army Research Office W911NF-14-1-0401.