Coulomb Blockade of Stochastic Permeation in Biological Ion Channels

Total Page:16

File Type:pdf, Size:1020Kb

Coulomb Blockade of Stochastic Permeation in Biological Ion Channels Outline Coulomb Blockade of Stochastic Permeation in Biological Ion Channels W.A.T.Gibby1 I.Kh. Kaufman1 D.G. Luchinsky1,2 P.V.E. McClintock1 R.S. Eisenberg3 1Department of Physics, Lancaster University, UK 2Mission Critical Technologies Inc., El Segundo, CA USA 3Molecular Biophysics, Rush University, Chicago, USA UPoN Barcelona – 13 July 2015 Gibby, Kaufman, Luchinsky, McClintock, Eisenberg Coulomb Blockade Dynamics Outline Outline 1 The unsolved problem 10 (a) Background 5 ) T 0 B Electrostatic model k ( −5 E Effect of the fixed charge −10 −15 2 Solving the problem? 0 1 20 40 Modelling – Brownian dynamics 2 −40 −20 0 Qf (e) x(A˚) Ionic Coulomb blockade 10 ) (b) T Mechanisms of permeation B 0 k ( 3 Conclusions E −10 −30 −20 −10 0 10 20 30 Prospects x(A)˚ Summary How are ions transported selectively through Ca2+ and Na+ channels? Gibby, Kaufman, Luchinsky, McClintock, Eisenberg Coulomb Blockade Dynamics The unsolved problem Background Solving the problem? Electrostatic model Conclusions Effect of the fixed charge Ion channels Cellular membrane has numerous ion channels (+ pumps and transporters). Ion channel is a natural nanotube through the membrane. Allows ion exchange between inside and outside of cell. Essential to physiology – bacteria to humans. Highly selective for particular ions. Gibby, Kaufman, Luchinsky, McClintock, Eisenberg Coulomb Blockade Dynamics The unsolved problem Background Solving the problem? Electrostatic model Conclusions Effect of the fixed charge Gating Channels spontaneously “gating” open/shut. A stochastic process influenced by e.g. – Voltage Chemicals We are interested only in open channels. Gibby, Kaufman, Luchinsky, McClintock, Eisenberg Coulomb Blockade Dynamics The unsolved problem Background Solving the problem? Electrostatic model Conclusions Effect of the fixed charge Puzzles 1. Selectivity? E.g. Calcium channel favours Ca2+ over Na+ by up to 1000:1, even though they are the same size – example of valence selectivity. Also alike charge selectivity, e.g. potassium channel strongly disfavours sodium, even though Na+ is smaller K+. 2. Fast permeation? Almost at the rate of free diffusion (open hole). 3. AMFE? Na+ goes easily through a calcium channel but is blocked by tiny traces of Ca2+. Gibby, Kaufman, Luchinsky, McClintock, Eisenberg Coulomb Blockade Dynamics The unsolved problem Background Solving the problem? Electrostatic model Conclusions Effect of the fixed charge Puzzles 4. Function of the fixed charge at the SF Ion channels have narrow “selectivity filters” with fixed negative charge... somehow associated with selectivity. What does the charge do, and how does it determine selectivity? 5. Mutations Mutations that alter the fixed charge (alone) can – (a) Destroy the channel (so it no longer conducts), or (b) Change the channel selectivity, e.g. Ca2+ to Na+ or vice versa. 6. Gating Why/how does the channel continually open and close? Gibby, Kaufman, Luchinsky, McClintock, Eisenberg Coulomb Blockade Dynamics The unsolved problem Background Solving the problem? Electrostatic model Conclusions Effect of the fixed charge Atomic structure of KcsA Potassium ion channel So-called “crystal structure” of a bacterial ion channel. Very complicated. Knowledge of the structure does not immediately explain the function – the famous “structure-function problem”. Which features are important? Selectivity filter? Generic features of structure? Need to pick out aspects important for modelling the permeation process. Gibby, Kaufman, Luchinsky, McClintock, Eisenberg Coulomb Blockade Dynamics The unsolved problem Background Solving the problem? Electrostatic model Conclusions Effect of the fixed charge Minimal model of calcium/sodium ion channel Na + Ca 2+ L Qf Ra X R Left bath, Right bath, CL>0 CR=0 Membrane protein Awater-filled, cylindrical hole, radius R = 3Å and length L = 16Å through the protein hub in the cellular membrane. Water and protein described as continuous media with dielectric constants εw = 80 (water) and εp = 2 (protein) The selectivity filter (charged residues) represented by a rigid ring of negative charge Qf = 0 − 6.5e. Gibby, Kaufman, Luchinsky, McClintock, Eisenberg Coulomb Blockade Dynamics The unsolved problem Background Solving the problem? Electrostatic model Conclusions Effect of the fixed charge Calcium and sodium channels We focus on voltage-gated calcium and sodium ion channels – which control muscle contraction, neurotransmitter secretion, transmission of action potentials) – because: They are very similar in structure, but have differing SF loci and hence different Qf . The many mutation experiments to change Qf (destroying the channel or changing its selectivity) are potentially revealing but not yet properly understood. But results and conclusions are more generally applicable. Gibby, Kaufman, Luchinsky, McClintock, Eisenberg Coulomb Blockade Dynamics The unsolved problem Background Solving the problem? Electrostatic model Conclusions Effect of the fixed charge Permeation through an ion channel Typically, the fixed charge Qf provides a strong binding site at the selectivity filter. Dimensions & electrostatics enforce single-file motion. Arriving ions captured at binding site, then fluctuational escape occurs over the potential barrier Eb. Gibby, Kaufman, Luchinsky, McClintock, Eisenberg Coulomb Blockade Dynamics The unsolved problem Modelling – Brownian dynamics Solving the problem? Ionic Coulomb blockade Conclusions Mechanisms of permeation Electrostatics The electrostatic field is derived by self-consistent numerical solution of Poisson’s equation: −∇(εε0∇u)= ezi ni X where ε0 is the dielectric permittivity of vacuum, ε is the dielectric permittivity of the medium (water or protein), u is the electric potential, e is the elementary charge, zi is the valence, and ni is the number density of ions. This equation accounts for both ion-ion interaction and self-action for all ions in their current positions. One result of the calculations is the axial potential energy profile = the Potential of the Mean Force (PMF). Gibby, Kaufman, Luchinsky, McClintock, Eisenberg Coulomb Blockade Dynamics The unsolved problem Modelling – Brownian dynamics Solving the problem? Ionic Coulomb blockade Conclusions Mechanisms of permeation Brownian dynamics The BD simulations use numerical solution of the 1-D overdamped, time-discretized, Langevin equation for the i-th ion: dx ∂u i = −D z + 2D ξ(t) dt i i ∂x i i p where xi stands for the ion’s position, Di is its diffusion coefficient, zi is the valence, u is the self-consistent potential in kBT /e units, and ξ(t) is normalized white noise. Numerical solution is implemented with the Euler forward scheme. We use an ion injection scheme that allows us to avoid simulations in the bulk. The arrival rate jarr is connected to the bulk concentration C through the Smoluchowski diffusion rate: jarr = 2πDRC. Gibby, Kaufman, Luchinsky, McClintock, Eisenberg Coulomb Blockade Dynamics The unsolved problem Modelling – Brownian dynamics Solving the problem? Ionic Coulomb blockade Conclusions Mechanisms of permeation Calcium conduction and occupancy bands Current J and 1 occupation P as a function of charge Qf at J (a.u.) 0 selectivity filter. 50 + 5 6 2 0 0 1 2 3 4 For different Ca [Ca] (mM) Q (e) f concentrations [Ca]. M0 M1 M2 1 We find – 20mM J (a.u.) 40mM 80mM Pattern of narrow 0 3 conduction bands 2 and stop bands. P 20mM 1 40mM Conduction bands 80mM 0 occur at transitions 0 1 2 3 4 5 6 7 Q (e) in channel f occupancy P. Gibby, Kaufman, Luchinsky, McClintock, Eisenberg Coulomb Blockade Dynamics The unsolved problem Modelling – Brownian dynamics Solving the problem? Ionic Coulomb blockade Conclusions Mechanisms of permeation Electrostatic exclusion principle In absence of fixed charge Qf , self-energy barrier Us prevents entry of ion to SF. But Qf compensates Us and allows cation to enter. This effectively restores the impermeable Us for 2nd ion at channel mouth. So for this Qf only one ion can occupy the SF. Implications 1. The SF’s forbidden multi-occupancy is an electrostatic exclusion principle. 2. Like the Pauli exclusion principle in quantum mechanics, it implies a Fermi-Dirac occupancy distribution. 3. For larger Qf similar arguments apply for occupancies of 2,3... Gibby, Kaufman, Luchinsky, McClintock, Eisenberg Coulomb Blockade Dynamics The unsolved problem Modelling – Brownian dynamics Solving the problem? Ionic Coulomb blockade Conclusions Mechanisms of permeation Coulomb blockade: channels vs. quantum dots + Ca2 channel Quantum dot 15 /s) 7 10 (a) M M M 0 1 2 (10 5 Ca J 20mM 0 40mM 4 80mM (b) Z Ca 2 P 2 Z Z 1 3 0 20 (c) U n T) U B M M M G 10 0 1 2 U/(k n=0 n=1 n=2 n=3 0 0 1 2 3 4 5 6 |Q /e| f Ion(s) trapped at SF ⇔ Electron(s) trapped in quantum dot Periodic conduction bands ⇔ Coulomb blockade oscillations Steps in occupation number ⇔ Coulomb staircase Classical mechanics for ion ⇔ Quantum mechanics for electron Stochastic permeation by ion ⇔ Quantum tunnelling by electron Gibby, Kaufman, Luchinsky, McClintock, Eisenberg Coulomb Blockade Dynamics The unsolved problem Modelling – Brownian dynamics Solving the problem? Ionic Coulomb blockade Conclusions Mechanisms of permeation Singular points in ionic Coulomb blockade 15 /s) 7 10 (a) M M M 0 1 2 (10 5 Ca J 20mM 0 40mM 4 80mM (b) Z Ca 2 P 2 Z Z 1 3 0 20 (c) U n T) U B M M M G 10 0 1 2 U/(k n=0 n=1 n=2 n=3 0 0 1 2 3 4 5 6 |Q /e| f Zn = zen ± δZn, Coulomb blockade Mn = ze(n + 1/2) ± δMn Resonant conduction Gibby, Kaufman, Luchinsky, McClintock, Eisenberg Coulomb Blockade Dynamics The unsolved problem Modelling – Brownian
Recommended publications
  • UC San Diego UC San Diego Electronic Theses and Dissertations
    UC San Diego UC San Diego Electronic Theses and Dissertations Title Particle dynamics in nanopore systems Permalink https://escholarship.org/uc/item/035178cn Author Krems, Matthew A. Publication Date 2011 Peer reviewed|Thesis/dissertation eScholarship.org Powered by the California Digital Library University of California UNIVERSITY OF CALIFORNIA, SAN DIEGO Particle Dynamics in Nanopore Systems A dissertation submitted in partial satisfaction of the requirements for the degree Doctor of Philosophy in Physics by Matthew A. Krems Committee in charge: Professor Massimiliano Di Ventra, Chair Professor Prabhakar Bandaru Professor Olga Dudko Professor Clifford Kubiak Professor Oleg Shpyrko 2011 Copyright Matthew A. Krems, 2011 All rights reserved. The dissertation of Matthew A. Krems is approved, and it is acceptable in quality and form for publication on microfilm and electronically: Chair University of California, San Diego 2011 iii DEDICATION To my mom iv EPIGRAPH Live as if you were to die tomorrow, learn as if you were to live forever. |Gandhi v TABLE OF CONTENTS Signature Page.................................. iii Dedication..................................... iv Epigraph.....................................v Table of Contents................................. vi List of Figures.................................. viii Acknowledgements................................x Vita and Publications.............................. xi Abstract of the Dissertation........................... xii Chapter 1 Introduction............................1 Chapter
    [Show full text]
  • Coulomb Blockade Model of Permeation and Selectivity in Biological Ion Channels
    Home Search Collections Journals About Contact us My IOPscience Coulomb blockade model of permeation and selectivity in biological ion channels This content has been downloaded from IOPscience. Please scroll down to see the full text. 2015 New J. Phys. 17 083021 (http://iopscience.iop.org/1367-2630/17/8/083021) View the table of contents for this issue, or go to the journal homepage for more Download details: IP Address: 162.247.145.43 This content was downloaded on 08/12/2015 at 06:57 Please note that terms and conditions apply. New J. Phys. 17 (2015) 083021 doi:10.1088/1367-2630/17/8/083021 PAPER Coulomb blockade model of permeation and selectivity in biological OPEN ACCESS ion channels RECEIVED 12 March 2015 I Kh Kaufman1, P V E McClintock1,3 and R S Eisenberg2 REVISED 1 Department of Physics, Lancaster University, Lancaster LA1 4YB, UK 22 June 2015 2 Department of Molecular Biophysics and Physiology, Rush Medical College, 1750 West Harrison, Chicago, IL 60612, USA ACCEPTED FOR PUBLICATION 3 Author to whom any correspondence should be addressed. 2 July 2015 PUBLISHED E-mail: [email protected], [email protected] and [email protected] 11 August 2015 Keywords: ion channel, Coulomb blockade, electrostatics, nonlinear, stochastic, fluctuational Content from this work may be used under the terms of the Creative Abstract Commons Attribution 3.0 licence. Biological ion channels are protein nanotubes embedded in, and passing through, the bilipid Any further distribution of membranes of cells. Physiologically, they are of crucial importance in that they allow ions to pass into this work must maintain fi attribution to the and out of cells, fast and ef ciently, though in a highly selective way.
    [Show full text]
  • Coulomb Blockade Oscillations in Biological Ion Channels
    Coulomb Blockade Oscillations in Biological Ion Channels I. Kh. Kaufman∗, W. Gibby∗, D. G. Luchinsky∗†, P. V. E. McClintock∗, R. S. Eisenberg‡ ∗Department of Physics, Lancaster University, Lancaster LA1 4YB, UK Email: [email protected] †Mission Critical Technologies Inc., 2041 Rosecrans Ave. Suite 225 El Segundo, CA 90245, USA ‡Department of Molecular Biophysics, Rush University, Chicago IL 60612, USA Abstract—The conduction and selectivity of calcium/sodium ion channels are described in terms of ionic Coulomb blockade, a phenomenon based on charge discreteness, an electrostatic exclusion principle, and stochastic ion motion through the channel. This novel approach provides a unified explanation of numerous observed and modelled conductance and selectivity phenomena, including the anomalous mole fraction effect and discrete conduction bands. Ionic Coulomb blockade and resonant conduction are similar to electronic Coulomb blockade and resonant tunnelling in quantum dots. The model is equally applicable to other nanopores. Biological ion channels are natural nanopores providing for the fast and highly selective permeation of physiologically important ions (e.g. Na+, K+ and Ca2+) through cellular membranes. [1]. The conduction and selectivity of e.g. voltage- gated Ca2+ [2] and Na+ channels [3] are defined by the ions’ stochastic movements and interactions inside a short, Fig. 1. (Color online) (a) Electrostatic model of a Ca2+ or Na+ channel. Ions narrow selectivity filter (SF) lined with negatively-charged 2+ 1 move in single file along the channel axis. (b) Energetics of moving Ca protein residues providing a net fixed charge Qf . Permeation ion for fixed charge Qf = −1e. The dielectric self-energy barrier Us (full through the SF sometimes involves the concerted motion of blue line) is balanced by the site attraction Ua (dashed green line) resulting more than one ion [4], [5].
    [Show full text]
  • Frontiers of Quantum and Mesoscopic Thermodynamics (Online Conference) 18 - 24 July 2021, Prague, Czech Republic
    Frontiers of Quantum and Mesoscopic Thermodynamics (online conference) 18 - 24 July 2021, Prague, Czech Republic Under the auspicies of Prof. RNDr. Eva Zažímalová, CSc. President of the Czech Academy of Sciences Supported by • Committee on Education, Science, Culture, Human Rights and Petitions of the Senate of the Parliament of the Czech Republic • Institute of Physics, the Czech Academy of Sciences • Institute for Theoretical Physics, University of Amsterdam, The Netherlands • College of Engineering and Science, University of Detroit Mercy, USA • Department of Physics, Texas A&M University, USA • Quantum Optics Lab at the BRIC, Baylor University, USA • Institut de Physique Théorique, CEA/CNRS Saclay, France Topics • Non-equilibrium quantum phenomena • Dissipation, dephasing and noise • Quantum statistical physics and thermodynamics • Foundations of quantum physics • Quantum measurement, entanglement and coherence • Many body physics, quantum field theory • Light matter interactions, quantum optics • Physics of quantum information and computing • Topological states of quantum matter, quantum phase transitions • Macroscopic quantum behavior • Atomic physics, cold atoms and molecules • Mesoscopic, nano-electromechanical and nano-optical systemss • Molecular motors, quantum heat engines • Biological systems • Cosmology, gravitation and astrophysics Scientific Committee Chair: Václav Špickaˇ (Institute of Physics, Czech Academy of Sciences, Prague) Co-Chair: Theo M. Nieuwenhuizen (University of Amsterdam) Raymond Dean Astumian (University
    [Show full text]
  • Origin and Control of Ionic Hydration Patterns in Nanopores ✉ Miraslau L
    ARTICLE https://doi.org/10.1038/s43246-021-00162-x OPEN Origin and control of ionic hydration patterns in nanopores ✉ Miraslau L. Barabash 1 , William A. T. Gibby 1, Carlo Guardiani 1,3, Alex Smolyanitsky 2, Dmitry G. Luchinsky 1,4 & Peter V. E. McClintock 1 In order to permeate a nanopore, an ion must overcome a dehydration energy barrier caused by the redistribution of surrounding water molecules. The redistribution is inhomogeneous, anisotropic and strongly position-dependent, resulting in complex patterns that are routinely observed in molecular dynamics simulations. Here, we study the physical origin of these patterns and of how they can be predicted and controlled. We introduce an analytic model 1234567890():,; able to predict the patterns in a graphene nanopore in terms of experimentally accessible radial distribution functions, giving results that agree well with molecular dynamics simula- tions. The patterns are attributable to a complex interplay of ionic hydration shells with water layers adjacent to the graphene membrane and with the hydration cloud of the nanopore rim atoms, and we discuss ways of controlling them. Our findings pave the way to designing required transport properties into nanoionic devices by optimising the structure of the hydration patterns. 1 Department of Physics, Lancaster University, Lancaster, UK. 2 Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, CO, USA. 3Present address: Department of Mechanical and Aerospace Engineering, Sapienza University, Rome, Italy. 4Present address: KBR Inc., ✉ Ames Research Center, Moffett Field, CA, USA. email: [email protected] COMMUNICATIONS MATERIALS | (2021) 2:65 | https://doi.org/10.1038/s43246-021-00162-x | www.nature.com/commsmat 1 ARTICLE COMMUNICATIONS MATERIALS | https://doi.org/10.1038/s43246-021-00162-x onic permeation through narrow water-filled channels and Results pores is of profound importance in both biophysics and Model.
    [Show full text]
  • Selective Molecular Transport in Nanopore Systems
    SELECTIVE MOLECULAR TRANSPORT IN NANOPORE SYSTEMS BY KARL STEVEN DECKER DISSERTATION Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Physics in the Graduate College of the University of Illinois at Urbana-Champaign, 2017 Urbana, Illinois Doctoral Committee: Professor Paul Selvin, Chair Professor Aleksei Aksimentiev, Director of Research Professor Emad Tajkhorshid Professor John Stack ABSTRACT A detailed characterization of the physics of novel nanopore systems has the potential to revolutionize water filtration, nanofluidics, and biomolecule detection technologies. I give my characterizations of five nanopore systems as my dissertation. First, I present my study of nanopores in polyethylene terephthalate (PET), revealing the mechanism for variance in current rectification based on cation species. Second, I demonstrate the mechanism of selective probe capture in bacterial toxin protein α-hemolysin (aHL) using dielectrophoresis. Third, I introduce the first simulation of molecular artificial water channel pillar[5]arene (PAP) and uncover the mechanics of its water transport and self-aggregation properties. Fourth, I characterize the water permeability and ion rejection of truncated human membrane protein aquaporin-1 (AQP) in simulation. Finally, I present MD simulation of truncated AQP as a voltage-gated ionic diode and as the functional element of a double-membrane ionic pump. ii TABLE OF CONTENTS CHAPTER 1: INTRODUCTION ....................................................................................................1
    [Show full text]
  • Coulomb Blockade Oscillations in Biological Ion Channels
    Paper 1570067345 at 23rd International Conference on Noise and Fluctuations (ICNF 2015) Coulomb Blockade Oscillations in Biological Ion Channels I. Kh. Kaufman*, W. Gibby*, D. G. Luchinsky*t, P. V. E. McClintock*, R. S. Eisenberg+ *Department of Physics, Lancaster University, Lancaster LAI 4YB, UK tMission Critical Technologies Inc., 2041 Rosecrans Ave. Suite 225 El Segundo, CA 90245, USA +Department of Molecular Biophysics, Rush University, Chicago IL 60612, USA Email: [email protected] Abstract-The conduction and selectivity of calcium/sodium 1.2 ion channels are described in terms of ionic Coulomb blockade, ii =80 a phenomenon based on charge discreteness, an electrostatic 0.6 0 • exclusion principle, and stochastic ion motion through the 00 channel. This novel approach provides a unified explanation of E 0 numerous observed and modelled conductance and selectivity ~ 0 phenomena, including the anomalous mole fraction effect and -0.6 discrete conduction bands. Ionic Coulomb blockade and resonant • conduction are similar to electronic Coulomb blockade and -1.2'--------'----'---....__ _...__ _._ __ .1__ _ __L__ _J resonant tunnelling in quantum dots. The model is equally applicable to other nanopores. Biological ion channels are natural nanopores providing for the fast and highly selective permeation of physiologically J important ions (e.g. Na+, K+ and Ca2 +) through cellular ~ :}f :: ~-~--~--~-~--~--~-~--~g __ ttt:-t membranes. [1]. The conduction and selectivity of e.g. voltage­ -1 .6 -0.8 0 0.8 1.6 x(nm) gated Ca2 + [2] and Na+ channels [3] are defined by the ions' stochastic movements and interactions inside a short, Fig.
    [Show full text]
  • Ionic Coulomb Blockade in Nanopores
    Ionic Coulomb Blockade in Nanopores Matt Krems and Massimiliano Di Ventra Department of Physics, University of California, San Diego, La Jolla, CA 92093 (Dated: December 6, 2018) Understanding the dynamics of ions in nanopores is essential for applications ranging from single- molecule detection to DNA sequencing. We show both analytically and by means of molecular dynamics simulations that under specific conditions ion-ion interactions in nanopores lead to the phenomenon of ionic Coulomb blockade, namely the build-up of ions inside a nanopore with specific capacitance impeding the flow of additional ions due to Coulomb repulsion. This is the counterpart of electronic Coulomb blockade observed in mesoscopic systems. We discuss the analogies and differences with the electronic case as well as experimental situations in which this phenomenon could be detected. Recently, there has been a surge of interest in takes the same ion to propagate unhindered by the pore. nanopores due to their potential in several technologi- In other words, the “contact” ionic resistance of the pore cal applications, the most notable being DNA sequencing with at least one reservoir has to be much larger than and detection [1–7]. In addition, nanopores offer unprece- the resistance offered to ion flow without the pore. From dented opportunities to study several fundamental physi- previous work [8, 9] we know that such a situation may cal processes related to ionic transport in confined geome- be realized when the pore opening is in the “quantum” tries. For example, due to the tightly bound hydration regime, namely when the pore radius is small enough that layers surrounding an ion, it may be possible to observe an ion has to shed part of its hydration layers to cross steps in the ionic conductance due to the shedding of the pore aperture.
    [Show full text]