Domain and interdomain energetics underlying gating in Shaker-type Kv channels

Closed Formula for Transport across Constrictions

In the last decade, the Fick-Jacobs approximation has been exploited to capture transport across constrictions. Here, we review the derivation of the Fick-Jacobs equation with particular emphasis on its linear response regime. We show that, for fore-aft symmetric channels, the flux of noninteracting systems is fully captured by its linear response regime. For this case, we derive a very simple formula that captures the correct trends and can be exploited as a simple tool to design experiments or simulations. Lastly, we show that higher-order corrections in the flux may appear for nonsymmetric channels.

Computational methods and theory for ion channel research

Ion channels are fundamental biological devices that act as gates in order to ensure selective ion transport across cellular membranes; their operation constitutes the molecular mechanism through which basic biological functions, such as nerve signal transmission and muscle contraction, are carried out. Here, we review recent results in the field of computational research on ion channels, covering theoretical advances, state-of-the-art simulation approaches, and frontline modeling techniques. We also report on few selected applications of continuum and atomistic methods to characterize the mechanisms of permeation, selectivity, and gating in biological and model channels.

Transport of neutral and charged nanorods across varying-section channels

We study the dynamics of neutral and charged rods embedded in varying-section channels. By means of systematic approximations, we derive the dependence of the local diffusion coefficient on both the geometry and charge of the rods. This microscopic insight allows us to provide predictions for the permeability of varying-section channels to rods with diverse lengths, aspect ratios and charge. Our analysis shows that the dynamics of charged rods is sensitive to the geometry of the channel and that their transport can be controlled by tuning both the shape of the confining walls and the charge of the rod. Interestingly, we find that the channel permeability does not depend monotonically on the charge of the rod. This opens the possibility of a novel mechanism to separate charged rods.

Reversible heat production during electric double layer buildup depends sensitively on the electrolyte and its reservoir

Several modern technologies for energy storage and conversion are based on the screening of electric charge on the surface of porous electrodes by ions in an adjacent electrolyte. This so-called electric double layer (EDL) exhibits an intricate interplay with the electrolyte's temperature that was the focus of several recent studies. In one of them, Janssen et al. [Phys. Rev. Lett. 119, 166002 (2017)] experimentally determined the ratio Q_rev/W_el of reversible heat flowing into a supercapacitor during an isothermal charging process and the electric work applied therein. To rationalize that data, here, we determine Q_rev/W_el within different models of the EDL using theoretical approaches such as density functional theory (DFT) as well as molecular dynamics simulations. Applying mainly the restricted primitive model, we find quantitative support for a speculation of Janssen et al. that steric ion interactions are key to the ratio Q_rev/W_el. Here, we identified the entropic contribution of certain DFT functionals, which grants direct access to the reversible heat. We further demonstrate how Q_rev/W_el changes when calculated in different thermodynamic ensembles and processes. We show that the experiments of Janssen et al. are explained best by a charging process at fixed bulk density or in a "semi-canonical" system. Finally, we find that Q_rev/W_el significantly depends on parameters such as pore and ion size, salt concentration, and valencies of the cations and anions of the electrolyte. Our findings can guide further heat production measurements and can be applied in studies on, for instance, nervous conduction, where reversible heat is a key element.

Multiscale modeling shows that dielectric differences make NaV channels faster than KV channels

The generation of action potentials in excitable cells requires different activation kinetics of voltage-gated Na (Na_V) and K (K_V) channels. Na_V channels activate much faster and allow the initial Na⁺ influx that generates the depolarizing phase and propagates the signal. Recent experimental results suggest that the molecular basis for this kinetic difference is an amino acid side chain located in the gating pore of the voltage sensor domain, which is a highly conserved isoleucine in K_V channels but an equally highly conserved threonine in Na_V channels. Mutagenesis suggests that the hydrophobicity of this side chain in Shaker K_V channels regulates the energetic barrier that gating charges cross as they move through the gating pore and control the rate of channel opening. We use a multiscale modeling approach to test this hypothesis. We use high-resolution molecular dynamics to study the effect of the mutation on polarization charge within the gating pore. We then incorporate these results in a lower-resolution model of voltage gating to predict the effect of the mutation on the movement of gating charges. The predictions of our hierarchical model are fully consistent with the tested hypothesis, thus suggesting that the faster activation kinetics of Na_V channels comes from a stronger dielectric polarization by threonine (Na_V channel) produced as the first gating charge enters the gating pore compared with isoleucine (K_V channel).

Bubble gating in biological ion channels: A density functional theory study

We study the behavior of a waterlike liquid inside the gate of a biological ion channel following the basic geometry of the well studied potassium channel KcsA. We calculate the three-dimensional density distribution ρ(r) of the liquid within the framework of classical density functional theory and observe the formation of a low density region (bubble) when the gate is narrow. This observation corresponds to a finite-size form of capillary evaporation and supports the so-called bubble-gate theory. From the density profile we also compute the energy landscape of the gate and the energy required to change the gate from a closed (narrow) to an open (wide) state and vice versa.