A new box iterative method for a class of nonlinear interface problems with application in solving Poisson-Boltzmann equation

In this paper, a new box iterative method for solving a class of nonlinear interface problems is proposed by intermixing linear and nonlinear boundary value problems based on a special seven-overlapped-boxes partition. It is then applied to the construction of a new finite element and finite difference hybrid scheme for solving the Poisson-Boltzmann equation (PBE) - a second order nonlinear elliptic interface problem for computing electrostatics of an ionic solvated protein. Furthermore, a modified Newton minimization algorithm accelerated by a multigrid preconditioned conjugate gradient method is presented to efficiently solve each involved nonlinear boundary value problem. In addition, the analytical solution of a Poisson dielectric test model with a spherical solute region containing multiple charges is expressed in a simple series of Legendre polynomials, resulting in a new PBE test model that works for a large number of point charges. The new PBE hybrid solver is programmed as a software package, and numerically validated on the new PBE test model with 892 point charges. It is also compared to a commonly used finite difference scheme in the accuracy of computing solution and electrostatic free energy for three proteins with up to 2124 atomic charges. Numerical results on six proteins demonstrate its high performance in comparison to the PBE finite element program package reported in Xie (2014). Present a new box iterative method for a class of nonlinear interface problems.Obtain a new hybrid solver of Poisson-Boltzmann Equation (PBE) as application.Develop a new Newton-PCG-MG scheme for nonlinear boundary value problems.Obtain a simple series solution for Poisson ball test model with multiple charges.Validate the new PBE hybrid software and demonstrate its high performance.

[1]  Ray Luo,et al.  Exploring accurate Poisson-Boltzmann methods for biomolecular simulations. , 2013, Computational & theoretical chemistry.

[2]  D. A. Dunnett.  Classical Electrodynamics , 2020, Nature.

[3]  Nathan A. Baker,et al.  Biomolecular electrostatics and solvation: a computational perspective , 2012, Quarterly Reviews of Biophysics.

[4]  Tao Lin,et al.  Partially Penalized Immersed Finite Element Methods For Elliptic Interface Problems , 2015, SIAM J. Numer. Anal..

[5]  Barry Honig,et al.  Focusing of electric fields in the active site of Cu‐Zn superoxide dismutase: Effects of ionic strength and amino‐acid modification , 1986, Proteins.

[6]  A. Brigo,et al.  The Poisson–Boltzmann equation for biomolecular electrostatics: a tool for structural biology , 2002, Journal of molecular recognition : JMR.

[7]  J. Jackson.  Classical Electrodynamics, 3rd Edition , 1998 .

[8]  Ronald Fedkiw,et al.  The immersed interface method. Numerical solutions of PDEs involving interfaces and irregular domains , 2007, Math. Comput..

[9]  Jinyong Ying,et al.  A new finite element and finite difference hybrid method for computing electrostatics of ionic solvated biomolecule , 2015, J. Comput. Phys..

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

[11]  J. Kirkwood,et al.  Theory of Solutions of Molecules Containing Widely Separated Charges with Special Application to Zwitterions , 1934 .

[12]  Guo-Wei Wei,et al.  MIB Galerkin method for elliptic interface problems , 2014, J. Comput. Appl. Math..

[13]  William J. Rider,et al.  A Multigrid Preconditioned Newton-Krylov Method , 1999, SIAM J. Sci. Comput..

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

[15]  Tarek P. Mathew,et al.  Domain Decomposition Methods for the Numerical Solution of Partial Differential Equations , 2008, Lecture Notes in Computational Science and Engineering.

[16]  Andrea Toselli,et al.  Domain decomposition methods : algorithms and theory , 2005 .

[17]  Zi-Cai Li,et al.  Schwarz Alternating Method , 1998 .

[18]  D. Brandt,et al.  Multi-level adaptive solutions to boundary-value problems math comptr , 1977 .

[19]  Maria Teresa Neves-Petersen,et al.  Protein electrostatics: a review of the equations and methods used to model electrostatic equations in biomolecules--applications in biotechnology. , 2003, Biotechnology annual review.

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

[21]  M. Born.  Volumen und Hydratationswärme der Ionen , 1920 .

[22]  Jinchao Xu,et al.  Some Nonoverlapping Domain Decomposition Methods , 1998, SIAM Rev..

[23]  Barry Smith,et al.  Domain Decomposition Methods for Partial Differential Equations , 1997 .

[24]  M. B. Pinto,et al.  Optimized δ expansion for relativistic nuclear models , 1997, nucl-th/9709049.

[25]  A. Toselli.  Domain Decomposition Methods , 2004 .

[26]  Jinyong Ying,et al.  A Poisson-Boltzmann Equation Test Model for Protein in Spherical Solute Region and its Applications , 2014 .

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

[28]  Barbara I. Wohlmuth,et al.  Mortar Finite Elements for Interface Problems , 2004, Computing.

[29]  Zhilin Li,et al.  The Immersed Interface Method: Numerical Solutions of PDEs Involving Interfaces and Irregular Domains (Frontiers in Applied Mathematics) , 2006 .

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

[31]  Shan Zhao,et al.  WEAK GALERKIN METHODS FOR SECOND ORDER ELLIPTIC INTERFACE PROBLEMS. , 2013, Journal of computational physics.

[32]  Sunhwan Jo,et al.  PBEQ-Solver for online visualization of electrostatic potential of biomolecules , 2008, Nucleic Acids Res..

[33]  Eftychios Sifakis,et al.  A second order virtual node method for elliptic problems with interfaces and irregular domains in three dimensions , 2012, J. Comput. Phys..

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

[35]  B. Honig,et al.  Classical electrostatics in biology and chemistry. , 1995, Science.

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

[37]  Guo-Wei Wei,et al.  A Galerkin formulation of the MIB method for three dimensional elliptic interface problems , 2014, Comput. Math. Appl..

[38]  Ray Luo,et al.  Recent progress in adapting Poisson–Boltzmann methods to molecular simulations , 2014 .

[39]  P. Deuflhard.  Newton Methods for Nonlinear Problems: Affine Invariance and Adaptive Algorithms , 2011 .

[40]  Eftychios Sifakis,et al.  ' s personal copy A second order virtual node method for elliptic problems with interfaces and irregular domains , 2010 .

[41]  Maxim A. Olshanskii,et al.  Iterative Methods for Linear Systems - Theory and Applications , 2014 .

[42]  H. Schwarz.  Ueber einige Abbildungsaufgaben. , 1869 .

[43]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[44]  Barry Honig,et al.  Extending the Applicability of the Nonlinear Poisson−Boltzmann Equation: Multiple Dielectric Constants and Multivalent Ions† , 2001 .

[45]  L. R. Scott,et al.  The Mathematical Theory of Finite Element Methods , 1994 .

[46]  Dexuan Xie,et al.  New solution decomposition and minimization schemes for Poisson-Boltzmann equation in calculation of biomolecular electrostatics , 2014, J. Comput. Phys..

[47]  Yunrong Zhu,et al.  ADAPTIVE FINITE ELEMENT MODELING TECHNIQUES FOR THE POISSON-BOLTZMANN EQUATION. , 2010, Communications in computational physics.

[48]  C. Peskin.  The immersed boundary method , 2002, Acta Numerica.

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