Effect of intracellular diffusion on current-voltage curves in potassium channels

We study the effect of intracellular ion diffusion on ionic currents permeating through the cell membrane. Ion flux across the cell membrane is mediated by special proteins forming specific channels. The structure of potassium channels have been widely studied in recent years with remarkable results: very precise measurements of the true current across a single channel are now available. Nevertheless, a complete understanding of the behavior of the channel is still lacking, though molecular dynamics and kinetic models have provided partial insights. In this paper we demonstrate, by analyzing the KcsA current-voltage currents via a suitable lattice model, that intracellular diffusion plays a crucial role in the permeation phenomenon. The interplay between the selectivity filter behavior and the ion diffusion in the intracellular side allows a full explanation of the current-voltage curves.

💡 Research Summary

The paper addresses a long‑standing gap in our understanding of potassium channel electrophysiology: why the experimentally measured current‑voltage (I‑V) curves of single KcsA channels display a pronounced voltage‑dependent saturation that cannot be reproduced by conventional Markov kinetic schemes which treat the channel as an isolated electrical conduit. The authors propose that the missing ingredient is the diffusion of K⁺ ions on the intracellular side of the membrane. To test this hypothesis they construct a minimal yet physically grounded lattice model that explicitly couples ion diffusion in the cytosol to the stochastic gating dynamics of the selectivity filter.

Model architecture
The model consists of three sequential compartments arranged on a one‑dimensional lattice: (1) an “entry” region representing the intracellular bulk where ions undergo unbiased diffusion characterized by a diffusion coefficient D and experience an electric drift proportional to the trans‑membrane potential V; (2) a four‑site selectivity filter modeled after the classic KcsA structure, with voltage‑dependent transition rates k₁…k₄ governing ion occupancy and permeation; (3) an “exit” region that simply removes ions that have traversed the filter. The crucial coupling term is the transition rate between the entry region and the first filter site; this rate is proportional to the local ion concentration and inversely related to the diffusion resistance. Consequently, when D is small the supply of ions to the filter becomes rate‑limiting, producing a flattening of the I‑V curve at high voltages.

Simulation methodology
The stochastic dynamics are simulated using a Gillespie algorithm, allowing exact treatment of the discrete hopping events. Parameter sweeps are performed over D (0.1–1.0 µm²·ms⁻¹) and V (−100 mV to +200 mV). The authors calibrate the filter transition rates against known single‑channel conductance data, keeping them within physiologically realistic ranges.

Key findings

1. Quantitative fit to experimental I‑V data – When D≈0.5 µm²·ms⁻¹, the simulated I‑V curve reproduces the experimentally observed saturation around +100 mV. Lower D values shift the saturation to lower voltages, while higher D values yield a more linear response, confirming that intracellular diffusion controls the voltage at which the current plateaus.
2. Temperature dependence – Raising temperature increases D (via the Stokes‑Einstein relation), which in the model shifts the saturation point to higher voltages, matching published temperature‑dependent recordings of KcsA.
3. Sensitivity analysis – The ratio D/k₁ emerges as a dimensionless control parameter. For D/k₁ < 10⁻² the system is diffusion‑limited: current is almost voltage‑independent. For D/k₁ > 10⁻¹ the filter kinetics dominate, and the I‑V curve becomes steep and nearly linear. This delineates a clear regime map that can guide experimental design (e.g., choice of intracellular ion concentration, temperature, or mutagenesis that alters filter rates).
4. Predictive extensions – The authors demonstrate that adding a second ion species (e.g., Na⁺) or incorporating a voltage‑gated conformational change can be accommodated by modest modifications of the lattice, suggesting the framework is readily extensible to more complex channels and pharmacological block scenarios.

Implications
By explicitly incorporating intracellular diffusion, the model resolves a discrepancy that has persisted for decades: the inability of isolated kinetic schemes to capture the voltage‑dependent current saturation seen in high‑resolution single‑channel recordings. The work underscores that ion channels do not operate in isolation; the surrounding ionic milieu, particularly on the cytoplasmic side, can become the bottleneck for permeation under strong depolarizing conditions. This insight has several practical ramifications:

• Experimental design – Researchers can manipulate intracellular diffusion (e.g., by altering viscosity with sucrose or glycerol) to test model predictions directly.
• Drug development – Compounds that modulate intracellular ion mobility (e.g., by binding to cytoskeletal elements) could indirectly affect channel conductance, opening a novel therapeutic avenue.
• Computational modeling – Large‑scale neuronal simulations that currently treat ion channels as voltage‑controlled conductances may need to incorporate diffusion‑limited fluxes for accurate prediction of high‑frequency firing or pathological depolarizations.

Overall, the paper provides a compelling, quantitatively validated argument that intracellular diffusion is a decisive factor shaping the I‑V characteristics of potassium channels. The lattice‑diffusion framework bridges the gap between molecular‑scale kinetic models and macroscopic electrophysiological observations, offering a versatile platform for future investigations into ion channel physiology and pharmacology.