Abstract:Potassium (K+) channels are widely distributed in many types of organisms. Potassium channels combine high efficiency (K+ ions permeation rates ~100 pS) and exquisite K+/Na+ selectivity by a conserved selectivity filter (SF) with four adjacent potassium-binding sites. Molecular Dynamics (MD) simulations can, in principle, provide a detailed, atomistic mechanism of this sophisticated ion permeation. However, currently there are clear inconsistencies between computational and experimental predictions. 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 from single-channel electrophysiology. This discrepancy is in part because the force fields used in MD simulations of potassium channels do not account for explicit electronic polarization. One of the proposed solutions is the Electronic Continuum Correction (ECC), a force field modification that scales down formal atomic charges, to introduce the polarization in a mean-field way. In this work, we apply the ECC to MD simulations of K+ channels under applied voltage. When the ECC is used in conjunction with the Charmm36m force field, the simulated K+ conductance increases 13-fold, whereas applying the ECC to Amber14 seemingly has no major effect. Following the analysis of ion occupancy states using Hamiltonian Replica Exchange (HRE) simulations, we propose a new parameter set for Amber14sb, that 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 tested, for three K+ channels. In all simulated systems, the direct knock-on permeation mechanism with near-full K+ ion occupancy of the SF is consistently observed. In general, our ECC-enabled simulations are in excellent agreement with experimental data, in terms of ion occupancy, conductance, current-voltage response, and K+/Na+ selectivity. These simulations thus form a robust basis to study structure-function relations in K+ channels on the atomistic level with unprecedented accuracy, as exemplified by deciphering the rectification mechanism of the MthK channel.

The Ion Permeation and Selectivity Mechanisms of Ryanodine Receptors-Molecular Dynamics Study with a New Multi-Site Ca2+ Model

Atomistic Details of Charge/Space Competition in the Ca2+ Selectivity of Ryanodine Receptors.

The Ca2+ permeation mechanism of the ryanodine receptor revealed by a multi-site ion model

Simulation of the Effect of Rogue Ryanodine Receptors on a Calcium Wave in Ventricular Myocytes with Heart Failure.

Isoform-Dependent Formation of Heteromeric Ca2+ Release Channels (ryanodine Receptors)

Deciphering Ca2+ Permeation and Valence Selectivity in Cav1: Molecular Dynamics Simulations Reveal the Three-Ion Knock-on Mechanism

Calcium Binding and Permeation in TRPV Channels: Insights from Molecular Dynamics Simulations

Atomistic modeling and molecular dynamics analysis of human CaV1.2 channel using external electric field and ion pulling simulations

Chemistry of cation hydration and conduction in a skeletal muscle ryanodine receptor

K(+) and Na(+) conduction in selective and nonselective ion channels via molecular dynamics simulations

Structural insights into Ca 2+ -activated long-range allosteric channel gating of RyR1

Molecular Dynamics Study of Ion Transport Through an Open Model of Voltage-Gated Sodium Channel

Charge Scaling in Potassium Channel Simulations: Conductance, Ion Occupancy, Voltage Response, and Selectivity

Structure-based mechanism of RyR channel operation by calcium and magnesium ions

The Role of Hydration and Amino Acid Interactions on the Ion Permeation Mechanism in the hNaV 1.5 Channel

Local control models of cardiac excitation-contraction coupling. A possible role for allosteric interactions between ryanodine receptors.

The Molecular Mechanism of Ion Selectivity in Nanopores

Metadynamics simulations reveal mechanisms of Na+ and Ca2+ transport in two open states of the channelrhodopsin chimera, C1C2

Ion permeation mechanism of the potassium channel

Decisive structural elements in water and ion permeation through mechanosensitive channels of large conductance: insights from molecular dynamics simulation