Solving the Linearized Poisson-Boltzmann Equation on GPUs Using CUDA

In this work an implementation of a linearized Poisson-Boltzmann equation solver based on a Finite Differences scheme on the GPU architecture is presented. The algorithm exploits the checkerboard structure of the discretized Laplace operator and follows the footprints of a popular solver called DelPhi, which is widely used in the Computational Biology community. The algorithm has been implemented using CUDA. This implementation has then been integrated with the DelPhi solver and tested over a few representative cases of biological interest. Details of the implementation as well as performance test results are illustrated.

[1]  Matthew G. Knepley,et al.  Biomolecular electrostatics using a fast multipole BEM on up to 512 gpus and a billion unknowns , 2010, Comput. Phys. Commun..

[2]  G Neshich,et al.  Electrostatic potential calculation for biomolecules--creating a database of pre-calculated values reported on a per residue basis for all PDB protein structures. , 2007, Genetics and molecular research : GMR.

[3]  高等学校計算数学学報編輯委員会編.  高等学校計算数学学報 = Numerical mathematics , 1979 .

[4]  M. Spagnuolo,et al.  Between algorithm and model: different Molecular Surface definitions for the Poisson-Boltzmann based electrostatic characterization of biomolecules in solution. , 2013, Communications in computational physics.

[5]  Arnaldo J. Montagner,et al.  JavaProtein Dossier: a novel web-based data visualization tool for comprehensive analysis of protein structure , 2004, Nucleic Acids Res..

[6]  Chuan Li,et al.  Highly efficient and exact method for parallelization of grid‐based algorithms and its implementation in DelPhi , 2012, J. Comput. Chem..

[7]  Thomas B Woolf,et al.  Discrimination of native loop conformations in membrane proteins: Decoy library design and evaluation of effective energy scoring functions , 2003, Proteins.

[8]  Aaas News,et al.  Book Reviews , 1893, Buffalo Medical and Surgical Journal.

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

[10]  J. Trylska,et al.  Continuum molecular electrostatics, salt effects, and counterion binding—A review of the Poisson–Boltzmann theory and its modifications , 2008, Biopolymers.

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

[12]  Walter Rocchia,et al.  Poisson-boltzmann equation boundary conditions for biological applications , 2005, Math. Comput. Model..