Computational Electrophysiology from a Single Molecular Dynamics Simulation and the Electrodiffusion Model

The availability of high-resolution structures of ion channels opens the doors to reliable computations of electrophysiological properties, such as the dependence of ionic currents and selectivities on applied voltage. We develop two theoretical approaches for calculating these properties from molecular dynamics simulations at a single voltage, or even in the absence of voltage, combined with the electrodiffusion model in which ion motion in the channel is represented as one-dimensional diffusion in the potential of mean force exerted by other components of the system and the applied electric field. No knowledge of diffusivity or ion densities at other voltages is needed. Instead, in one approach, one-sided ion fluxes and density profiles are used to determine the free energy profile. In the other approach, committor probabilities for ions transported at the selected voltage are used for this purpose. Both approaches have been validated in an example of a simple ion channel formed by trichotoxin. The potentials of mean force calculated by way of the proposed approaches and obtained from traditional methods are in excellent agreement. Furthermore, the current-voltage dependence agrees very well with results obtained by way of computationally more demanding methods. We also have readily calculated the reversal potential, a computationally challenging electrophysiological property. The key assumptions of the electrodiffusion model, such as the independence of crossing events or the insensitivity of the potential of mean force to applied voltage, have been found to be satisfied. We also show that the voltage changes linearly in the hydrophobic core of the membrane and is constant elsewhere.