We present a comparison of a classical and a quantum mechanical calculation of the motion of K+ ions in the highly conserved KcsA selectivity filter motive of voltage gated ion channels. We first show that the de Broglie wavelength of thermal ions is not much smaller than the periodic structure of Coulomb potentials in the nano-pore model of the selectivity filter. This implies that an ion may no longer be viewed to be at one exact position at a given time but can better be described by a quantum mechanical wave function. Based on first principle methods, we demonstrate solutions of a non-linear Schrödinger model that provide insight into the role of short-lived (~1 ps) coherent ion transition states and attribute an important role to subsequent decoherence and the associated quantum to classical transition for permeating ions. It is found that short coherences are not just beneficial but also necessary to explain the fast-directed permeation of ions through the potential barriers of the filter. Certain aspects of quantum dynamics and non-local effects appear to be indispensable to resolve the discrepancy between potential barrier height, as reported from classical thermodynamics, and experimentally observed transition rates of ions through channel proteins.
We simulate the transmission of K+ ions through the KcsA potassium ion channel filter region at physiological temperatures, employing classical molecular dynamics (MD) at the atomic scale together with a quantum mechanical version of MD simulation (QMD), treating single ions as quantum wave packets. We provide a direct comparison between both concepts, embedding the simulations into identical force fields and thermal fluctuations. The quantum simulations permit the estimation of coherence times and wave packet dispersions of a K+ ion over a range of 0.5 nm (a range that covers almost 50% of the filter domains longitudinal extension). We find that this observed extension of particle delocalization changes the mean orientation of the coordinating carbonyl oxygen atoms significantly, transiently suppressing their ‘caging action’ responsible for selective ion coordination. Compared to classical MD simulations, this particular quantum effect allows the K+ ions to ‘escape’ more easily from temporary binding sites provided by the surrounding filter atoms. To further elucidate the role of this observation for ion conduction rates, we compare the temporal pattern of single conduction events between classical MD and quantum QMD simulations at a femto-sec time scale. A finding from both approaches is that ion permeation follows a very irregular time pattern, involving flushes of permeation interrupted by non-conductive time intervals. However, as compared with classical behavior, the QMD simulation shortens non-conductive time by more than a half. As a consequence, and given the same force-fields, the QMD-simulated ion current appears to be considerably stronger as compared with the classical current. To bring this result in line with experimentally observed ion currents and the predictions based on Nernst–Planck theories, the conclusion is that a transient short time quantum behavior of permeating ions can successfully compromise high conduction rates with ion selectivity in the filter of channel proteins.
We provide a high resolution, all-atom, femto-second molecular dynamics (MD) simulation of the passage of K+ ions and H2O molecules through the selectivity filter of the KcsA potassium ion channel, based on first principle physical methods. Our results show that a change in the length of the selectivity filter of as little as 3%, regardless of whether the filter is made longer or shorter, will reduce the K+ ion current by around 50%. In addition, further squeezing or stretching by about 9% can effectively stop the current. Our results demonstrate optimized conformational dynamics that associate an increased mobility of parts in the filter linings with a standard configuration, leading to maximized conduction rates that are highly sensitive to geometrical distortions. We discuss this latter aspect in relation to lateral membrane effects on the filter region of ion channels and the ‘force from lipids’ hypothesis.
We investigate the possible influence of poly unsaturated fatty acids (PUFAs) in the lipid bilayer of cell membranes on the conduction properties of the selectivity filter of voltage gated ion channels. Any change of this conductivity can have major consequences on electrical signalling in brain and other cells. At the microscopic level a change in the concentration of PUFAs can cause a change in the pressure distribution within the lipid bilayer, and this in turn can lead to a change in the length of the selectivity filter. In order to estimate the consequences this might entail for the ion channel's conductivity we performed high resolution molecular dynamics (MD) simulations of the passage of K+ ions and H2O molecules through the selectivity filter of the KcsA potassium ion channel. Our results show that a change in length of the selectivity filter of as little as 4%, independent of whether the filter is made longer or shorter, will reduce the K+ ion current by around 50%. And further squeezing or stretching by about 10% can effectively stop the current. For instructive purposes and to facilitate further interpretation of the results, we give details of the force models employed in the MD simulations. We finally discuss the present results in the context of possible effects of membrane lipid compositions on Kv ion permeation.