A parallel build-up algorithm for global energy minimizations of molecular clusters using effective energy simulated annealing

This work studies the build-up method for the global minimization problem for molecular conformation, especially protein folding. The problem is hard to solve for large molecules using general minimization approaches because of the enormous amount of required computation. We therefore propose a build-up process to systematically “construct” the optimal molecular structures. A prototype algorithm is designed using the anisotropic effective energy simulated annealing method at each build-up stage. The algorithm has been implemented on the Intel iPSC/860 parallel computer, and tested with the Lennard-Jones microcluster conformation problem. The experiments showed that the algorithm was effective for relatively large test problems, and also very suitable for massively parallel computation. In particular, for the 72-atom Lennard-Jones microcluster, the algorithm found a structure whose energy is lower than any others found in previous studies.

[1]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[2]  Shang‐keng Ma Introduction to the renormalization group , 1973 .

[3]  K. Wilson The renormalization group: Critical phenomena and the Kondo problem , 1975 .

[4]  W. Murray Numerical Methods for Unconstrained Optimization , 1975 .

[5]  Philip E. Gill,et al.  Practical optimization , 1981 .

[6]  M. Pincus,et al.  Calculation of the three-dimensional structure of the membrane-bound portion of melittin from its amino acid sequence. , 1982, Proceedings of the National Academy of Sciences of the United States of America.

[7]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[8]  John E. Dennis,et al.  Numerical methods for unconstrained optimization and nonlinear equations , 1983, Prentice Hall series in computational mathematics.

[9]  Steven G. Louie,et al.  A Monte carlo simulated annealing approach to optimization over continuous variables , 1984 .

[10]  J. Northby Structure and binding of Lennard‐Jones clusters: 13≤N≤147 , 1987 .

[11]  L. Wille Minimum-energy configurations of atomic clusters: new results obtained by simulated annealing , 1987 .

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

[13]  M. Karplus,et al.  Proteins: A Theoretical Perspective of Dynamics, Structure, and Thermodynamics , 1988 .

[14]  T. Hayden,et al.  Prediction of DNA structure from sequence: A build‐up technique , 1989, Biopolymers.

[15]  H. Scheraga,et al.  On the multiple-minima problem in the conformational analysis of molecules: deformation of the potential energy hypersurface by the diffusion equation method , 1989 .

[16]  B. Hendrickson The Molecular Problem: Determining Conformation from Pairwise Distances , 1990 .

[17]  T. Ackermann C. L. Brooks III, M. Karplus, B. M. Pettitt. Proteins: A Theoretical Perspective of Dynamics, Structure and Thermodynamics, Volume LXXI, in: Advances in Chemical Physics, John Wiley & Sons, New York 1988. 259 Seiten, Preis: US $ 65.25 , 1990 .

[18]  H. Scheraga,et al.  Performance of the diffusion equation method in searches for optimum structures of clusters of Lennard-Jones atoms , 1991 .

[19]  Richard H. Byrd,et al.  Parallel global optimization: numerical methods, dynamic scheduling methods, and application to molecular configuration , 1994 .

[20]  Jae Kwang Shin,et al.  High directional Monte Carlo procedure coupled with the temperature heating and annealing as a method to obtain the global energy minimum structure of polypeptides and proteins , 1991, Biopolymers.

[21]  David Shalloway,et al.  Packet annealing: a deterministic method for global minimization , 1992 .

[22]  David Shalloway,et al.  Application of the renormalization group to deterministic global minimization of molecular conformation energy functions , 1992, J. Glob. Optim..

[23]  Thomas F. Coleman,et al.  Isotropic effective energy simulated annealing searches for low energy molecular cluster states , 1993, Comput. Optim. Appl..

[24]  Guoliang Xue Improvement on the northby algorithm for molecular conformation: Better solutions , 1994, J. Glob. Optim..