The determination of a stable molecular structure can often be formulated in terms of calculating the global (or approximate global) minimum of a potential energy function. Computing the global minimum of this function is very difficult because it typically has a very large number of local minima which may grow exponentially with molecule size. The optimization method presented involves collecting a large number of conformers, each attained by finding a local minimum of the potential energy function from a random starting point. The information from these conformers is then used to form a convex quadratic global underestimating function for the potential energy of all known conformers. This underestimator is an L 1 approximation to all known local minima, and is obtained by a linear programming formulation and solution. The minimum of this underestimator is used to predict the global minimum for the function, allowing a localized conformer search to be performed based on the predicted minimum. The new set of conformers generated by the localized search serves as the basis for another quadratic underestimation step in an iterative algorithm. This algorithm has been used to determine the structures of homopolymers of lengthn ≤ 30 with no sidechains. While it is estimated that there areO(3n) local minima for a chain of length n, this method requires O(n4) computing time on average. It is also shown that the global minimum potential energy values lie on a concave quadratic curve for n ≤ 30. This important property permits estimation of the minimum energy for larger molecules, and also can be used to accelerate the global minimization algorithm.
[1]
A. Lehninger.
Biochemistry: The Molecular Basis of Cell Structure and Function
,
1970
.
[2]
David A. Case,et al.
Computer simulations of protein dynamics and thermodynamics
,
1993,
Computer.
[3]
K. Dill.
Dominant forces in protein folding.
,
1990,
Biochemistry.
[4]
Scott Wilson.
The Protein Folding Problem and Tertiary Structure Prediction
,
1994,
Birkhäuser Boston.
[5]
H. Bradford Thompson,et al.
Calculation of Cartesian Coordinates and Their Derivatives from Internal Molecular Coordinates
,
1967
.
[6]
R. Abagyan.
Towards protein folding by global energy optimization
,
1993,
FEBS letters.
[7]
S A Benner,et al.
Predicting the conformation of proteins man versus machine
,
1993,
FEBS letters.
[8]
L. Greengard.
The Rapid Evaluation of Potential Fields in Particle Systems
,
1988
.
[9]
Thomas G. Metzger,et al.
Conformational searches for the global minimum of protein models
,
1994,
J. Glob. Optim..
[10]
S. Sun,et al.
Reduced representation model of protein structure prediction: Statistical potential and genetic algorithms
,
1993,
Protein science : a publication of the Protein Society.