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.

Computational study of a calcium release-activated calcium channel

Functional Study Of Cbs Domain Interaction During Common Gating Of Clc-0 Chloride Channel

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

Homology Modeling Identifies C-Terminal Residues That Contribute to the Ca2+ Sensitivity of A Bkca Channel

Exploring the Membrane Potential of Simple Dual-Membrane Systems As Models for Gap-Junction Channels.

Pore waters regulate ion permeation in a calcium release-activated calcium channel

Molecular Dynamics Study of the KcsA Potassium Channel

Energetics of discrete selectivity bands and mutation-induced transitions in the calcium-sodium ion channels family

Conservation of cardiac L-type Ca2+ channels and their regulation in Drosophila: A novel genetically-pliable channelopathic model

Molecular Dynamics Simulation of Transmembrane Transport of Chloride Ions in Mutants of Channelrhodopsin

Simulations on Simple Models of Connexin Hemichannels Indicate That Ca2+ Blocking Is Not a Pure Electrostatic Effect

Structure and assembly of calcium homeostasis modulator proteins

Calcium Permeation across the AMPA Receptor Channel: a Combined MD and QM/MM Study

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

Mechanism of Calcium Ion-Selective Channel Opening in the ChR2_L132C Mutant: A Molecular Dynamics Simulation

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

Challenges and advances in atomistic simulations of potassium and sodium Ion channel gating and permeation

An all-atom model of the human cardiac sodium channel in a lipid bilayer

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

Cryo-EM structure of the heptameric calcium homeostasis modulator 1 channel

Revealing a hidden conducting state by manipulating the intracellular domains in K 10.1 exposes the coupling between two gating mechanisms.