Determining minimum energy conformations of polypeptides by dynamic programming

A combinatorial optimization approach is used for solving the multiple‐minima problem when determining the low‐energy conformations of short polypeptides. Each residue is represented by a finite number of discrete states corresponding to single residue local minima of the energy function. These precomputed values constitute a search table and define the conformational space for discrete minimization by a generalized dynamic programming algorithm that significantly limits the number of intermediate conformations to be generated during the search. Since dynamic programming involves stagewise decisions, it results in buildup‐type procedures implemented in two different forms. The first procedure predicts a number of conformations by a completely discrete search and these are subsequently refined by local minimization. The second involves limited continuous local minimization within the combinatorial algorithm, generally restricted to two dihedral angles in a buildup step. Both procedures are tested on 17 short peptides previously studied by other global minimization methods but involving the same potential energy function. The discrete method is extremely fast, but proves to be successful only in 14 of the 17 test problems. The version with limited local minimization finds, however, conformations in all the 17 examples that are close to the ones previously presented in the literature or have lower energies. In addition, resutls are almost independent of the cutoff energy, the most important parameter governing the search. Although the limited local minimization increases the number of energy evaluations, the method still offers substantial advantages in speed.

[1]  H. Bremermann.  A method of unconstrained global optimization , 1970 .

[2]  A. Nijenhuis.  Combinatorial algorithms , 1975 .

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

[4]  H. Scheraga,et al.  Investigation of the conformations of four tetrapeptides by nuclear magnetic resonance and circular dichroism spectroscopy, and conformational energy calculations. , 1975, Macromolecules.

[5]  H. Scheraga,et al.  Conformational analysis of the 20 naturally occurring amino acid residues using ECEPP. , 1977, Macromolecules.

[6]  H. Scheraga,et al.  Enkephalin: conformational analysis by means of empirical energy calculations. , 1977, Proceedings of the National Academy of Sciences of the United States of America.

[7]  H. Scheraga,et al.  Conformational Energy Calculations on the Contraceptive Tetrapeptide H-Thr-Pro-Arg-Lys-OH , 1978 .

[8]  Harold A. Scheraga,et al.  Conformational Energy Calculations of the Effects of Sequence Variations on the Conformations of Two Tetrapeptides , 1978 .

[9]  H. Scheraga,et al.  Conformational Energy Study of Tuftsin , 1978 .

[10]  Ellis Horowitz,et al.  Fundamentals of Computer Algorithms , 1978 .

[11]  H. Zimmermann.  Towards global optimization 2: L.C.W. DIXON and G.P. SZEGÖ (eds.) North-Holland, Amsterdam, 1978, viii + 364 pages, US $ 44.50, Dfl. 100,-. , 1979 .

[12]  G. Bekey,et al.  A global optimization algorithm using adaptive random search , 1980 .

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

[14]  D. Eisenberg,et al.  The structure of melittin in the form I crystals and its implication for melittin's lytic and surface activities. , 1982, Biophysical journal.

[15]  H. Scheraga,et al.  Computed conformational states of the 20 naturally occurring amino acid residues and of the prototype residue α-aminobutyric acid , 1983 .

[16]  H. Scheraga,et al.  Energy parameters in polypeptides. 9. Updating of geometrical parameters, nonbonded interactions, and hydrogen bond interactions for the naturally occurring amino acids , 1983 .

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

[18]  M S Waterman,et al.  Sequence alignments in the neighborhood of the optimum with general application to dynamic programming. , 1983, Proceedings of the National Academy of Sciences of the United States of America.

[19]  Eric Walter,et al.  A general-purpose global optimizer: implementation and applications , 1984 .

[20]  H A Scheraga,et al.  β‐Bend conformation of CH3CO‐Pro‐Pro‐Gly‐Pro‐NHCH3: Implications for posttranslational proline hydroxylation in collagen , 1984, Biopolymers.

[21]  A. V. Levy,et al.  The Tunneling Algorithm for the Global Minimization of Functions , 1985 .

[22]  H. Scheraga,et al.  Prediction of the native conformation of a polypeptide by a statistical‐mechanical procedure. I. Backbone structure of enkephalin , 1985, Biopolymers.

[23]  H. Scheraga,et al.  Use of buildup and energy‐minimization procedures to compute low‐energy structures of the backbone of enkephalin , 1985, Biopolymers.

[24]  J. Moult,et al.  An algorithm for determining the conformation of polypeptide segments in proteins by systematic search , 1986, Proteins.

[25]  H A Scheraga,et al.  An approach to the multiple-minima problem by relaxing dimensionality. , 1986, Proceedings of the National Academy of Sciences of the United States of America.

[26]  M. E. Johnson,et al.  Generalized simulated annealing for function optimization , 1986 .

[27]  P. Schiller,et al.  Theoretical conformational analysis of a μ‐selective cyclic opioid peptide analog , 1987, Biopolymers.

[28]  C. Levinthal,et al.  Predicting antibody hypervariable loop conformation. I. Ensembles of random conformations for ringlike structures , 1987, Biopolymers.

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

[30]  H. Scheraga,et al.  Revised algorithms for the build‐up procedure for predicting protein conformations by energy minimization , 1987 .

[31]  H. Scheraga,et al.  Proline‐induced constraints in α‐helices , 1987, Biopolymers.

[32]  H A Scheraga,et al.  Prediction of the native conformation of a polypeptide by a statistical‐mechanical procedure. III. Probable and average conformations of enkephalin , 1987, Biopolymers.

[33]  M. Karplus,et al.  Prediction of the folding of short polypeptide segments by uniform conformational sampling , 1987, Biopolymers.

[34]  H A Scheraga,et al.  An approach to the multiple-minima problem in protein folding by relaxing dimensionality. Tests on enkephalin. , 1987, Journal of molecular biology.

[35]  H. Scheraga,et al.  Low‐energy conformations of two lysine‐containing tetrapeptides of collagen: Implications for posttranslational lysine hydroxylation , 1987, Biopolymers.

[36]  W. C. Still,et al.  The multiple minimum problem in molecular modeling. Tree searching internal coordinate conformational space , 1988 .

[37]  H. Scheraga,et al.  Calculation of protein conformation by the build-up procedure. Application to bovine pancreatic trypsin inhibitor using limited simulated nuclear magnetic resonance data. , 1988, Journal of biomolecular structure & dynamics.

[38]  Stephen R. Wilson,et al.  Conformational Analysis of Flexible Molecules: Location of the Global Minimum Energy Conformation by the Simulated Annealing Method , 1988 .

[39]  H. Mizuno,et al.  Prediction of the conformation of ice‐nucleation protein by conformational energy calculation , 1989, Proteins.

[40]  William H. Press,et al.  Numerical recipes , 1990 .