MIBPB: A software package for electrostatic analysis
Tóm tắt
The Poisson–Boltzmann equation (PBE) is an established model for the electrostatic analysis of biomolecules. The development of advanced computational techniques for the solution of the PBE has been an important topic in the past two decades. This article presents a matched interface and boundary (MIB)‐based PBE software package, the MIBPB solver, for electrostatic analysis. The MIBPB has a unique feature that it is the first interface technique‐based PBE solver that rigorously enforces the solution and flux continuity conditions at the dielectric interface between the biomolecule and the solvent. For protein molecular surfaces, which may possess troublesome geometrical singularities, the MIB scheme makes the MIBPB by far the only existing PBE solver that is able to deliver the second‐order convergence, that is, the accuracy increases four times when the mesh size is halved. The MIBPB method is also equipped with a Dirichlet‐to‐Neumann mapping technique that builds a Green's function approach to analytically resolve the singular charge distribution in biomolecules in order to obtain reliable solutions at meshes as coarse as 1 Å — whereas it usually takes other traditional PB solvers 0.25 Å to reach similar level of reliability. This work further accelerates the rate of convergence of linear equation systems resulting from the MIBPB by using the Krylov subspace (KS) techniques. Condition numbers of the MIBPB matrices are significantly reduced by using appropriate KS solver and preconditioner combinations. Both linear and nonlinear PBE solvers in the MIBPB package are tested by protein–solvent solvation energy calculations and analysis of salt effects on protein–protein binding energies, respectively. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011
Từ khóa
Tài liệu tham khảo
Baker N. A., 2006, Implicit solvent electrostatics in biomolecular simulation, 10.1007/3-540-31618-3_15
Lamm G., 2003, The Poisson–Boltzmann Equation, 147
Baker N. A., 2003, Electrostatic interactions, 427
Kirkwood J. G., 1934, J Comput Phys, 7, 351
Jo S, 2008, Nucl Acids Res, 270, 36
Yu S. N.Matched interface and boundary (MIB) method for geometric singularities and its application to molecular biology and structural analysis. Dissertation of Michigan State University: Michigan State University 2007.
Holst M. J., 1993, Multilevel Methods for the Poisson‐Boltzmann Equation. University of Illinois
Ortega J. M., 1999, Numerical Analysis: A Second Course