An efficient algorithm for searching low-energy conformers of cyclic and acyclic molecules

A set of strategies for exhaustively finding low-energy conformers of cyclic, acyclic or alicyclic molecules is presented. Starting from any conformation, local perturbation is systematically applied to all the flexible portions of the molecule in question to produce candidates of new conformation. The perturbations consist of flapping and/or flipping for endocyclic bonds and stepwise rotation for acyclic bonds. The conformations they produced are believed to lie close to the initial geometry in the conformational space. The global energy minimum (GEM) structure of the starting domain of conformational space can be quickly reached by always choosing the most stable of the conformers produced in the last perturbation cycle as the next initial structure. Once GEM of the domain is reached, the local perturbations direct the search gradually to higher and higher energy regions while exhaustively finding all the low-energy conformers therein. The variable search-limit strategy allows one to use unstable conformers as the initial structure for perturbation to ensure the exhaustiveness of the search in the low-energy region. By further increasing the search-limit, new domains of conformational space may be found. A program CONFLEX3 containing several additional strategies for improved performance has been tested for n-alkanes up to decane and cycloalkanes up to cyclododecane.

[1]  C. Swahn,et al.  Exploratory Calculations of Medium and Large Rings. Part 1. Conformational Minima of Cycloalkanes. , 1973 .

[2]  H. Hosoya,et al.  Investigation of intramolecular interactions in n-alkanes. Cooperative energy increments associated with GG and GTG' [G = gauche, T = trans] sequences , 1991 .

[3]  Robert S. H. Liu,et al.  A bioorganic view of the chemistry of vision: H.T.-n and B.P.-m,n mechanisms for reactions of confined, anchored polyenes , 1986 .

[4]  S. Ōmura,et al.  (+)-Trienomycins A, B, and C : relative and absolute stereochemistry , 1990 .

[5]  E. Corey,et al.  Studies on the conformational mobility of arachidonic acid. Facile macrolactonization of 20-hydroxyarachidonic acid. , 1983 .

[6]  Martin Saunders,et al.  Searching for conformers of nine‐ to twelve‐ring hydrocarbons on the MM2 and MM3 energy surfaces: Stochastic search for interconversion pathways , 1991 .

[7]  Peter A. Kollman,et al.  Conformational sampling and ensemble generation by molecular dynamics simulations: 18‐Crown‐6 as a test case , 1992 .

[8]  D J Barlow,et al.  RAMBLE: a conformational search program. , 1990, Journal of molecular graphics.

[9]  Scott H. Northrup,et al.  Molecular mechanics of kink formation in lipid monolayers , 1984 .

[10]  J. G. Vinter,et al.  DMC: A multifunctional hybrid dynamics/Monte Carlo simulation algorithm for the evaluation of conformational space , 1992 .

[11]  N. Nakajima,et al.  Stereoselective Total Synthesis of 16-Membered Macrolide Aglycons, Leuconolides and Maridonolides. Macrocylic Stereocontrol Based on Conformational Analysis of the 16-Membered Macrolide Ring , 1991 .

[12]  P. Kollman,et al.  Comparison of ab initio, semiempirical, and molecular mechanics calculations for the conformational analysis of ring systems , 1992 .

[13]  Carol A. Venanzi,et al.  A comparison of distance geometry and molecular dynamics simulation techniques for conformational analysis of β‐cyclodextrin , 1992 .

[14]  Peter A. Kollman,et al.  Computer modeling of the interactions of complex molecules , 1990 .

[15]  W. Clark Still,et al.  An unbounded systematic search of conformational space , 1991 .

[16]  G. Chang,et al.  An internal-coordinate Monte Carlo method for searching conformational space , 1989 .

[17]  E. Ōsawa,et al.  Application of molecular mechanics to natural product chemistry , 1989 .

[18]  E. Ōsawa,et al.  On the possibility of determining stereochemistry in acyclic polyhydroxylated compounds by the combined vicinal coupling constant/molecular mechanics method. A test with alditol peracetates , 1991 .

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

[20]  I. Shamovsky,et al.  Comparison of global-minimum-energy conformations of phospholipids , 1991 .

[21]  Eiji Ōsawa,et al.  Further developments in the algorithm for generating cyclic conformers. Test with cycloheptadecane , 1992 .

[22]  Catherine E. Peishoff,et al.  Improvements to the distance geometry algorithm for conformational sampling of cyclic structures , 1992 .

[23]  Jay W. Ponder,et al.  Stereochemistry of the hygrolidins , 1984 .

[24]  Gordon M. Crippen,et al.  Exploring the conformation space of cycloalkanes by linearized embedding , 1992 .

[25]  A E Howard,et al.  An analysis of current methodologies for conformational searching of complex molecules. , 1988, Journal of medicinal chemistry.

[26]  Yi-der Chen,et al.  Monte Carlo study of freely jointed ring polymers. I. Generation of ring polymers by dimerization method , 1981 .

[27]  H. Goto.  A revised nomenclature for the ring conformation and a note on the conformational distance in cyclododecane , 1992 .

[28]  E. Clementi,et al.  Global scientific and engineering simulations on scalar, vector and parallel LCAP-type supercomputers , 1988, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[29]  Harold A. Scheraga,et al.  Standard‐geometry chains fitted to X‐ray derived structures: Validation of the rigid‐geometry approximation. II. Systematic searches for short loops in proteins: Applications to bovine pancreatic ribonuclease A and human lysozyme , 1992 .

[30]  R. Elber,et al.  Reaction path study of conformational transitions in flexible systems: Applications to peptides , 1990 .

[31]  J. Uiterwijk,et al.  The number of ideal rings on the diamond lattice; application to crown ethers , 1983 .

[32]  W. Guida,et al.  Probing the conformational space available to inhibitors in the thermolysin active site using Monte Carlo/energy minimization techniques , 1992 .

[33]  Eiji Osawa,et al.  Corner flapping: a simple and fast algorithm for exhaustive generation of ring conformations , 1989 .

[34]  Martin Saunders,et al.  Conformations of cycloheptadecane. A comparison of methods for conformational searching , 1990 .

[35]  E. Clementi,et al.  Supercomputing and supercomputers for science and engineering in general and for chemistry and biosciences in particular , 1989 .

[36]  D. Ferguson,et al.  Molecular mechanics conformational analysis of cyclononane using the RIPS method and comparison with quantum‐mechanical calculations , 1989 .

[37]  M. Volkenstein,et al.  Statistical mechanics of chain molecules , 1969 .

[38]  Werner Braun,et al.  Efficient search for all low energy conformations of polypeptides by Monte Carlo methods , 1991 .

[39]  K. Hassler.  Zum Redundanzproblem substituierter Ringe , 1979 .

[40]  M. Snow.  Powerful simulated‐annealing algorithm locates global minimum of protein‐folding potentials from multiple starting conformations , 1992 .

[41]  G. Chang,et al.  Macromodel—an integrated software system for modeling organic and bioorganic molecules using molecular mechanics , 1990 .

[42]  J. Snyder,et al.  Syntheses of Macrocyclic Enzyme Models. 8. ∠ Conformational Mobility and Molecular Recognition by the Internal Cage of Kyuphane , 1991 .

[43]  E. Burlakova,et al.  Influence of fatty acid composition on the structure and stability of fatty acid complexes with vitamin E , 1992 .

[44]  P. Camps,et al.  Application of empirical potential energy calculations to organic chemistry. Part 19. Conformational preference in 2,4-dimethoxybicyclo[3.3.1]nonan-9-one and related molecules. Analysis of vicinal NMR coupling constants in multiple rotor system by combined molecular mechanics and generalized Karplus , 1983 .

[45]  Applications of the simulated annealing method to intermolecular interactions , 1991 .