Fast and Accurate Langevin Simulations of Stochastic Hodgkin-Huxley Dynamics

DOI10.1162/NECO_A_01312zbMATH Open1453.92060arXiv2005.10598OpenAlexW3049668116WikidataQ98469098 ScholiaQ98469098MaRDI QIDQ3386401FDOQ3386401

Authors: Shusen Pu, Peter J. Thomas 

Publication date: 4 January 2021

Published in: Neural Computation (Search for Journal in Brave)

Abstract: Fox and Lu introduced a Langevin framework for discrete-time stochastic models of randomly gated ion channels such as the Hodgkin-Huxley (HH) system. They derived a Fokker-Planck equation with state-dependent diffusion tensor $D$ and suggested a Langevin formulation with noise coefficient matrix $S$ such that $S{S}^{i}ntercal=D$. Subsequently, several authors introduced a variety of Langevin equations for the HH system. In this paper, we present a natural 14-dimensional dynamics for the HH system in which each emph{directed} edge in the ion channel state transition graph acts as an independent noise source, leading to a $14imes28$ noise coefficient matrix $S$. We show that (i) the corresponding 14D system of ordinary differential ev{equations} is consistent with the classical 4D representation of the HH system; (ii) the 14D representation leads to a noise coefficient matrix $S$ that can be obtained cheaply on each timestep, without requiring a matrix decomposition; (iii) sample trajectories of the 14D representation are pathwise equivalent to trajectories of Fox and Lu's system, as well as trajectories of several existing Langevin models; (iv) our 14D representation (and those equivalent to it) give the most accurate interspike-interval distribution, not only with respect to moments but under both the $L_1$ and $L_{i}nfty$ metric-space norms; and (v) the 14D representation gives an approximation to exact Markov chain simulations that are as fast and as efficient as all equivalent models. Our approach goes beyond existing models, in that it supports a stochastic shielding decomposition that dramatically simplifies $S$ with minimal loss of accuracy under both voltage- and current-clamp conditions.

Full work available at URL: https://arxiv.org/abs/2005.10598

Recommendations

• Multithreaded stochastic PDES for reactions and diffusions in neurons
• Stochastic representations of ion channel kinetics and exact stochastic simulation of neuronal dynamics
• An investigation of the stochastic Hodgkin-Huxley models under noisy rate functions
• Stochastic simulation of neurons, axons, and action potentials
• Fast and exact simulation methods applied on a broad range of neuron models
• Efficient calculation of heterogeneous non-equilibrium statistics in coupled firing-rate models
• An improved parameter estimation method for Hodgkin-Huxley models
• High-order accurate difference schemes for the Hodgkin-Huxley equations
• Stochastic Hodgkin-Huxley equations with colored noise terms in the conductances
• Hybrid stochastic simulations of intracellular reaction-diffusion systems

Mathematics Subject Classification ID

Neural biology (92C20) Applications of stochastic analysis (to PDEs, etc.) (60H30)

Cites Work

• Empirical Processes with Applications to Statistics
• Asymptotic Minimax Character of the Sample Distribution Function and of the Classical Multinomial Estimator
• Title not available (Why is that?)
• Table for Estimating the Goodness of Fit of Empirical Distributions
• Title not available (Why is that?)
• Mathematical physiology
• A stochastic analysis of first-order reaction networks
• Title not available (Why is that?)
• Construction of Equivalent Stochastic Differential Equation Models
• Title not available (Why is that?)
• Title not available (Why is that?)
• Stochastic analysis of biochemical systems
• Global Asymptotic Stability of Solutions of Nonautonomous Master Equations
• Invariant Manifolds of Binomial-Like Nonautonomous Master Equations
• Title not available (Why is that?)
• Invariant manifold reductions for Markovian ion channel dynamics
• Stochastic representations of ion channel kinetics and exact stochastic simulation of neuronal dynamics
• An investigation of the stochastic Hodgkin-Huxley models under noisy rate functions
• Stochastic modeling of excitable dynamics: improved Langevin model for mesoscopic channel noise
• Dynamics of gamma bursts in local field potentials
• Measuring edge importance: a quantitative analysis of the stochastic shielding approximation for random processes on graphs
• Stochastic Hodgkin-Huxley equations with colored noise terms in the conductances
• Exact reductions of Markovian dynamics for ion channels with a single permissive state
• Spontaneous action potentials and neural coding in unmyelinated axons

Cited In (11)

• Stochastic representations of ion channel kinetics and exact stochastic simulation of neuronal dynamics
• Measuring edge importance: a quantitative analysis of the stochastic shielding approximation for random processes on graphs
• Fluctuation-response relations for integrate-and-fire models with an absolute refractory period
• Efficient Implicit Runge-Kutta Methods for Fast-Responding Ligand-Gated Neuroreceptor Kinetic Models
• Mean-return-time phase of a stochastic oscillator provides an approximate renewal description for the associated point process
• Dynamical consequences of sensory feedback in a half-center oscillator coupled to a simple motor system
• Stochastic many-particle model for LFP electrodes
• Stochastic modeling of excitable dynamics: improved Langevin model for mesoscopic channel noise
• Title not available (Why is that?)
• Heteroclinic cycling and extinction in May-Leonard models with demographic stochasticity
• Resolving molecular contributions of ion channel noise to interspike interval variability through stochastic shielding

This page was built for publication: Fast and Accurate Langevin Simulations of Stochastic Hodgkin-Huxley Dynamics

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q3386401)