Efficient parallelization of the energy, surface, and derivative calculations for internal coordinate mechanics

An efficient algorithm for parallelization of a molecular mechanics program operating in the space of internal coordinates such as dihedral angles, bond angles, and bond lengths is described. The iterative procedure to calculate analytical energy derivatives with respect to the internal coordinates was modified to allow parallelization. Computationally intensive modules that calculate energy and its derivatives, solvent‐accessible surface, electrostatic polarization energy and that update lists of interactions were parallelized with nearly 100% efficiency. The proposed strategy for the shared‐memory computer architecture is easily scalable and requires minimum changes in a program code. The overall speedup for a realistic calculation minimizing the energy of a myoglobin reaches a factor of 3 for 4 processors. © 1994 by John Wiley & Sons, Inc.

[1]  Peter A. Kollman,et al.  AMBERCUBE MD, parallelization of Amber's molecular dynamics module for distributed‐memory hypercube computers , 1993, J. Comput. Chem..

[2]  M J Sternberg,et al.  New algorithm to model protein-protein recognition based on surface complementarity. Applications to antibody-antigen docking. , 1992, Journal of molecular biology.

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

[4]  P Argos,et al.  A method to configure protein side-chains from the main-chain trace in homology modelling. , 1993, Journal of molecular biology.

[5]  S. Phillips,et al.  Structure of oxymyoglobin , 1978, Nature.

[6]  H. Scheraga,et al.  Energy parameters in polypeptides. VII. Geometric parameters, partial atomic charges, nonbonded interactions, hydrogen bond interactions, and intrinsic torsional potentials for the naturally occurring amino acids , 1975 .

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

[8]  James Andrew McCammon,et al.  Parallelization of a molecular dynamics non-bonded force algorithm for MIMD architecture , 1990, Comput. Chem..

[9]  Julia M. Goodfellow,et al.  Parallelization strategies for molecular simulation using the Monte Carlo algorithm , 1993, J. Comput. Chem..

[10]  Protein calculations on parallel processors. II. Parallel algorithm for the forces and molecular dynamics , 1992 .

[11]  Douglas J. Tobias,et al.  Vector and parallel algorithms for the molecular dynamics simulation of macromolecules on shared‐memory computers , 1991 .

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

[13]  Ruben Abagyan,et al.  ICM—A new method for protein modeling and design: Applications to docking and structure prediction from the distorted native conformation , 1994, J. Comput. Chem..

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

[15]  I. Kuntz,et al.  Protein docking and complementarity. , 1991, Journal of molecular biology.

[16]  J. W. Causey,et al.  Accelerated molecular dynamics simulation with the parallel fast multipole algorithm , 1992 .

[17]  Robert C. Schweitzer,et al.  Application of parallel processing techniques to improving the efficiency of MM2 molecular mechanics calculations , 1993, J. Comput. Chem..

[18]  Hiroshi Wako,et al.  Algorithm for rapid calculation of Hessian of conformational energy function of proteins by supercomputer , 1987 .

[19]  R. Abagyan,et al.  Biased probability Monte Carlo conformational searches and electrostatic calculations for peptides and proteins. , 1994, Journal of molecular biology.

[20]  W. Braun,et al.  Rapid calculation of first and second derivatives of conformational energy with respect to dihedral angles for proteins general recurrent equations , 1984, Comput. Chem..

[21]  R Abagyan,et al.  Design, creation, and characterization of a stable, monomeric triosephosphate isomerase. , 1994, Proceedings of the National Academy of Sciences of the United States of America.

[22]  J. F. Janak,et al.  Protein calculations on parallel processors. I. Parallel algorithm for the potential energy , 1992 .

[23]  H. Scheraga,et al.  Monte Carlo-minimization approach to the multiple-minima problem in protein folding. , 1987, Proceedings of the National Academy of Sciences of the United States of America.

[24]  W F van Gunsteren,et al.  Can the stability of protein mutants be predicted by free energy calculations? , 1993, Protein engineering.