Numerical Methods for Biomolecule ElectrostaticsJaydeep P. Bardhan, Michael D. Altman, Jacob K. White & Bruce TidorIntroductionElectrostatic interactions play important roles in biomolecule structure and function. In contrast to hydrogen bonding and van der Waals interactions, however, electrostatic interactions act over much longer distances. This longrange nature significantly complicates computational modeling. Since computational expense prohibits explicitly modeling the many solvent molecules surrounding a biomolecule of interest, approximate models have been introduced that capture solvent effects in an average way [1]. In these models, macroscopic laws of electrostatics are assumed to hold in the molecular interior and in the solvent region exterior; the charge distribution in the molecule is assumed to be a set of discrete point charges located at the atom centers, and the potential satisfies the Poisson equation in this region. In the solvent exterior, mobile water molecules and salt ions screen charges and the PoissonBoltzmann equation governs the electrostatic potential. Even these simplified models, however, require large quantities of computer time and memory. In this project we focus on designing new, more efficient numerical techniques for computing the electrostatic interactions between biomolecules and surrounding solvent. Previous WorkFinite difference methods (FDM) are one of the most popular approaches for molecular electrostatics [1]; in these methods, a grid is laid down over the problem domain, and one finds a set of electrostatic potentials at the grid points that approximately solve the underlying partial differential equations. The finite difference approach is straightforward to understand and implement, but suffers from several sources of error. First, the irregular moleculesolvent interface cannot be accurately resolved on a regular grid. Second, the point charges that represent the molecular charge distribution must be spread to nearby grid points. Finally, the mathematical model has an infinite domain, which the finite difference method is forced to truncate. The boundary element method (BEM) is another technique that can be used to model molecular electrostatics [2]. The partial differential equations governing the potential throughout space are converted to an integral equation, or set of integral equations, over the molecular surface. The transformation therefore reduces the dimensionality of unknowns from three to two. In addition, the integral equation suffers from none of the gridbased inaccuracies inherent to finite difference methods: boundary conditions at infinity are treated exactly, point charges can be treated exactly, and the molecular surface can be represented much more accurately. ApproachThe linear systems that arise in the BEM are dense. Solving these systems by Gaussian elimination therefore has memory requirements that grow quadratically in time with respect to the number of unknowns, and time requirements grow cubically. For large systems, such cost is prohibitive, and we use Krylov subspace iterative methods such as GMRES instead. Using Krylov methods allows us to avoid the cubic time cost, but the dense matrixvector multiply still limits performance. We therefore use techniques that compute the matrixvector product approximately but in much less time and memory. We have shown, for instance, that the precorrectedFFT algorithm [3] can be used to solve the BEM formulation of the biomolecule electrostatics problem [2]. To reduce memory requirements, we have developed a new algorithm, which we call FFTSVD, for biomolecule electrostatics and bioMEMS simulations [4]. The algorithm is based on compressing longrange interactions using an approximate singular value decomposition (SVD), and then quickly computing the interactions using the FFT. SummaryThe boundary element method holds significant promise for improving the accuracy and speed of biomolecule electrostatics simulations. We have already demonstrated the viability of a fast solver approach [2]. We believe that the FFTSVD algorithm and future techniques will allow us to dramatically reduce the memory and time requirements for these calculations. References[1] M. K. Gilson, K. A. Sharp, and B. H. Honig. Calculating the electrostatic potential of molecules in solution: Method and error assessment. Journal of Computational Chemistry, 9:327335, 1987. [2] S. S. Kuo, M. D. Altman, J. P. Bardhan, B. Tidor, and J. K. White. Fast Methods for Biomolecule Electrostatics. International Conference on ComputerAided Design (ICCAD), 2002. [3] J. Phillips and J. K. White. A PrecorrectedFFT Method for Electrostatic Analysis of Complicated 3D Structures. IEEE Transactions on ComputerAided Design of Integrated Circuits and Systems, 16:10591072, 1997. [4] M. D. Altman, J. P. Bardhan, B. Tidor, and J. K. White. FFTSVD: A Fast MultiScale BEM Algorithm Suitable for BioMEMS Design. submitted to IEEE Transactions on ComputerAided Design of Integrated Circuits and Systems. 

