GHK flux equation

Not to be confused with Goldman equation or Goldman–Hodgkin–Katz voltage equation.

The Goldman–Hodgkin–Katz flux equation (or GHK flux equation or GHK current density equation) describes the ionic flux across a cell membrane as a function of the transmembrane potential and the concentrations of the ion inside and outside of the cell. Since both the voltage and the concentration gradients influence the movement of ions, this process is a simplified version of electrodiffusion. Electrodiffusion is most accurately defined by the Nernst–Planck equation and the GHK flux equation is a solution to the Nernst–Planck equation with the assumptions listed below.

Origin

The American David E. Goldman of Columbia University, and the English Nobel laureates Alan Lloyd Hodgkin and Bernard Katz derived this equation.

Assumptions

Several assumptions are made in deriving the GHK flux equation (Hille 2001, p. 445) :

The membrane is a homogeneous substance
The electrical field is constant so that the transmembrane potential varies linearly across the membrane
The ions access the membrane instantaneously from the intra- and extracellular solutions
The permeant ions do not interact
The movement of ions is affected by both concentration and voltage differences

Equation

The GHK flux equation for an ion S (Hille 2001, p. 445):

\Phi_{S} = P_{S}z_{S}^2\frac{V_{m}F^{2}}{RT}\frac{[\mbox{S}]_{i} - [\mbox{S}]_{o}\exp(-z_{S}V_{m}F/RT)}{1 - \exp(-z_{S}V_{m}F/RT)}

where

$\Phi$ _S is the current density (flux) across the membrane carried by ion S, measured in amperes per square meter (A·m⁻²)
P_S is the permeability of the membrane for ion S measured in m·s⁻¹
z_S is the valence of ion S
V_m is the transmembrane potential in volts
F is the Faraday constant, equal to 96,485 C·mol⁻¹ or J·V⁻¹·mol⁻¹
R is the gas constant, equal to 8.314 J·K⁻¹·mol⁻¹
T is the absolute temperature, measured in Kelvin (= degrees Celsius + 273.15)
[S]_i is the intracellular concentration of ion S, measured in mol·m⁻³ or mmol·l⁻¹
[S]_o is the extracellular concentration of ion S, measured in mol·m⁻³

Implicit definition of reversal potential

The reversal potential is shown to be contained in the GHK flux equation (Flax 2008). The proof is replicated from the reference (Flax 2008) here.

We wish to show that when the flux is zero, the transmembrane potential is not zero. Formally it is written $\lim_{\Phi_{S}\rightarrow0} V_{m}\ne0$ which is equivalent to writing $\lim_{V_{m}\rightarrow0} \Phi_{S}\ne0$ , which states that when the transmembrane potential is zero, the flux is not zero.

However due to the form of the GHK flux equation when $V_{m}=0$ , $\Phi_{S}=\frac{0}{0}$ . This is a problem as the value of $\frac{0}{0}$ is indeterminate.

We turn to l'Hôpital's rule to find the solution for the limit:

\lim_{V_{m}\rightarrow0} \Phi_{S} = P_{S}\frac{z_{S}^2F^{2}}{RT}\frac{[V_{m}([\mbox{S}]_{i} - [\mbox{S}]_{o}\exp(-z_{S}V_{m}F/RT))]'}{[1 - \exp(-z_{S}V_{m}F/RT)]'}

where $[f]'$ represents the differential of f and the result is :

\lim_{V_{m}\rightarrow0} \Phi_{S} = P_{S}z_{S}F([\mbox{S}]_{i} - [\mbox{S}]_{o})

It is evident from the previous equation that when $V_{m}=0$ , $\Phi_{S}\ne0$ if $([\mbox{S}]_{i} - [\mbox{S}]_{o})\ne0$ and thus

\lim_{\Phi_{S}\rightarrow0}V_{m}\ne0

which is the definition of the reversal potential.

By setting $\Phi_{S}=0$ we can also obtain the reversal potential :

\Phi_{S}=0=P_{S}\frac{z_{S}^2F^{2}}{RT}\frac{V_{m}([\mbox{S}]_{i} - [\mbox{S}]_{o}\exp(-z_{S}V_{m}F/RT))}{1 - \exp(-z_{S}V_{m}F/RT)}

which reduces to :

[\mbox{S}]_{i} - [\mbox{S}]_{o}\exp(-z_{S}V_{m}F/RT)=0

and produces the Nernst equation :

V_{m}=-\frac{RT}{z_{S}F}\ln\left (\frac{[\mbox{S}]_{i}}{[\mbox{S}]_{o}}\right )

Rectification

Since one of the assumptions of the GHK flux equation is that the ions move independently of each other, the total flow of ions across the membrane is simply equal to the sum of two oppositely directed fluxes. Each flux approaches an asymptotic value as the membrane potential diverges from zero. These asymptotes are

\Phi_{S|i\to o} = P_{S}z_{S}^2 \frac{V_{m}F^{2}}{RT}[\mbox{S}]_{i}\ \mbox{for}\ V_{m} \gg \; 0

and

\Phi_{S|o\to i} = P_{S}z_{S}^2 \frac{V_{m}F^{2}}{RT}[\mbox{S}]_{o}\ \mbox{for}\ V_{m} \ll \; 0

where subscripts 'i' and 'o' denote the intra- and extracellular compartments, respectively. Keeping all terms except V_m constant, the equation yields a straight line when plotting $\Phi$ _S against V_m. It is evident that the ratio between the two asymptotes is merely the ratio between the two concentrations of S, [S]_i and [S]_o. Thus, if the two concentrations are identical, the slope will be identical (and constant) throughout the voltage range (corresponding to Ohm's law scaled by the surface area). As the ratio between the two concentrations increases, so does the difference between the two slopes, meaning that the current is larger in one direction than the other, given an equal driving force of opposite signs. This is contrary to the result obtained if using Ohm's law scaled by the surface area, and the effect is called rectification.

The GHK flux equation is mostly used by electrophysiologists when the ratio between [S]_i and [S]_o is large and/or when one or both of the concentrations change considerably during an action potential. The most common example is probably intracellular calcium, [Ca²⁺]_i, which during a cardiac action potential cycle can change 100-fold or more, and the ratio between [Ca²⁺]_o and [Ca²⁺]_i can reach 20,000 or more.

References

Hille, Bertil (2001) Ion channels of excitable membranes, 3rd ed., Sinauer Associates, Sunderland, Massachusetts. ISBN 978-0-87893-321-1
Flax, Matt R. and Holmes, W.Harvey (2008) Goldman-Hodgkin-Katz Cochlear Hair Cell Models - a Foundation for Nonlinear Cochlear Mechanics, Conference proceedings: Interspeech 2008