The Hodgkin-Huxley equations are a foundational set of nonlinear differential equations that model the electrical activity of excitable cells, particularly neurons. Developed by Alan Hodgkin and Andrew Huxley in 1952, this model earned them the Nobel Prize for its accuracy in describing how action potentials are initiated and propagated through ion channel dynamics.
Formulation
The basic Hodgkin-Huxley model includes the following core components:
- Membrane voltage equation: Describes how membrane potential varies over time due to ionic currents.
- Ion channel equations: Governed by gating variables that determine the open probability of sodium and potassium channels.
Differential Equations
The membrane potential \( V \) evolves according to:
\[ C_m \frac{dV}{dt} = \mminus I_{\text{Na}} \mminus I_{\text{K}} \mminus I_{\text{L}} + I_{\text{ext}} \]
Where:
- \( C_m \) is the membrane capacitance.
- \( I_{\text{Na}} \), \( I_{\text{K}} \), and \( I_{\text{L}} \) are sodium, potassium, and leak currents.
- \( I_{\text{ext}} \) is the externally applied current.
Each ion current is modeled as:
\[ I_x = g_x \cdot m^p h^q (V \mminus E_x) \]
Where:
- \( g_x \): Maximum conductance for ion type \( x \).
- \( m \), \( h \): Gating variables (with exponents \( p \), \( q \)) controlling channel state.
- \( E_x \): Reversal (Nernst) potential for ion \( x \).
Gating Variables
Each gating variable \( x \in \{m, h, n\} \) follows a voltage-dependent first-order kinetic equation:
\[ \frac{dx}{dt} = \alpha_x(V)(1 \mminus x) \mminus \beta_x(V)x \]
Jacobian and Stability Analysis
The Jacobian matrix is constructed from the partial derivatives of the system’s differential equations with respect to each variable. Evaluating this matrix at equilibrium points helps assess the local stability and can guide numerical integration strategies.
Equilibrium Points
To find equilibria, we solve:
\[ \frac{dV}{dt} = 0,\quad \frac{dm}{dt} = 0,\quad \frac{dh}{dt} = 0,\quad \frac{dn}{dt} = 0 \]
Simultaneously solving these equations gives the resting state(s) of the neuron.
Maple Implementation
In Maple, these equations can be encoded symbolically for simulation, bifurcation analysis, or Jacobian evaluation. This facilitates exploration of action potential dynamics and parameter sensitivities.
Gating Summary
- \( m \): Sodium activation gate
- \( h \): Sodium inactivation gate
- \( n \): Potassium activation gate
Voltage-Dependent Rate Constants
-
Sodium activation (\( m \))
\[ \alpha_m(V) = \frac{0.1 (V + 40)}{1 \mminus e^{\mminus 0.1 (V + 40)}},\quad \beta_m(V) = 4 e^{\mminus 0.0556 (V + 65)} \] -
Sodium inactivation (\( h \))
\[ \alpha_h(V) = 0.07 e^{\mminus 0.05 (V + 65)},\quad \beta_h(V) = \frac{1}{1 + e^{\mminus 0.1 (V + 35)}} \] -
Potassium activation (\( n \))
\[ \alpha_n(V) = \frac{0.01 (V + 55)}{1 \mminus e^{\mminus 0.1 (V + 55)}},\quad \beta_n(V) = 0.125 e^{\mminus 0.0125 (V + 65)} \]
These empirically derived rate constants enable the model to replicate the voltage-sensitive gating behavior of ion channels critical to action potential generation.
