Improved strategy in analytic surface calculation for molecular systems: Handling of singularities and computational efficiency

Computer methods for analytic surface calculations of molecular systems suffer from numerical instabilities and are CPU time consuming. In this article, we present proposals toward the solution of both problems. Singularities arise when nearly collinear triples of neighboring atoms or multiple vertices are encountered during the calculation. Topological decisions in analytic surface calculation algorithms (accessibility of vertices and arcs) are based upon the comparison of distances or angles. If two such numbers are nearly equal, then currently used computer programs may not resolve this ambiguity correctly and can subsequently fail. In this article, modifications in the analytic surface calculation algorithm are described that recognize singularities automatically and treat them appropriately without restarting parts of the computation. The computing time required to execute these alterations is minimal. The basic modification consists in defining an accuracy limit within which two values may be assumed as equal. The search algorithm has been reformulated to reduce the computational effort. A new set of formulas makes it possible to avoid mostly the extraction of square roots. Tests for small‐and medium‐sized intersection circles and for pairs of vertices with small vertex height help recognize fully buried circles and vertex pairs at an early stage. The new program can compute the complete topology of the surface and accessible surface area of the protein crambin in 1.50–4.29 s (on a single R3000 processor of an SGI 4D/480) depending on the compactness of the conformation where the limits correspond to the fully extended or fully folded chain, respectively. The algorithm, implemented in a computer program, will be made available on request. © John Wiley & Sons, Inc.

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

[2]  A. Shrake,et al.  Environment and exposure to solvent of protein atoms. Lysozyme and insulin. , 1973, Journal of molecular biology.

[3]  B. Fedorov,et al.  Improved technique for calculating X‐ray scattering intensity of biopolymers in solution: Evaluation of the form, volume, and surface of a particle , 1983 .

[4]  K. D. Gibson,et al.  Surface area of the intersection of three spheres with unequal radii A simplified analytical formula , 1988 .

[5]  R. Abagyan,et al.  New methodology for computer-aided modelling of biomolecular structure and dynamics. 1. Non-cyclic structures. , 1989, Journal of biomolecular structure & dynamics.

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

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

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

[9]  Cyrus Levinthal,et al.  A vectorized algorithm for calculating the accessible surface area of macromolecules , 1991 .

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

[11]  Manfredo P. do Carmo,et al.  Differential geometry of curves and surfaces , 1976 .

[12]  P Argos,et al.  Optimal protocol and trajectory visualization for conformational searches of peptides and proteins. , 1992, Journal of molecular biology.

[13]  B Honig,et al.  Extracting hydrophobic free energies from experimental data: relationship to protein folding and theoretical models. , 1991, Biochemistry.

[14]  Iñaki Tuñón,et al.  GEPOL: An improved description of molecular surfaces II. Computing the molecular area and volume , 1991 .

[15]  M. L. Connolly Analytical molecular surface calculation , 1983 .

[16]  A. Lipanov,et al.  Charge grouping approaches to calculation of electrostatic forces in molecular dynamics of macromolecules. , 1991, Journal of biomolecular structure & dynamics.

[17]  Michael L. Connolly,et al.  Molecular surface Triangulation , 1985 .

[18]  A. D. McLachlan,et al.  Solvation energy in protein folding and binding , 1986, Nature.

[19]  R M Stroud,et al.  Inclusion of solvation free energy with molecular mechanics energy: Alanyl dipeptide as a test case , 1992, Protein science : a publication of the Protein Society.

[20]  F M Richards,et al.  Areas, volumes, packing and protein structure. , 1977, Annual review of biophysics and bioengineering.

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

[22]  M. Karplus,et al.  An analysis of incorrectly folded protein models. Implications for structure predictions. , 1984, Journal of molecular biology.

[23]  C Sander,et al.  Polarity as a criterion in protein design. , 1989, Protein engineering.

[24]  Kenneth M. Merz,et al.  Rapid approximation to molecular surface area via the use of Boolean logic and look‐up tables , 1993, J. Comput. Chem..

[25]  I. Tuñón,et al.  Molecular surface area and hydrophobic effect. , 1992, Protein engineering.

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

[27]  B Honig,et al.  Reconciling the magnitude of the microscopic and macroscopic hydrophobic effects. , 1991, Science.

[28]  W. Kauzmann Some factors in the interpretation of protein denaturation. , 1959, Advances in protein chemistry.

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

[30]  R A Abagyan,et al.  New methodology for computer-aided modelling of biomolecular structure and dynamics. 2. Local deformations and cycles. , 1989, Journal of biomolecular structure & dynamics.

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

[32]  H. Scheraga,et al.  Empirical solvation models can be used to differentiate native from near‐native conformations of bovine pancreatic trypsin inhibitor , 1991, Proteins.

[33]  R. Fausto,et al.  The CαC internal rotation in α‐alkyl substituted carbonyls and thiocarbonyls: CH(CH3)2C(X)YH (X, Y O or S) , 1991 .