On a novel rate theory for transport in narrow ion channels and its application to the study of flux optimization via geometric effects

Calculating tracer currents through narrow ion channels: beyond the independent particle model

Discrete state models of single-file ion permeation through a narrow ion channel pore are employed to analyze the ratio of forward to backward tracer current. Conditions under which the well-known Ussing formula for this ratio hold are explored in systems where ions do not move independently through the channel. Building detailed balance into the rate constants for the model in such a way that under equilibrium conditions (equal rate of forward versus backward permeation events) the Nernst equation is satisfied, it is found that in a model where only one ion can occupy the channel at a time, the Ussing formula is always obeyed for any number of binding sites, reservoir concentrations of the ions and electric potential difference across the membrane which the ion channel spans, independent of the internal details of the permeation pathway. However, numerical analysis demonstrates that when multiple ions can occupy the channel at once, the nonequilibrium forward/backward tracer flux ratio deviates from the prediction of the Ussing model. Assuming an appropriate effective potential experienced by ions in the channel, we provide explicit formulae for the rate constants in these models.

Tight-binding derivation of a discrete-continuous description of mechanochemical coupling in a molecular motor

We present a tight-binding derivation of a discrete-continuous description of mechanochemical coupling in a molecular motor. Our derivation is based on the continuous diffusion equation for overdamped Brownian motion on a time-independent tilted periodic potential in two dimensions. The free-energy potential is nonseparable to allow coupling between the chemical and mechanical degrees of freedom. We formally discretize the chemical coordinate by expanding in Wannier states that are localized along the chemical coordinate and parametrized along the mechanical coordinate. A discrete-continuous equation is derived that is valid for anisotropic systems with weak mechanochemical coupling and deep potential wells along the chemical coordinate. The discrete-continuous description is consistent with established theoretical models of molecular motors with discrete chemical states but is constrained by the underlying continuous two-dimensional potential. In particular, we derive analytic expressions for the effective potential along the mechanical coordinate and for the rate of thermal hopping between chemical states. We determine the thermodynamic efficiency of energy conversion and find that, for a molecular motor with one chemical state per cycle, the derived discrete-continuous equation can accurately describe mechanochemical coupling but cannot describe energy conversion.

Local discretization method for overdamped Brownian motion on a potential with multiple deep wells

We present a general method for transforming the continuous diffusion equation describing overdamped Brownian motion on a time-independent potential with multiple deep wells to a discrete master equation. The method is based on an expansion in localized basis states of local metastable potentials that match the full potential in the region of each potential well. Unlike previous basis methods for discretizing Brownian motion on a potential, this approach is valid for periodic potentials with varying multiple deep wells per period and can also be applied to nonperiodic systems. We apply the method to a range of potentials and find that potential wells that are deep compared to five times the thermal energy can be associated with a discrete localized state while shallower wells are better incorporated into the local metastable potentials of neighboring deep potential wells.

Model reduction for slow-fast stochastic systems with metastable behaviour

The quasi-steady-state approximation (or stochastic averaging principle) is a useful tool in the study of multiscale stochastic systems, giving a practical method by which to reduce the number of degrees of freedom in a model. The method is extended here to slow-fast systems in which the fast variables exhibit metastable behaviour. The key parameter that determines the form of the reduced model is the ratio of the timescale for the switching of the fast variables between metastable states to the timescale for the evolution of the slow variables. The method is illustrated with two examples: one from biochemistry (a fast-species-mediated chemical switch coupled to a slower varying species), and one from ecology (a predator-prey system). Numerical simulations of each model reduction are compared with those of the full system.

Energetics of Multi-Ion Conduction Pathways in Potassium Ion Channels

Potassium ion channels form pores in cell membranes, allowing potassium ions through while preventing the passage of sodium ions. Despite numerous high-resolution structures, it is not yet possible to relate their structure to their single molecule function other than at a qualitative level. Over the past decade, there has been a concerted effort using molecular dynamics to capture the thermodynamics and kinetics of conduction by calculating potentials of mean force (PMF). These can be used, in conjunction with the electro-diffusion theory, to predict the conductance of a specific ion channel. Here, we calculate seven independent PMFs, thereby studying the differences between two potassium ion channels, the effect of the CHARMM CMAP forcefield correction, and the sensitivity and reproducibility of the method. Thermodynamically stable ion-water configurations of the selectivity filter can be identified from all the free energy landscapes, but the heights of the kinetic barriers for potassium ions to move through the selectivity filter are, in nearly all cases, too high to predict conductances in line with experiment. This implies it is not currently feasible to predict the conductance of potassium ion channels, but other simpler channels may be more tractable.

A permeation theory for single-file ion channels: one- and two-step models

How many steps are required to model permeation through ion channels? This question is investigated by comparing one- and two-step models of permeation with experiment and MD simulation for the first time. In recent MD simulations, the observed permeation mechanism was identified as resembling a Hodgkin and Keynes knock-on mechanism with one voltage-dependent rate-determining step [Jensen et al., PNAS 107, 5833 (2010)]. These previously published simulation data are fitted to a one-step knock-on model that successfully explains the highly non-Ohmic current-voltage curve observed in the simulation. However, these predictions (and the simulations upon which they are based) are not representative of real channel behavior, which is typically Ohmic at low voltages. A two-step association/dissociation (A/D) model is then compared with experiment for the first time. This two-parameter model is shown to be remarkably consistent with previously published permeation experiments through the MaxiK potassium channel over a wide range of concentrations and positive voltages. The A/D model also provides a first-order explanation of permeation through the Shaker potassium channel, but it does not explain the asymmetry observed experimentally. To address this, a new asymmetric variant of the A/D model is developed using the present theoretical framework. It includes a third parameter that represents the value of the "permeation coordinate" (fractional electric potential energy) corresponding to the triply occupied state n of the channel. This asymmetric A/D model is fitted to published permeation data through the Shaker potassium channel at physiological concentrations, and it successfully predicts qualitative changes in the negative current-voltage data (including a transition to super-Ohmic behavior) based solely on a fit to positive-voltage data (that appear linear). The A/D model appears to be qualitatively consistent with a large group of published MD simulations, but no quantitative comparison has yet been made. The A/D model makes a network of predictions for how the elementary steps and the channel occupancy vary with both concentration and voltage. In addition, the proposed theoretical framework suggests a new way of plotting the energetics of the simulated system using a one-dimensional permeation coordinate that uses electric potential energy as a metric for the net fractional progress through the permeation mechanism. This approach has the potential to provide a quantitative connection between atomistic simulations and permeation experiments for the first time.