Approximate atomic surfaces from linear combinations of pairwise overlaps (LCPO)

A fast analytical formula was derived for the calculation of approximate atomic and molecular van der Waals (vdWSA), and solvent‐accessible surface areas (SASAs), as well as the first and second derivatives of these quantities with respect to atomic coordinates. This method makes use of linear combinations of terms composed from pairwise overlaps of hard spheres; therefore, we term this the LCPO method for linear combination of pairwise overlaps. For higher performance, neighbor‐list reduction (NLR) was applied as a preprocessing step. Eighteen compounds of different sizes (8–2366 atoms) and classes (organic, proteins, DNA, and various complexes) were chosen as representative test cases. LCPO/NLR computed the SASA and first derivatives of penicillopepsin, a protein with 2366 atoms, in 0.87 s (0.22 s for the creation of the neighbor list, 0.35 s for NLR, and 0.30 s for SASA and first derivatives) on an SGI R10000/194 Mhz processor. This appears comparable to or better than timings reported previously for other algorithms. The vdWSAs were in good agreement with the numerical results: relative errors for total molecular surface areas ranged from 0.1 to 2.0% and average absolute atomic surface area deviations from 0.3 to 0.7 Å2. For SASAs without NLR, the LCPO method exhibited relative errors in the range of 0.4–9.2% for total molecular surface areas and average absolute atomic surface area deviations of 2.0–2.7 Å2; with NLR the relative molecular errors ranged from 0.1 to 7.8% and the average absolute atomic surface area deviation from 1.6 to 3.0 Å2. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 217–230, 1999

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

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

[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]  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.

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

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

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

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

[9]  H. Scheraga,et al.  Accessible surface areas as a measure of the thermodynamic parameters of hydration of peptides. , 1987, Proceedings of the National Academy of Sciences of the United States of America.

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

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

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

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

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

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

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

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

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

[19]  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.

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

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

[22]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

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

[24]  W. C. Still,et al.  Semianalytical treatment of solvation for molecular mechanics and dynamics , 1990 .

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

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

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