applied sciences
Article Quantum Mechanical Coherence of K+ Ion Wave Packets Increases Conduction in the KcsA Ion Channel
Johann Summhammer 1,*, Georg Sulyok 1 and Gustav Bernroider 2
1 Atominstitut, Technische Universität Wien, Stadionallee 1, 1020 Vienna, Austria; [email protected] 2 Department of Biosciences, Universität Salzburg (retired), 5020 Salzburg, Austria (retired); [email protected] * Correspondence: [email protected]
Received: 3 June 2020; Accepted: 18 June 2020; Published: 21 June 2020
Abstract: 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.
Keywords: potassium ion channel; KcsA channel; selectivity filter; quantum wave packet; coherence time; biological quantum coherence
1. Introduction Electrical signals generated in nerve cells are generated by membrane proteins that catalyze a selective electro-diffusion of ions across plasma membranes. The proteins are embedded in the lipid bilayer of membranes and provide currents at the pico-ampere level, combining fast conduction with a high selectivity for particular ions. In voltage-gated potassium channels, this selective control of ion conduction is organized by a specific arrangement of amino-acids within a narrow ‘selectivity filter’ domain of the protein [1,2]. The most frequently studied filter model is based on the bacterial-derived KcsA channel motif (KcsA, potassium crystallographically sited activation channel) where the filter
Appl. Sci. 2020, 10, 4250; doi:10.3390/app10124250 www.mdpi.com/journal/applsci Appl. Sci. 2019, 9, x FOR PEER REVIEW 3 of 17 for their carbonyl groups oriented towards the pore axis to change their vibrational motion. Each Appl. Sci. 2020, 10, 4250 2 of 17 oxygen atom of a carbonyl (oxygens are shown as red spheres) has a constant distance to its carbon, but it can vibrate around its mean position in both angular directions. The carbon atoms have a positiveconstriction partial is formed charge, by and four the subunits, oxygen each atoms composed have a bynegative the same partial amino charge acid sequence of equal TVGYG magnitude. [3–5]. WhereasIon translocation the charge is largely of a coordinatedcarbon atom by is amino centered acid carbonylat this location, groups pointingthe charge towards on the the oxygen axis of the is dispersedconstriction, representing with oxygen oxygen lone pair lone electrons pair electrons. creating a seriesNon-shared of four potentialoxygen lone valleys. pair By charges convention, are consideredand as shown to play in Figure an important1, these transientrole in coordinating cages are labeled and guiding as sites K S1+ ions to S4 and as seenwater from dipoles outside during to diffusioninside of thethrough membrane the narrow with coordinating SF constriction interactions atoms. We mimicking initially the located hydration these shell charges of ions along in bulk an extensionwater [6]. of The the relevant C=O binding caging axis forces and are implemen Coulombicted attractionscarbonyl angle-bending between ion and modes electrons, in dependence together ofwith ion repulsive location. e ffects, if lone pair electrons and the inner shell of the K+ ion start to overlap [2,6]. Further,The themobile diff usionspecies of in K +theions SF seems are K+ to ions occur (green in a hoppingspheres) and manner water between molecules potential (pink minima,oxygen atomsand this and in turnattached is driven grey byhydrogen the transmembrane atoms in Figure field 1). resulting Ionic motion from the can charge occur concentration bidirectionally gradient along theacross SF thecentral cell membrane.axis. Water Additionalmolecules tocarry K+ ions,partia alsol charges water that molecules compensate diffuse in through sum but the carry channel, an electricexperiencing dipole similar moment carbonyl-derived with two opposing caging axis forces. orientations and arbitrary rotational orientation.
+ Figure 1.1. AA ribbonribbon model model K Kchannel+ channel embedded embedded inthe in lipidthe lipid bilayer bilayer (left), ( withleft), two with of thetwo four of the subunits four subunitsand the extracellularand the extracellular side on top side (middle on top). The(middle four). bindingThe four sites binding (S1–S4) sites with (S1–S4) four dehydrated with four dehydratedions in locations ions S1in locations and S3, and S1 oneand partiallyS3, and one dehydrated partially iondehydrated in the neighboring ion in the cavityneighboring location cavity Scav, are demonstrated (right). Ions are kept in green color, water is shown in pink, carbonyl oxygens in location Scav, are demonstrated (right) . Ions are kept in green color, water is shown in pink, carbonyl oxygensred, and in backbone red, and carbon backbone atoms carbon in brown atoms color. in brow Lonen paircolor. electrons Lone pair attached electrons to oxygens attached are to indicatedoxygens areby blueindicated dots. by The blue close-up dots. The view close-up on the view right on side the sketches right side the sketches relative the sizes relative of particles, sizes ofwhich particles, are + whichassumed are to assumed be as follows: to be as ion follows: radius ofion K radius=138 pm;of K covalent+ =138 pm; radius covalent of carbon radius (when of carbon in double (when bond in doublewith O) bond= 64.1 with pm; O)= covalent 64.1 pm; radius covalent of oxygen radius (when of oxygen in double (when bond in withdouble C) bond= 58.9 with pm; C) center–center = 58.9 pm; distance of C and O-atoms in carbonyls = 123 pm; H O: covalent radius of O-atoms 73 pm; bond radius center–center distance of C and O-atoms in carbonyls2 = 123 pm; H2O: covalent radius of O-atoms 73 of H-atoms = 38 pm; center–center distance between O and H= 95.84 pm. pm; bond radius of H-atoms = 38 pm; center–center distance between O and H= 95.84 pm. The small interaction distances of temporary caging sites together with the underlying conduction 2.2. Classical MD Simulations strategy, as hopping from one site to the next, stimulate the question of whether quantum effects in these transitionsFrom a could geometrical play a functional perspective, role. the In variation a recent study of ion we coordination have provided is strongly a first assessment confined ofto thisthe axialquestion extension by simulating of the afilter single domain site situation (z-axis), and and simulating it can be a Ksafely+ ion asassumed a quantum that mechanical lateral motion particle. is negligible.We found that,This dueaspect to thealso slow reduces dispersion computational of the initially efforts well-localizedsignificantly. Furthermore, wave packet, to the achieve enclosing the presentcarbonyl intention oxygen atomsto compare tend to a beclassical pushed with some a quan distancetum aside, mechanical and this interpretation leads to a higher of ion probability motion, the for classicalthe particle environment to escape from had this to sitebe designed within a givenin a wa timey [that7]. Insubsequently the present work, allowed we extendthe extension this study to toa quantuminclude the algorithmic complete TVGYG implementation. sequence and This permit view an necessitates arbitrary number a special of K +implementationions and water molecules of force fieldsto permeate. and in Again,our case we a findcomputational a considerable design enhancement based on ofVerlet’s a K+ ion’s algorithm probability [14]. toWe propagate calculated to the the forcenext sitebetween due to any its quantumtwo particles mechanical A and B, nature consisting as a wave of a Coulomb packet. Further, attractive the and present of a studyrepulsive confirms part, asour follows: previous finding about an optimal decoherence time of the particles’ wave-function regarding a directional ‘guidance’ through subsequent potential barriers provided by the surrounding electrons. ( , ) = − ̂ if | − | ≥ (1) Previous studies were generally ( based ) on classical MD simulations. Whereas early studies within + this frameand considered an alternating sequence of K and water molecules (e.g., KwKw) passing through the constriction [4], subsequent MD simulations also revealed other conduction mechanisms with ( , ) = − ̂ if | − | <