Grid positioning independence and the reduction of self‐energy in the solution of the Poisson—Boltzmann equation

A common problem in the solution of the Poisson–Boltzmann equation using finite difference methods is the self‐energy of the system, also known as the grid energy. Because atoms are typically modeled as a point charge, the infinite self‐energy of a point charge is likewise modeled. In this article, a simple, alternate treatment of atomic charge is described where each atom is represented as a sphere of uniform charge. Unlike the point charge model, this method converges as the grid spacing is reduced. The uniform charge model generates the same electrostatic field outside the atoms. In addition, the use of fine grids reduces the variations in the potential due to variations in the position of atoms relative to the grid. Calculations of Born ion solvation energies, small‐molecule solvation energies, and the electrostatic field of superoxide dismutase are used to demonstrate that this method yields the same results as the point charge model. © John Wiley & Sons, Inc.

[1]  W. C. Still,et al.  Semianalytical treatment of solvation for molecular mechanics and dynamics , 1990 .

[2]  Malcolm E. Davis,et al.  Electrostatics in biomolecular structure and dynamics , 1990 .

[3]  Ralph G. Pearson,et al.  Ionization potentials and electron affinities in aqueous solution , 1986 .

[4]  F. Momany,et al.  Ab initio studies of molecular geometries. 27. Optimized molecular structures and conformational analysis of N.alpha.-acetyl-N-methylalaninamide and comparison with peptide crystal data and empirical calculations , 1983 .

[5]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[6]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[7]  J. Warwicker,et al.  Calculation of the electric potential in the active site cleft due to alpha-helix dipoles. , 1982, Journal of molecular biology.

[8]  Anna Tempczyk,et al.  Electrostatic contributions to solvation energies: comparison of free energy perturbation and continuum calculations , 1991 .

[9]  B. Honig,et al.  A rapid finite difference algorithm, utilizing successive over‐relaxation to solve the Poisson–Boltzmann equation , 1991 .

[10]  J. Richardson,et al.  Determination and analysis of the 2 A-structure of copper, zinc superoxide dismutase. , 1980, Journal of molecular biology.

[11]  Vincenzo Mollica,et al.  Group contributions to the thermodynamic properties of non-ionic organic solutes in dilute aqueous solution , 1981 .

[12]  M J Sternberg,et al.  Electrostatic interactions in globular proteins. Different dielectric models applied to the packing of alpha-helices. , 1984, Journal of molecular biology.

[13]  J. A. McCammon,et al.  Solving the finite difference linearized Poisson‐Boltzmann equation: A comparison of relaxation and conjugate gradient methods , 1989 .

[14]  Robert E. Bruccoleri,et al.  Application of Systematic Conformational Search to Protein Modeling , 1993 .