A new FFT-based algorithm to compute Born radii in the generalized Born theory of biomolecule solvation
- Department of Mathematics and Statistics, University of North Carolina at Charlotte, Charlotte, NC 28223 (United States)
- Department of Physics and Optical Science, University of North Carolina at Charlotte, Charlotte, NC 28223 (United States)
In this paper, a new method for calculating effective atomic radii within the generalized Born (GB) model of implicit solvation is proposed, for use in computer simulations of biomolecules. First, a new formulation for the GB radii is developed, in which smooth kernels are used to eliminate the divergence in volume integrals intrinsic in the model. Next, the fast Fourier transform (FFT) algorithm is applied to integrate smoothed functions, taking advantage of the rapid spectral decay provided by the smoothing. The total cost of the proposed algorithm scales as O(N{sup 3}logN+M) where M is the number of atoms comprised in a molecule and N is the number of FFT grid points in one dimension, which depends only on the geometry of the molecule and the spectral decay of the smooth kernel but not on M. To validate our algorithm, numerical tests are performed for three solute models: one spherical object for which exact solutions exist and two protein molecules of differing size. The tests show that our algorithm is able to reach the accuracy of other existing GB implementations, while offering much lower computational cost.
- OSTI ID:
- 21167743
- Journal Information:
- Journal of Computational Physics, Vol. 227, Issue 24; Other Information: DOI: 10.1016/j.jcp.2008.08.015; PII: S0021-9991(08)00448-8; Copyright (c) 2008 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9991
- Country of Publication:
- United States
- Language:
- English
Similar Records
Treecode-based generalized Born method
Exploiting graphical processing units to enable quantum chemistry calculation of large solvated molecules with conductor-like polarizable continuum models