Calculate ion equilibrium potentials and membrane potential using the Nernst and GHK equations. Real-time simulation of action potential waveforms and threshold effects with a simplified Hodgkin-Huxley model.
What exactly is a "membrane potential"? I see the simulator shows a voltage, but what's physically happening?
🎓
Basically, it's the voltage difference across a cell's membrane, like a tiny biological battery. It's created because there are different concentrations of ions (like potassium K⁺ and sodium Na⁺) inside versus outside the cell. Try moving the "[K]_in" and "[K]_out" sliders above—you'll see the voltage change instantly as you change the concentration gradient.
🙋
Wait, really? So if the concentrations are equal, the voltage is zero? But the simulator shows a negative voltage even when I set them equal. Why?
🎓
Great observation! That's because a real cell membrane isn't equally permeable to all ions. In practice, it's much more leaky to potassium (K⁺) than sodium (Na⁺) at rest. That's what the "Permeability ratio P_Na/P_K" control models. Set it to 1 (equal permeability) and you'll see the voltage shift dramatically. The resting potential is a weighted average of all ion forces.
🙋
Okay, that makes sense for a steady voltage. But what causes the sudden spike—the "action potential"—when I click "Stimulate"?
🎓
That's the magic of voltage-gated ion channels! When your stimulus current (controlled by the "I_stim" slider) pushes the membrane past a threshold, special sodium channels snap open. Sodium rushes in, causing the sharp upward spike. Then potassium channels open to repolarize it. The simulator models this dynamic sequence, letting you see how changing ion concentrations affects the spike's shape and timing.
Physical Model & Key Equations
The Nernst Equation calculates the equilibrium potential for a single ion species. This is the voltage at which the electrical force exactly balances the diffusion force from the concentration gradient.
Where $E_{ion}$ is the Nernst potential (Volts), $R$ is the gas constant, $T$ is absolute Temperature (Kelvin), $z$ is the ion's charge (e.g., +1 for K⁺), $F$ is Faraday's constant, and $[ion]_{in/out}$ are the intracellular and extracellular concentrations.
Since the membrane is permeable to multiple ions at once, the Goldman-Hodgkin-Katz (GHK) voltage equation gives the resting membrane potential $V_m$. It's a weighted average of the Nernst potentials, with each ion's permeability $P$ as the weight.
Here, $P_K$, $P_{Na}$, and $P_{Cl}$ are the absolute permeabilities. In the simulator, you control the ratios $P_{Na}/P_K$ and $P_{Cl}/P_K$. This equation is foundational for predicting how changes in permeability—like during an action potential—alter the voltage.
Frequently Asked Questions
The Nernst equation calculates the equilibrium potential for a single ion, while the GHK equation accounts for the permeability of multiple ions. Since the actual membrane potential is influenced by multiple ion species, the GHK equation is more physiologically accurate. Try adjusting the permeability ratios (e.g., P_K:P_Na) and compare the results.
The threshold depends on the initial state of the membrane and ion channel parameters (e.g., maximum conductance, time constants). For example, increasing the external K⁺ concentration depolarizes the resting membrane potential, making it easier to reach the threshold. Observe waveform changes by adjusting each parameter with the sliders.
Both the Nernst and GHK equations include the absolute temperature T. As temperature increases, the RT/F term becomes larger, increasing the absolute value of the equilibrium potential for the same concentration gradient. Additionally, in the Hodgkin-Huxley model, the temperature coefficient Q10 alters channel gating kinetics, affecting the action potential waveform.
Cl⁻ is an anion with a charge of -1, so -1 is substituted for z in the Nernst equation. This inverts the sign of the logarithmic term, meaning that a higher external concentration results in a more negative equilibrium potential. If the valence is entered incorrectly, the sign of the reversal potential will be reversed, so ensure the correct valence when selecting the ion species.
Real-World Applications
Deep Brain Stimulation (DBS) Parameter Design: DBS treats Parkinson's disease by delivering electrical pulses to specific brain nuclei. Engineers use these exact equations in computer models to simulate how the stimulus current affects the membrane potential of neurons around the electrode. This helps optimize pulse shape, frequency, and amplitude to maximize therapeutic effect while minimizing side effects.
Cochlear and Retinal Implants: These neuroprosthetics convert sound or light into electrical signals to stimulate auditory or optic nerves. The models predict the required stimulus current to reliably trigger action potentials in the target neurons, which is critical for designing the implant's electrode array and signal processing strategy to create clear sensory perception.
Cardiac Electrophysiology & Arrhythmia Modeling: Abnormal heart rhythms (arrhythmias) arise from disruptions in the coordinated action potentials of cardiac cells. Biomedical engineers create detailed "in-silico" heart models using these principles to locate the source of arrhythmias and to plan surgical or ablation procedures, all within a virtual simulation before touching a patient.
BCI Neural Interfaces & Patch-Clamp Validation: In Brain-Computer Interfaces (BCIs), it's vital to interpret the electrical signals from neurons. These models help decode neural activity. Furthermore, the simulator acts as a virtual patch-clamp experiment, allowing researchers to validate their experimental setups and hypotheses about ion channel behavior before conducting costly and time-consuming lab work.
Common Misconceptions and Points to Note
When you start using this simulator, there are a few common pitfalls. First is memorizing that "the resting membrane potential is always more positive than the K⁺ equilibrium potential". That's true for typical neurons. But try increasing the P_Cl/P_K value in the simulator. When chloride permeability is high, the membrane potential is pulled toward E_Cl (usually around -70mV). In some cases, the calculated resting potential can actually become more negative than E_K. It's important to play with the tool and break your preconceptions of "how things should be".
Next, the "sensitivity" of the membrane potential to slider changes is not uniform. For example, at rest (P_Na/P_K=0.01), raising [K]_out from 5mM to 10mM shifts the membrane potential significantly, from about -80mV to -65mV. However, making the same change to [Na]_out barely moves it. This is because the membrane is selective for K⁺. When adjusting parameters, always be aware of which ion is dominant.
Finally, the misconception that "the action potential waveform does not change with stimulus strength". If you actually try it, you'll see that setting I_stim just at threshold delays spike generation and results in a more gradual waveform. Conversely, applying a strong current pulse causes almost instantaneous firing and a sharp spike. Although the model in the simulator is simplified, the same phenomenon occurs in real neurons. Use the tool to experience how "threshold" is not an absolute value but an "event" dependent on how the stimulus is applied.