Skip to content


  • Poster presentation
  • Open Access

Fast and accurate representations of stochastic ion channel fluctuations

  • 1,
  • 2,
  • 3,
  • 3,
  • 4, 5,
  • 6,
  • 7 and
  • 6Email author
BMC Neuroscience201516 (Suppl 1) :P258

  • Published:


  • Stochastic Differential Equation
  • Channel Noise
  • Channel Gating
  • Isochron
  • Stochastic Simulation Algorithm

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 [1] 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 [2]. 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 [3] 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 [4] 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.



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.

Authors’ Affiliations

Department of Mathematics, University of Wisconsin, Madison, WI 53706, USA
Department of Mathematics, University of Pittsburgh, Pittsburgh, PA 15260, USA
Department of Neurosciences, Case Western Reserve University, Cleveland, OH 44106, USA
Bernstein Center for Computational Neuroscience, 10115 Berlin, Germany
Department of Physics, Humboldt University, 12489 Berlin, Germany
Department of Mathematics, Applied Mathematics and Statistics, Case Western Reserve University, Cleveland, OH 44106, USA
Department of Mathematics and Statistics, University of Nevada, Reno, NV 89557, USA


  1. 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.PubMedView ArticleGoogle Scholar
  2. 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-PubMedView ArticleGoogle Scholar
  3. 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-PubMedPubMed CentralView ArticleGoogle Scholar
  4. Schwabedal Justus TC, Arkady Pikovsky: Phase description of stochastic oscillations. Phys rev lett. 2013, 110 (20): 204102-PubMedView ArticleGoogle Scholar
  5. Thomas Peter J, Benjamin Lindner: Asymptotic Phase for Stochastic Oscillators. Phys rev lett. 2014, 113 (25): 254101-PubMedView ArticleGoogle Scholar
  6. Potoyan Davit A, Peter G Wolynes: On the dephasing of genetic oscillators. Proc Natl Acad Sci. 2011, 111 (6): 2391-2396.View ArticleGoogle Scholar


© Anderson et al. 2015

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.