ADAPTIVE FINITE ELEMENT MODELING TECHNIQUES FOR THE POISSON-BOLTZMANN EQUATION.

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.

[1]  Christian Kreuzer,et al.  Quasi-Optimal Convergence Rate for an Adaptive Finite Element Method , 2008, SIAM J. Numer. Anal..

[2]  James H. Bramble,et al.  A finite element method for interface problems in domains with smooth boundaries and interfaces , 1996, Adv. Comput. Math..

[3]  Michael J. Holst,et al.  Optimality of Multilevel Preconditioners for Local Mesh Refinement in Three Dimensions , 2006, SIAM J. Numer. Anal..

[4]  B. Lee,et al.  The interpretation of protein structures: estimation of static accessibility. , 1971, Journal of molecular biology.

[5]  Shan Zhao,et al.  High order matched interface and boundary method for elliptic equations with discontinuous coefficients and singular sources , 2006, J. Comput. Phys..

[6]  Michael J. Holst,et al.  Local multilevel preconditioners for elliptic equations with jump coefficients on bisection grids , 2010, Comput. Vis. Sci..

[7]  Michael J. Holst,et al.  The Finite Element Approximation of the Nonlinear Poisson-Boltzmann Equation , 2007, SIAM J. Numer. Anal..

[8]  James F. Blinn,et al.  A Generalization of Algebraic Surface Drawing , 1982, TOGS.

[9]  Jinchao Xu Two-grid Discretization Techniques for Linear and Nonlinear PDEs , 1996 .

[10]  E. M.,et al.  Statistical Mechanics , 2021, Manual for Theoretical Chemistry.

[11]  Benzhuo Lu,et al.  Order N algorithm for computation of electrostatic interactions in biomolecular systems , 2006, Proceedings of the National Academy of Sciences.

[12]  A. Olson,et al.  Shape analysis of molecular surfaces , 1993, Biopolymers.

[13]  Ricardo H. Nochetto,et al.  Convergence of Adaptive Finite Element Methods , 2002, SIAM Rev..

[14]  Charles Tanford,et al.  Physical Chemistry of Macromolecules , 1961 .

[15]  I-Liang Chern,et al.  Accurate Evaluation of Electrostatics for Macromolecules in Solution , 2003 .

[16]  Benzhuo Lu,et al.  Computation of electrostatic forces between solvated molecules determined by the Poisson-Boltzmann equation using a boundary element method. , 2005, The Journal of chemical physics.

[17]  Michael Levitt,et al.  Finite‐difference solution of the Poisson–Boltzmann equation: Complete elimination of self‐energy , 1996, J. Comput. Chem..

[18]  Burak Aksoylu,et al.  Local Refinement and Multilevel Preconditioning: Implementation and Numerical Experiments , 2010 .

[19]  Klaus Gärtner,et al.  Meshing Piecewise Linear Complexes by Constrained Delaunay Tetrahedralizations , 2005, IMR.

[20]  Benzhuo Lu,et al.  Calculation of the Maxwell stress tensor and the Poisson-Boltzmann force on a solvated molecular surface using hypersingular boundary integrals. , 2005, The Journal of chemical physics.

[21]  Chandrajit L. Bajaj,et al.  Quality meshing of implicit solvation models of biomolecular structures , 2006, Comput. Aided Geom. Des..

[22]  Zhilin Li AN OVERVIEW OF THE IMMERSED INTERFACE METHOD AND ITS APPLICATIONS , 2003 .

[23]  Giuseppe Savaré,et al.  Regularity Results for Elliptic Equations in Lipschitz Domains , 1998 .

[24]  Jitendra Malik,et al.  Scale-Space and Edge Detection Using Anisotropic Diffusion , 1990, IEEE Trans. Pattern Anal. Mach. Intell..

[25]  Kun Huang,et al.  Enhancing the activity of insulin at the receptor interface: crystal structure and photo-cross-linking of A8 analogues. , 2004, Biochemistry.

[26]  Michael Oevermann,et al.  A Cartesian grid finite volume method for elliptic equations with variable coefficients and embedded interfaces , 2006, J. Comput. Phys..

[27]  Michael Holst,et al.  Local Convergence of Adaptive Methods for Nonlinear Partial Differential Equations , 2010, 1001.1382.

[28]  Rüdiger Verfürth,et al.  A posteriori error estimates for nonlinear problems , 1994 .

[29]  J Andrew McCammon,et al.  Feature-preserving adaptive mesh generation for molecular shape modeling and simulation. , 2008, Journal of molecular graphics & modelling.

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

[31]  Michael J. Holst,et al.  Adaptive Numerical Treatment of Elliptic Systems on Manifolds , 2001, Adv. Comput. Math..

[32]  Weihua Geng,et al.  Treatment of charge singularities in implicit solvent models. , 2007, The Journal of chemical physics.

[33]  Ricardo H. Nochetto,et al.  Convergence of Adaptive Finite Element Methods for General Second Order Linear Elliptic PDEs , 2005, SIAM J. Numer. Anal..

[34]  Kazufumi Ito,et al.  Maximum Principle Preserving Schemes for Interface Problems with Discontinuous Coefficients , 2001, SIAM J. Sci. Comput..

[35]  Guo-Wei Wei,et al.  Highly accurate biomolecular electrostatics in continuum dielectric environments , 2008, J. Comput. Chem..

[36]  William E. Lorensen,et al.  Marching cubes: A high resolution 3D surface construction algorithm , 1987, SIGGRAPH.

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

[38]  Joachim Weickert,et al.  Anisotropic diffusion in image processing , 1996 .

[39]  Michael Holst,et al.  The Poisson-Boltzmann equation: Analysis and multilevel numerical solution , 1994 .

[40]  Ivo Babuska,et al.  The finite element method for elliptic equations with discontinuous coefficients , 1970, Computing.

[41]  J. Andrew McCammon,et al.  Computation of electrostatic forces on solvated molecules using the Poisson-Boltzmann equation , 1993 .

[42]  Michael Holst,et al.  Green's Functions and Boundary Value Problems: Stakgold/Green's Functions , 2011 .

[43]  J. Zou,et al.  Finite element methods and their convergence for elliptic and parabolic interface problems , 1998 .

[44]  M. L. Connolly Analytical molecular surface calculation , 1983 .

[45]  Michael J. Holst,et al.  An Odyssey into Local Refinement and Multilevel Preconditioning III: Implementation and Numerical Experiments , 2003, SIAM J. Sci. Comput..

[46]  Marcia O. Fenley,et al.  Hybrid boundary element and finite difference method for solving the nonlinear Poisson–Boltzmann equation , 2004, J. Comput. Chem..

[47]  Khamron Mekchay,et al.  Convergence of Adaptive Finite Element Methods , 2005 .

[48]  J. A. McCammon,et al.  R EVIEW A RTICLE Recent Progress inNumericalMethods for the Poisson-Boltzmann Equation in Biophysical Applications , 2008 .

[49]  José Jesús Fernández,et al.  An improved algorithm for anisotropic nonlinear diffusion for denoising cryo-tomograms. , 2003, Journal of structural biology.

[50]  Zeyun Yu,et al.  A segmentation-free approach for skeletonization of gray-scale images via anisotropic vector diffusion , 2004, CVPR 2004.

[51]  Zeyun Yu,et al.  High-Fidelity Geometric Modelling for Biomedical Applications , 2008 .

[52]  F M Richards,et al.  Areas, volumes, packing and protein structure. , 1977, Annual review of biophysics and bioengineering.

[53]  M. Sanner,et al.  Reduced surface: an efficient way to compute molecular surfaces. , 1996, Biopolymers.

[54]  G. Burton Sobolev Spaces , 2013 .

[55]  Rob P. Stevenson,et al.  An Optimal Adaptive Finite Element Method , 2004, SIAM J. Numer. Anal..

[56]  Marcia O. Fenley,et al.  A new outer boundary formulation and energy corrections for the nonlinear Poisson–Boltzmann equation , 2007, J. Comput. Chem..

[57]  W. Dörfler A convergent adaptive algorithm for Poisson's equation , 1996 .

[58]  Wei-Cheng Wang A Jump Condition Capturing Finite Difference Scheme for Elliptic Interface Problems , 2004, SIAM J. Sci. Comput..

[59]  Carl Ollivier-Gooch,et al.  Tetrahedral mesh improvement using swapping and smoothing , 1997 .

[60]  Kenji Shimada,et al.  An Angle-Based Approach to Two-Dimensional Mesh Smoothing , 2000, IMR.

[61]  Zeyun Yu,et al.  Computational Approaches for Automatic Structural Analysis of Large Biomolecular Complexes , 2008, IEEE/ACM Transactions on Computational Biology and Bioinformatics.

[62]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[63]  J. A. Grant,et al.  A Gaussian Description of Molecular Shape , 1995 .

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

[65]  Kuo-Young Cheng,et al.  A sharpness dependent filter for mesh smoothing , 2005, Comput. Aided Geom. Des..

[66]  R Abagyan,et al.  The contour-buildup algorithm to calculate the analytical molecular surface. , 1996, Journal of structural biology.

[67]  Timothy S. Newman,et al.  2D FE Quad Mesh Smoothing via Angle-Based Optimization , 2005, International Conference on Computational Science.

[68]  Dexuan Xie,et al.  A new minimization protocol for solving nonlinear Poisson–Boltzmann mortar finite element equation , 2007 .

[69]  Nathan A. Baker,et al.  PDB2PQR: an automated pipeline for the setup of Poisson-Boltzmann electrostatics calculations , 2004, Nucleic Acids Res..