Optimization of battery strengths in the Hodgkin-Huxley model

Abstract

Simulations show that neurons are capable of operating over a much broader range of values of ionic reversal potentials than what is actually observed. Since the reversal potentials, which depend on the ionic concentration gradients across the membrane, have strong effects on the signaling properties and metabolic energy consumption rates of neurons, it is natural to hypothesize that the actual values are optimal for some functional property of the neuron involving energy and information, much as the ion channel densities and leak conductance appear to be [1,2]. Here we consider the Hodgkin-Huxley model of the squid giant axon and systematically investigate whether the reversal potentials for sodium and potassium, ENa and EK, whose experimental values are respectively around +50 and -80 mV, minimize or maximize any simple combination of: (a) metabolic energy consumption during an action potential; (b) metabolic energy consumption during quiescence; (c) action potential velocity; (d) maximum firing frequency; and (e) the “energy efficiency” of the action potentials, which is the proportion of the metabolic energy associated with them that actually contributes to net inward or outward currents. None of these quantities, nor any algebraically simple combination of them, has a maximum or minimum over the range of {ENa, EK} parameter space that we investigated. However, there are many functions for which the potassium reversal potential alone has an optimal value. Further analysis reveals that there are three basic underlying optimizations involving EK. The velocity of action potentials is maximized for a value of EK in the experimental range (Fig. ​(Fig.1),1), and the product of the metabolic energy consumption during an action potential and at rest is minimized. The third optimization is for the ratio of action potential energy and maximum firing frequency, which is maximized for biological EK. The interpretation of this optimization may be clear only in the context of the larger squid escape-jet system. Figure 1 The velocity of action potentials (vAP) in the Hodgkin-Huxley model of the squid giant axon as a function of the potassium reversal potential (EK), with the sodium reversal potential (ENa) held fixed at various values. In all cases, vAP attains a maximum ...