Next Article in Journal
Universal Dimensions of Meaning Derived from Semantic Relations among Words and Senses: Mereological Completeness vs. Ontological Generality
Previous Article in Journal
A 3-D Model of a Perennial Ryegrass Primary Cell Wall and Its Enzymatic Degradation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Can the Thermodynamic Hodgkin-Huxley Model of Voltage-Dependent Conductance Extrapolate for Temperature?

1
Department of Computer Science, University of Warwick, Coventry, CV4 7AL, UK
2
Institute for Adaptive and Neural Computation, School of Informatics, University of Edinburgh, Edinburgh, EH1 2QL, UK 
Submission received: 25 February 2014 / Revised: 12 April 2014 / Accepted: 4 May 2014 / Published: 14 May 2014
(This article belongs to the Section Computational Biology)

Abstract

:
Hodgkin and Huxley (H-H) fitted their model of voltage-dependent conductances to experimental data using empirical functions of voltage. The thermodynamic H-H model of voltage dependent conductances is more physically plausible, as it constrains and parameterises its empirical fit by assuming that ion channel transition rates depend exponentially on a free energy barrier that in turn, linearly or non-linearly, depends on voltage. The original H-H model contains no explicit temperature terms and requires Q10 factors to describe data at different temperatures. The thermodynamic H-H model does have explicit terms for temperature. Do these endow the model with extrapolation for temperature? We utilised voltage clamp data for a voltage-gated K+ current, recorded at three different temperatures. The thermodynamic H-H model’s free parameters were fitted (Marquardt-Levenberg algorithm) to a data set recorded at one (or more) temperature(s). Then we assessed whether it could describe another data set, recorded at a different temperature, with these same free parameter values and its temperature terms set to the new temperature. We found that it could not.

Graphical Abstract

1. Introduction

Hodgkin and Huxley (H-H) quantitatively characterised the voltage-dependence of membrane currents in the giant squid axon, and showed how they can generate and propagate action potentials [1]. Their model can describe ion channels in the membrane, opening and closing in a voltage-dependent manner [2]. However, membrane currents are not just voltage-dependent, but temperature dependant also, and the H-H description is without any explicit temperature term. To describe data over a range of temperatures, the H-H model requires an empirical multiplicative constant, Q10, to be applied to a number of its parameters [1].
Hodgkin and Huxley fitted their model of voltage-dependent conductances to experimental data using empirical functions of voltage. Thermodynamic Hodgkin-Huxley models describe voltage-dependence empirically as well, but can be construed to be more physically plausible, as they constrain and parameterise their fit with thermodynamic principles of transition state theory [3,4,5,6,7]. They consider that the rate of transition between channel states depends exponentially on the free energy barrier that separates them. The free energy barrier, in turn, is voltage dependent. It can vary linearly or non-linearly with voltage: linear and non-linear thermodynamic H-H models respectively [6,7].
In contrast to the original H-H model, thermodynamic H-H models do have explicit terms for temperature, on account of their incorporation of thermodynamic theory. We hypothesise that these temperature terms permit these models to describe voltage-gated currents at different temperatures i.e.,describe how currents change with voltage and temperature.
Voltage clamp data for a K+ current was taken from [8]. The data was grouped at three temperatures: T15 (15 °C), T25 (25 °C) and T35 (35 °C). The thermodynamic H-H model (linear and non-linear) was fitted to each temperature set separately; the Marquardt-Levenberg algorithm [9] optimised the model’s free parameters to find a best fit in each case. This was a success, with close fits being produced. We then investigated if the model, with its free parameters fitted to one temperature data set, could replicate the data of a different temperature set, when the temperature terms of the model were changed to the new temperature. We found that it could not.
We repeated this process but, in the fitting phase, we fitted the thermodynamic H-H model (linear and non-linear) to the combined data sets of two different temperatures: (T15 + T25) or (T15 + T35) or (T25 + T35). We then investigated if the model, with its free parameters fitted to two temperature data sets, could replicate the data of a third temperature set, when the temperature terms of the model were changed to the new temperature. We found that it could not.
We find that the thermodynamic H-H model (linear or non-linear) cannot accurately predict data at temperatures that it is not specifically fitted to i.e., it cannot extrapolate for temperature.

2. Experimental Section

2.1. K+ Current Data

Tiwari and Sikdar studied non-inactivating K+ currents, in a gonadotroph cell line, across a temperature range [8]. They assembled whole cell, voltage clamp data at a number of depolarising potentials, grouped at three temperatures: T15 (15 °C, 288 K), T25 (25 °C, 298 K), T35 (35 °C, 308 K). We obtained this data via personal communication and modified it, as shown in Figure 1.
Figure 1. K+ current in a gonadotroph cell, when the holding voltage is stepped from −10 mV to +116 mV using a voltage clamp. We modified this data before we used it ourselves, removing the labelled capacitance spikes and time lag. We repeated this action for all other current data used. The x-axis presents time (mS); the y-axis presents K+ current.
Figure 1. K+ current in a gonadotroph cell, when the holding voltage is stepped from −10 mV to +116 mV using a voltage clamp. We modified this data before we used it ourselves, removing the labelled capacitance spikes and time lag. We repeated this action for all other current data used. The x-axis presents time (mS); the y-axis presents K+ current.
Computation 02 00047 g001

2.2. Thermodynamic H-H Models

α and β are forward and backward rate constants between the Closed (C) and Open (O) channel states.
Computation 02 00047 i001
When the channel switches between these states, it must pass through a high energy intermediate: the transition state (not shown). α and β are voltage (V) and temperature (T) dependent:

2.2.1. Linear Variant

Computation 02 00047 i002
Computation 02 00047 i003
R is the gas constant, V is the membrane voltage. α0, a1, b1 are parameters describing the free energy barrier between the Closed and Transition state. β0, a2, b2 are parameters describing the free energy barrier between the Open and Transition state. T is temperature; so this model has terms for temperature, unlike the original Hodgkin-Huxley formulation [1].

2.2.2. Non-Linear Variant (Quadratic)

Computation 02 00047 i004
Computation 02 00047 i005
R is the gas constant, V is the membrane voltage. α0, a1, b1, c1 are parameters describing the free energy barrier between the Closed and Transition state. β0, a2, b2, c2 are parameters describing the free energy barrier between the Open and Transition state. T is temperature; so this model has terms for temperature, unlike the original Hodgkin-Huxley formulation [1].

2.3. Thermodynamic H-H Model of a Non-Inactivating K+ Current (IK), Recorded in a Gonadotroph Cell Line

IK = gK(VEK)
gK = gK · n
Computation 02 00047 i006
Computation 02 00047 i007
Computation 02 00047 i008
IK is the K+ current, gK is the K+ conductance, gK is the maximal K+ conductance, V is the membrane potential, t is time and EK is the reversal/Nernst potential for K+ (set to −9.8 mV; [8]). To produce a linear thermodynamic H-H model, αn and βn are set by Equations (2) and (3) (respectively). To produce a non-linear thermodynamic H-H model, αn and βn are instead set by Equations (4) and (5) (respectively).
The maximal conductance (gK) for the T15, T25 and T35 membrane patches was calculated by deriving a maximal current (IK) value for the largest depolarising potential explored at each temperature (Figure 1D; [8]). Then:
Computation 02 00047 i009

2.4. Fitting the Thermodynamic H-H Model to the K+ Current (IK) Data

The free parameters in the thermodynamic H-H model, adjusted in order to fit the model to the data, are α0, β0, a1, b1, a2, b2 for the linear variant (refer to Equations (2) and (3)); with the addition of c1 and c2 for the non-linear variant (quadratic) (refer to Equations (4) and (5)). The fitting was performed using the criterion of least squares minimisation, implemented by the Marquardt-Levenberg algorithm [9] in Matlab (MathWorks Inc., Natick, MA, USA). Initial free parameter values were chosen arbitrarily. Repeated runs with different initial free parameter values checked that our fits were globally, and not just locally, optimal. No constraints were set for free parameter values.

2.5. Temperature Extrapolation

The thermodynamic H-H model (linear and non-linear) was fitted to each temperature set separately: T15, T25 and T35. We then investigated if the model, with its free parameters fitted to one temperature data set, could replicate the data of a different temperature set, when the temperature terms of the model (in Equations (2)–(5)) were changed to the new temperature.
We repeated this process but, in the fitting phase, we fitted the thermodynamic H-H model (linear and non-linear) to the combined data sets of two different temperatures: (T15 + T25), (T25 + T35), (T15 + T35). At each temperature, the temperature setting in the thermodynamic model equations was set appropriately and the model’s free parameters were optimised with experimental data spanning more than one temperature. We then investigated if the model, with its free parameters fitted to two temperature data sets, could replicate the data of a third temperature set, when the temperature terms of the model were changed to the new temperature.
So, to elaborate further on what was done for both the linear and non-linear models: where p denotes parameters and d is data: (T15p on T25d), (T15p on T35d), (T25p on T15d), (T25p on T35d), (T35p on T15d) and (T35p on T25d). Combined data sets were also used to “train” the free parameter values: (T15+T25p on T35d), (T15+T35p on T25d), (T25+T35p on T15d).
Figure 2. Thermodynamic Hodgkin-Huxley models, linear and non-linear, (red lines) fitted to K+ current data (black circles). (A) Linear thermodynamic H-H model fitted to the T15 data set; (B) Linear thermodynamic H-H model fitted to the T25 data set; (C) Linear thermodynamic H-H model fitted to the T35 data set; (D) Non-linear thermodynamic H-H model fitted to the T15 data set; (E) Non-linear thermodynamic H-H model fitted to the T25 data set; (F) Non-linear thermodynamic H-H model fitted to the T35 data set. For all panels: the x-axis presents time (ms), the y-axis presents K+ current (µA).
Figure 2. Thermodynamic Hodgkin-Huxley models, linear and non-linear, (red lines) fitted to K+ current data (black circles). (A) Linear thermodynamic H-H model fitted to the T15 data set; (B) Linear thermodynamic H-H model fitted to the T25 data set; (C) Linear thermodynamic H-H model fitted to the T35 data set; (D) Non-linear thermodynamic H-H model fitted to the T15 data set; (E) Non-linear thermodynamic H-H model fitted to the T25 data set; (F) Non-linear thermodynamic H-H model fitted to the T35 data set. For all panels: the x-axis presents time (ms), the y-axis presents K+ current (µA).
Computation 02 00047 g002

3. Results and Discussion

3.1. Curve Fitting

The thermodynamic H-H model (linear and non-linear) was fitted to each temperature data set separately. Figure 2 presents the linear and non-linear models fitted separately to the T15, T25 and T35 data sets. Each panel shows K+ current in a gonadotroph cell (black circles; every tenth data point plotted) when the holding voltage is stepped from −10 mV to +57 mV, from −10 mV to +77 mV, from −10 mV to +96 mV and from −10 mV to +116 mV [8]. So, there are 4 data plots at each temperature. Greater depolarisation produces greater K+ current. Red lines show the thermodynamic Hodgkin-Huxley (H-H) model—linear or non-linear—fitted to the data; its free parameters optimised to the data by the Marquardt-Levenberg algorithm [9]. This was a success, with close fits being produced. For example, the non-linear model fit to the T15 data set has a logged reduced chi-square metric: log10 Computation 02 00047 i010 = −0.05. The worst fit was still good in absolute terms: the linear model fit to the T35 data set has a logged reduced chi-square metric: log10 Computation 02 00047 i010 = 1.69. Temperature and goodness of fit were inversely related, possibly due to data variance being a function of temperature. The non-linear model provided a slightly better fit than the linear model.
The thermodynamic H-H model (linear and non-linear) was then fitted to the combined data sets of two different temperatures: (T15 + T25), (T15 + T35), (T25 + T35). At each temperature, the temperature setting in the thermodynamic model equations was set appropriately and the model’s free parameters were optimised with experimental data spanning more than one temperature. Figure 3 presents the linear and non-linear models fitted separately to the (T15 + T25), (T15 + T35) and (T25 + T35) data sets. This was a reasonable success, with reasonably close fits being produced. But, the fit to an amalgamation data set (recorded at 2 temperatures) is inferior to the fit of any of its individual component data sets (recorded at 1 temperature). Again, temperature and goodness of fit were inversely related. The non-linear model provided a slightly better fit than the linear model.
To conclude, the thermodynamic H-H formalisms can describe K+ current data recorded at a certain temperature(s) well, if their free parameters are specifically tuned to describe that data.
Figure 3. Thermodynamic Hodgkin-Huxley models, linear and non-linear, (red lines) fitted to K+ current data (black circles). (A) Linear thermodynamic H-H model fitted to the T15 and T25 data set combined (T15, T25); (B) Linear thermodynamic H-H model fitted to the T15 and T35 data set combined (T15, T35); (C) Linear thermodynamic H-H model fitted to the T25 and T35 data set combined (T25, T35); (D) Non-linear thermodynamic H-H model fitted to the T15 and T25 data set combined (T15, T25); (E) Non-linear thermodynamic H-H model fitted to the T15 and T35 data set combined (T15, T35); (F) Non-linear thermodynamic H-H model fitted to the T25 and T35 data set combined (T25, T35). For all panels: the x-axis presents time (ms), the y-axis presents K+ current (µA).
Figure 3. Thermodynamic Hodgkin-Huxley models, linear and non-linear, (red lines) fitted to K+ current data (black circles). (A) Linear thermodynamic H-H model fitted to the T15 and T25 data set combined (T15, T25); (B) Linear thermodynamic H-H model fitted to the T15 and T35 data set combined (T15, T35); (C) Linear thermodynamic H-H model fitted to the T25 and T35 data set combined (T25, T35); (D) Non-linear thermodynamic H-H model fitted to the T15 and T25 data set combined (T15, T25); (E) Non-linear thermodynamic H-H model fitted to the T15 and T35 data set combined (T15, T35); (F) Non-linear thermodynamic H-H model fitted to the T25 and T35 data set combined (T25, T35). For all panels: the x-axis presents time (ms), the y-axis presents K+ current (µA).
Computation 02 00047 g003
Figure 4. The LINEAR thermodynamic Hodgkin-Huxley (H-H) model, with its free parameters fitted to one temperature data set, cannot replicate the data of a different temperature set, when its temperature terms are changed to the new temperature. So, the model’s temperature terms do not permit it to replicate data at temperatures that it has not been fitted to. (A) Experimental K+ current data recorded at 15 °C (black circles). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 15 °C, and its temperature terms set to 15 °C. The model describes the data well; (B) Experimental K+ current data recorded at 15 °C (black circles). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 25 °C, and its temperature terms set to 15 °C. The model does not describe the data well; (C) Experimental K+ current data recorded at 15 °C (black circles). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 35 °C, and its temperature terms set to 15 °C. The model does not describe the data well; (D) Experimental K+ current data recorded at 25 °C (black). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 15 °C, and its temperature terms set to 25 °C. The model does not describe the data well; (E) Experimental K+ current data recorded at 25 °C (black circles). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 25 °C, and its temperature terms set to 25 °C. The model describes the data well; (F) Experimental K+ current data recorded at 25 °C (black circles). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 35 °C, and its temperature terms set to 25 °C. The model does not describe the data well; (G) Experimental K+ current data recorded at 35 °C (black circles). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 15 °C, and its temperature terms set to 35 °C. The model does not describe the data well; (H) Experimental K+ current data recorded at 35 °C (black circles). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 25 °C, and its temperature terms set to 35 °C. The model does not describe the data well; (I) Experimental K+ current data recorded at 35 °C (black circles). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 35 °C, and its temperature terms set to 35 °C. The model describes the data well. For all panels: the x-axis presents time (mS), the y-axis presents K+ current (µA).
Figure 4. The LINEAR thermodynamic Hodgkin-Huxley (H-H) model, with its free parameters fitted to one temperature data set, cannot replicate the data of a different temperature set, when its temperature terms are changed to the new temperature. So, the model’s temperature terms do not permit it to replicate data at temperatures that it has not been fitted to. (A) Experimental K+ current data recorded at 15 °C (black circles). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 15 °C, and its temperature terms set to 15 °C. The model describes the data well; (B) Experimental K+ current data recorded at 15 °C (black circles). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 25 °C, and its temperature terms set to 15 °C. The model does not describe the data well; (C) Experimental K+ current data recorded at 15 °C (black circles). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 35 °C, and its temperature terms set to 15 °C. The model does not describe the data well; (D) Experimental K+ current data recorded at 25 °C (black). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 15 °C, and its temperature terms set to 25 °C. The model does not describe the data well; (E) Experimental K+ current data recorded at 25 °C (black circles). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 25 °C, and its temperature terms set to 25 °C. The model describes the data well; (F) Experimental K+ current data recorded at 25 °C (black circles). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 35 °C, and its temperature terms set to 25 °C. The model does not describe the data well; (G) Experimental K+ current data recorded at 35 °C (black circles). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 15 °C, and its temperature terms set to 35 °C. The model does not describe the data well; (H) Experimental K+ current data recorded at 35 °C (black circles). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 25 °C, and its temperature terms set to 35 °C. The model does not describe the data well; (I) Experimental K+ current data recorded at 35 °C (black circles). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 35 °C, and its temperature terms set to 35 °C. The model describes the data well. For all panels: the x-axis presents time (mS), the y-axis presents K+ current (µA).
Computation 02 00047 g004

3.2. Temperature Extrapolation of the Model

The ability of a model to describe data not used in determining its parameters is an independent test of how well that model approximates reality. Can the thermodynamic H-H model, with its free parameters fitted to data at one temperature, replicate data at a different temperature upon the temperature terms of the model (in Equations (2)–(5)) being changed to the new temperature? We found that it could not.
Figure 4 shows linear model functions against experimental data. Figure 5 shows non- linear model functions against experimental data. In each panel, the plotted function (red line) utilises best-fit free parameter values derived at one temperature to predict the data set recorded at a second temperature (black circles). Figure 6 shows plotted thermodynamic functions (red lines), utilising best-fit free parameter values derived at two temperatures, to predict a data set at a third temperature (black circles). Figure 7 consolidates all the findings; it presents reduced chi-square values, which indicate the agreement between model and data [9]. The thermodynamic H-H model (linear and non-linear) can describe the data that it has been fitted to well; much better than other data, recorded at other temperatures. Indeed, relatively speaking, it cannot describe data at other temperatures very well at all. Ceteris paribus, the models describe data at higher temperatures worse, possibly due to data variance scaling with temperature. All other things being equal, the model’s performance can be improved by fitting it to data recorded at two temperatures, rather than just one; but this effect is slight. The non-linear model can describe the data that it has been fitted to better than the linear model; but it offers little superiority in describing other data sets, recorded at other temperatures.
The data suggests that the thermodynamic H-H model (linear or non-linear) cannot accurately predict data at temperatures that it is not specifically fitted to i.e., it cannot extrapolate for temperature.
Figure 5. The NON-LINEAR thermodynamic Hodgkin-Huxley (H-H) model, with its free parameters fitted to one temperature data set, cannot replicate the data of a different temperature set, when its temperature terms are changed to the new temperature. So, the model’s temperature terms do not permit it to replicate data at temperatures that it has not been fitted to. (A) Experimental K+ current data recorded at 15 °C (black circles). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 15 °C, and its temperature terms set to 15 °C. The model describes the data well; (B) Experimental K+ current data recorded at 15 °C (black circles). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 25 °C, and its temperature terms set to 15 °C. The model does not describe the data well; (C) Experimental K+ current data recorded at 15 °C (black circles). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 35 °C, and its temperature terms set to 15 °C. The model does not describe the data well; (D) Experimental K+ current data recorded at 25 °C (black). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 15 °C, and its temperature terms set to 25 °C. The model does not describe the data well; (E) Experimental K+ current data recorded at 25 °C (black circles). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 25 °C, and its temperature terms set to 25 °C. The model describes the data well; (F) Experimental K+ current data recorded at 25 °C (black circles). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 35 °C, and its temperature terms set to 25 °C. The model does not describe the data well; (G) Experimental K+ current data recorded at 35 °C (black circles). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 15 °C, and its temperature terms set to 35 °C. The model does not describe the data well; (H) Experimental K+ current data recorded at 35 °C (black circles). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 25 °C, and its temperature terms set to 35 °C. The model does not describe the data well; (I) Experimental K+ current data recorded at 35 °C (black circles). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 35 °C, and its temperature terms set to 35 °C. The model describes the data well. For all panels: the x-axis presents time (mS), the y-axis presents K+ current (µA).
Figure 5. The NON-LINEAR thermodynamic Hodgkin-Huxley (H-H) model, with its free parameters fitted to one temperature data set, cannot replicate the data of a different temperature set, when its temperature terms are changed to the new temperature. So, the model’s temperature terms do not permit it to replicate data at temperatures that it has not been fitted to. (A) Experimental K+ current data recorded at 15 °C (black circles). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 15 °C, and its temperature terms set to 15 °C. The model describes the data well; (B) Experimental K+ current data recorded at 15 °C (black circles). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 25 °C, and its temperature terms set to 15 °C. The model does not describe the data well; (C) Experimental K+ current data recorded at 15 °C (black circles). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 35 °C, and its temperature terms set to 15 °C. The model does not describe the data well; (D) Experimental K+ current data recorded at 25 °C (black). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 15 °C, and its temperature terms set to 25 °C. The model does not describe the data well; (E) Experimental K+ current data recorded at 25 °C (black circles). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 25 °C, and its temperature terms set to 25 °C. The model describes the data well; (F) Experimental K+ current data recorded at 25 °C (black circles). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 35 °C, and its temperature terms set to 25 °C. The model does not describe the data well; (G) Experimental K+ current data recorded at 35 °C (black circles). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 15 °C, and its temperature terms set to 35 °C. The model does not describe the data well; (H) Experimental K+ current data recorded at 35 °C (black circles). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 25 °C, and its temperature terms set to 35 °C. The model does not describe the data well; (I) Experimental K+ current data recorded at 35 °C (black circles). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 35 °C, and its temperature terms set to 35 °C. The model describes the data well. For all panels: the x-axis presents time (mS), the y-axis presents K+ current (µA).
Computation 02 00047 g005
Figure 6. The thermodynamic Hodgkin-Huxley (H-H) model (linear and non-linear), with its free parameters fitted to data recorded at two different temperatures, does not replicate data recorded at a third temperature, when its temperature terms are changed to this third temperature. So, the model’s temperature terms do not permit it to replicate data at temperatures that it has not been fitted to. (A) Experimental K+ current data recorded at 15 °C (black circles). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 15 °C and 25 °C, and its temperature terms set to 15 °C. The model does not describe the data well; (B) Experimental K+ current data recorded at 25 °C (black circles). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 15 °C and 35 °C, and its temperature terms set to 25 °C. The model does not describe the data well; (C) Experimental K+ current data recorded at 35 °C (black circles). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 15 °C and 25 °C, and its temperature terms set to 35 °C. The model does not describe the data well; (D) Experimental K+ current data recorded at 15 °C (black). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 25 °C and 35 °C, and its temperature terms set to 15 °C. The model does not describe the data well; (E) Experimental K+ current data recorded at 25 °C (black circles). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 15 °C and 25 °C, and its temperature terms set to 25 °C. The model does not describe the data well; (F) Experimental K+ current data recorded at 35 °C (black circles). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 15 °C and 25 °C, and its temperature terms set to 35 °C. The model does not describe the data well. For all panels: the x-axis presents time (mS), the y-axis presents K+ current (µA).
Figure 6. The thermodynamic Hodgkin-Huxley (H-H) model (linear and non-linear), with its free parameters fitted to data recorded at two different temperatures, does not replicate data recorded at a third temperature, when its temperature terms are changed to this third temperature. So, the model’s temperature terms do not permit it to replicate data at temperatures that it has not been fitted to. (A) Experimental K+ current data recorded at 15 °C (black circles). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 15 °C and 25 °C, and its temperature terms set to 15 °C. The model does not describe the data well; (B) Experimental K+ current data recorded at 25 °C (black circles). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 15 °C and 35 °C, and its temperature terms set to 25 °C. The model does not describe the data well; (C) Experimental K+ current data recorded at 35 °C (black circles). Plots of the linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 15 °C and 25 °C, and its temperature terms set to 35 °C. The model does not describe the data well; (D) Experimental K+ current data recorded at 15 °C (black). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 25 °C and 35 °C, and its temperature terms set to 15 °C. The model does not describe the data well; (E) Experimental K+ current data recorded at 25 °C (black circles). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 15 °C and 25 °C, and its temperature terms set to 25 °C. The model does not describe the data well; (F) Experimental K+ current data recorded at 35 °C (black circles). Plots of the non-linear thermodynamic H-H model (red lines), with its free parameters having been previously fitted to data recorded at 15 °C and 25 °C, and its temperature terms set to 35 °C. The model does not describe the data well. For all panels: the x-axis presents time (mS), the y-axis presents K+ current (µA).
Computation 02 00047 g006
Figure 7. Logged (log10), reduced chi-square values indicating if the linear and non-linear thermodynamic H-H models—with their parameters optimised to data recorded at one (or more) temperature(s)—can describe different data recorded at a different temperature, if their temperature terms are set to this new temperature. (A) Bars show how well the linear variant describes K+ current data recorded at 15 °C, if its free parameters have been fitted to this data (black bar) or if they have been fitted to different data recorded at 25 °C (turquoise bar), or 35 °C (blue bar) or (25 °C and 35 °C) (olive bar); (B) Bars show how well the linear variant describes K+ current data recorded at 25 °C, if its free parameters have been fitted to this data (turquoise bar) or if they have been fitted to different data recorded at 15 °C (black bar), or 35 °C (blue bar) or (15 °C and 35 °C) (wine bar); (C) Bars show how well the linear variant describes K+ current data recorded at 35 °C, if its free parameters have been fitted to this data (blue bar) or if they have been fitted to different data recorded at 15 °C (black bar), or 25 °C (turquoise bar) or (15 °C and 25 °C) (purple bar); (D) Bars show how well the non-linear variant describes K+ current data recorded at 15 °C, if its free parameters have been fitted to this data (black bar) or if they have been fitted to different data recorded at 25 °C (turquoise bar), or 35 °C (blue bar) or (25 °C and 35 °C) (olive bar); (E) Bars show how well the non-linear variant describes K+ current data recorded at 25 °C, if its free parameters have been fitted to this data (turquoise bar) or if they have been fitted to different data recorded at 15 °C (black bar), or 35 °C (blue bar) or (15 °C and 35 °C) (wine bar); (F) Bars show how well the non-linear variant describes K+ current data recorded at 35 °C, if its free parameters have been fitted to this data (blue bar) or if they have been fitted to different data recorded at 15 °C (black bar), or 25 °C (turquoise bar) or (15 °C and 25 °C) (purple bar). All bars are scaled by the presented y-scale bar.
Figure 7. Logged (log10), reduced chi-square values indicating if the linear and non-linear thermodynamic H-H models—with their parameters optimised to data recorded at one (or more) temperature(s)—can describe different data recorded at a different temperature, if their temperature terms are set to this new temperature. (A) Bars show how well the linear variant describes K+ current data recorded at 15 °C, if its free parameters have been fitted to this data (black bar) or if they have been fitted to different data recorded at 25 °C (turquoise bar), or 35 °C (blue bar) or (25 °C and 35 °C) (olive bar); (B) Bars show how well the linear variant describes K+ current data recorded at 25 °C, if its free parameters have been fitted to this data (turquoise bar) or if they have been fitted to different data recorded at 15 °C (black bar), or 35 °C (blue bar) or (15 °C and 35 °C) (wine bar); (C) Bars show how well the linear variant describes K+ current data recorded at 35 °C, if its free parameters have been fitted to this data (blue bar) or if they have been fitted to different data recorded at 15 °C (black bar), or 25 °C (turquoise bar) or (15 °C and 25 °C) (purple bar); (D) Bars show how well the non-linear variant describes K+ current data recorded at 15 °C, if its free parameters have been fitted to this data (black bar) or if they have been fitted to different data recorded at 25 °C (turquoise bar), or 35 °C (blue bar) or (25 °C and 35 °C) (olive bar); (E) Bars show how well the non-linear variant describes K+ current data recorded at 25 °C, if its free parameters have been fitted to this data (turquoise bar) or if they have been fitted to different data recorded at 15 °C (black bar), or 35 °C (blue bar) or (15 °C and 35 °C) (wine bar); (F) Bars show how well the non-linear variant describes K+ current data recorded at 35 °C, if its free parameters have been fitted to this data (blue bar) or if they have been fitted to different data recorded at 15 °C (black bar), or 25 °C (turquoise bar) or (15 °C and 25 °C) (purple bar). All bars are scaled by the presented y-scale bar.
Computation 02 00047 g007

4. Conclusions

The Hodgkin-Huxley (H-H) model of voltage-dependent conductances, for which Hogkin and Huxley were awarded the Nobel Prize in 1963, was and is a triumph. However, their model cannot intrinsically account for temperature. A modified variant of their model has been proposed, the thermodynamic H-H model, which suggests that channel transition rates depend on a free energy barrier by analogy with reaction rates. In this study, we investigate if this variant can intrinsically account for temperature unlike the original H-H description, which can’t. We hypothesise that it might because this variant—unlike the original H-H form—has temperature terms in its equations.
We conclude that the temperature terms of the thermodynamic H-H model (linear or non-linear) do not permit it to describe voltage-gated K+ current data, over a range of temperatures, with a single set of free parameter values. The model cannot extrapolate for temperature. Why this failure?
Eyring rate theory, which underpins the thermodynamic H-H formalisms, assumes that temperature has little or no effect on the relative energies of the states (Closed, Transition, Open). Its tenet is that temperature changes the probability of moving between states, it does not significantly change the states themselves [4]. However, an ion channel is likely to breach this assumption. Protein structure is temperature dependent [10,11,12] and membrane lipid structure is temperature dependent. The latter is relevant because it can change ion channel structure [13,14] and also because voltage-gated channels are gated by “voltage sensor paddles”, which float and move in the lipid medium [15].
Our study only investigates one K+ current data set and so the conclusions we draw are somewhat provisional. We hope that others will build on this work and follow our methodology with other, further data sets; because our hypothesis is an important issue to address.
The Hodgkin-Huxley model, in all its forms, is an abstraction. The reality is better approximated by a Markov model [2]. It can represent different voltage-gated channel states, and how these states change with voltage and time. They can be accurate and powerful but are computationally expensive to simulate (e.g., [16,17]). So, neuron modelling studies typically use the H-H model as opposed to Markov descriptions (e.g., [18,19]). Hence the H-H model is still relevant and important [20]. Presently, thermodynamic H-H variants are not typically to be found in neuron modelling studies. However, if they can be found to provide intrinsic temperature tenability, then their incorporation will increase the predictive power of neuron models immensely.

Acknowledgments

The author conducted this work in 2003, as a component of his Master degree at the University of Edinburgh [21]. He was funded on this course by the Medical Research Council (MRC) of the UK. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. This work was reviewed by the author, and written up for publication, in 2012 (with no funding appropriated for this task). The author would like to thank David C. Sterratt and Andrew Gillies at the University of Edinburgh who conceived this investigation and who provided much help and guidance; they were the author’s project supervisors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hodgkin, A.L.; Huxley, A.F. A quantitative description of membrane current and its application to conduction and excitation in nerve. J. Physiol. (Lond.) 1952, 117, 500–544. [Google Scholar]
  2. Hille, B. Ion Channels of Excitable Membranes, 3rd ed.; Sinauer Associates: Sunderland, MA, USA, 2001. [Google Scholar]
  3. Eyring, H. The activated complex in chemical reactions. J. Chem. Phys. 1935, 3, 107–115. [Google Scholar] [CrossRef]
  4. Tsien, R.W.; Noble, D. A transition state theory approach to kinetics of conductance changes in excitable membranes. J. Membr. Biol. 1969, 1, 248–273. [Google Scholar] [CrossRef]
  5. Hill, T.L.; Chen, Y.-D. On the theory of ion transport across the nerve membrane. VI. Free energy and activation free energies of conformational change. Proc. Natl. Acad. Sci. USA 1972, 69, 1723–1726. [Google Scholar]
  6. Destexhe, A.; Huguenard, J.R. Which formalism to use for modelling voltage-dependent conductances? In Computational Neuroscience: Realistic modelling for experimentalists; CRC Press: Boca Raton, FL, USA, 2000. [Google Scholar]
  7. Destexhe, A.; Huguenard, J.R. Nonlinear thermodynamic models of voltage dependent currents. J. Comput. Neurosci. 2000, 9, 259–270. [Google Scholar] [CrossRef]
  8. Tiwari, J.K.; Sikdar, S.K. Temperature dependent conformational changes in a voltage-gated potassium channel. Eur. Biophys. J. 1999, 28, 338–345. [Google Scholar] [CrossRef]
  9. Press, W.H.; Flannery, B.P.; Teukolsky, S.A.; Vetterling, W.T. Numerical Recipes in C: The Art of Scientific Computing; Cambridge University Press: Cambridge, UK, 1992. [Google Scholar]
  10. Creighton, T.E. Proteins; WH Freeman and Company: New York, NY, USA, 1984. [Google Scholar]
  11. Alberts, B.; Bray, D.; Lewis, J.; Raff, M.; Roberts, K.; Watson, J.D. Molecular Biology of the Cell; Garland Science: New York, NY, USA, 1994. [Google Scholar]
  12. Stryer, L. Biochemistry; WH Freeman and Company: New York, NY, USA, 1995. [Google Scholar]
  13. Chapman, D. Phase transitions and fluidity characteristics of lipids and cell membranes. Q. Rev. Biophys. 1975, 8, 185–235. [Google Scholar] [CrossRef]
  14. Rosen, A.D. Nonlinear temperature modulation of sodium channel kinetics in GH3 cells. Biochim. Biophys. Acta 2001, 1511, 391–396. [Google Scholar] [CrossRef]
  15. Jiang, Y.; Ruta, V.; Chen, J.; Lee, A.; MacKinnon, R. The principle of gating charge movement in a voltage-dependent K+ channel. Nature 2003, 423, 42–48. [Google Scholar] [CrossRef]
  16. Vandenberg, C.A.; Bezanilla, F. A model of sodium channel gating based on single channel, macroscopic ionic, and gating currents in the squid giant axon. Biophys. J. 1991, 60, 1511–1533. [Google Scholar] [CrossRef]
  17. Perozo, E.; Bezanilla, F. Phosphorylation affects voltage gating of the delayed rectifier K+ channel by electrostatic interactions. Neuron 1990, 5, 685–690. [Google Scholar] [CrossRef]
  18. Forrest, M.D.; Wall, M.J.; Press, D.A.; Feng, J. The Sodium-Potassium Pump Controls the Intrinsic Firing of the Cerebellar Purkinje Neuron. PLoS One 2012, 7, e51169. [Google Scholar] [CrossRef]
  19. Forrest, M.D. Mathematical Model of Bursting in Dissociated Purkinje Neurons. PLoS One 2013, 8, e68765. [Google Scholar] [CrossRef]
  20. Vandenberg, J.I.; Waxman, S.G. Hodgkin and Huxley and the basis for electrical signalling: A remarkable legacy still going strong. J. Physiol. 2012, 590, 2569–2570. [Google Scholar] [CrossRef]
  21. Forrest, M.D. Can Thermodynamic Models of Voltage-Dependent Conductances Extrapolate for Temperature? Master’s Thesis, September 2003. [Google Scholar]

Share and Cite

MDPI and ACS Style

Forrest, M.D. Can the Thermodynamic Hodgkin-Huxley Model of Voltage-Dependent Conductance Extrapolate for Temperature? Computation 2014, 2, 47-60. https://0-doi-org.brum.beds.ac.uk/10.3390/computation2020047

AMA Style

Forrest MD. Can the Thermodynamic Hodgkin-Huxley Model of Voltage-Dependent Conductance Extrapolate for Temperature? Computation. 2014; 2(2):47-60. https://0-doi-org.brum.beds.ac.uk/10.3390/computation2020047

Chicago/Turabian Style

Forrest, Michael D. 2014. "Can the Thermodynamic Hodgkin-Huxley Model of Voltage-Dependent Conductance Extrapolate for Temperature?" Computation 2, no. 2: 47-60. https://0-doi-org.brum.beds.ac.uk/10.3390/computation2020047

Article Metrics

Back to TopTop