1
|
Becker M, Loche P, Rezaei M, Wolde-Kidan A, Uematsu Y, Netz RR, Bonthuis DJ. Multiscale Modeling of Aqueous Electric Double Layers. Chem Rev 2024; 124:1-26. [PMID: 38118062 PMCID: PMC10785765 DOI: 10.1021/acs.chemrev.3c00307] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2023] [Revised: 11/17/2023] [Accepted: 11/30/2023] [Indexed: 12/22/2023]
Abstract
From the stability of colloidal suspensions to the charging of electrodes, electric double layers play a pivotal role in aqueous systems. The interactions between interfaces, water molecules, ions and other solutes making up the electrical double layer span length scales from Ångströms to micrometers and are notoriously complex. Therefore, explaining experimental observations in terms of the double layer's molecular structure has been a long-standing challenge in physical chemistry, yet recent advances in simulations techniques and computational power have led to tremendous progress. In particular, the past decades have seen the development of a multiscale theoretical framework based on the combination of quantum density functional theory, force-field based simulations and continuum theory. In this Review, we discuss these theoretical developments and make quantitative comparisons to experimental results from, among other techniques, sum-frequency generation, atomic-force microscopy, and electrokinetics. Starting from the vapor/water interface, we treat a range of qualitatively different types of surfaces, varying from soft to solid, from hydrophilic to hydrophobic, and from charged to uncharged.
Collapse
Affiliation(s)
| | - Philip Loche
- Fachbereich
Physik, Freie Universität Berlin, 14195 Berlin, Germany
- Laboratory
of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Majid Rezaei
- Fachbereich
Physik, Freie Universität Berlin, 14195 Berlin, Germany
- Institute
of Theoretical Chemistry, Ulm University, 89081 Ulm, Germany
| | | | - Yuki Uematsu
- Department
of Physics and Information Technology, Kyushu
Institute of Technology, 820-8502 Iizuka, Japan
- PRESTO,
Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
| | - Roland R. Netz
- Fachbereich
Physik, Freie Universität Berlin, 14195 Berlin, Germany
| | - Douwe Jan Bonthuis
- Institute
of Theoretical and Computational Physics, Graz University of Technology, 8010 Graz, Austria
| |
Collapse
|
2
|
|
3
|
Hennequin T, Manghi M, Palmeri J. Competition between Born solvation, dielectric exclusion, and Coulomb attraction in spherical nanopores. Phys Rev E 2021; 104:044601. [PMID: 34781526 DOI: 10.1103/physreve.104.044601] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2021] [Accepted: 09/15/2021] [Indexed: 11/07/2022]
Abstract
The recent measurement of a very low dielectric constant, ε, of water confined in nanometric slit pores leads us to reconsider the physical basis of ion partitioning into nanopores. For confined ions in chemical equilibrium with a bulk of dielectric constant ε_{b}>ε, three physical mechanisms, at the origin of ion exclusion in nanopores, are expected to be modified due to this dielectric mismatch: dielectric exclusion at the water-pore interface (with membrane dielectric constant, ε_{m}<ε), the solvation energy related to the difference in Debye-Hückel screening parameters in the pore, κ, and in the bulk κ_{b}, and the classical Born solvation self-energy proportional to ε^{-1}-ε_{b}^{-1}. Our goal is to clarify the interplay between these three mechanisms and investigate the role played by the Born contribution in ionic liquid-vapor (LV) phase separation in confined geometries. We first compute analytically the potential of mean force (PMF) of an ion of radius R_{i} located at the center of a nanometric spherical pore of radius R. Computing the variational grand potential for a solution of confined ions, we then deduce the partition coefficients of ions in the pore versus R and the bulk electrolyte concentration ρ_{b}. We show how the ionic LV transition, directly induced by the abrupt change of the dielectric contribution of the PMF with κ, is favored by the Born self-energy and explore the decrease of the concentration in the pore with ε both in the vapor and liquid states. Phase diagrams are established for various parameter values and we show that a signature of this phase transition can be detected by monitoring the total osmotic pressure as a function of R. For charged nanopores, these exclusion effects compete with the electrostatic attraction that imposes the entry of counterions into the pore to enforce electroneutrality. This study will therefore help in deciphering the respective roles of the Born self-energy and dielectric mismatch in experiments and simulations of ionic transport through nanopores.
Collapse
Affiliation(s)
- Théo Hennequin
- Laboratoire de Physique Théorique, Université de Toulouse, CNRS, UPS, 31062 Toulouse cedex 4, France
| | - Manoel Manghi
- Laboratoire de Physique Théorique, Université de Toulouse, CNRS, UPS, 31062 Toulouse cedex 4, France
| | - John Palmeri
- Laboratoire Charles Coulomb (L2C), Université de Montpellier, CNRS, 34095 Montpellier cedex 5, France
| |
Collapse
|
4
|
Liu X, Lu B. Incorporating Born solvation energy into the three-dimensional Poisson-Nernst-Planck model to study ion selectivity in KcsA K^{+} channels. Phys Rev E 2017; 96:062416. [PMID: 29347452 DOI: 10.1103/physreve.96.062416] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2017] [Indexed: 06/07/2023]
Abstract
Potassium channels are much more permeable to potassium than sodium ions, although potassium ions are larger and both carry the same positive charge. This puzzle cannot be solved based on the traditional Poisson-Nernst-Planck (PNP) theory of electrodiffusion because the PNP model treats all ions as point charges, does not incorporate ion size information, and therefore cannot discriminate potassium from sodium ions. The PNP model can qualitatively capture some macroscopic properties of certain channel systems such as current-voltage characteristics, conductance rectification, and inverse membrane potential. However, the traditional PNP model is a continuum mean-field model and has no or underestimates the discrete ion effects, in particular the ion solvation or self-energy (which can be described by Born model). It is known that the dehydration effect (closely related to ion size) is crucial to selective permeation in potassium channels. Therefore, we incorporated Born solvation energy into the PNP model to account for ion hydration and dehydration effects when passing through inhomogeneous dielectric channel environments. A variational approach was adopted to derive a Born-energy-modified PNP (BPNP) model. The model was applied to study a cylindrical nanopore and a realistic KcsA channel, and three-dimensional finite element simulations were performed. The BPNP model can distinguish different ion species by ion radius and predict selectivity for K^{+} over Na^{+} in KcsA channels. Furthermore, ion current rectification in the KcsA channel was observed by both the PNP and BPNP models. The I-V curve of the BPNP model for the KcsA channel indicated an inward rectifier effect for K^{+} (rectification ratio of ∼3/2) but indicated an outward rectifier effect for Na^{+} (rectification ratio of ∼1/6).
Collapse
Affiliation(s)
- Xuejiao Liu
- State Key Laboratory of Scientific and Engineering Computing, National Center for Mathematics and Interdisciplinary Sciences, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 100190, China and School of Mathematical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Benzhuo Lu
- State Key Laboratory of Scientific and Engineering Computing, National Center for Mathematics and Interdisciplinary Sciences, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 100190, China and School of Mathematical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
5
|
Fahrenberger F, Holm C. Computing the Coulomb interaction in inhomogeneous dielectric media via a local electrostatics lattice algorithm. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 90:063304. [PMID: 25615219 DOI: 10.1103/physreve.90.063304] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/18/2013] [Indexed: 06/04/2023]
Abstract
The local approach to computing electrostatic interactions proposed by Maggs and adapted by Rottler and Pasichnyk for molecular-dynamics simulations is extended to situations where the dielectric background medium is inhomogeneous. We furthermore correct a problem of the original algorithm related to the correct treatment of the global dipole moment, provide an error estimate for the accuracy of the algorithm, and suggest a different form of the treatment of the self-energy problem. Our implementation is highly scalable on many cores, and we have validated and compared its performance against theoretical predictions and simulation data obtained by other algorithmic approaches.
Collapse
Affiliation(s)
- F Fahrenberger
- Institut für Computerphysik, Universität Stuttgart, Allmandring 3, 70569 Stuttgart, Germany
| | - C Holm
- Institut für Computerphysik, Universität Stuttgart, Allmandring 3, 70569 Stuttgart, Germany
| |
Collapse
|
6
|
Boda D. Monte Carlo Simulation of Electrolyte Solutions in Biology. ANNUAL REPORTS IN COMPUTATIONAL CHEMISTRY 2014. [DOI: 10.1016/b978-0-444-63378-1.00005-7] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
7
|
Bonthuis DJ, Netz RR. Beyond the Continuum: How Molecular Solvent Structure Affects Electrostatics and Hydrodynamics at Solid–Electrolyte Interfaces. J Phys Chem B 2013; 117:11397-413. [DOI: 10.1021/jp402482q] [Citation(s) in RCA: 94] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Douwe Jan Bonthuis
- Rudolf
Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX1 3NP, United Kingdom
| | - Roland R. Netz
- Fachbereich
Physik, Freie Universität Berlin, 14195 Berlin, Germany
| |
Collapse
|
8
|
Dyrka W, Bartuzel MM, Kotulska M. Optimization of 3D Poisson-Nernst-Planck model for fast evaluation of diverse protein channels. Proteins 2013; 81:1802-22. [PMID: 23720356 DOI: 10.1002/prot.24326] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2012] [Revised: 05/02/2013] [Accepted: 05/09/2013] [Indexed: 12/12/2022]
Abstract
We show the accuracy and applicability of our fast algorithmic implementation of a three-dimensional Poisson-Nernst-Planck (3D-PNP) flow model for characterizing different protein channels. Due to its high computational efficiency, our model can predict the full current-voltage characteristics of a channel within minutes, based on the experimental 3D structure of the channel or its computational model structure. Compared with other methods, such as Brownian dynamics, which currently needs a few weeks of the computational time, or even much more demanding molecular dynamics modeling, 3D-PNP is the only available method for a function-based evaluation of very numerous tentative structural channel models. Flow model tests of our algorithm and its optimal parametrization are provided for five native channels whose experimental structures are available in the protein data bank (PDB) in an open conductive state, and whose experimental current-voltage characteristics have been published. The channels represent very different geometric and structural properties, which makes it the widest test to date of the accuracy of 3D-PNP on real channels. We test whether the channel conductance, rectification, and charge selectivity obtained from the flow model, could be sufficiently sensitive to single-point mutations, related to unsignificant changes in the channel structure. Our results show that the classical 3D-PNP model, under proper parametrization, is able to achieve a qualitative agreement with experimental data for a majority of the tested characteristics and channels, including channels with narrow and irregular conductivity pores. We propose that although the standard PNP model cannot provide insight into complex physical phenomena due to its intrinsic limitations, its semiquantitative agreement is achievable for rectification and selectivity at a level sufficient for the bioinformatical purpose of selecting the best structural models with a great advantage of a very short computational time.
Collapse
Affiliation(s)
- Witold Dyrka
- Group of Bioinformatics and Biophysics of Nanopores, Institute of Biomedical Engineering and Instrumentation, Wroclaw University of Technology, Wroclaw, Poland
| | | | | |
Collapse
|
9
|
Boda D, Henderson D, Gillespie D. The role of solvation in the binding selectivity of the L-type calcium channel. J Chem Phys 2013; 139:055103. [DOI: 10.1063/1.4817205] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
|
10
|
Finnerty JJ, Eisenberg R, Carloni P. Localizing the Charged Side Chains of Ion Channels within the Crowded Charge Models. J Chem Theory Comput 2012; 9:766-73. [DOI: 10.1021/ct300768j] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Justin J. Finnerty
- Computational Biophysics, German
Research School for Simulation Sciences, 52425 Jülich, Germany
| | - Robert Eisenberg
- Department of Molecular Biophysics
and Physiology, Rush University, Chicago, Illinois 60612, United States
| | - Paolo Carloni
- Computational Biophysics, German
Research School for Simulation Sciences, 52425 Jülich, Germany
and Institute for Advanced Simulation IAS-5, Computational Biomedicine,
Forschungszentrum Jülich, 52425 Jülich, Germany
| |
Collapse
|
11
|
Horng TL, Lin TC, Liu C, Eisenberg B. PNP Equations with Steric Effects: A Model of Ion Flow through Channels. J Phys Chem B 2012; 116:11422-41. [DOI: 10.1021/jp305273n] [Citation(s) in RCA: 118] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Affiliation(s)
- Tzyy-Leng Horng
- Department of Applied Mathematics, Feng Chia University, 100 Wen-Hwa Road, Taichung, Taiwan
40724
| | - Tai-Chia Lin
- Department of Mathematics, Taida Institute for Mathematical
Sciences (TIMS), No. 1, Sec. 4, National Taiwan University, Roosevelt Road, Taipei 106, Taiwan
| | - Chun Liu
- Department of Mathematics, Pennsylvania State University University Park, Pennsylvania 16802,
United States
| | - Bob Eisenberg
- Department of Molecular Biophysics and Physiology, Rush University, 1653 West Congress Parkway, Chicago,
Illinois 60612, United States
| |
Collapse
|
12
|
Csányi E, Boda D, Gillespie D, Kristóf T. Current and selectivity in a model sodium channel under physiological conditions: Dynamic Monte Carlo simulations. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2011; 1818:592-600. [PMID: 22080102 DOI: 10.1016/j.bbamem.2011.10.029] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/02/2011] [Revised: 10/06/2011] [Accepted: 10/26/2011] [Indexed: 11/15/2022]
Abstract
A reduced model of a sodium channel is analyzed using Dynamic Monte Carlo simulations. These include the first simulations of ionic current under approximately physiological ionic conditions through a model sodium channel and an analysis of how mutations of the sodium channel's DEKA selectivity filter motif transform the channel from being Na(+) selective to being Ca(2+) selective. Even though the model of the pore, amino acids, and permeant ions is simplified, the model reproduces the fundamental properties of a sodium channel (e.g., 10 to 1 Na(+) over K(+) selectivity, Ca(2+) exclusion, and Ca(2+) selectivity after several point mutations). In this model pore, ions move through the pore one at a time by simple diffusion and Na(+) versus K(+) selectivity is due to both the larger K(+) not fitting well into the selectivity filter that contains amino acid terminal groups and K(+) moving more slowly (compared to Na(+)) when it is in the selectivity filter.
Collapse
Affiliation(s)
- Eva Csányi
- Department of Physical Chemistry, University of Pannonia, Veszprém, Hungary
| | | | | | | |
Collapse
|
13
|
Nagy T, Henderson D, Boda D. Simulation of an electrical double layer model with a low dielectric layer between the electrode and the electrolyte. J Phys Chem B 2011; 115:11409-19. [PMID: 21848262 DOI: 10.1021/jp2063244] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We report Monte Carlo simulation results for double layers of 1:1 and 2:1 electrolytes near an electrode with an inner layer that has a dielectric constant, ε(2), smaller than that of the electrolyte, ε(3). The electrolyte is modeled in the implicit solvent framework (primitive model), while the electrode is a metal electrode in this study (ε(1) → ∞). The charged hard sphere ions are not allowed to enter into the inner layer. We show that the capacitance of the inner layer is C(δ) = ε(0)(ε(2) + ε(3))/2δ, where δ is the thickness of the inner layer. This result is different from that obtained from solutions of the Poisson-Boltzmann equation (ε(0)ε(2)/δ), indicating that interpretation of experimental data with a fitted ε(2) dielectric constant of the inner layer must be done using a different equation. We also show that the properties of the diffuse layer are not independent of the value of ε(2), which is a usual assumption of the Poisson-Boltzmann theory. This is mainly because the repulsive image charges repel both the counterions and the co-ions, while the electrode charge attracts the counterions and repels the co-ions.
Collapse
Affiliation(s)
- Tímea Nagy
- Department of Physical Chemistry, University of Pannonia, P.O. Box 158, H-8201 Veszprém, Hungary
| | | | | |
Collapse
|