Hostname: page-component-76fb5796d-2lccl Total loading time: 0 Render date: 2024-04-25T22:55:50.029Z Has data issue: false hasContentIssue false
  • English
  • Français

Mathematical and Numerical Aspects of the Adaptive Fast Multipole Poisson-Boltzmann Solver

Published online by Cambridge University Press:  03 June 2015

Bo Zhang ,
Benzhuo Lu ,
Xiaolin Cheng ,
Jingfang Huang ,
Nikos P. Pitsianis ,
Xiaobai Sun  and
J. Andrew McCammon
Bo Zhang*
Affiliation:
Department of Computer Science, Duke University, NC 27708, USA
Benzhuo Lu*
Affiliation:
Institute of Computational Mathematics and Scientific/Engineering Computing, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 100910, China
Xiaolin Cheng*
Affiliation:
Center for Molecular Biophysics, Oak Ridge National Laboratory, TN 37831, USA
Jingfang Huang*
Affiliation:
Department of Mathematics, University of North Carolina, Chapel Hill, NC 27599, USA
Nikos P. Pitsianis*
Affiliation:
Department of Computer Science, Duke University, NC 27708, USA Department of Electrical and Computer Engineering, Aristotle University, Thessaloniki, 54124, Greece
Xiaobai Sun*
Affiliation:
Department of Computer Science, Duke University, NC 27708, USA
J. Andrew McCammon*
Affiliation:
Department of Chemistry & Biochemistry, Center for Theoretical Biological Physics, Department of Pharmacology, Howard Hughes Medical Institute, University of California, San Diego, CA 92093, USA
*
Email:zhangb@cs.duke.edu
Email:bzlu@lsec.cc.ac.cn
Email:chengx@ornl.gov
Corresponding author.Email:huang@email.unc.edu
Email:Nikos.P.Pitsianis@duke.edu
Email:xiaobai@cs.duke.edu
Email:jmccammo@ucsd.edu
Get access

Abstract

This paper summarizes the mathematical and numerical theories and computational elements of the adaptive fast multipole Poisson-Boltzmann (AFMPB) solver. We introduce and discuss the following components in order: the Poisson-Boltzmann model, boundary integral equation reformulation, surface mesh generation, the nodepatch discretization approach, Krylov iterative methods, the new version of fast multipole methods (FMMs), and a dynamic prioritization technique for scheduling parallel operations. For each component, we also remark on feasible approaches for further improvements in efficiency, accuracy and applicability of the AFMPB solver to large-scale long-time molecular dynamics simulations. The potential of the solver is demonstrated with preliminary numerical results.

Keywords

45B0565Y0568W1090B1092C0592C40Biomolecular systemelectrostaticsPoisson-Boltzmann equationfast multipole methodsmesh generationdirected acyclic graphdynamic prioritizationparallelization
Type
Research Article
Information
Communications in Computational Physics , Volume 13 , Issue 1 , January 2013 , pp. 107 - 128
Copyright
Copyright © Global Science Press Limited 2013

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]http://wiki.c2b2.columbia.edu/honiglab_public/index.php/Software:DelPhi.Google Scholar
[2]http://wiki.c2b2.columbia.edu/honiglab_public/index.php/Software:GRASP.Google Scholar
[3]http://mccammon.ucsd.edu/uhbd.html.Google Scholar
[4]http://www.poissonboltzmann.org/apbs.Google Scholar
[5]http://mgltools.scripps.edu/packages/MSMS/.Google Scholar
[6]http://fetk.org/codes/gamer/.Google Scholar
[7]http://lsec.cc.ac.cn/~lubz/Meshing.html.Google Scholar
[8]http://www-users.cs.umn.edu/~saad/software/SPARSKIT/index.html.Google Scholar
[9]http://www.fastmultipole.org/.Google Scholar
[10]Abramowitz, M. and Stegun, I. A.Handbook of Mathematical Functions. Dover, 1970.Google Scholar
[11]Appel, A. W.An efficient program for many-body simulation. SIAM Journal on Scientific and Statistical Computing, 6: 85–103,1985.CrossRefGoogle Scholar
[12]Barnes, J. and Hut, P.A hierarchical O(NlogN) force-calculation algorithm. Nature, 324: 446–449, 1986.Google Scholar
[13]Bashford, D.An object-oriented programming suite for electrostatic effects in biological molecules. Lecture Notes in Computer Science, 1343: 233–240, 1997.Google Scholar
[14]Bates, P. W., Wei, G. W., and Zhao, S.Minimal molecular surfaces and their applications. Journal of Computational Chemistry, 29: 380–391, 2008.Google Scholar
[15]Bondi, A.van der Waals volumes and radii. The Journal of Physical Chemistry, 68: 441–451, 1964.Google Scholar
[16]Boschitsch, A. H. and Fenley, M. O.Hybrid boundary element and finite difference method for solving the nonlinear Poisson-Boltzmann equation. Journal of Computational Chemistry, 25: 935–955, 2004.Google Scholar
[17]Boschitsch, A. H., Fenley, M. O., and Zhou, H. X.Fast boundary element method for the linear Poisson-Boltzmann equation. The Journal of Physical Chemistry B, 106: 2741–2754, 2002.Google Scholar
[18]Case, D. A., Cheatham, T. E. III, Darden, T., Gohlke, H., Luo, R., Merz, K. M., Onufriev, A., Simmerling, C., Wang, B., and Woods, R. J.The Amber biomolecular simulation programs. Journal of Computational Chemistry, 26: 1668–1688, 2005.Google Scholar
[19]Chen, M. X. and Lu, B. Z.TMSmesh: A robust method for molecular surface mesh generation using a trace technique. Journal of Chemical Theory and Computation, 7(1): 203–212, 2011.Google Scholar
[20]Cheng, H., Greengard, L., and Rokhlin, V.A fast adaptive multipole algorithm in three dimensions. Journal of Computational Physics, 155: 468–498, 1999.CrossRefGoogle Scholar
[21]D’Agostino, D., Clematis, A., Merelli, I., Milanesi, L., and Coloberti, M.A grid service based parallel molecular surface reconstruction system. In Proceedings of 16th Euromicro Conference on Parallel, Distributed and Network-Based Processing, pages 455–462, 2008.Google Scholar
[22]Darden, T., York, D., and Pedersen, L.Particle mesh Ewald: an NlogN method for Ewald sums in large systems. Journal of Chemical Physics, 98: 10089–10092, 1993.Google Scholar
[23]Dzubiella, J., Swanson, J. M. J., and McCammon, J. A.Coupling hydrophobicity, dispersion, and electrostatics in continuum solvent models. Physics Review Letters, 96: 087802, 2006.CrossRefGoogle ScholarPubMed
[24]Fogolari, F., Brigo, A., and Molinari, H.The Poisson-Boltzmann equation for biomolecular electrostatics: a tool for structural biology. Journal of Molecular Recognition, 15:377-392, 2002.Google Scholar
[25]Fritsch, S., Ivanov, I., Wang, H. L., and Cheng, X. L.Ion selectivity mechanism in a bacterial pentameric ligand-gated ion channel. Biophysical Journal, 100: 390–398, 2011.Google Scholar
[26]Gouaillard, A., Gelas, A., and Megason, S.Triangular meshes Delaunay conforming filter. Insight Journal, July-December, 2008.Google Scholar
[27]Grant, J. A., Pickup, B. T., and Nicholls, A.A smooth permittivity function for Poisson-Boltzmann solvation methods. Journal of Computational Chemistry, 22: 608–640, 2001.Google Scholar
[28]Greengard, L. and Gropp, W.A parallel version of the fast multipole method. Computers & Mathematics with Applications, 20: 63–71, 1990.Google Scholar
[29]Greengard, L. and Huang, J.A new version of the fast multipole method for screened Coulomb interactions in three dimensions. Journal of Computational Physics, 180: 642–658, 2002.CrossRefGoogle Scholar
[30]Greengard, L. and Rokhlin, V.A fast algorithm for particle simulations. Journal of Computational Physics, 73: 325–348, 1987.Google Scholar
[31]Greengard, L. and Rokhlin, V.A new version of the fast multipole method for the Laplace equation in three dimensions. Acta Numerica, 6: 229–269, 1997.Google Scholar
[32]Hsieh, M. J. and Luo, R.Physical scoring function based on AMBER force field and Poisson-Boltzmann implicit solvent for protein structure prediction. Proteins, 56: 475–486, 2004.Google Scholar
[33]Hu, T. C.Parallel sequencing and assembly line problems. Operations Research, 9: 841–848, 1961.Google Scholar
[34]Huang, J., Jia, J., and Zhang, B.FMM-Yukawa: an adaptive fast multipole method for screened Coulomb interactions. Computer Physics Communications, 180: 2331–2338, 2009.Google Scholar
[35]Im, W., Beglov, D., and Roux, B.Continuum solvation model: computation of electrostatic forces from numerical solutions to the Poisson-Boltzmann equation. Computer Physics Communications, 111: 59–75, 1998.Google Scholar
[36]Jo, S., Vargyas, M., Vasko-Szedlar, J., Roux, B., and Im, V.PBEQ-solver for online visualization of electrostatic potential of biomolecules. Nucleic Acids Research, 36: 270–275, 2008.Google Scholar
[37]Juffer, A. H., Botta, E. F. F., Vankeulen, B. A. M., Vanderploeg, A., and Berendsen, H. J. C.The electric potential of a macromolecule in a solvent: a fundamental approach. Journal of Computational Physics, 97: 144–171, 1991.Google Scholar
[38]Kuo, S., Altman, M., Bardhan, J., Tidor, B., and White, J.Fast methods for biomolecule charge optimization. In Proceedings of International Conference on Modeling and Simulation of Microsystems, 2002.Google Scholar
[39]Kuo, S. and White, J.A spectrally accurate integral equation solver for molecular surface electrostatics. In Proceedings of the IEEE Conference on Computer-Aided Design, 2006.Google Scholar
[40]Lee, B. and Richards, F.The interpretation of protein structures: estimation of static accessibility. Journal of Molecular Biology, 55: 379–400, 1971.Google Scholar
[41]Liang, J. and Subramaniam, S.Computation of molecular electrostatics with boundary element methods. Biophysical Journal, 73: 1830–1841, 1997.Google Scholar
[42]Lu, B., Cheng, X., Huang, J., and McCammon, J. A.AFMPB: an adaptive fast multipole Poisson-Boltzmann solver for calculating electrostatics in biomolecular systems. Computer Physics Communications, 181: 1150–1160, 2010.Google Scholar
[43]Lu, B., Cheng, X., and McCammon, J. A. “New-version-fast-multipolemethod” accelerated electrostatic calculations in biomolecular systems. Journal of Computational Physics, 226: 1348–1366, 2007.Google Scholar
[44]Lu, B. and McCammon, J. A.Improved boundary element methods for Poisson-Boltzmann electrostatic potential and force calculations. Journal of Chemical Theory and Computation, 3: 1134–1142, 2007.Google Scholar
[45]Lu, B., Zhou, Y. C., Huber, G. A., Bond, S. D., Holst, M. J., and McCammon, J. A.Electrodiffusion: a continuum modeling framework for biomolecular systems with realistic spatiotem-poral resolution. Journal of Chemical Physics, 127: 135102, 2007.Google Scholar
[46]Madura, J. D., Briggs, J. M., Wade, R. C., Davis, M. E., Luty, B. A., Ilin, A., Antosiewicz, J., Gilson, M. K., Bagheri, B., Scott, L. R., and McCammon, J. A.Electrostatics and diffusion of molecules in solution – simulations with the University of Houston Brownian Dynamics Program. Computer Physics Communication, 91: 57–95, 1995.CrossRefGoogle Scholar
[47]Nicholls, A., Sharp, K. A., and Honig, B.Protein folding and association: insights from the interfacial and thermodynamic properties of hydrocarbons. Proteins: Structure, Function, and Bioinformatics, 11: 281–296, 1991.CrossRefGoogle ScholarPubMed
[48]Orozco, M. and Lugue, F.Theoretical methods for the description of the solvent effect in biomolecular systems. Chemical Reviews, 100: 4187–4226, 2000.Google Scholar
[49]Phillips, J. R. and White, J. K.A precorrected-FFT method for electrostatic analysis of complicated 3D structures. IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems, 16: 1059–1072, 1997.Google Scholar
[50]Richards, F. M.Areas, volumes, packing and protein structure. Annual Review in Biophysics and Bioengineering, 6: 151–176, 1977.Google Scholar
[51]Rokhlin, V.Solution of acoustic scattering problems by means of second kind integral equations. Wave Motion, 5: 257–272, 1983.CrossRefGoogle Scholar
[52]Roux, B. and Simonson, T.Implicit solvent models. Biophysical Chemistry, 78: 1–20, 1999.Google Scholar
[53]Salmon, J.Parallel Hierarchical N-Body Methods. PhD thesis, California Institute of Technology, 1990.Google Scholar
[54]Sevilgen, F. E. and Aluru, S.A unifying data structure for hierarchical methods. In Proceedings of the ACM/IEEE Conference on Supercomputing, 1999.Google Scholar
[55]Singh, J., Holt, C., Hennessy, J., and Gupta, A.A parallel adaptive fast multipole method. In Proceedings of the ACM/IEEE Conference on Supercomputing, 1993.Google Scholar
[56]Singh, J., Holt, C., Totsuka, T., Gupta, A., and Hennessy, J.Load balancing and data locality in adaptive hierarchical N-body methods: Barnes-Hut, fast multipole, and radiosity. Journal of Parallel and Distributed Computing, 27: 118–141, 1995.Google Scholar
[57]Sun, X. and Pitsianis, N. P.A matrix version of the fast multipole method. SIAM Review, 43: 289–300, 2001.Google Scholar
[58]Teng, S.Provably good partitioning and load balancing algorithms for parallel adaptive N-body simulation. SIAM Journal on Scientific Computing, 19: 635–656, 1998.Google Scholar
[59]Ullman, J. D.NP-complete scheduling problems. Journal of Computer and System Sciences, 10: 384–393, 1975.Google Scholar
[60]Warren, M. and Salmon, J.A parallel hashed oct-tree N-body algorithm. In Proceedings of the ACM/IEEE Conference on Supercomputing, 1993.Google Scholar
[61]Weiser, J., Shenkin, P. S., and Still, W. C.Optimization of Gaussian surface calculations and extension to solvent-accessible surface areas. Journal of Computational Chemistry, 20: 688–703, 1999.Google Scholar
[62]Zhou, Y. C. and Feig, M. and Wei, G. W.Highly accurate biomolecular electrostatics in continuum dielectric environments. Journal of Computational Chemistry, 29: 87–97, 2007.Google Scholar
[63]Zhang, B., Huang, J., Pitsianis, N. P., and Sun, X.Dynamic prioritization for parallel traversal of irregularly structured spatio-temporal graphs. In Proceedings of the 3rd USENIX Workshop on Hot Topics in Parallelism, 2011.Google Scholar
[64]Zhang, Y., Xu, G., and Bajaj, C.Quality meshing of implicit solvation models of biomolecular structures. The Special Issue of Computer Aided Geometric Design on Applications of Geometric Modeling in the Life Sciences, 23: 510–530, 2006.CrossRefGoogle ScholarPubMed