Premium
Stochastic molecular optimization using generalized simulated annealing
Author(s) -
Moret M. A.,
Pascutti P. G.,
Bisch P. M.,
Mundim K. C.
Publication year - 1998
Publication title -
journal of computational chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.907
H-Index - 188
eISSN - 1096-987X
pISSN - 0192-8651
DOI - 10.1002/(sici)1096-987x(19980430)19:6<647::aid-jcc6>3.0.co;2-r
Subject(s) - maxima and minima , simulated annealing , molecular dynamics , force field (fiction) , statistical physics , potential energy , molecule , energy minimization , computational chemistry , chemistry , computer science , physics , algorithm , mathematics , classical mechanics , quantum mechanics , mathematical analysis
We propose a stochastic optimization technique based on a generalized simulated annealing (GSA) method for mapping minima points of molecular conformational energy surfaces. The energy maps are obtained by coupling a classical molecular force field (THOR package) with a GSA procedure. Unlike the usual molecular dynamics (MD) method, the method proposed in this study is force independent; that is, we obtain the optimized conformation without calculating the force, and only potential energy is involved. Therefore, we do not need to know the conformational energy gradient to arrive at equilibrium conformations. Its utility in molecular mechanics is illustrated by applying it to examples of simple molecules (H 2 O and H 2 O 3 ) and to polypeptides. The results obtained for H 2 O and H 2 O 3 using Tsallis thermostatistics suggest that the GSA approach is faster than the other two conventional methods (Boltzmann and Cauchy machines). The results for polypeptides show that pentalanine does not form a stable α‐helix structure, probably because the number of hydrogen bonds is insufficient to maintain the helical array. On the contrary, the icoalanine molecule forms an α‐helix structure. We obtain this structure simulating all Φ, Ψ pairs using only a few steps, as compared with conventional methods. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 647–657, 1998