The Molecular Mechanics Poisson-Boltzmann Surface Area (MM/PBSA) methods have been widely used in calculating protein-ligand binding affinities. However, the efficiency is still limited due to the high dimensionality of typical biomolecular systems while numerous snapshots need to be processed. Here we report our efforts to develop a highly efficient MM/PBSA program on graphics processing units (GPUs) utilizing Nvidia<sup>®</sup> CUDA. First, we implemented several linear Poisson-Boltzmann equation (PBE) solvers which are the bottlenecks of PBSA calculation onto GPUs. The preconditioned conjugate gradient (CG) and geometric multigrid (MG) solvers were found to be among the most efficient. Second, we rewrote most of the computation-intensive setup routines of PBSA into GPU kernels. The level set density calculation, the surface area non-bonded list determination, and the reaction field energy computation were the most time-consuming steps. Finally, we implemented the full GPU-supported MM/PBSA program by hooking the GPU-accelerated molecular mechanics (MM) energy calculation with the PBSA CUDA code. In addition, different matrix formats such as DIA or CSR, the unified memory, and the CUSP library were also explored to optimize the code efficiency. The second-order IIM-based interface algorithm was also implemented on GPUs using the biconjugate gradient (BiCG) method. Impressive efficiency speedups of MM/PBSA calculation with a range of biomolecular simulations were achieved. The new GPU-accelerated MM/PBSA code is to be released via the Amber package, and will significantly improve the computational throughput of the protein-ligand binding affinities.

Accelerating MM/PBSA calculation of protein-ligand binding on graphics processing units

2019 Fall National Meeting

The official archive and repository for all ACS Meetings & Expositions abstracts and select posters or presentations.