Optimization of Gaussian surface calculations and extension to solvent‐accessible surface areas

We explored the use of several breadth‐first and depth‐first algorithms for the computation of Gaussian atomic and molecular surface areas. Our results for whole‐molecule van der Waals surface areas (vdWSAs) were 10 times more accurate in relative error, relative to actual hard‐sphere areas, than those reported by earlier workers. We were also able to extend the method to the computation of solvent‐accessible surface areas (SASAs). This was made possible by an appropriate combination of algorithms, parameters, and preprocessing steps. For united‐atom 3app, a 2366‐atom protein, we obtained an average absolute atomic error of 1.16 Å2 with respect to the hard‐sphere atomic SASA results in 7 s of CPU time on an R10000/194 MHz processor. Speed and accuracy were both optimized for SASA by the use of neighbor‐list reduction (NLR), buried‐atom elimination (BAE), and a depth‐first search of the tree of atomic intersections. Accuracy was further optimized by the application of atom type specific parameters to the raw Gaussian results. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 688–703, 1999

[1]  G. Chang,et al.  Macromodel—an integrated software system for modeling organic and bioorganic molecules using molecular mechanics , 1990 .

[2]  J. Janin,et al.  Analytical approximation to the accessible surface area of proteins. , 1980, Proceedings of the National Academy of Sciences of the United States of America.

[3]  Doros N. Theodorou,et al.  Analytical treatment of the volume and surface area of molecules formed by an arbitrary collection of unequal spheres intersected by planes , 1991 .

[4]  T. Richmond,et al.  Solvent accessible surface area and excluded volume in proteins. Analytical equations for overlapping spheres and implications for the hydrophobic effect. , 1984, Journal of molecular biology.

[5]  B. Lee,et al.  Hydrophobic potential by pairwise surface area sum. , 1995, Protein engineering.

[6]  B. Lee,et al.  The interpretation of protein structures: estimation of static accessibility. , 1971, Journal of molecular biology.

[7]  W. C. Still,et al.  A rapid approximation to the solvent accessible surface areas of atoms , 1988 .

[8]  Akbar Nayeem,et al.  MSEED: A program for the rapid analytical determination of accessible surface areas and their derivatives , 1992 .

[9]  Jörg Weiser,et al.  Fast, approximate algorithm for detection of solvent‐inaccessible atoms , 1999 .

[10]  Armin A. Weiser,et al.  Neighbor‐list reduction: Optimization for computation of molecular van der Waals and solvent‐accessible surface areas , 1998 .

[11]  J. Weiser,et al.  HUNTER: A conformational search program for acyclic to polycyclic molecules with special emphasis on stereochemistry , 1997, J. Comput. Chem..

[12]  Y. Cui,et al.  Protein folding simulation with genetic algorithm and supersecondary structure constraints , 1998, Proteins.

[13]  W. Braun,et al.  Folding of Globular Proteins by Energy Minimization and Monte Carlo Simulations with Hydrophobic Surface Area Potentials , 1995 .

[14]  S. Sridharan,et al.  A rapid method for calculating derivatives of solvent accessible surface areas of molecules , 1995, J. Comput. Chem..

[15]  Roberto Cammi,et al.  Analytical first derivatives of molecular surfaces with respect to nuclear coordinates , 1996, J. Comput. Chem..

[16]  Donald G. Truhlar,et al.  Improved methods for semiempirical solvation models , 1995, J. Comput. Chem..

[17]  J. Ponder,et al.  Protein structure prediction using a combination of sequence homology and global energy minimization: II. Energy functions , 1998 .

[18]  J. A. Grant,et al.  A fast method of molecular shape comparison: A simple application of a Gaussian description of molecular shape , 1996, J. Comput. Chem..

[19]  S L Mayo,et al.  Pairwise calculation of protein solvent-accessible surface areas. , 1998, Folding & design.

[20]  Frank Eisenhaber,et al.  Improved strategy in analytic surface calculation for molecular systems: Handling of singularities and computational efficiency , 1993, J. Comput. Chem..

[21]  J. A. Grant,et al.  A Gaussian Description of Molecular Shape , 1995 .

[22]  Werner Braun,et al.  Minimization of empirical energy functions in proteins including hydrophobic surface area effects , 1993, J. Comput. Chem..

[23]  G J Williams,et al.  The Protein Data Bank: a computer-based archival file for macromolecular structures. , 1977, Journal of molecular biology.

[24]  C. Sander,et al.  An effective solvation term based on atomic occupancies for use in protein simulations , 1993 .

[25]  Harold A. Scheraga,et al.  An efficient, differentiable hydration potential for peptides and proteins , 1996, J. Comput. Chem..

[26]  K. D. Gibson,et al.  Exact calculation of the volume and surface area of fused hard-sphere molecules with unequal atomic radii , 1987 .

[27]  D. Eisenberg,et al.  Atomic solvation parameters applied to molecular dynamics of proteins in solution , 1992, Protein science : a publication of the Protein Society.

[28]  W. C. Still,et al.  Approximate atomic surfaces from linear combinations of pairwise overlaps (LCPO) , 1999 .

[29]  M. Sanner,et al.  Reduced surface: an efficient way to compute molecular surfaces. , 1996, Biopolymers.

[30]  R C Wade,et al.  Analytically defined surfaces to analyze molecular interaction properties. , 1996, Journal of molecular graphics.

[31]  Werner Braun,et al.  Exact and efficient analytical calculation of the accessible surface areas and their gradients for macromolecules , 1998 .

[32]  Robert B. Hermann,et al.  Theory of hydrophobic bonding. II. Correlation of hydrocarbon solubility in water with solvent cavity surface area , 1972 .