Paper 1570067345 at 23rd International Conference on Noise and Fluctuations (ICNF 2015) 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 . 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. 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 Q f = - l e. 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]. Mutation studies [6], [7], [8], [3] in a barrier-less energy profile Ub (red solid line). See text for details. and simulations [9], [ 10], [ 11] show that Q f is a determinant of the channel's conductivity and valence selectivity. It has general physical picture of the phenomena remained unclear. recently been shown that nanopores can exhibit ionic Coulomb Here, we reinterpret and generalize the electrostatic anal­ blockade (CB) [12], a quantized electrostatic phenomenon ysis of the multi-ion energetics of conduction bands [18] equivalent to electronic CB in mesoscopic systems [13], [14], by introducing a novel ionic CB model of conduction and [15]. selectivity in biological ion channels thereby bringing them An electrostatic theory of ionic transport in water-filled into the context of mesoscopic phenomena. We show that periodically-charged nanopores revealed the phenomenon of the experimentally-observed valence selectivity phenomena in ion-exchange through low-barrier phase transitions as the ion Ca2 + /Na+ channels, including AMFE and mutation-induced concentration and fixed charge Qf were varied [16]. Compa­ transformations of selectivity, and the simulated conduction rable transitions in Brownian dynamics (BD) simulations of bands vs. Qf, can be well-described in terms of CB conduc­ Ca2 + channels result in discrete conduction and selectivity tance oscillations while the occupancy represents a Coulomb bands as functions of Qf [17], [18] consistent with earlier staircase with Fermi-Dirac (FD) step shapes, so that a ion speculations [19], [14] and explaining both the anomalous channel can be thought as a discrete electrostatic device. mole fraction effect (AMFE) [2] and some of the puzzling In what follows with SI units e is the proton charge, z the mutation-induced transformations of selectivity in Ca2+ /Na+ ionic valence, Eo the vacuum permittivity, T the temperature, channels [6], [3]. We have already connected the periodicity of and kB Boltzmann's constant. the pattern of conduction bands with single/multi-ion barrier­ We consider the generic electrostatic model of the SF of less conduction and sequential neutralisation of the SF [18], a Ca2 + /Na+ ion channel shown in Fig. 1. It is described but the shapes of the occupancy and conduction bands and the as an axisymmetric, water-filled, cylindrical pore of radius R = 0.3nm and length L = l.6nm through a protein hub in 1 The protein residues are amino acids, of which aspartate (D) and glutamate (E) have negatively charged side chains while lysine (K) and arganine (R) have the cellular membrane. A centrally-placed, uniform, rigid ring positively charged side chains; we also mention neutral alanine (A). of negative charge Q f in the range 0 :::; IQf I :::; 7 e is embedded in the wall at RQ = R. The left-hand bath, modeling the 10 extracellular space, contains non-zero concentrations of Ca2 + ,.._.!!!. 0 5 and/or Na+ ions. We take both the water and the protein to ro ...,u be homogeneous continua with relative permittivities Ew = 80 0 and Ep = 2, respectively, but describe the ions as discrete 3 charges qi = ze within the framework of the implicit hydration c'.:l 2 model moving in single file within the channel, with bulk 0.. values of coefficients Di. We take no account of 0 negative counterions inside the SF, which will be few on ac­ j'.:' count of repulsion by the negative Q . The implicit hydration 20 1 ~"' model (10] works well for small Ca2+ and Na+ ions both of 3 10 n=O which have ionic radii Ri :=:;j O.lnm. 00 -1 -2 -3 -4 -5 -7 BD simulations solve the coupled 3D axisymmetrical Pois­ 0/e son electrostatic equation and lD overdamped Langevin stochastic equation numerically and self-consistently at each Fig. 2. (Color online) Brownian dynamics simulations of multi-ion conduction simulation step [9] , [4] , (20]. The model obviously represents and occupancy in a Ca2 + channel model vs the effective fixed charge Q f; (a),(b) are reworked from [17]. (a) Plots of the Ca2 + current J for pure a considerable simplification of the actual and Ca2 + baths of different concentration (20, 40 and 80mM as indicated). (b) dynamics of the ions and water molecules moving within The occupancy P. (c) The excess self-energy Un and ground state energy Uc the narrow SF (21]. However, reduced models successfully vs Qf for channels with n = 0, 1, 2 and 3 Ca2+ ions inside. The conduction bands Mnand the blockade/neutralisation points Zn are discussed in the text. reproduce significant properties of real biological channels (10], [9] thereby illustrating their applicability. Simple models can be also applied to artificial nanopores (16], (12], (21]. We can readily account for the pattern of bands in terms Details of the model, its range of validity and its limitations, of CB oscillations (13]. The discreteness of the ionic charge have been presented and discussed elsewhere [17], [18]. allows to us to introduce exclusive "eigenstates" { n} of the Fig. 1(b) shows the potential energy profiles for a single channel with fixed integer numbers of ions inside its SF having Ca2 + ion for the case Qf = - l e: a barrier-less effective en­ total electrostatic energy Un . The transition { n } ---+ { n+ 1} ergy profile appears due to the balance between the dielectric corresponds to the entry of a new ion, whereas { n } ---+ { n-1} self-energy barrier Us (see e.g. (16]) and and the site attraction corresponds to the escape of a trapped ion. These { n }-states Ua (18]: form a discrete exclusive set (23] :

(1) n = {O , 1, 2, .. . } Len= 1; Pc= L nen, (2) n n 2 Here, Cs = 41rE0 Ew R / L stands for the geometry-depend where en is the occupancy of the state { n} and Pc is the channel self-capacitance. The captured ion neutralizes Q1 and average SF occupancy. In equilibrium em is defined by the restores the effective barrier, thereby preventing the next ion Boltzmann factor en ex exp(-Un/(kBT)). The exact distri­ from entering the channel. When Us » kBT this amounts bution for en and Pc (which is FD (24], (25]) will be derived to a single-occupancy requirement or electrostatic exclusion below. principle. It plays a significant role in ion permeation and leads The total energy Un for a channel in state { n} can be to the CB phenomenon, causing the ion channel to behave as expressed as: single-charge device (14], (15]. (3) The single- and multi-ion conduction bands found in the BD Un = Un,s + Un,uttr + Un,int simulations [17], (18] are shown in Fig. 2(a),(b) which plot where Un,s is the self-energy, Un,attr is the energy of attrac­ the Ca2 + current J and channel occupancy P for pure baths tion, and Un ,int is the ions' mutual interaction energy. of different concentration. Fig. 2(a) shows narrow conduction We approximate Un as the dielectric self-energy Un ,s of the bands M 0 , M 1, M 2 separated by stop-bands of almost zero­ excess charge Qn based on the assumption that both the ions conductance centred on the blockade points Z1 , Z2 , Z3 . and Qf are located within the central part of the SF, so that Fig. 2(b) shows that the M n peaks in J correspond to tran­ (1) gives the standard CB equation: sition regions in channel occupancy, where P jumps from one (4) integer value to the next, and that the stop-bands correspond to saturated regions with integer P = 1, 2, 3 .... Here, Qn represents the excess charge at the SF for then ions Band Mo corresponds to single-ion barrier-less conduction as function of Q1. Binomial expansion of Q~ in (4) gives first

(see Fig. l(b)). M 1 corresponds to the double-ion knock­ approximations for Un ,s , Un ,attr and Un,int consistent with 1- on conduction, which is well-established for L-type Ca2 + D Coulomb gas model (16] and with the energetics analysis channels [2]; and M 2 corresponds to triple-ion conduction in (18]. which can be connected with Ryanodine receptor calcium CB appears in low-capacitance systems from quantization channels (22] (see Table I and related discussion). of the quadratic energy form (4) at a grid of discrete states TABLE I IDENTIFICATION OF C A2+ RESONANT CONDUCTION POINTS Mn WITH SOME KNOWN CHANNELS AND MUTANTS (EXTENDED FROM [17], [18])).

Channel Locus Nominal QJ Band I Jump Modelled QJ Non-selective OmpF [8] (K)RRRDE - l e Mo {O -+ l} - l e Ca- selective Nav mutant [6] EEEA -3e M1 {l -+ 2} -3e Ca L-type channel [2] , Ca-selective OmpF mutant [8] EEEE -4e M1 {l -+ 2} -3e Calcium RyR [22] DDDD(ED) - 6e M2 {2 -+ 3} -5e Calcium-selective mutants of NavAB [3] and NavChBac [7] EEEE+DDDD -Se M3 {3 -+ 4} - 7e

(2), providing a Coulomb energy gap D.Un = Us large enough Ca2 + /Na+ channels family, thereby confirming the validity of (D.Un » k 8 T) to block transitions between neighbouring {n} the CB-based band model. We speculate that the discrepancy states. We calculate Un as a function of Q1 for n = 0, 1, 2, 3 between Q1 = M1 = - 3e derived from CB model and the and focus on the ground state energy Uc as a function of Qr nominal Qf = -4e for the EEEE loci of L-type calcium channels may be associated with protonation (26]. (5) Derivation of the distribution of Pc in the vicinity of the Mn points (and hence calculation of the shapes of Pc(Uc) or Fig. 2(c) plots Un and Uc as functions of Q1. We see a Pc( Qf)) follows standard CB theory. For divalent Ca2+ ions, periodic pattern with two kinds of Uc singular points, marked the Coulomb gap D.Un » k T and (2) reduces to a simple as M n and Zn. The minima of Uc (and the blockade regions) 8 2-state exclusive set: appear around the neutralisation points Zn = - zen where Qc = 0 and the occupancy Pc is saturated at an integer value m = {n, n + l}; Bn+ Bn+1 = l; Pc= n + Bn+l· (7) 2 [16], [18]. For divalent Ca + ions D.Un :=:;j 20ksT and hence CB is strong. The crossover points Mn (Un = Un+1) allow The strong electrostatic exclusion principle (7) plays the same role as the Pauli exclusion principle plays in quantum mechan­ barrier-less { n} '.:::::+ { n + 1} transitions; they correspond to the Pc transition regions and to the conduction peaks in J (18]. ics (25], [27]. Hence the standard derivation via a partition The positions of the singular Qf points in Fig. 2(c) can be function, taking account of (7) leads (24] to FD statistics for written as: Bn+l and to an excess (fractional) occupancy P; = Pc mod 1: 1 Zn = -zen ± i5Zn, Coulomb blockade P; = (1 + pb- l exp(Uc/ksT))- , Fermi-Dirac (8) (6) Mn = -ze(n + 1/ 2) ± i5Mn Resonant conduction where Uc = Un+l - Un, and Pb is a reference occupancy where i5Zn, i5Mn are possible corrections for the singular part related to the bulk concentration. Note, that (8) is equivalent to of the affinity and ion-ion interactions, not accounted for here. the Langmuir isotherm [24] or to Michaelis-Menten saturation. Equation (6) is exactly the same as its counterpart in electronic A similar Fermi function was obtained earlier (16] for the CB. We may therefore interpret the ED-simulated conduction variation of Pc with concentration. bands (Fig. 2(a)) as CB oscillations of conductance (13] A self-consistent calculation of the conductance can be appearing as Qf increases and the corresponding occupancy effected via linear response theory, leading to the standard steps (Fig. 2(b)) as a Coulomb staircase [15]. CB "classical" approximation [13]: 1 The positions of the Mn and Zn points in the theory (6 ) l e/ lmax = (Uc/ksT) sinh- (Uc/ksT) (9) and BD simulations (Fig. 2) are consistent with an energetics analysis [18], supporting above interpretation; the deviations where l max is the barrier-less diffusive current. The current in the precise positions of Mn and Zn can be attributed to Jc from (9) exhibits a resonant peak coinciding with the field leaks and model simplifications. maximum in the derivative of Pc, dP/ dUc (Fig. 3), similar That is our main result - that the ionic CB model to that of the tunneling current in a quantum dot [13]. The predicts a universal, valence-dependent, periodic pattern of FD function (8) predicts that u; / ksT = - ln(P; / (1 - P;)) stop/conduction bands similar to the electronic CB oscillations of conductance in quantum dots (13], (15]. It allows us _ J/Jma,' "Classical" to identify conduction bands with real calcium-conductive 1 x - . - . Pc' Fermi-Dirac channels/mutants. ~ Table I shows putative identifcations of the ·' ~~u 0.5 ,' bands/singularities of the CB model with real channels [l_ u in the Ca2 + /Na+ family, extended from [ 17], [ 18]. We identify wild-type calcium channels and Ca2+ -selective mutants of Na-selective (Nav) channels with the Ca2 + resonant points M1 (L-type) and M2 (RyR). We also identify Ca2+ -selective mutants of bacterial sodium channels Fig. 3. (Color online) Resonant conduction peak shapes vs Uc (or equivalently Q f ). The occupancy Pc (blue, dash-dot) shows the Fermi - Dirac transition NaChBac [7] and NavAB [3] with point M 3 . In a similar way, from Pc = 0 to Pc = 1. The theoretical current Jc/ Jmax (red,solid) exhibits we can account for numerous mutation transformations in a resonant peak in the transitional region. [4] B. Roux, T. Allen, S. Berneche, and W. Im, 'Theoretical and computa­ tional models of biological ion channels," Quart. Rev. Biophys. , vol. 37, no. 1, pp. 15-103, 2004. [5] V. N. Kharkyanen, S. 0. Yesylevskyy, and N. M. Berezetskaya, "Approx­ imation of super-ions for single-file diffusion of multiple ions through narrow pores," Phys. Rev. E, vol. 82, p. 051103, 2010. [6] S. H. Heinemann, H. Teriau, W. Stuhmer, K. Imoto, and S. Numa, .~;:r•••••·· ~ ··•••••••• N •••••••••••• S ••'"'••j "Calcium-channel characteristics conferred on the sodium-channel by single mutations," Nature, vol. 356, no. 6368, pp. 441-443, 1992. [7] L. X. Yue, B. Navarro, D. J. Ren, A. Ramos, and D. E. Clapham, 'The cation selectivity filter of the bacterial sodium channel, NaChBac," J. Gen. Physiol. , vol. 120, no. 6, pp. 845- 853, 2002. [8] H. Miedema, A. Meter-Arkema, J. Wierenga, J. Tang, B. Eisenberg, (r0 !\-1 -2l ,f -3 \ -4 .. -5 J\ -6 ~ ' -7 W. Nonner, H. Hektor, D. Gillespie, and W. Meijberg, "Permeation o;e properties of an engineered bacterial OmpF porin containing the EEEE­ locus of Ca2 + channels," Biophys. J., vol. 87, no. 5, pp. 3137-3147, 2004. Fig. 4. (Color online) Comparisons of the ionic Coulomb blockade model [9] B. Corry, T. Vora, and S. H. Chung, "Electrostatic basis of valence with Brownian dynamics simulation results as Qf is varied. (a) The effective selectivity in cationic channels," Biochim. Biophys. Acta-Biomembranes, well depth U~ (blue point-down triangles) fitted by Fermi-Dirac function (full vol. 1711 , no. 1, pp. 72-86, 2005. red lines). (b) The J peaks in the"classical" approximation (red solid lines) [10] D. Boda, "Monte Carlo simulation of electrolyte solutions in biology: are compared with simulation results (blue, point-up triangles). in and out of equilibrium," in Ann. Rep. Comp. Chem., ser. Ann. Rep. Comp. Chem., D. Gillespie and N. Baker, Eds. Elsevier, 2014, vol. 10, ch. 5, pp. 127-164. should be linear in Uc and also (due to the relatively narrow [11] B. Corry, "Na+ /Ca2+ selectivity in the bacterial voltage-gated sodium channel NavAb," Peer], vol. 1, no. DOil0.7717/peerj.16, p. el6 transition region) in Q f. Fig. 4(a) plots the predicted sawtooth (DOil0.7717/peerj.16), 2013. dependence of u; on Qf, confirming that the P; transitions [12] M. Krems and M. Di Ventra, "Ionic Coulomb blockade in nanopores," obey the FD function (8) of Uc. Fig. 4(b) compares the CB J. Phys. Condens. Matter, vol. 25, p. 065101 , 2013. [13] C. W. J. Beenakker, "Theory of Coulomb-blockade oscillations in the 2 , model with the the ED-simulated Ca + conduction bands M 0 conductance of a quantum dot," Phys. Rev. B, vol. 44, no. 4, pp. 1646- M 1 , M 2 [17]. Theoretical CB oscillations (9) fit the BD­ 1656, 1991. simulated peak shapes and positions reasonably well, given the [14] M. Stopa, "Rectifying behavior in Coulomb blockades: charging recti­ fiers," Phys. Rev. Lett., vol. 88, no. 14, p. 146802, 2002. model simplifications. Although an ion moving inside a chan­ [15] K. K. Likharev, "SET: Coulomb blockade devices," Nano et Micro Tech. , nel or is a classical system described by Newtonian vol. 3, no. 1-2, pp. 71-114, 2003. dynamics, it exhibits some quantum-like mesoscopic features [16] J. Zhang, A. Kamenev, and B. I. Shklovskii, "Ion exchange phase transitions in water-filled channels with charged walls," Phys. Rev. E, [12]. We attribute such behavior to charge discreteness, the vol. 73, p. 051205, 2006. electrostatic exclusion principle, and confinement effects. [17] I. K. Kaufman, D. G. Luchinsky, R. Tindjong, P. V. E. McClintock, In conclusion, we have shown that Ca2+ channel permeation and R. S. Eisenberg, "Multi-ion conduction bands in a simple model of calcium ion channels," Phys. Biol., vol. 10, no. 2, p. 026007, 2013. is analogous to mesoscopic transport in quantum dots: the [18] --, "Energetics of discrete selectivity bands and mutation-induced electrostatic exclusion principle leads to an FD distribution of transitions in the calcium-sodium ion channels family," Phys. Rev. E, channel occupancy; the stop-bands correspond to blockade; the vol. 88, no. 5, p. 052712, 2013. [19] R. S. Eisenberg, "Atomic biology, electrostatics and ionic channels," barrier-less conduction peaks are similar to those in resonant in New Developments in Theoretical Studies of Proteins, R. Elber, Ed. tunneling and can be described by standard CB formul~. Singapore: World Scientific, 1996, pp. 269-357. The ionic CB model provides a good account of both the [20] C. Berti, S. Furini, D. Gillespie, D. Boda, R. S. Eisenberg, E. Sangiorgi, and C. Fiegna, "Three-dimensional Brownian dynamics simulator for experimental (AMFE and valence selectivity) and the sim­ the study of ion permeation through membrane pores," J. Chem. Theo r. ulated (discrete multi-ion conduction and occupancy bands) Comp. dx.doi.org/10.1021/ct4011008, 2014. phenomena observed in model Ca2+ channels. The results [21] A. Laio and V. Torre, "Physical origin of selectivity in ionic channels of biological membranes," Biophys. J., vol. 76, no. 1, pp. 129- 148, 1999. should be applicable to other ion channels and to biomimetic [22] D. Gillespie, "Energetics of divalent selectivity in a calcium channel: The nanopores with charged walls. ryanodine receptor case study," Biophys. J. , vol. 94, no. 4, pp. 1169- We are grateful to Aneta Stefanovska, Boris Shklovskii 1184, 2008. [23] P. H. Nelson, "A permeation theory for single-file ion channels: One­ and Max Di Ventra for helpful comments and discussions. and two-step models," J. Chem. Phys. , vol. 134, pp. 165 102-165 114, The research was supported by the Engineering and Physical 2011. Sciences Research Council UK (grant No. EP/M015831/1). [24] R. H. Fowler, "A statistical derivation of Langmuir's adsorption isotherm," in Mathematical Proceedings of the Cambridge Philosophical No new data were created during this study. Society, vol. 31. Cambridge: CUP, 1935, pp. 260- 264. [25] G. Kaniadak:is and P. Quarati, "Kinetic equation for classical particles REFERENCES obeying an exclusion principle," Phys. Rev. E, vol. 48, no. 6, pp. 4263- 4270, 1993. [l] B. Hille, Ion Channels Of Excitable Membranes, 3rd ed. Sunderland, [26] S. Furini, P. Barbini, and C. Domene, "Effects of the protonation state MA: Sinauer Associates, 2001. of the EEEE motif of a bacterial Na+ -channel on conduction and pore [2] W. A. Sather and E. W. McCleskey, "Permeation and selectivity in structure," Biophys. J. , vol. 106, pp. 2175-2183, 2014. calcium channels," Ann. Rev. Physiol. , vol. 65, pp. 133-159, 2003. [27] J.-L. Liu and R. Eisenberg, "Correlated ions in a calcium channel model: [3] L. Tang, T. M. G. El-Din, J. Payandeh, G. Q. Martinez, T. M. Heard, A Poisson-Fermi theory," J. Phys. Chem. B, vol. 117, pp. 12051-12 058, T. Scheuer, N. Zheng, and W. A. Catterall, "Structural basis for Ca2 + 2013. selectivity of a voltage-gated calcium channel," Nature, vol. 505, no. 7481, pp. 56- 61 , 2014.