Effective Polarization in Potassium Channel Simulations: Ion Conductance, Occupancy, Voltage Response, and Selectivity

Potassium (K ⁺ ) channels are widely distributed in many types of organisms. They combine high efficiency (∼100 pS) and K ⁺ /Na ⁺ selectivity by a conserved selectivity filter (SF). Molecular Dynamics (MD) simulations can provide detailed, atomistic mechanisms of this sophisticated ion permeation. However, currently there are clear inconsistencies between computational predictions and experimental results. Firstly, the ion occupancy of the SF in simulations is lower than expected (∼2.5 in MD compared to ∼4 in X-ray crystallography). Secondly, in many reported MD simulations of K ⁺ channels, K ⁺ conductance is typically an order of magnitude lower than experimental values. This discrepancy is in part because the force fields used in MD simulations of potassium channels do not account for polarization. One of the proposed solutions is the Electronic Continuum Correction (ECC), a force field modification that scales down formal charges, to introduce the polarization in a mean-field way. When the ECC is used in conjunction with the Charmm36m force field, the simulated K ⁺ conductance increases 13-fold. Following the analysis of ion occupancy states using Hamiltonian Replica Exchange (HRE) simulations, we propose a new parameter set for Amber14sb, that also leads to a similar increase in conductance. These two force fields are then used to compute, for the first time, the full current-voltage (I-V) curves from MD simulations, approaching quantitative agreement with experiments at all voltages. In general, the ECC-enabled simulations are in excellent agreement with experiment, in terms of ion occupancy, conductance, current-voltage response, and K ⁺ /Na ⁺ selectivity.