GENESIS: Documentation

Related Documentation:

Further Details of the Hodgkin-Huxley Model

In their model, Hodgkin and Huxley represented the rate of change of the potassium activation variable n as a rate equation for a first-order kinetic process:

---=  αn(V ) (1 - n) - βn(V ) n.

The modern interpretation of this equation is that the opening of a K channel involves the motion of four physical gates between “permissive” and “non-permissive” states. If all gates are in the permissive state, then the channel is open, and ions are allowed to flow. In a population of channels, n is the fraction of gates that are open. Thus, the conductance of these of channels is proportional to the fourth power of n. Here, αn is the voltage-dependent rate constant for transitions to the permissive state, and βn is the rate constant for transitions to the non-permissive state. By comparison with Eq. 4, we can see that

αn(V ) =  τ (V)

βn(V ) =    τ (V)

Likewise, the Na activation and inactivation variables obey the equations

----= αm(V  ) (1 - m) - βm(V ) m

--- = αh(V ) (1 - h) - βh(V ) h

with similar relationships between the rate constants and the activation variables and their time constants.

From their voltage clamp measurements of the activation variables and their time constants, Hodgkin and Huxley were able to make the empirical fits:

          0.01(10 - V )
αn(V  ) = ----10-V------
          exp( 10  ) - 1

βn(V ) = 0.125 exp( - V ∕80)

αm(V  ) = exp(25-V-) - 1

βm(V ) = 4 exp(- V∕18)

αh(V ) = 0.07 exp(- V ∕20)

β (V) =  ------1-------
 h       exp(30-V-) + 1

Here voltages are measured with respect to a resting potential of zero (rather than the actual value of about -70 mV), and are given in mV . With these expressions for the rate constants, one can numerically solve Eqs. 69, and 10 for the time dependence of the activation variables, and use these in Eqs. 2 and 5 to solve for the conductances. These can then be used in Eq. 1 to solve for the time behavior of the membrane potential.

As the Traub hippocampal pyramidal cell model shows, there are many other varieties of voltage-activated conductances found in neurons, in addition to those found in the squid giant axon. Most of these can be modeled using variations of the Hodgkin-Huxley model, with a similar notation, but with different values of the parameters and different values of the exponents for activation and inactivation. Some modelers use equations similar to Eqs. 1116, and others use expressions for the steady-state activation and time constant, which can be related to these by Eqs. 7 and 8.

Further Reading

It is well worth reading the original paper:

Hodgkin, A. and Huxley, A. (1952). A quantitative description of membrane current and its application to conduction and excitation in nerve, J. Physiol. (London) 117: 500–544.

For the details of implementing models of ionic conductances see (UNDER CONSTRUCTION).