Adaptive multilevel finite element solution of the Poisson-Boltzmann equation II. Refinement at solvent-accessible surfaces in biomolecular systems

We apply the adaptive multilevel finite element techniques described in [20] to the nonlinear Poisson-Boltzmann equation (PBE) in the context of biomolecules. Fast and accurate numerical solution of the PBE in this setting is usually difficult to accomplish due to presence of discontinuous coefficients, delta functions, three spatial dimensions, unbounded domains, and rapid (exponential) nonlinearity. However, these adaptive techniques have shown substantial improvement in solution time over conventional uniform-mesh finite difference methods. One important aspect of the adaptive multilevel finite element method is the robust a posteriori error estimators necessary to drive the adaptive refinement routines. This paper discusses the choice of solvent accessibility for a posteriori error estimation of PBE solutions and the implementation of such routines in the “Adaptive Poisson-Boltzmann Solver” (APBS) software package based on the “Manifold Code” (MC) libraries. Results are shown for the application of this method to several biomolecular

[1]  R W Hockney,et al.  Computer Simulation Using Particles , 1966 .

[2]  B. Quentrec,et al.  New method for searching for neighbors in molecular dynamics computations , 1973 .

[3]  I. Babuska,et al.  A‐posteriori error estimates for the finite element method , 1978 .

[4]  W. Rheinboldt,et al.  Error Estimates for Adaptive Finite Element Computations , 1978 .

[5]  Randolph E. Bank,et al.  An optimal order process for solving finite element equations , 1981 .

[6]  M. Rivara.  Algorithms for refining triangular grids suitable for adaptive and multigrid techniques , 1984 .

[7]  Wolfgang Hackbusch,et al.  Multi-grid methods and applications , 1985, Springer series in computational mathematics.

[8]  K. Sharp,et al.  Calculating the electrostatic potential of molecules in solution: Method and error assessment , 1988 .

[9]  C. Brooks.  Computer simulation of liquids , 1989 .

[10]  B Jayaram,et al.  The electrostatic potential of B‐DNA , 1989, Biopolymers.

[11]  Andrzej Wierzbicki,et al.  Solutions of the full Poisson-Boltzmann equation with application to diffusion-controlled reactions , 1989 .

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

[13]  B. J. Yoon,et al.  A boundary element method for molecular electrostatics with electrolyte effects , 1990 .

[14]  K. Sharp,et al.  Electrostatic interactions in macromolecules: theory and applications. , 1990, Annual review of biophysics and biophysical chemistry.

[15]  Barry Honig,et al.  Calculating total electrostatic energies with the nonlinear Poisson-Boltzmann equation , 1990 .

[16]  Alexander A. Rashin,et al.  New method for the computation of ionic distribution around rod‐like polyelectrolytes with the helical distribution of charges. I. General approach and a nonlinearized Poisson–Boltzmann equation , 1991 .

[17]  M. Rivara,et al.  Local modification of meshes for adaptive and/or multigrid finite-element methods , 1991 .

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

[19]  H. Berendsen,et al.  The electric potential of a macromolecule in a solvent: A fundamental approach , 1991 .

[20]  J. Andrew McCammon,et al.  Solving the finite‐difference non‐linear Poisson–Boltzmann equation , 1992 .

[21]  James M. Briggs,et al.  Computation Unravels Mysteries of Molecular Biophysics , 1992 .

[22]  Jinchao Xu,et al.  Iterative Methods by Space Decomposition and Subspace Correction , 1992, SIAM Rev..

[23]  Klaus Schulten,et al.  Molecular Dynamics Simulations in Heterogeneous Dielectrica and Debye-Hückel Media - Application to the Protein Bovine Pancreatic Trypsin Inhibitor , 1992 .

[24]  R. Verfürth.  A posteriori error estimates for nonlinear problems: finite element discretizations of elliptic equations , 1994 .

[25]  Michael J. Holst,et al.  Numerical solution of the nonlinear Poisson–Boltzmann equation: Developing more robust and efficient methods , 1995, J. Comput. Chem..

[26]  Tony You,et al.  An analytical algorithm for the rapid determination of the solvent accessibility of points in a three‐dimensional lattice around a solute molecule , 1995, J. Comput. Chem..

[27]  Rüdiger Verführt,et al.  A review of a posteriori error estimation and adaptive mesh-refinement techniques , 1996, Advances in numerical mathematics.

[28]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[29]  Bowen,et al.  Adaptive Finite-Element Solution of the Nonlinear Poisson-Boltzmann Equation: A Charged Spherical Particle at Various Distances from a Charged Cylindrical Pore in a Charged Planar Surface , 1997, Journal of colloid and interface science.

[30]  Richard A. Friesner,et al.  An automatic three-dimensional finite element mesh generation system for the Poisson-Boltzmann equation , 1997, J. Comput. Chem..

[31]  P Taylor,et al.  Inhibition of mouse acetylcholinesterase by fasciculin: crystal structure of the complex and mutagenesis of fasciculin. , 1998, Toxicon : official journal of the International Society on Toxinology.

[32]  Randolph E. Bank,et al.  PLTMG - a software package for solving elliptic partial differential equations: users' guide 8.0 , 1998, Software, environments, tools.

[33]  Michael J. Holst,et al.  Adaptive multilevel finite element solution of the Poisson–Boltzmann equation I. Algorithms and examples , 2001 .

[34]  M. Holst,et al.  Adaptive multilevel finite element solution of the Poisson–Boltzmann equation I. Algorithms and examples , 2001 .