Generating Triangulated Macromolecular Surfaces by Euclidean Distance Transform

Macromolecular surfaces are fundamental representations of their three-dimensional geometric shape. Accurate calculation of protein surfaces is of critical importance in the protein structural and functional studies including ligand-protein docking and virtual screening. In contrast to analytical or parametric representation of macromolecular surfaces, triangulated mesh surfaces have been proved to be easy to describe, visualize and manipulate by computer programs. Here, we develop a new algorithm of EDTSurf for generating three major macromolecular surfaces of van der Waals surface, solvent-accessible surface and molecular surface, using the technique of fast Euclidean Distance Transform (EDT). The triangulated surfaces are constructed directly from volumetric solids by a Vertex-Connected Marching Cube algorithm that forms triangles from grid points. Compared to the analytical result, the relative error of the surface calculations by EDTSurf is <2–4% depending on the grid resolution, which is 1.5–4 times lower than the methods in the literature; and yet, the algorithm is faster and costs less computer memory than the comparative methods. The improvements in both accuracy and speed of the macromolecular surface determination should make EDTSurf a useful tool for the detailed study of protein docking and structure predictions. Both source code and the executable program of EDTSurf are freely available at http://zhang.bioinformatics.ku.edu/EDTSurf.

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

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

[3]  B. Bush,et al.  Macromolecular shape and surface maps by solvent exclusion. , 1978, Proceedings of the National Academy of Sciences of the United States of America.

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

[5]  M. L. Connolly Solvent-accessible surfaces of proteins and nucleic acids. , 1983, Science.

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

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

[8]  William E. Lorensen,et al.  Marching cubes: A high resolution 3D surface construction algorithm , 1987, SIGGRAPH.

[9]  Q. Ye The signed Euclidean distance transform and its applications , 1988, [1988 Proceedings] 9th International Conference on Pattern Recognition.

[10]  Luc Vincent,et al.  Exact Euclidean distance function by chain propagations , 1991, Proceedings. 1991 IEEE Computer Society Conference on Computer Vision and Pattern Recognition.

[11]  Herbert Edelsbrunner,et al.  Three-dimensional alpha shapes , 1992, VVS.

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

[13]  M L Connolly,et al.  The molecular surface package. , 1993, Journal of molecular graphics.

[14]  A. Olson,et al.  Approximation and characterization of molecular surfaces , 1993, Biopolymers.

[15]  Xiang Zhexin,et al.  Calculating the electric potential of macromolecules: A simple method for molecular surface triangulation , 1995, J. Comput. Chem..

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

[17]  S Subramaniam,et al.  Computation of molecular electrostatics with boundary element methods. , 1997, Biophysical journal.

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

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

[20]  André H. Juffer,et al.  A flexible triangulation method to describe the solvent-accessible surface of biopolymers , 1998, J. Comput. Aided Mol. Des..

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

[22]  Maria I. Zavodszky,et al.  Distilling the essential features of a protein surface for improving protein-ligand docking, scoring, and virtual screening , 2002, J. Comput. Aided Mol. Des..

[23]  Ruth Nussinov,et al.  Taking geometry to its edge: Fast unbound rigid (and hinge‐bent) docking , 2003, Proteins.

[24]  Conrad C. Huang,et al.  UCSF Chimera—A visualization system for exploratory research and analysis , 2004, J. Comput. Chem..

[25]  Shura Hayryan,et al.  A new analytical method for computing solvent‐accessible surface area of macromolecules and its gradients , 2005, J. Comput. Chem..

[26]  Yuan-Fang Wang,et al.  Efficient molecular surface generation using level-set methods. , 2006, Journal of molecular graphics & modelling.

[27]  Chandrajit L. Bajaj,et al.  Quality meshing of implicit solvation models of biomolecular structures , 2006, Comput. Aided Geom. Des..

[28]  Dong Xu,et al.  Euclidean Distance Transform of Digital Images in Arbitrary Dimensions , 2006, PCM.

[29]  Michael Petukhov,et al.  Joint neighbors approximation of macromolecular solvent accessible surface area , 2007, J. Comput. Chem..

[30]  Kengo Kinoshita,et al.  Docking of protein molecular surfaces with evolutionary trace analysis , 2007, Proteins.

[31]  Yang Zhang Progress and challenges in protein structure prediction. , 2008, Current opinion in structural biology.

[32]  Sitao Wu,et al.  MUSTER: Improving protein sequence profile–profile alignments by using multiple sources of structure information , 2008, Proteins.

[33]  Helen M Berman,et al.  Signatures of protein-DNA recognition in free DNA binding sites. , 2009, Journal of molecular biology.