Comparison of algorithms for the simulation of action potentials with stochastic sodium channels.
This paper presents a comparison of computational algorithms to simulate action potentials using stochastic sodium channels. Four algorithms are compared in single-node models: Strassberg and DeFelice (1993) (SD), Rubinstein (1995) (R), Chow and White (1996) (CW), and Fox (1997) (F). Neural responses are simulated to a simple and a preconditioned monophasic current pulse. Three exact algorithms implementing Markov jumping processes (SD, R, CW) resulted in similar responses, while the approximation algorithm using Langevin's equation (F) showed quite different responses from those in the exact algorithms. The computational time was measured as well: 1(F), 7(CW), 32(SD), 39(R) relative to that of the F algorithm. Furthermore, it is shown that as the sampling step for integration of the transmembrane potential increases, neural responses in all algorithms tended to be different from those in dense sampling steps, however, the CW algorithm was robust even at a sparse sampling step. It is concluded that the most computationally efficient algorithm having appropriate properties of neural excitability is the CW algorithm.
Mino
Hiroyuki
H
Department of Otolaryngology, Head and Neck Surgery, The University of Iowa College of Medicine, Iowa City 52242, USA. hiromino@earpower.oto.uiowa.edu
Rubinstein
Jay T
JT
White
John A
JA
D000200:Action Potentials / Q000502:physiology*; D000465:Algorithms; D000818:Animals; D003198:Computer Simulation; D004558:Electric Stimulation; D008322:Mammals; D008390:Markov Chains; D008564:Membrane Potentials / Q000502:physiology; D008959:Models, Neurological; D011901:Ranvier's Nodes / Q000502:physiology*; D015203:Reproducibility of Results; D012680:Sensitivity and Specificity; D015222:Sodium Channels / Q000502:physiology*; D013269:Stochastic Processes
D015222:Sodium Channels