Hierarchical molecular interfaces and solvation electrostatics

Electrostatic interactions play a significant role in determining the binding affinity of molecules and drugs. While significant effort has been devoted to the accurate computation of biomolecular electrostatics based on an all-atomic solution of the Poisson-Boltzmann (PB) equation for smaller proteins and nucleic acids, relatively little has been done to optimize the efficiency of electrostatic energetics and force computations of macromolecules at varying resolutions (also called coarse-graining). We have developed an efficient and comprehensive framework for computing coarse-grained PB electrostatic potentials, polarization energetics and forces for smooth multi-resolution representations of almost all molecular structures, available in the PDB. Important aspects of our framework include the use of variational methods for generating C2-smooth and multi-resolution molecular surfaces (as dielectric interfaces), a parameterization and discretization of the PB equation using an algebraic spline boundary element method, and the rapid estimation of the electrostatic energetics and forces using a kernel independent fast multipole method. We present details of our implementation, as well as several performance results on a number of examples.