Skip to content


  • Poster presentation
  • Open Access

Effect of power-law ionic conductances in the Hodgkin and Huxley model

BMC Neuroscience201516 (Suppl 1) :P250

  • Published:


  • Firing Rate
  • Fractional Derivative
  • Spike Activity
  • Past Activity
  • Spike Rate

An increasing number of results show that the voltage and spiking activity of a neuron follows scale free adaptation. Under such conditions, the voltage, or firing rate, of a neuron cannot be characterized by a unique time constant, instead these processes are characterized by power-laws. Power-law behaviors suggest that the components responsible for the voltage are strongly interacting across temporal scales, slow processes affect the rates of fast processes and vice versa. We recently introduced the fractional leaky integrate-and-fire model, which we have used to replicate the firing rate activity of adapting cortical neurons [1]. Our results have shown that spike rate adaptation can be modeled by a fractional derivative of low order. Thus suggesting, a strong interaction across conductances in cortical neurons.

this work we decided to study the effects of power-law behavior in a biophysical model of spiking activity, the Hodgkin-Huxley model. The classical Hodgkin-Huxley model is described by
C d V d t = ( g m ( V E r e s t ) + i = 1 N g i ( V E i ) ) + I .
where C is the capacitance; V, voltage; gm, the membrane conductance; Erest, the resting potential; gi={gNa, gK } the sodium or potassium conductances with corresponding reversal potentials (Ei={ ENa, EK }); and I, the input curret. The conductances are g K = g K ¯ n 4 and g N a = g N a ¯ m 3 h , with g K ¯ and g N a ¯ the maximum conductances. The gating variables n, m, and h are defined by the general equation
d x d t = α x V , t 1 - x - β x ( V , t ) x
where x={n, m, h}. To implement power-law behavior in either gating variable we substitute the classical derivative by the fractional derivative of order η using the Caputo definition
d η f d t η = 1 Γ n - η 0 t f ( n ) t t - u η + 1 - n d u

where is the Gamma function, and n is the nearest integer larger than η. In our case, n = 1 because this models how the interaction of previous activity slows down the activation of the gating variables. The fractional derivative value is the result of integrating the activity of the function over all past activities weighted by a function that follows a power-law. The weighted values are called the memory trace. Thus, the fractional derivative provides information over all past activity, in contrast with the classical derivative that only takes into account the value of the function in the immediate previous time point. We have recently developed efficient ways to computationally solve these equations [2].

Our results show the emergence of a wide range of spiking behaviors (bursting, spiking and sub-threshold oscillations) in response to constant stimulation as a function of the fractional order in the different activation/inactivation variables. In the case of n, the neuron shows reduction if spiking response and emergence of sub-threshold oscillations. While fractional h results in bursting activity. This emergent richness in spiking activity, while only modeling two conductances, allows study of the overall effects of power-law behavior in neuronal activity. At the biophysical level, our results suggest that voltage dependent ion channels that deviate from a classical Markov process could increase the amount of spiking patterns in single cells, thus increasing their information content capacity.



NSF DBI 1451032 (BRAIN-EAGER) and NSF EF 1137897

Authors’ Affiliations

Wodimu TekaUTSA Neurosciences Institute, University of Texas at San Antonio, San Antonio, TX 78249, USA


  1. Teka W, Marinov TM, Santamaria F: Neuronal Spike Timing Adaptation Described with a Fractional Leaky Integrate-and-Fire Model. PLoS Comput Biol. 2014, 10 (3): e1003526-PubMedPubMed CentralView ArticleGoogle Scholar
  2. Marinov T, Ramirez N, Santamaria F: Fractional integration toolbox. fcaa. 2013, 16 (3): 670-681.View ArticleGoogle Scholar


© Santamaria 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.