Atomic Biology, Electrostatics, and Ionic Channels

I believe an atomic biology is needed to supplement present day molecular biology, if we are to design and understand proteins, as well as define, make, and use them. Topics in the paper are molecular biology and atomic biology. Electrodiffusion in the open channel. Electrodiffusion in mixed electrolytes. Models of permeation. State Models of Permeation are Inconsistent with the Electric Field. Making models in atomic biology. Molecular dynamics. Temporal Limitations; Spatial Limitations; Periodic boundary conditions. Hierarchy of models of the open channel. Stochastic Motion of the Channel. Langevin Dynamics. Simulations of the Reaction Path: the Permion. Chemical reactions. What was wrong? Back to the hierarchy: Occam’s razor can slit your throat. Poisson-Nernst-Planck PNP Models Flux Ratios; Pumping by Field Coupling. Gating in channels of one conformation. Gating by Field Switching; Gating Current; Gating in Branched Channels; Blocking. Back to the hierarchy: Linking levels. Is there a theory? At what level will the adaptation be found? Simplicity, evolution, and natural function.

💡 Research Summary

The manuscript puts forward a bold thesis: contemporary molecular biology, while highly successful in describing sequences, structures, and macroscopic functions of biomolecules, lacks the resolution needed to fully understand systems whose operation is governed by atomic‑scale electrostatics. To fill this gap the author proposes “atomic biology,” a discipline that treats proteins, and especially ion channels, as assemblies of charged atoms whose positions, electronic clouds, and mutual electrostatic interactions are modeled explicitly.

The paper begins with a concise review of electrodiffusion in open channels. When a membrane potential is applied, a simultaneous electric field and concentration gradient drive ion flux. In a single‑electrolyte solution the situation reduces to the classic Nernst‑Planck description, but real biological membranes contain mixed electrolytes (Na⁺, K⁺, Ca²⁺, Cl⁻, etc.) with disparate mobilities and valences. The author derives the coupled Poisson‑Nernst‑Planck equations for this mixed case, showing how charge shielding, selectivity, and field‑induced barrier modulation emerge naturally from the mathematics.

Next, the author critiques traditional state‑model approaches to permeation. Those models treat a channel as a set of discrete conformational states linked by voltage‑independent rate constants. Because they ignore the explicit electric field, they cannot reproduce the experimentally observed voltage‑dependent flux ratios or the non‑Ohmic I‑V curves of many channels. This inconsistency motivates a shift toward models that embed the field directly into the transition dynamics.

The core of the manuscript is a systematic comparison of three modeling hierarchies:

1. Molecular Dynamics (MD) – All atoms are represented with classical force fields; electrostatics are computed via Ewald summation or particle‑mesh methods. MD provides the most faithful picture of atomic charge distributions, water structuring, and ion‑protein interactions. However, its temporal window (typically up to microseconds) and spatial limits (tens of nanometers) make it unsuitable for studying long‑time gating cycles or macroscopic fluxes. Periodic boundary conditions, while convenient for bulk simulations, obscure the asymmetric trans‑membrane potential that drives real channels.

2. Langevin Dynamics (LD) – By coarse‑graining the solvent into a friction term and a stochastic noise term, LD extends the accessible timescale to milliseconds while retaining a realistic description of ion motion under an external field. The Langevin equation incorporates the electric force directly, allowing the simulation of field‑driven barrier crossing and stochastic gating events. The author demonstrates how LD can reproduce the “permion” reaction path—a hypothesized transition state for ion permeation—without the prohibitive cost of full MD.

3. Poisson‑Nernst‑Planck (PNP) Continuum Theory – At the highest level of abstraction, ion concentrations and electric potential are treated as continuous fields governed by coupled differential equations. PNP excels at predicting steady‑state flux ratios, reversal potentials, and the behavior of ion pumps that couple transport to an imposed field (field‑coupled pumping). The paper shows that, when appropriate boundary conditions are imposed, PNP can capture the essential physics of electrodiffusion in mixed electrolytes and even predict gating currents arising from charge redistribution within the channel protein.

The author weaves these three levels together using the principle of Occam’s razor: a model should be as simple as possible but no simpler. Over‑parameterized atomistic simulations waste computational resources, while overly simplified continuum models miss critical charge‑redistribution effects that underlie selectivity and gating. The recommended workflow is hierarchical: start with MD to calibrate force‑field parameters and identify key charged residues, translate those insights into LD simulations to explore stochastic gating pathways, and finally employ PNP to extrapolate the results to physiological scales.

A substantial portion of the manuscript is devoted to the physics of gating. The author argues that a single protein conformation can support multiple functional states purely through electric‑field‑induced reorganization of internal charges—a “field‑switching” mechanism. In this view, the gating current measured experimentally corresponds to the movement of gating charges as the field reshapes the electrostatic landscape, not necessarily to large‑scale protein domain motions. The paper also discusses blocking phenomena, showing how external ions or drugs can alter the local field and raise the effective energy barrier, thereby reducing conductance.

Finally, the manuscript asks whether a unifying theory of ion‑channel function exists. By linking atomic‑scale electrostatics to macroscopic observables through the hierarchical modeling framework, the author suggests that such a theory is attainable, provided that simplicity (as dictated by evolutionary pressure) is kept central. The paper concludes that atomic biology offers the quantitative tools needed to test evolutionary hypotheses about channel design, to predict the effects of mutations, and ultimately to guide the rational engineering of synthetic channels with desired conductance and selectivity properties.

In sum, the work is a comprehensive call to integrate atomic‑level electrostatic calculations with stochastic dynamics and continuum electrodiffusion models, thereby creating a multi‑scale theoretical edifice capable of explaining ion‑channel permeation, gating, and pumping in a physically rigorous yet computationally tractable manner.