Modeling Microtubule Counterion Distributions and Conductivity Using the Poisson-Boltzmann Equation.
Ontology highlight
ABSTRACT:
SUBMITTER: Eakins BB
PROVIDER: S-EPMC8027483 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:We consider the design of an effective and reliable adaptive finite element method (AFEM) for the nonlinear Poisson-Boltzmann equation (PBE). We first examine the two-term regularization technique for the continuous problem recently proposed by Chen, Holst, and Xu based on the removal of the singular electrostatic potential inside biomolecules; this technique made possible the development of the first complete solution and approximation theory for the Poisson-Boltzmann equation, the first provably convergent discretization, and also allowed for the development of a provably convergent AFEM. However, in practical implementation, this two-term regularization exhibits numerical instability. Therefore, we examine a variation of this regularization technique which can be shown to be less susceptible to such instability. We establish a priori estimates and other basic results for the continuous regularized problem, as well as for Galerkin finite element approximations. We show that the new approach produces regularized continuous and discrete problems with the same mathematical advantages of the original regularization. We then design an AFEM scheme for the new regularized problem, and show that the resulting AFEM scheme is accurate and reliable, by proving a contraction result for the error. This result, which is one of the first results of this type for nonlinear elliptic problems, is based on using continuous and discrete a priori L(∞) estimates to establish quasi-orthogonality. To provide a high-quality geometric model as input to the AFEM algorithm, we also describe a class of feature-preserving adaptive mesh generation algorithms designed specifically for constructing meshes of biomolecular structures, based on the intrinsic local structure tensor of the molecular surface. All of the algorithms described in the article are implemented in the Finite Element Toolkit (FETK), developed and maintained at UCSD. The stability advantages of the new regularization scheme are demonstrated with FETK through comparisons with the original regularization approach for a model problem. The convergence and accuracy of the overall AFEM algorithm is also illustrated by numerical approximation of electrostatic solvation energy for an insulin protein.
| S-EPMC3178275 | biostudies-literature
Project description:Macromolecular interactions are essential for understanding numerous biological processes and are typically characterized by the binding free energy. Important component of the binding free energy is the electrostatics, which is frequently modeled via the solutions of the Poisson-Boltzmann Equations (PBE). However, numerous works have shown that the electrostatic component (ΔΔGelec ) of binding free energy is very sensitive to the parameters used and modeling protocol. This prompted some researchers to question the robustness of PBE in predicting ΔΔGelec . We argue that the sensitivity of the absolute ΔΔGelec calculated with PBE using different input parameters and definitions does not indicate PBE deficiency, rather this is what should be expected. We show how the apparent sensitivity should be interpreted in terms of the underlying changes in several numerous and physical parameters. We demonstrate that PBE approach is robust within each considered force field (CHARMM-27, AMBER-94, and OPLS-AA) once the corresponding structures are energy minimized. This observation holds despite of using two different molecular surface definitions, pointing again that PBE delivers consistent results within particular force field. The fact that PBE delivered ΔΔGelec values may differ if calculated with different modeling protocols is not a deficiency of PBE, but natural results of the differences of the force field parameters and potential functions for energy minimization. In addition, while the absolute ΔΔGelec values calculated with different force field differ, their ordering remains practically the same allowing for consistent ranking despite of the force field used. © 2016 Wiley Periodicals, Inc.
| S-EPMC5030180 | biostudies-literature
Project description:Experimental results have demonstrated that the numbers of counterions surrounding nucleic acids differ from those predicted by the nonlinear Poisson-Boltzmann equation, NLPBE. Some studies have fit these data against the ion size in the size-modified Poisson-Boltzmann equation, SMPBE, but the present study demonstrates that other parameters, such as the Stern layer thickness and the molecular surface definition, can change the number of bound ions by amounts comparable to varying the ion size. These parameters will therefore have to be fit simultaneously against experimental data. In addition, the data presented here demonstrate that the derivative, SK, of the electrostatic binding free energy, ΔGel, with respect to the logarithm of the salt concentration is sensitive to these parameters, and experimental measurements of SK could be used to parameterize the model. However, although better values for the Stern layer thickness and ion size and better molecular surface definitions could improve the model's predictions of the numbers of ions around biomolecules and SK, ΔGel itself is more sensitive to parameters, such as the interior dielectric constant, which in turn do not significantly affect the distributions of ions around biomolecules. Therefore, improved estimates of the ion size and Stern layer thickness to use in the SMPBE will not necessarily improve the model's predictions of ΔGel.
| S-EPMC3977794 | biostudies-literature
Project description:Grid-based solvers of the Poisson-Boltzmann, PB, equation are routinely used to estimate electrostatic binding, ??Gel, and solvation, ?Gel, free energies. The accuracies of such estimates are subject to grid discretization errors from the finite difference approximation to the PB equation. Here, we show that the grid discretization errors in ??Gel are more significant than those in ?Gel, and can be divided into two parts: (i) errors associated with the relative positioning of the grid and (ii) systematic errors associated with grid spacing. The systematic error in particular is significant for methods, such as the molecular mechanics PB surface area, MM-PBSA, approach that predict electrostatic binding free energies by averaging over an ensemble of molecular conformations. Although averaging over multiple conformations can control for the error associated with grid placement, it will not eliminate the systematic error, which can only be controlled by reducing grid spacing. The present study indicates that the widely-used grid spacing of 0.5 Å produces unacceptable errors in ??Gel, even though its predictions of ?Gel are adequate for the cases considered here. Although both grid discretization errors generally increase with grid spacing, the relative sizes of these errors differ according to the solute-solvent dielectric boundary definition. The grid discretization errors are generally smaller on the Gaussian surface used in the present study than on either the solvent-excluded or van der Waals surfaces, which both contain more surface discontinuities (e.g., sharp edges and cusps). Additionally, all three molecular surfaces converge to very different estimates of ??Gel.
| S-EPMC3756143 | biostudies-literature
Project description:Accurate calculation of electrostatic potential and gradient on the molecular surface is highly desirable for the continuum and hybrid modeling of large scale deformation of biomolecules in solvent. In this article a new numerical method is proposed to calculate these quantities on the dielectric interface from the numerical solutions of the Poisson-Boltzmann equation. Our method reconstructs a potential field locally in the least square sense on the polynomial basis enriched with Green's functions, the latter characterize the Coulomb potential induced by charges near the position of reconstruction. This enrichment resembles the decomposition of electrostatic potential into singular Coulomb component and the regular reaction field in the Generalized Born methods. Numerical experiments demonstrate that the enrichment recovery produces drastically more accurate and stable potential gradients on molecular surfaces compared to classical recovery techniques.
| S-EPMC7461612 | biostudies-literature
Project description:The Poisson-Boltzmann (PB) formalism is among the most popular approaches to modeling the solvation of molecules. It assumes a continuum model for water, leading to a dielectric permittivity that only depends on position in space. In contrast, the dipolar Poisson-Boltzmann-Langevin (DPBL) formalism represents the solvent as a collection of orientable dipoles with nonuniform concentration; this leads to a nonlinear permittivity function that depends both on the position and on the local electric field at that position. The differences in the assumptions underlying these two models lead to significant differences in the equations they generate. The PB equation is a second order, elliptic, nonlinear partial differential equation (PDE). Its response coefficients correspond to the dielectric permittivity and are therefore constant within each subdomain of the system considered (i.e., inside and outside of the molecules considered). While the DPBL equation is also a second order, elliptic, nonlinear PDE, its response coefficients are nonlinear functions of the electrostatic potential. Many solvers have been developed for the PB equation; to our knowledge, none of these can be directly applied to the DPBL equation. The methods they use may adapt to the difference; their implementations however are PBE specific. We adapted the PBE solver originally developed by Holst and Saied [J. Comput. Chem. 16, 337 (1995)] to the problem of solving the DPBL equation. This solver uses a truncated Newton method with a multigrid preconditioner. Numerical evidences suggest that it converges for the DPBL equation and that the convergence is superlinear. It is found however to be slow and greedy in memory requirement for problems commonly encountered in computational biology and computational chemistry. To circumvent these problems, we propose two variants, a quasi-Newton solver based on a simplified, inexact Jacobian and an iterative self-consistent solver that is based directly on the PBE solver. While both methods are not guaranteed to converge, numerical evidences suggest that they do and that their convergence is also superlinear. Both variants are significantly faster than the solver based on the exact Jacobian, with a much smaller memory footprint. All three methods have been implemented in a new code named AQUASOL, which is freely available.
| S-EPMC2833186 | biostudies-literature
Project description:The ion size-modified Poisson Boltzmann equation (SMPBE) is applied to the simple model problem of a low-dielectric spherical cavity containing a central charge, in an aqueous salt solution to investigate the finite ion size effect upon the electrostatic free energy and its sensitivity to changes in salt concentration. The SMPBE is shown to predict a very different electrostatic free energy than the nonlinear Poisson-Boltzmann equation (NLPBE) due to the additional entropic cost of placing ions in solution. Although the energy predictions of the SMPBE can be reproduced by fitting an appropriatelysized Stern layer, or ion-exclusion layer to the NLPBE calculations, the size of the Stern layer is difficult to estimate a priori. The SMPBE also produces a saturation layer when the central charge becomes sufficiently large. Ion-competition effects on various integrated quantities such the total number of ions predicted by the SMPBE are qualitatively similar to those given by the NLPBE and those found in available experimental results.
| S-EPMC3378335 | biostudies-literature
Project description:Although knowledge of the pKa values and charge states of individual residues is critical to understanding the role of electrostatic effects in protein structure and function, calculating these quantities is challenging because of the sensitivity of these parameters to the position and distribution of charges. Values for many different proteins which agree well with experimental results have been obtained with modified Tanford-Kirkwood theory in which the protein is modeled as a sphere (reviewed in Ref. 1); however, convergence is more difficult to achieve with finite difference methods, in which the protein is mapped onto a grid and derivatives of the potential function are calculated as differences between the values of the function at grid points (reviewed in Ref. 6). Multigrid methods, in which the size of the grid is varied from fine to coarse in several cycles, decrease computational time, increase rates of convergence, and improve agreement with experiment. Both the accuracy and computational advantage of the multigrid approach increase with grid size, because the time required to achieve a solution increases slowly with grid size. We have implemented a multigrid procedure for solving the nonlinear Poisson-Boltzmann equation, and, using lysozyme as a test case, compared calculations for several crystal forms, different refinement procedures, and different charge assignment schemes. The root mean square difference between calculated and experimental pKa values for the crystal structure which yields best agreement with experiment (1LZT) is 1.1 pH units, with the differences in calculated and experimental pK values being less than 0.6 pH units for 16 out of 21 residues. The calculated titration curves of several residues are biphasic.
| S-EPMC1225699 | biostudies-literature
Project description:The research of adsorption theory has recently gained renewed attention due to its critical relevance to a number of trending industrial applications, hydrogen storage and shale gas exploration for instance. The existing theoretical foundation, laid mostly in the early twentieth century, was largely based on simple heuristic molecular interaction models and static interaction potential which, although being insightful in illuminating the fundamental mechanisms, are insufficient for computations with realistic adsorbent structure and adsorbate hydrodynamics, both critical for real-life applications. Here we present and validate a novel lattice Boltzmann model incorporating both adsorbate-adsorbate and adsorbate-adsorbent interactions with hydrodynamics which, for the first time, allows adsorption to be computed with real-life details. Connection with the classic Ono-Kondo lattice theory is established and various adsorption isotherms, both within and beyond the IUPAC classification are observed as a pseudo-potential is varied. This new approach not only enables an important physical to be simulated for real-life applications, but also provides an enabling theoretical framework within which the fundamentals of adsorption can be studied.
| S-EPMC4891696 | biostudies-other
Project description:We present a boundary-element method (BEM) implementation for accurately solving problems in biomolecular electrostatics using the linearized Poisson-Boltzmann equation. Motivating this implementation is the desire to create a solver capable of precisely describing the geometries and topologies prevalent in continuum models of biological molecules. This implementation is enabled by the synthesis of four technologies developed or implemented specifically for this work. First, molecular and accessible surfaces used to describe dielectric and ion-exclusion boundaries were discretized with curved boundary elements that faithfully reproduce molecular geometries. Second, we avoided explicitly forming the dense BEM matrices and instead solved the linear systems with a preconditioned iterative method (GMRES), using a matrix compression algorithm (FFTSVD) to accelerate matrix-vector multiplication. Third, robust numerical integration methods were employed to accurately evaluate singular and near-singular integrals over the curved boundary elements. Fourth, we present a general boundary-integral approach capable of modeling an arbitrary number of embedded homogeneous dielectric regions with differing dielectric constants, possible salt treatment, and point charges. A comparison of the presented BEM implementation and standard finite-difference techniques demonstrates that for certain classes of electrostatic calculations, such as determining absolute electrostatic solvation and rigid-binding free energies, the improved convergence properties of the BEM approach can have a significant impact on computed energetics. We also demonstrate that the improved accuracy offered by the curved-element BEM is important when more sophisticated techniques, such as nonrigid-binding models, are used to compute the relative electrostatic effects of molecular modifications. In addition, we show that electrostatic calculations requiring multiple solves using the same molecular geometry, such as charge optimization or component analysis, can be computed to high accuracy using the presented BEM approach, in compute times comparable to traditional finite-difference methods.
| S-EPMC3465726 | biostudies-literature