1
|
Zhao S, Ijaodoro I, McGowan M, Alexov E. Calculation of electrostatic free energy for the nonlinear Poisson-Boltzmann model based on the dimensionless potential. JOURNAL OF COMPUTATIONAL PHYSICS 2024; 497:112634. [PMID: 38045553 PMCID: PMC10688429 DOI: 10.1016/j.jcp.2023.112634] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/05/2023]
Abstract
The Poisson-Boltzmann (PB) equation governing the electrostatic potential with a unit is often transformed to a normalized form for a dimensionless potential in numerical studies. To calculate the electrostatic free energy (EFE) of biological interests, a unit conversion has to be conducted, because the existing PB energy functionals are all described in terms of the original potential. To bypass this conversion, this paper proposes energy functionals in terms of the dimensionless potential for the first time in the literature, so that the normalized PB equation can be directly derived by using the Euler-Lagrange variational analysis. Moreover, alternative energy forms have been rigorously derived to avoid approximating the gradient of singular functions in the electrostatic stress term. A systematic study has been carried out to examine the surface integrals involved in alternative energy forms and their dependence on finite domain size and mesh step size, which leads to a recommendation on the EFE forms for efficient computation of protein systems. The calculation of the EFE in the regularization formulation, which is an analytical approach for treating singular charge sources of the PB equation, has also been studied. The proposed energy forms have been validated by considering smooth dielectric settings, such as diffuse interface and super-Gaussian, for which the EFE of the nonlinear PB model is found to be significantly different from that of the linearized PB model. All proposed energy functionals and EFE forms are designed such that the dimensionless potential can be simply plugged in to compute the EFE in the unit of kcal/mol, and they can also be applied in the classical sharp interface PB model.
Collapse
Affiliation(s)
- Shan Zhao
- Department of Mathematics, University of Alabama, Tuscaloosa, AL 35487, USA
| | - Idowu Ijaodoro
- Department of Mathematics, University of Alabama, Tuscaloosa, AL 35487, USA
| | - Mark McGowan
- Department of Mathematics, University of Alabama, Tuscaloosa, AL 35487, USA
| | - Emil Alexov
- Department of Physics and Astronomy, Clemson University, Clemson, SC 29634, USA
| |
Collapse
|
2
|
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
|
3
|
Beckinghausen M, Spakowitz AJ. Interplay of Polymer Structure, Solvent Ordering, and Charge Fluctuations in Polyelectrolyte Solution Thermodynamics. Macromolecules 2022. [DOI: 10.1021/acs.macromol.2c01826] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Michael Beckinghausen
- Department of Chemical Engineering, Stanford University, Stanford, California94305, United States
| | - Andrew J. Spakowitz
- Department of Chemical Engineering, Stanford University, Stanford, California94305, United States
- Department of Materials Science and Engineering, Stanford University, Stanford, California94305, United States
- Department of Applied Physics, Stanford University, Stanford, California94305, United States
- Biophysics Program, Stanford University, Stanford, California94305, United States
| |
Collapse
|
4
|
Xie D, Audi SH, Dash RK. A size-modified poisson-boltzmann ion channel model in a solvent of multiple ionic species: Application to voltage-dependent anion channel. J Comput Chem 2020; 41:218-230. [PMID: 31845398 PMCID: PMC8189662 DOI: 10.1002/jcc.26091] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2019] [Revised: 08/30/2019] [Accepted: 09/24/2019] [Indexed: 12/30/2022]
Abstract
We present a new size-modified Poisson-Boltzmann ion channel (SMPBIC) model and use it to calculate the electrostatic potential, ionic concentrations, and electrostatic solvation free energy for a voltage-dependent anion channel (VDAC) on a biological membrane in a solution mixture of multiple ionic species. In particular, the new SMPBIC model adopts a membrane surface charge density and a natural Neumann boundary condition to reflect the charge effect of the membrane on the electrostatics of VDAC. To avoid the singularity difficulties caused by the atomic charges of VDAC, the new SMPBIC model is split into three submodels such that the solution of one of the submodels is obtained analytically and contains all the singularity points of the SMPBIC model. The other two submodels are then solved numerically much more efficiently than the original SMPBIC model. As an application of this SMPBIC submodel partitioning scheme, we derive a new formula for computing the electrostatic solvation free energy. Numerical results for a human VDAC isoform 1 (hVDAC1) in three different salt solutions, each with up to five different ionic species, confirm the significant effects of membrane surface charges on both the electrostatics and ionic concentrations. The results also show that the new SMPBIC model can describe well the anion selectivity property of hVDAC1, and that the new electrostatic solvation free energy formula can significantly improve the accuracy of the currently used formula. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Dexuan Xie
- Department of Mathematical Sciences, University of Wisconsin-Milwaukee, Milwaukee, Wisconsin, 53201
| | - Said H Audi
- Department of Biomedical Engineering, Marquette University, Milwaukee, Wisconsin, 53233
| | - Ranjan K Dash
- Department of Biomedical Engineering, Medical College of Wisconsin, Milwaukee, Wisconsin, 53226
| |
Collapse
|
5
|
Sugimoto Y, Kitazumi Y, Shirai O, Nishikawa K, Higuchi Y, Yamamoto M, Kano K. Electrostatic roles in electron transfer from [NiFe] hydrogenase to cytochrome c 3 from Desulfovibrio vulgaris Miyazaki F. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2017; 1865:481-487. [DOI: 10.1016/j.bbapap.2017.02.009] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/06/2017] [Revised: 02/08/2017] [Accepted: 02/08/2017] [Indexed: 10/20/2022]
|
6
|
Xie Y, Ying J, Xie D. SMPBS: Web server for computing biomolecular electrostatics using finite element solvers of size modified Poisson-Boltzmann equation. J Comput Chem 2017; 38:541-552. [DOI: 10.1002/jcc.24703] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2016] [Revised: 11/21/2016] [Accepted: 11/27/2016] [Indexed: 11/11/2022]
Affiliation(s)
- Yang Xie
- Department of Computer Science; University of Wisconsin-Milwaukee; Milwaukee Wisconsin 53201
| | - Jinyong Ying
- Department of Mathematical Sciences; University of Wisconsin-Milwaukee; Milwaukee Wisconsin 53201
| | - Dexuan Xie
- Department of Mathematical Sciences; University of Wisconsin-Milwaukee; Milwaukee Wisconsin 53201
| |
Collapse
|
7
|
Buyukdagli S, Blossey R. Beyond Poisson-Boltzmann: fluctuations and fluid structure in a self-consistent theory. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2016; 28:343001. [PMID: 27357125 DOI: 10.1088/0953-8984/28/34/343001] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Poisson-Boltzmann (PB) theory is the classic approach to soft matter electrostatics and has been applied to numerous physical chemistry and biophysics problems. Its essential limitations are in its neglect of correlation effects and fluid structure. Recently, several theoretical insights have allowed the formulation of approaches that go beyond PB theory in a systematic way. In this topical review, we provide an update on the developments achieved in the self-consistent formulations of correlation-corrected Poisson-Boltzmann theory. We introduce a corresponding system of coupled non-linear equations for both continuum electrostatics with a uniform dielectric constant, and a structured solvent-a dipolar Coulomb fluid-including non-local effects. While the approach is only approximate and also limited to corrections in the so-called weak fluctuation regime, it allows us to include physically relevant effects, as we show for a range of applications of these equations.
Collapse
Affiliation(s)
- S Buyukdagli
- Department of Physics, Bilkent University, Ankara 06800, Turkey
| | | |
Collapse
|
8
|
Poitevin F, Delarue M, Orland H. Beyond Poisson-Boltzmann: Numerical Sampling of Charge Density Fluctuations. J Phys Chem B 2016; 120:6270-7. [PMID: 27075231 DOI: 10.1021/acs.jpcb.6b02650] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present a method aimed at sampling charge density fluctuations in Coulomb systems. The derivation follows from a functional integral representation of the partition function in terms of charge density fluctuations. Starting from the mean-field solution given by the Poisson-Boltzmann equation, an original approach is proposed to numerically sample fluctuations around it, through the propagation of a Langevin-like stochastic partial differential equation (SPDE). The diffusion tensor of the SPDE can be chosen so as to avoid the numerical complexity linked to long-range Coulomb interactions, effectively rendering the theory completely local. A finite-volume implementation of the SPDE is described, and the approach is illustrated with preliminary results on the study of a system made of two like-charge ions immersed in a bath of counterions.
Collapse
Affiliation(s)
- Frédéric Poitevin
- Institut de Physique Théorique, Université Paris Saclay, CEA, UMR3681 du CNRS, F-91191 Gif-sur-Yvette, France
| | - Marc Delarue
- Unit of Structural Dynamics of Macromolecules, UMR 3528 du CNRS, Institut Pasteur, 75015 Paris, France
| | - Henri Orland
- Institut de Physique Théorique, Université Paris Saclay, CEA, UMR3681 du CNRS, F-91191 Gif-sur-Yvette, France.,Beijing Computational Science Research Center, Haidian District, Beijing 100094, China
| |
Collapse
|
9
|
Cumberworth A, Bui JM, Gsponer J. Free energies of solvation in the context of protein folding: Implications for implicit and explicit solvent models. J Comput Chem 2015; 37:629-40. [DOI: 10.1002/jcc.24235] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2015] [Revised: 09/25/2015] [Accepted: 10/06/2015] [Indexed: 12/13/2022]
Affiliation(s)
| | | | - Jörg Gsponer
- Center for High-Throughput Biology, UBC; Vancouver Canada
| |
Collapse
|
10
|
Alexov E, Luo R. Editorial. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2014. [DOI: 10.1142/s0219633614020015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Affiliation(s)
- Emil Alexov
- Computational Biophysics and Bioinformatics, Department of Physics, Clemson University, Clemson, SC 29634, USA
| | - Ray Luo
- Department of Biomedical Engineering, Department of Molecular Biology and Biochemistry, University of California, Irvine, CA 92697, USA
| |
Collapse
|