Novel computer program for fast exact calculation of accessible and molecular surface areas and average surface curvature

New computer programs, SurfRace and FastSurf, perform fast calculations of the solvent accessible and molecular (solvent excluded) surface areas of macromolecules. Program SurfRace also calculates the areas of cavities inaccessible from the outside. We introduce the definition of average curvature of molecular surface and calculate average molecular surface curvatures for each atom in a structure. All surface area and curvature calculations are analytic and therefore yield exact values of these quantities. High calculation speed of this software is achieved primarily by avoiding computationally expensive mathematical procedures wherever possible and by efficient handling of surface data structures. The programs are written initially in the language C for PCs running Windows 2000/98/NT, but their code is portable to other platforms with only minor changes in input‐output procedures. The algorithm is robust and does not ignore either multiplicity or degeneracy of atomic overlaps. Fast, memory‐efficient and robust execution make this software attractive for applications both in computationally expensive energy minimization algorithms, such as docking or molecular dynamics simulations, and in stand‐alone surface area and curvature calculations. © 2002 Wiley Periodicals, Inc. J Comput Chem 23: 600–609, 2002

[1]  K S Wilson,et al.  Crystal structure of thermitase at 1.4 A resolution. , 1990, Journal of molecular biology.

[2]  N. Allewell,et al.  Functionally linked hydration changes in Escherichia coli aspartate transcarbamylase and its catalytic subunit. , 1997, Biochemistry.

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

[4]  H. Bartunik,et al.  Accuracy and precision in protein structure analysis: restrained least-squares refinement of the structure of poplar plastocyanin at 1.33 A resolution. , 1992, Acta crystallographica. Section B, Structural science.

[5]  M Levitt,et al.  Quantification of the hydrophobic interaction by simulations of the aggregation of small hydrophobic solutes in water , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[6]  Frederic M. Richards,et al.  Packing of α-helices: Geometrical constraints and contact areas☆ , 1978 .

[7]  Jacopo Tomasi,et al.  Free energy and entropy for inserting cavities in water: Comparison of Monte Carlo simulation and scaled particle theory results , 1997 .

[8]  H Edelsbrunner,et al.  Analytical shape computation of macromolecules: II. Inaccessible cavities in proteins , 1998, Proteins.

[9]  D. M. Blow,et al.  Structure of crystalline -chymotrypsin. V. The atomic structure of tosyl- -chymotrypsin at 2 A resolution. , 1972, Journal of molecular biology.

[10]  H. Meirovitch,et al.  Optimization of solvation models for predicting the structure of surface loops in proteins , 2001, Proteins.

[11]  E. Baker,et al.  Crystallographic refinement of the structure of actinidin at 1.7 Å resolution by fast Fourier least‐squares methods , 1980 .

[12]  Benno P. Schoenborn,et al.  Neutron diffraction reveals oxygen–histidine hydrogen bond in oxymyoglobin , 1981, Nature.

[13]  W. Steigemann,et al.  Structure of erythrocruorin in different ligand states refined at 1.4 A resolution. , 1979, Journal of molecular biology.

[14]  S Subramaniam,et al.  Analytical shape computation of macromolecules: I. molecular area and volume through alpha shape , 1998, Proteins.

[15]  R. Huber,et al.  The Geometry of the Reactive Site and of the Peptide Groups in Trypsin, Trypsinogen and its Complexes with Inhibitors , 1983 .

[16]  R. Dickerson,et al.  Redox conformation changes in refined tuna cytochrome c. , 1980, Proceedings of the National Academy of Sciences of the United States of America.

[17]  W. Royer,et al.  High-resolution crystallographic analysis of a co-operative dimeric hemoglobin. , 1994, Journal of molecular biology.

[18]  G. Brayer,et al.  High-resolution refinement of yeast iso-1-cytochrome c and comparisons with other eukaryotic cytochromes c. , 1990, Journal of molecular biology.

[19]  S V Evans,et al.  Refinement of recombinant oncomodulin at 1.30 A resolution. , 1993, Journal of molecular biology.

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

[21]  M. Record,et al.  Vapor pressure osmometry studies of osmolyte-protein interactions: implications for the action of osmoprotectants in vivo and for the interpretation of "osmotic stress" experiments in vitro. , 2000, Biochemistry.

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

[23]  D Tsernoglou,et al.  Structure and function of snake venom curarimimetic neurotoxins. , 1981, Molecular pharmacology.

[24]  Djordje Musil,et al.  The high-resolution X-ray crystal structure of the complex formed between subtilisin Carlsberg and eglin c, an elastase inhibitor from the leech Hirudo medicinalis Structural analysis, subtilisin structure and interface geometry , 1987 .

[25]  J C Fontecilla-Camps,et al.  Structure of variant-3 scorpion neurotoxin from Centruroides sculpturatus Ewing, refined at 1.8 A resolution. , 1983, Journal of molecular biology.

[26]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[27]  R. S. Spolar,et al.  Coupling of local folding to site-specific binding of proteins to DNA. , 1994, Science.

[28]  H. Piaggio.  Differential Geometry of Curves and Surfaces , 1952, Nature.

[29]  W. Kabsch,et al.  Refined crystal structure of the triphosphate conformation of H‐ras p21 at 1.35 A resolution: implications for the mechanism of GTP hydrolysis. , 1990, The EMBO journal.

[30]  D S Moss,et al.  Segmented anisotropic refinement of bovine ribonuclease A by the application of the rigid-body TLS model. , 1989, Acta crystallographica. Section A, Foundations of crystallography.

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

[32]  S E Ealick,et al.  Structure of scorpion toxin variant-3 at 1.2 A resolution. , 1992, Journal of molecular biology.

[33]  Ken A. Dill,et al.  The Mechanism of Hydrophobic Solvation Depends on Solute Radius , 2000 .

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

[35]  R. S. Spolar,et al.  Contribution to the thermodynamics of protein folding from the reduction in water-accessible nonpolar surface area. , 1991, Biochemistry.

[36]  M Hirose,et al.  The primary structure and structural characteristics of Achromobacter lyticus protease I, a lysine-specific serine protease. , 1993, The Journal of biological chemistry.

[37]  Robert Huber,et al.  On the disordered activation domain in trypsinogen: chemical labelling and low‐temperature crystallography , 1982 .

[38]  David S. Moss,et al.  The refined structure of ribonuclease-A at 1.45 Å resolution , 1984 .

[39]  B Honig,et al.  Internal cavities and buried waters in globular proteins. , 1986, Biochemistry.

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

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

[42]  M Gerstein,et al.  Packing at the protein-water interface. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[43]  M J Sternberg,et al.  A continuum model for protein-protein interactions: application to the docking problem. , 1995, Journal of molecular biology.

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

[45]  F. Richards,et al.  Crystallographic structures of ribonuclease S variants with nonpolar substitution at position 13: packing and cavities. , 1993, Biochemistry.

[46]  R Diamond,et al.  Real-space refinement of the structure of hen egg-white lysozyme. , 1977, Journal of molecular biology.

[47]  J C Sacchettini,et al.  Escherichia coli-derived rat intestinal fatty acid binding protein with bound myristate at 1.5 A resolution and I-FABPArg106-->Gln with bound oleate at 1.74 A resolution. , 1994, The Journal of biological chemistry.