A spectrally accurate integral equation solver for molecular surface electrostatics
Pages 899 - 906
Abstract
Electrostatic analysis of complicated molecular surfaces arises in a number of nanotechnology applications including: biomolecule design, carbon nanotube simulation, and molecular electron transport. Molecular surfaces are typically smooth, without the corners common in electrical interconnect problems, and are candidates for methods with higher order convergence than that of the commonly used flat panel methods. In this paper we describe and demonstrate a spectrally accurate approach for analyzing molecular surfaces described by a collection of surface points. The method is a synthesis of several techniques, and starts by using least squares to fit a high order spherical harmonic surface representation to the given points. Then this analytic representation is used to construct a differentiable map from the molecular suface to a cube, an orthogonal basis is generated on the rectangular cube surfaces, and a change of variables is used to desingularize the required integrals of products of basis functions and Green's function. Finally, an efficient method for solving the discretized system using a matrix-implicit scheme is described. The combined method is demonstrated on an analytically solvable sphere problem, capacitance calculation of complicated molecular surface, and a coupled Poisson/Poisson-Boltzmann problem associated with a biomolecule. The results demonstrate that for a tolerance of 10-3 this new approach requires one to two orders of magnitude fewer unknowns than a flat panel method.
References
[1]
A. Baxansky and N. Kiryati. Calculating geometric properties of three-dimensional objects from the spherical harmonic representation. Technical Report VIA-2005-6-1, Tel Aviv University, June 2005.
[2]
H. M. Berman, J. Westbrook, Z. Feng, G. Gilliland, T. N. Bhat, H. Weissig, I. N. Shindyalov, and P. E. Bourne. The Protein Data Bank. Nucleic Acids Research, 28(1):235--242, 2000.
[3]
J.-P. Berrut and L. N. Trefethen. Barycentric lagrange interpolation. SIAM Review, 46(3):501--517, 2004.
[4]
G. Beylkin, R. Coifman, and V. Rokhlin. Fast wavelet transforms and numerical algorithms I. Communications on Pure and Applied Mathematics, 44:141--183, 1991.
[5]
R. Bharadwaj, A. Windemuth, S. Sridharan, B. Honig, and A. Nicholls. The fast multipole boundary element method for molecular electrostatics: An optimal approach for large systems. Journal of Computational Chemistry, 16:898--913, 1995.
[6]
A. Boschitsch, M. Fenley, and H.-X. Zhou. Fast boundary element method for the linear poisson-boltzmann equation. Journal of Physical Chemistry B, 106(10):2741--2754, 2002.
[7]
C. Brechbühler, G. Gerig, and O. Kübler. Parametrization of closed surfaces for 3-D shape description. Computer Vision and Image Understanding, 61(2):154--170, 1995.
[8]
O. Bruno and L. Kunyansky. A fast, high-order algorithm for the solution of surface scattering problems: basic implementation, tests, and applications. Journal of Computational Physics, 169(1):80--110, 2001.
[9]
B. Büchmann. Accuracy and stability of a set of free-surface time-domain boundary element models based on B-splines. International Journal for Numerical Methods in Fluids, 33(1):125--155, 2000.
[10]
M. L. Connolly. Analytical molecular surface calculation. Journal of Applied Crystallography, 16:548--558, 1983.
[11]
B. S. Duncan and A. J. Olson. Approximation and characterization of molecular surfaces. Biopolymers, 33(2):219--229, Feb. 1993.
[12]
E. J. Garboczi. Three-dimensional mathematical analysis of particle shape using x-ray tomography and spherical harmonics: Application to aggregates used in concrete. Cement and Concrete Research, 32(10):1621--1638, Oct. 2002.
[13]
L. Greengard. The Rapid Evaluation of Potential Fields in Particle Systems. MIT Press, 1988.
[14]
W. Hackbusch and Z. P. Nowak. On the fast matrix multiplication in the boundary element method by panel clustering. Numerische Mathematik, 54:463--491, 1989.
[15]
R. F. Harrington. Field Computation by Moment Methods. MacMillan, 1968.
[16]
J. L. Hess and A. M. O. Smith. Calculation of nonlifting potential flow about arbitrary three-dimensional bodies. Journal of Ship Research, 8(2):22--44, 1964.
[17]
B. Honig and A. Nicholls. Classical electrostatics in biology and chemistry. Science, 268:1144--1149, 1995.
[18]
A. H. Juffer, E. F. F. Botta, B. A. M. V. Keulen, A. V. D. Ploeg, and H. J. C. Berendsen. The electric potential of a macromolecule in a solvent: A fundamental approach. Journal of Computational Physics, 97:144--171, 1991.
[19]
S. Kapur and D. Long. High-order Nyström schemes for efficient 3-D capacitance extraction. In Proceedings of the International Conference on Computer-Aided Design, pages 178--185, 1998.
[20]
S. Kapur and D. E. Long. IES3: A fast integral equation solver for efficient 3-dimensional extraction. In Proceedings of the International Conference on Computer-Aided Design, pages 448--455, San Jose, CA, 1997.
[21]
S. S. Kuo, M. D. Altman, J. P. Bardhan, B. Tidor, and J. K. White. Fast methods for simulation of biomolecule electrostatics. In Proceedings of the International Conference on Computer-Aided Design, San Jose, CA, Nov. 2002.
[22]
A. Y. Lee, P. A. Karplus, B. Ganem, and J. Clardy. Atomic structure of the buried catalytic pocket of escherichia coli chorismate mutase. Journal of the American Chemical Society, 117(12):3627--3628, 1995.
[23]
S. Leicester, J. Finney, and R. Bywater. A quantitative representation of molecular surface shape. i: Theory and development of the method. Journal of Mathematical Chemistry, 16:315--341, 1994.
[24]
S. E. Leicester, J. L. Finney, and R. P. Bywater. Description of molecular surface shape using fourier descriptors. Journal of Molecular Graphics, 6(2):104--108, 1988.
[25]
J. D. Madura, J. M. Briggs, R. C. Wade, M. E. Davis, B. A. Luty, A. Ilin, J. Antosiewicz, M. K. Gilson, B. Bagheri, L. Ridgway-Scott, and J. A. McCammon. Electrostatics and diffusion of molecules in solution: Simulations with the University of Houston Brownian Dynamics program. Computer Physics Communications, 91:57--95, 1995.
[26]
H. D. Manier. A Three Dimensional Higher Order Panel Method based on B-splines. PhD thesis, Massachusetts Institute of Technology, Cambridge, MA, 1995.
[27]
N. L. Max and E. D. Getzoff. Spherical harmonic molecular surfaces. IEEE Computer Graphics and Applications, 8(4):42--50, 1988.
[28]
R. J. Morris, R. J. Najmanovich, A. Kahraman, and J. M. Thornton. Real spherical harmonic expansion coefficients as 3D shape descriptors for protein binding pocket and ligand comparisons. Bioinformatics, 21(10):2347--2355, Feb. 2005.
[29]
K. Nabors, F. T. Korsmeyer, F. T. Leighton, and J. White. Preconditioned, adaptive, multipole-accelerated iterative methods for three-dimensional first-kind integral equations of potential theory. SIAM Journal on Scientific Computing, 15(3):713--735, 1994.
[30]
K. Nabors and J. White. FastCap: A multipole accelerated 3-D capacitance extraction program. IEEE Transactions on Computer-Aided Design, 10(11):1447--1459, Nov. 1991.
[31]
A. T. Patera. A spectral element method for fluid dynamics: Laminar flow in a channel expansion. Journal of Computational Physics, 54:468--488, June 1984.
[32]
J. R. Phillips and J. K. White. A precorrected-FFT method for electrostatic analysis of complicated 3-D structures. IEEE Transactions on Computer-Aided Design, 16(10):1059--1072, Oct. 1997.
[33]
D. Ramaswamy, W. Ye, X. Wang, and J. White. Fast algorithms for 3-D simulation. Journal of Modeling and Simulation of Microsystems, 1(1):77--82, Dec. 1999.
[34]
F. M. Richards. Areas, volumes, packing, and protein structure. Annual Review of Biophysics and Bioengineering, 6:151--176, 1977.
[35]
V. Rokhlin. Rapid solution of integral equations of classical potential theory. Journal of Computational Physics, 60:187--207, 1985.
[36]
Y. Saad and M. Schultz. GMRES: A generalized minimal residual algorithm for solving nonsymmetric linear systems. SIAM Journal of Scientific and Statistical Computing, 7:856--869, 1986.
[37]
M. Sanner, A. J. Olson, and J. C. Spehner. Reduced surface: An efficient way to compute molecular surfaces. Biopolymers, 38:305--320, 1996.
[38]
K. A. Sharp and B. Honig. Electrostatic interactions in macromolecules: Theory and applications. Annual Review of Biophysics and Biophysical Chemistry, 19:301--332, 1990.
[39]
C. Tanford and J. G. Kirkwood. Theory of protein titration curves I. general equations for impenetrable spheres. Journal of the American Chemical Society, 59:5333--5339, 1957.
[40]
L. N. Trefethen and D. Bau, III. Numerical Linear Algebra. Society for Industrial & Applied Mathematics, 1997.
[41]
J. C. Vassberg. A fast surface-panel method capable of solving million-element problems. Aerospace Sciences Meeting and Exhibit, 35th, Reno, NV, Jan. 1997. AIAA Paper 97--0168.
[42]
J. Warwicker and H. C. Watson. Calculation of the electric potential in the active site cleft due to alpha-helix dipoles. Journal of Molecular Biology, 157:671--679, 1982.
[43]
B. J. Yoon and A. M. Lenhoff. A boundary element method for molecular electrostatics with electrolyte effects. Journal of Computational Chemistry, 11:1080--1086, 1990.
[44]
R. J. Zauhar and R. S. Morgan. The rigorous computation of the molecular electric potential. Journal of Computational Chemistry, 9:171--187, 1988.
[45]
R. J. Zauhar and A. Varnek. A fast and space-efficient boundary element method for computing electrostatic and hydration effects in large molecules. Journal of Computational Chemistry, 17:864--877, 1996.
[46]
Z. Zhu, B. Song, and J. White. Algorithms in FastImp: A fast and wideband impedance extraction program for complicated 3-D geometries. In Proceedings of the Design Automation Conference, pages 712--717, Anaheim, CA, June 2003.
Index Terms
- A spectrally accurate integral equation solver for molecular surface electrostatics
Recommendations
A meshless, spectrally accurate, integral equation solver for molecular surface electrostatics
The need to determine electrostatic fields in domains bounded by molecular surfaces arises in a number of nanotechnology applications including: biomolecule design, carbon nanotube simulation, and molecular electron transport analysis. Molecular ...
Hierarchical molecular interfaces and solvation electrostatics
SPM '09: 2009 SIAM/ACM Joint Conference on Geometric and Physical ModelingElectrostatic 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 ...
A new boundary integral equation for molecular electrostatics with charges over whole space
AbstractIn this paper, a new integral equation of electrostatics is proposed as an integral form of a basic dielectric continuum model, which is traditionally represented in a form of Poisson differential equation. As an application in protein simulations,...
Comments
Please enable JavaScript to view thecomments powered by Disqus.Information & Contributors
Information
Published In
Copyright © 2006 ACM.
Permission to make digital or hard copies of all or part of this work for personal or classroom use is granted without fee provided that copies are not made or distributed for profit or commercial advantage and that copies bear this notice and the full citation on the first page. Copyrights for components of this work owned by others than ACM must be honored. Abstracting with credit is permitted. To copy otherwise, or republish, to post on servers or to redistribute to lists, requires prior specific permission and/or a fee. Request permissions from [email protected]
Sponsors
- SIGDA: ACM Special Interest Group on Design Automation
- IEEE-CAS: Circuits & Systems
- IEEE-CS: Computer Society
Publisher
Association for Computing Machinery
New York, NY, United States
Publication History
Published: 05 November 2006
Check for updates
Qualifiers
- Article
Conference
ICCAD06
Sponsor:
- SIGDA
- IEEE-CAS
- IEEE-CS
ICCAD06: The International Conference on Computer-Aided Design 2006
November 5 - 9, 2006
California, San Jose
Acceptance Rates
Overall Acceptance Rate 457 of 1,762 submissions, 26%
Contributors
Other Metrics
Bibliometrics & Citations
Bibliometrics
Article Metrics
- 0Total Citations
- 194Total Downloads
- Downloads (Last 12 months)0
- Downloads (Last 6 weeks)0
Reflects downloads up to 13 Nov 2024
Other Metrics
Citations
View Options
Get Access
Login options
Check if you have access through your login credentials or your institution to get full access on this article.
Sign in