Efficient algorithms for protein solvent accessible surface area

We present faster sequential and parallel algorithms for computing the solvent accessible surface area (ASA) of protein molecules. The ASA can be obtained by calculating the exposed surface area of the spheres obtained by increasing the van der Waals' radii of the atoms with the van der Waals' radius of the solvent. Using domain specific knowledge, we show that the number of sphere intersections is O(n) and present algorithms to compute the same in O(nlogn) sequential time and O(nlogn/p) parallel time, where n is the number of atoms and p is the number of processors. We also present a heuristic based on space-filling curves to improve performance in practice. These are significant improvements over previously known algorithms which take /spl Omega/(n/sup 2/) time sequentially and /spl Omega/(n/sup 2//p) time in parallel. While existing parallel algorithms achieve their run-time by dynamic load balancing, our algorithms are faster and do not need load balancing.

[1]  R. V. Shankar,et al.  Many-to-many Personalized Communication with Bounded Traac 1 Basic Communication Primitives , 2007 .

[2]  C A Johnson,et al.  Parallel computing in biomedical research. , 1994, Science.

[3]  Alok N. Choudhary,et al.  Parallel computation of solvent accessible surface area of protein molecules , 1993, [1993] Proceedings Seventh International Parallel Processing Symposium.

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

[5]  J. H. Lindsey,et al.  Sphere packing in R 3 , 1986 .

[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]  Robert L. Martino,et al.  Parallel Algorithms in Molecular Biology , 1997, HPCN Europe.

[8]  F. Richards,et al.  The protein folding problem. , 1991, Scientific American.

[9]  Sanjay Ranka,et al.  Many-to-many personalized communication with bounded traffic , 1995, Proceedings Frontiers '95. The Fifth Symposium on the Frontiers of Massively Parallel Computation.

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

[11]  Bhagirath Narahari,et al.  Dynamic load balancing schemes for computing accessible surface area of protein molecules , 1998, Proceedings. Fifth International Conference on High Performance Computing (Cat. No. 98EX238).