ABSTRACT
Biological molecules perform their functions surrounded by water and mobile ions, which strongly influence molecular structure and behavior. The electrostatic interactions between a molecule and solvent are particularly difficult to model theoretically, due to the forces' long range and the collective response of many thousands of solvent molecules. The dominant modeling approaches represent the two extremes of the trade-off between molecular realism and computational efficiency: all-atom molecular dynamics in explicit solvent, and macroscopic continuum theory (the Poisson or Poisson--Boltzmann equation). We present the first fast-solver implementation of an advanced nonlocal continuum theory that combines key advantages of both approaches. In particular, molecular realism is included by limiting solvent dielectric response on short length scales, using a model for nonlocal dielectric response allows the resulting problem (a linear integro-differential Poisson equation) to be reformulated as a system of coupled boundary-integral equations using double reciprocity. Whereas previous studies using the nonlocal theory had been limited to small model problems, owing to computational cost, our work opens the door to studying much larger problems including rational drug design, protein engineering, and nanofluidics.
- M. D. Altman, J. P. Bardhan, B. Tidor, and J. K. White. FFTSVD: A fast multiscale boundary-element method solver suitable for BioMEMS and biomolecule simulation. IEEE T. Comput.-Aid. D., 25:274--284, 2006. Google ScholarDigital Library
- M. D. Altman, J. P. Bardhan, J. K. White, and B. Tidor. Accurate solution of multi-region continuum electrostatic problems using the linearized Poisson--Boltzmann equation and curved boundary elements. J. Comput. Chem., 30:132--153, 2009.Google ScholarCross Ref
- J. P. Bardhan. Numerical solution of boundary-integral equations for molecular electrostatics. J. Chem. Phys., 130:094102, 2009.Google ScholarCross Ref
- M. V. Basilevsky and D. F. Parsons. An advanced continuum medium model for treating solvation effects: Nonlocal electrostatics with a cavity. J. Chem. Phys., 105(9):3734, Aug 1996.Google ScholarCross Ref
- R. Diamond. Real-space refinement of the structure of hen egg-white lysozyme. J. Mol. Biol., 82, 1974.Google Scholar
- L. Greengard and V. Rokhlin. A fast algorithm for particle simulations. J. Comput. Phys., 73:325--348, 1987. Google ScholarDigital Library
- J. L. Hess and A. M. O. Smith. Calculation of non-lifting potential flow about arbitrary three-dimensional bodies. J. Ship Res., 8(2):22--44, 1962.Google Scholar
- A. Hildebrandt, R. Blossey, S. Rjasanow, O. Kohlbacher, and H.-P. Lenhof. Novel formulation of nonlocal electrostatics. Phys. Rev. Lett., 93:108104, 2004.Google ScholarCross Ref
- A. Hildebrandt, R. Blossey, S. Rjasanow, O. Kohlbacher, and H.-P. Lenhof. Electrostatic potentials of proteins in water: a structured continuum approach. Bioinformatics, 23(2):e99--e103, Jan 2007. Google ScholarDigital Library
- J. C. Inkson. Dielectric response in semiconductors. J. Phys. Chem., 7:1571--1581, 1974.Google Scholar
- S. Kapur and D. E. Long. IES3: Efficient electrostatic and electromagnetic simulation. IEEE Comput. Sci. Eng., 5(4):60--7, 1998. Google ScholarDigital Library
- I. Klapper, R. Hagstrom, R. Fine, K. Sharp, and B. Honig. Focusing of electric fields in the active site of Cu-Zn superoxide dismutase: Effects of ionic strength and amino-acid modification. Proteins, 1:47--59, 1986.Google ScholarCross Ref
- A. A. Kornyshev, A. I. Rubinshtein, and M. A. Vorotyntsev. Model nonlocal electrostatics: I. J. Phys. C Solid State, 11:3307, Dec 1978.Google ScholarCross Ref
- K. Nabors and J. White. FASTCAP: A multipole accelerated 3-D capacitance extraction program. IEEE T. Comput.-Aid. D., 10(10):1447--1459, 1991.Google ScholarDigital Library
- S. Pennathur and J. G. Santiago. Electrokinetic transport in nanochannels. 1. theory. Anal. Chem., 77(21):6772--6781, 2005.Google ScholarCross Ref
- B. Roux and T. Simonson. Implicit solvent models. Biophys. Chem., 78:1--20, 1999.Google ScholarCross Ref
- A. Rubinstein and S. Sherman. Influence of the solvent structure on the electrostatic interactions in proteins. Biophys. J., 87(3):1544--1557, Sep 2004.Google ScholarCross Ref
- S. Rush, A. H. Turner, and A. H. Cherin. Computer solution for time-invariant electric fields. J. Appl. Phys., 37(6):2211--2217, 1966.Google ScholarCross Ref
- Y. Saad and M. Schultz. GMRES: A generalized minimal residual algorithm for solving nonsymmetric linear systems. SIAM J. Sci. Stat. Comput., 7:856--869, 1986. Google ScholarDigital Library
- M. Sanner, A. J. Olson, and J. C. Spehner. Reduced surface: An efficient way to compute molecular surfaces. Biopolymers, 38:305--320, 1996.Google ScholarCross Ref
- M. F. Sanner. Molecular surface computation home page. http://www.scripps.edu/sanner/html/msms_home.html, 1996.Google Scholar
- K. A. Sharp and B. Honig. Electrostatic interactions in macromolecules: Theory and applications. Annu. Rev. Biophys. Bio., 19:301--332, 1990.Google ScholarCross Ref
- D. Sitkoff, K. A. Sharp, and B. Honig. Accurate calculation of hydration free energies using macroscopic solvent models. J. Phys. Chem. B, 98:1978--1988, 1994.Google ScholarCross Ref
- J. Tausch, J. Wang, and J. White. Improved integral formulations for fast 3-D method-of-moment solvers. IEEE T. Comput.-Aid. D., 20(12):1398--1405, 2001. Google ScholarDigital Library
- B. J. Yoon and A. M. Lenhoff. A boundary element method for molecular electrostatics with electrolyte effects. J. Comput. Chem., 11(9):1080--1086, 1990.Google ScholarCross Ref
Index Terms
- A fast solver for nonlocal electrostatic theory in biomolecular science and engineering
Recommendations
Fast Analytical Methods for Macroscopic Electrostatic Models in Biomolecular Simulations
We review recent developments of fast analytical methods for macroscopic electrostatic calculations in biological applications, including the Poisson-Boltzmann (PB) and the generalized Born models for electrostatic solvation energy. The focus is on ...
"New-version-fast-multipole-method" accelerated electrostatic calculations in biomolecular systems
In this paper, we present an efficient and accurate numerical algorithm for calculating the electrostatic interactions in biomolecular systems. In our scheme, a boundary integral equation (BIE) approach is applied to discretize the linearized Poisson-...
Comments