Abstract
Voltage and charge distributions in cellular microdomains regulate communications, excitability, and signal transduction. We report here new electrical laws in a biological cell, which follow from a nonlinear electro-diffusion model. These newly discovered laws derive from the geometrical cell-membrane properties, such as membrane curvature, volume, and surface area. The electro-diffusion laws can now be used to predict and interpret voltage distribution in cellular microdomains such as synapses, dendritic spine, cilia and more.
Similar content being viewed by others
Introduction
Electro-diffusion in cellular microdmains remains difficult to study due to the absence of specific sensors and the theoretical hurdle of understanding the dynamics of charged particle in shaped domains. The diffusion of charged particles is largely influenced by the interaction of diffusing ions with the electrical field generated by all charges in the solution and possibly with external field.
The dielectric membrane of a charged biological cell also affects the electric field, because it creates image charges (see ref. 1 for an infinite plan). So far, only a few electro-diffusion systems are well understood. For example, although the electrical battery was invented more than 200 years ago, designing optimal configurations is still a challenge. On the other extreme, ionic flux and gating of voltage-channels2 is now well explained by the modern Poisson-Nernst-Planck (PNP) theory of electro-diffusion, because at the nanometer scale, cylindrical symmetry of a channel model reduces computations to a one-dimensional model for the electric field and charge densities in the channel pore3,4,5. However, cellular microdomains involve two- and three-dimensional neuronal geometry7,8,9,10,11, which makes the analysis of the PNP equations much more complicated than in the cylindrical geometry of a channel pore12.
We report here results about the distribution of charges and field, obtained from the analysis of the nonlinear PNP model of electro-diffusion in various geometries of microdomains in the absence of electro-neutrality. Although at this stage, we are still lacking experimental proofs of the nonelectro-neutrality at a scale of hundreds of nanometers, there is a clear unbalance in the concentrations of positive versus negative ions (see discussion below). This hypothesis should be considered as a prediction and is the basis of the present computations. In our model the entire boundary is impermeable to particles (ions) and the electric field satisfies the compatibility condition of Poisson’s equation. Phenomenological descriptions of the electro-diffusion, such as cable equations or the reduced electrical-engineering approximation by resistance, capacitance, and even electronic devices, are insufficient for the description of non-cylindrical geometry13, because they assume simple one-dimensional or reduced geometry.
Results
Electrostatic theory with no electro-neutrality
In the absence of electro-neutrality and with N charges distributed in a bounded domain Ω surrounded by a dielectric membrane, the PNP model for total charge Q = zeN, where e is the electron charge (1.6 · 10−19 C), z - the valence, and the charge density ρ(x, t) is given by ref. 14 for x ∈ Ω,
where D is the diffusion coefficient, kT represents the thermal energy and V(x, t) is the electric potential in Ω. It is the solution of the Poisson equation
with the boundary condition
where the surface charge density σ(x, t) is defined on the boundary ∂Ω using the the electrical permitivity εε 0 of the electrolyte solution. In the steady state,
Then (1) gives the density
hence (4) becomes
and (5) gives the boundary condition
which is the compatibility condition, obtained by integrating Poisson’s equation (4) over Ω.
Debye’s15 concept of charge screening, which makes the induced field decay exponentially fast away from a charge, does not apply when electroneutrality is broken and long-range correlations lead to a gradient of charges, as is the case, for example, in a ball without inward directed current. By solving (8) numerically and asymptotically, a new capacitance law was derived for an electrolytic solution in a ball of radius R 14, where the difference of potentials between the center C and any point of the spherical surface S, V(C) − V(S), increases with the total number of charges, first linearly and then logarithmically
The effect of the geometry on the voltage and the charge distribution for other cell shapes is described below.
Local boundary curvature affects field and charge distribution
Axons and dendrites are not perfect cylinders and the curvature of their surfaces has many local maxima16. It turns out that this local curvature can influence the local voltage significantly, as shown in numerical solutions (using Comsol classical packages) of the PNP equations (Fig. 1), which reveal that regions of high curvature correspond to local charge accumulation. This effect should be sufficient to influence the voltage by creating measurable local voltage increase of the order of a few millivolts. The voltage inside the cylinder and along curved surfaces can vary (Fig. 1A–D) and can also depend on changes in curvature (Fig. 1F–I) and in the total number of charges N (Fig. 1E–J). Curvature creates a narrow boundary layer in the voltage (Fig. 1H,I).
The previous computations assume one specie only and in practice, due to the presence of N other ions (i.e., λ = (ze)2 N/εε 0 kT), should only represent the excess of positive charges. Thus, we explore a large range for the parameter λ (a generalization of Bjerrum length l B = (ze)2/εε 0 kT, but has no physical units in dim 2 and is a length in dim 3). Note that λ is not simply given by the difference of positive and negative charges, due to the non-linear PNP equations. From the exploration of the graphs in Figs 1 and 2, we expect λ ≈ 1–500, which corresponds to a range N ≈ 102–5 · 104, leading to voltage fluctuations of few tens of mV. A more physiological range for the parameter λ is shown in Fig. 1E (inset). Consequently, the resting cross-membrane potential and the resting voltage across voltage-gated channels along a dendrite may vary, depending on curvature. This may affect the propagation and genesis of local depolarization or back-propagation action potential in dendrites of neuronal cells17.
A cusp-shaped funnel can influence charge distribution
The passage of ions and other particles between different cellular compartments is done through narrow passages that form cusp-shaped funnels, which have negative curvature (see the example described in Fig. 2A). We consider here the case of dendritic spines, whose geometry is approximated as a convex domain (head) connected by a cylindrical neck to the dendritic shaft. This geometrical model approximation was already validated for diffusion18. We observe that the cusp-shaped funnel prevents the entrance of charges, at least when their number does not exceed a given threshold (Fig. 2C–F). Therefore, there is a difference of steady state potential drop between the end of the funnel and the bulk of the domain.
Assume that for a dimensionless domain \({\rm{\Omega }}\subset {{\mathbb{R}}}^{2}\) with a cusp-shaped funnel F formed by two bounding circles of equal radii R c (see Fig. 2A), the opening of the funnel is \(\tilde{\varepsilon }\ll 1\) (in unit of micrometers). This radius is used in changing to dimensionless coordinates, so we may assume R c = 1. The surface area of the domain boundary is |∂Ω|. For an uncharged funnel domain, the condition (5) becomes ∂V(x)/∂n = 0 for x ∈ F. The potential drop between the center of Ω and the end of the cusp is found by solving (8) in a domain obtained by mapping Ω conformally by the Möbius transformation18
where \(\alpha =-1-\sqrt{\tilde{\varepsilon }/{R}_{c}}+O(\varepsilon )\) and R c is the (dimensional) radius of curvature at the cusp domain (Fig. 2A). In the variable w = Re iθ = X + iY, the dimensional solution in the image of the cusp is given by
Note that in the planar case λ has no dimensional units. It leads for \(\lambda \gg 1\), to the dimensional voltage drop
This difference should be compared with that between the north pole N and the center C, given in ref. 14, is given by
where λ D = 2πλR l /|∂Ω| (R l is the (dimensional) radius of the external domain Ω). For a three-dimensional symmetric domain with a cusp-funnel,
We conclude in the limit of \(\lambda \gg 1\) and \(\tilde{\varepsilon }\to 0\) that the difference of potential between the end of cusp S and the north pole N in the domain is obtained by adding (13) and (12) and we get
When Ω is charged, the boundary condition is still (9) and the potential drop is
The constant c depends only on the geometric center of mass C such as \(w(C)=c\sqrt{\frac{\tilde{\varepsilon }}{{R}_{c}}}\), where w is the Möbius transformation (10), R c is the curvature at the cusp and \(\tilde{\varepsilon }\) the width at the base of the cusp. The analytical solution in the charged cusp-funnel domain is (in the same conformal coordinates as in eq. 11)
which is in very good agreement with the numerical solution (Fig. 2C dotted orange). The present results are quite different from the Gouy-Chapman theory which concerns the exponential decay of voltage away from the double layer on a flat surface and derived for electro-neutral bulk, which is not the case here (see wiki “double layer”). It follows that the funnel curvature is an interesting free parameter for the design of optimized nano-pipettes to regulate the molecular and ionic fluxes. These pipettes were recently designed to record voltage in dendritic spines19. Thus considering the shape of the pipette as a new parameter to control the ion fluxes injected in small biological domain, could impact the design and the precision of future biological experiments.
Voltage distribution in an elongated ellipse
Numerical solutions of the PNP equations 8 (Fig. 3A) show the effect of local charge accumulation in an elongated ellipse, in particular the potential difference V(A) − V(B) is maximal on the ellipse for A and B the ends of the small and long axes (Fig. 3B,C). The voltage along each axis is shown in Fig. 3D–G. The elliptic domain can represent the cross section of an axon, and confirms the effect of curvature discussed in the previous sections. Specifically, that charges accumulate near the boundary of highest curvature. The equipotential contours are shown in Fig. 3G–L.
Discussion and Conclusions
As shown in the present study, inside a non-flat electrolyte domain characterized by significant changes in membrane curvature, the difference of voltage varies with the log of the curvature. Specifically, we derived various formula 16, 17 and 18 at a cusp-shaped funnel (curvature) R c and we showed that the changes are significant in the range R c ∈ [0.1–2] μm (Fig. 2E). These formula should be considered as predictions to study the effect of curvature inside neuronal microdomains, such as dendritic spines or part of dendrites.
We found here that local changes in geometry, non-electro-neutrality, and a dielectric boundary affect charge distribution in electro-diffusion of electrolytes, as previously suggested in refs 13, 14 and 20. This paper considers a single ionic specie, whereas in reality the solution contains multiple positive and negative ions. The present electro-diffusion model captures the effect of the excess of positive ions. Consequently, the large value of these voltage differences that we have found in the large charge limit λ → ∞ are probably much attenuated in a mixed positive and negative ionic solution, but the electro-neutrality remains broken. For large λ, the largest voltage differences occur near the boundary as shown in Figs 1B–D–I and 2B, and Fig. 3K, but this property will certainly remain true even if the difference between positive and negative ions decreases. The assumption of non-electroneutrality is based on the known cytoplasmic ionic concentrations given by Na + = 12 mM, K + = 155 mM and Cl − = 4.2 mM 6. There is a clear unbalance toward positive charges. Are negative charges missing? Probably there are molecules of various sizes with negative charges, but their diffusion coefficients are certainly much smaller than these of the ions. This difference of mobility is certainly a key feature in maintaining non-electro-neutrality and should be investigated both by ionic numerical simulations of electrolytes and experimentally. To conclude, the present analysis by considering only one positive specie, leads to an excessive potential difference, that will be diminished when adding negative charged proteins, but it provides a new qualitative insight, suggesting that the geometry of dendritic spines can influence significantly the voltage amplitude and shape of the excitatory-post-synaptic-potential, which depends on the membrane curvature. This prediction suggests that the neuronal code and communication is directly modulated by the geometry of nano-micro synaptic domains. This effect could have consequences in other electrical cellular domains, such as cardiac myocytes, glial protrusions, cilia and more. Indeed, local curvature is associated with a gradient of charge density that can affect the electrical properties of micro-compartments10, 11, 21,22,23,24. In all cases, charge accumulates near boundary points of locally maximal curvature.
These results and predictions can further be used to design nano-devices such as pipettes and to better understand voltage changes inside dendrites and axons. Finally, the non-electroneutrality could also be amplified during a synaptic current flow with an amplitude of 50 pA during 5 ms, leading to an entry of roughly 1.5106 (positive charges), which cannot be compensated with the membrane capacitance. Indeed, the capacitance effect could at most neutralize 1.6% of the positive charges (the value is computed from the formula \({N}_{m}=\frac{{\varepsilon }_{0}{\varepsilon }_{r}{S}_{urf}}{de}V\), with a dielectric membrane thickness d = 8 nm, the electric permittivity ε 0 ε r = 4.410−11 F/m, ε r = 5, the membrane potential V = 60 mV and the surface \({S}_{urf}\approx 4\pi {R}_{bulk}^{2}\)). Future analysis should reveal the charge distribution generated by a transient current.
References
Abrashkin, A. & Orland, H. “Dipolar Poisson-Boltzmann equation: ions and dipoles close to charge interfaces”. Physical Review Letters 99(7), 077801 (2007).
Bezanilla, F. “How membrane proteins sense voltage”. Nat. Rev. Mol. Cell Biol. 9, 323–32 (2008).
Eisenberg, R. S. “Ionic channels in biological membranes. Electrostatic analsis of a natural nanotube”. Contemp. Phys. 39(6), 447–466 (1998).
Eisenberg, R. S. “From structure to function in open ionic channels”. J. Membrane Biol. 171, 1–24 (1998).
Eisenberg, R. S., Klosek, M. M. & Schuss, Z. “Diffusion as a chemical reaction: Stochastic trajectories between fixed concentrations”. J. Chem. Phys. 102(4), 1767–1780 (1995).
Hille, B. Ion Channels of Excitable Membranes, Third Ed (2001).
Lee, P., Sobie, E. A. & Peskin, C. S. “Computer simulation of voltage sensitive calcium ion channels in a dendritic spine”. J. Theor. Biol. 338, 87–93 (2013).
Lopreore, C. L. et al. “Computational modeling of three-dimensional electrodiffusion in biological systems: application to the node of Ranvier”. Biophys. J. 15(95(6)), 2624–35 (2008).
Savtchenko, L. P., Kulahin, N., Korogod, S. M. & Rusakov, D. A. “Electric fields of synaptic currents could influence diffusion of charged neurotransmitter molecules”. Synapse 15(51(4)), 270–8 (2004).
Sylantyev, S. et al. “Electric fields due to synaptic currents sharpen excitatory transmission”. Science 28(319(5871)), 1845–9 (2008).
Sylantyev, S., Savtchenko, L. P., Ermolyuk, Y., Michaluk, P. & Rusakov, D. A. “Spike-driven glutamate electrodiffusion triggers synaptic potentiation via homer-dependent mGluR-NMDAR link”. Neuron 77, 528–541 (2013).
Qian, N. & Sejnowski, T. J. “An electro-diffusion model for computing membrane potentials and ionic concentrations in branching dendrites, spines and axons”. Biological Cybernetics 62.1, 1–15 (1989).
Holcman, D. & Yuste, R. “The new nanophysiology: regulation of ionic flow in neuronal subcompartments”. Nature Reviews Neuroscience 16, 685–692 (2015).
Cartailler, J., Schuss, Z. & Holcman, D. “Analysis of the PoissonNernstPlanck equation in a ball for modeling the voltagecurrent relation in neurobiological microdomains”. Physica D: Nonlinear Phenomena 339, 39–48 (2017).
Debye, P. & Hückel, E. “Zur Theorie der Elektrolyte. I. Gefrierpunktserniedrigung und verwandte Erscheinungen”. Physikalische Zeitschrift 24(9), 185–206 (1923).
Bourne, J. N. & Harris, K. M. “Balancing structure and function at hippocampal dendritic spines”. Annu Rev Neurosci. 31, 47–67 (2008).
Markram, H., Lubke, J., Frotscher, M. & Sakmann, B. “Regulation of synaptic efficacy by coincidence of postsynaptic APs and EPSPs”. Science 275(5297), 213–5 (1997).
Holcman, D. & Schuss, Z. Stochastic Narrow Escape in Molecular an Cellular Biology, Analysis and Applications, Springer Verlag, NY 2015 (2015).
Jayant, K. et al. “Targeted intracellular voltage recordings from dendritic spines using quantum-dot-coated nanopipettes”. Nat. Nanotechnol. (2016).
Mamonov, A., Coalson, R., Nitzan, A. & Kurnikova, M. “The role of the dielectric barrier in narrow biological channels: a novel composite approach to modeling single channel currents”. Biophys. J. 84, 3646–3661 (2003).
Yuste, R. Dendritic Spines, The MIT Press, Cambridge, MA (2010).
Araya, R., Jiang, J., Eisenthal, K. B. & Yuste, R. “The spine neck filters membrane potentials”. Proc. Natl. Acad. Sci. USA 103, 17961–17966 (2006).
Araya, R., Vogels, T. P., T. P. & Yuste, R. “Activity-dependent dendritic spine neck changes are correlated with synaptic strength. Proc. Natl. Acad. Sci. USA 111, 2895–904 (2014).
Lan, W. J. et al. Voltage-Rectified Current and Fluid Flow in Conical Nanopores. Accounts of Chemical Research 49(11), 2605–2613 (2016).
Acknowledgements
The research of DH research was partially supported by a Marie-Curie Award grant.
Author information
Authors and Affiliations
Contributions
J.C. and D.H. made to the figures and design the simulations and the computations. D.H., Z.S. and J.C. design the project, made the analytical derivation and wrote the main manuscript text and D.H. and J.C. prepared Figures 1–3. All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Competing Interests
The authors declare that they have no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Cartailler, J., Schuss, Z. & Holcman, D. Electrostatics of non-neutral biological microdomains. Sci Rep 7, 11269 (2017). https://doi.org/10.1038/s41598-017-11590-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-017-11590-6
This article is cited by
-
An Algorithm Based on a Cable-Nernst Planck Model Predicting Synaptic Activity throughout the Dendritic Arbor with Micron Specificity
Neuroinformatics (2023)
-
Modeling the voltage distribution in a non-locally but globally electroneutral confined electrolyte medium: applications for nanophysiology
Journal of Mathematical Biology (2021)
-
Reply to ‘Only negligible deviations from electroneutrality are expected in dendritic spines’
Nature Reviews Neuroscience (2020)
-
Electrodiffusion models of synaptic potentials in dendritic spines
Journal of Computational Neuroscience (2019)
-
Steady-state voltage distribution in three-dimensional cusp-shaped funnels modeled by PNP
Journal of Mathematical Biology (2019)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.