Analyzing the performance of conformational search programs on compound databases.

We have studied the sampling performance of conformational search programs using geometric and energetic criteria. Ideally, a conformational search algorithm should identify the largest possible number of low-energy structures (energy criterion) covering the widest possible range of molecular shapes (geometric criterion). Geometric analysis consisted in comparing the distribution of conformations within the generated ensembles by multidimensional scaling and by analysing the eigenvalue structure of the pairwise coordinate covariance matrices. The energetic comparison was carried out by assessing the energy distribution of conformers after minimizing them all using the same semi-empirical quantum mechanics optimization protocol. The present investigation focused on five conformational search programs: DGEOM, QXP, ROTATE, LMOD and OMEGA. We have applied these methodologies to a maximally diverse 604-compound subset of the LeadQuest library. The program LMOD performs best according to the energetic criterion, whereas a wider range of geometrically diverse conformations is sampled by the other programs, at the cost of higher median conformer energies. In terms of speed, OMEGA is fastest. We recommend the use of LMOD or OMEGA for high-quality conformational search applications.

[1]  Gordon M. Crippen,et al.  Distance Geometry and Molecular Conformation , 1988 .

[2]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[3]  D C Spellmeyer,et al.  Conformational analysis using distance geometry methods. , 1997, Journal of molecular graphics & modelling.

[4]  Johann Gasteiger,et al.  Handbook of Chemoinformatics , 2003 .

[5]  Colin McMartin,et al.  QXP: Powerful, rapid computer algorithms for structure-based drug design , 1997, J. Comput. Aided Mol. Des..

[6]  David A. Case,et al.  Modeling Unusual Nucleic Acid Structures , 1998 .

[7]  P. M. Kroonenberg,et al.  Consensus Molecular Alignment Based on Generalized Procrustes Analysis , 2003, J. Chem. Inf. Comput. Sci..

[8]  P. Kollman,et al.  Automatic atom type and bond type perception in molecular mechanical calculations. , 2006, Journal of molecular graphics & modelling.

[9]  Jonas Boström,et al.  Reproducing the conformations of protein-bound ligands: A critical evaluation of several popular conformational searching tools , 2001, J. Comput. Aided Mol. Des..

[10]  Nuria B. Centeno,et al.  New insights into the multiple minima problem , 1996 .

[11]  R. Shepard.  The analysis of proximities: Multidimensional scaling with an unknown distance function. I. , 1962 .

[12]  Cheuk-San Wang,et al.  Efficient algorithm for conformational search of macrocyclic molecules , 1997, J. Comput. Chem..

[13]  A R Leach,et al.  Exploring the conformational space of protein side chains using dead‐end elimination and the A* algorithm , 1998, Proteins.

[14]  J. Gower.  Generalized procrustes analysis , 1975 .

[15]  Robert D. Clark,et al.  OptiSim: An Extended Dissimilarity Selection Method for Finding Diverse Representative Subsets , 1997, J. Chem. Inf. Comput. Sci..

[16]  Robert Stanton,et al.  Conformational analysis by intersection: CONAN , 2003, J. Comput. Chem..

[17]  Richard M. Everson,et al.  Inferring the eigenvalues of covariance matrices from limited, noisy data , 2000, IEEE Trans. Signal Process..

[18]  Emelyn Smith,et al.  A comparison of the Low Mode and Monte Carlo conformational search methods. , 2002, Journal of molecular graphics & modelling.

[19]  Christof H. Schwab,et al.  Conformational Analysis and Searching , 2008 .

[20]  R. Shepard.  The analysis of proximities: Multidimensional scaling with an unknown distance function. II , 1962 .

[21]  Michael D. Miller,et al.  Comparison of Knowledge-Based and Distance Geometry Approaches for Generation of Molecular Conformations , 2001, J. Chem. Inf. Comput. Sci..

[22]  Andrew Smellie,et al.  Analysis of Conformational Coverage, 1. Validation and Estimation of Coverage , 1995, J. Chem. Inf. Comput. Sci..

[23]  István Kolossváry,et al.  Low-mode conformational search elucidated: Application to C39H80 and flexible docking of 9-deazaguanine inhibitors into PNP , 1999, J. Comput. Chem..