Conformational Sampling of Druglike Molecules with MOE and Catalyst: Implications for Pharmacophore Modeling and Virtual Screening

Computational conformational sampling is integral to small molecule pharmaceutical research, for detailed conformational analysis and high-throughput 3D library enumeration. These two regimes were tested in details for the general-purpose modeling program MOE, using its three conformational sampling methods, i.e. systematic search, stochastic search, and Conformation Import. The tests include i) identification of the global energy minimum, ii) reproduction of the bioactive conformation, iii) measures of conformational coverage with 3D descriptors, and iv) compute times. The bioactive conformers are from a new set of 256 diverse, druglike, protein-bound ligands compiled and analyzed with particular care. The MOE results are compared to those obtained from the established program Catalyst. Key parameters controlling the conformational coverage were varied systematically. Coverage and diversity of the conformational space were characterized with unique pharmacophore triplets or quadruplets. Overall, the protocols in both MOE and Catalyst performed well for their intended tasks. MOE performed at least as well as Catalyst for high-throughput library generation and detailed conformational modeling. This work provides a guide and specific recommendations regarding the usage of conformational sampling tools in MOE.

[1]  Paul M. Selzer,et al.  The Impact of Tautomer Forms on Pharmacophore-Based Virtual Screening , 2006, J. Chem. Inf. Model..

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

[3]  Andrew Smellie,et al.  Poling: Promoting conformational variation , 1995, J. Comput. Chem..

[4]  Tad Hurst,et al.  Flexible 3D searching: The directed tweak technique , 1994, J. Chem. Inf. Comput. Sci..

[5]  F. Momany,et al.  Validation of the general purpose QUANTA ®3.2/CHARMm® force field , 1992 .

[6]  F. Lombardo,et al.  Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. , 2001, Advanced drug delivery reviews.

[7]  Jonas Boström,et al.  Conformational energy penalties of protein-bound ligands , 1998, J. Comput. Aided Mol. Des..

[8]  Eric J. Martin,et al.  Conformational Sampling of Bioactive Molecules: A Comparative Study , 2007, J. Chem. Inf. Model..

[9]  Robert P. Sheridan,et al.  Using CONCORD to construct a large database of three-dimensional coordinates from connection tables , 1989, J. Chem. Inf. Comput. Sci..

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

[11]  J. Gasteiger,et al.  Automatic generation of 3D-atomic coordinates for organic molecules , 1990 .

[12]  C. Lemmen,et al.  FLEXS: a method for fast flexible ligand superposition. , 1998, Journal of medicinal chemistry.

[13]  D. E. Clark,et al.  Flexible docking using tabu search and an empirical estimate of binding affinity , 1998, Proteins.

[14]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[15]  Malcolm J. McGregor,et al.  Clustering of Large Databases of Compounds: Using the MDL "Keys" as Structural Descriptors , 1997, J. Chem. Inf. Comput. Sci..

[16]  Gordon M. Crippen,et al.  Note rapid calculation of coordinates from distance matrices , 1978 .

[17]  D. Ferguson,et al.  A new approach to probing conformational space with molecular mechanics: Random incremental pulse search , 1989 .

[18]  Gerhard Klebe,et al.  A fast and efficient method to generate biologically relevant conformations , 1994, J. Comput. Aided Mol. Des..

[19]  M C Nicklaus,et al.  Conformational changes of small molecules binding to proteins. , 1995, Bioorganic & medicinal chemistry.

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

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

[22]  Tudor I. Oprea,et al.  Is There a Difference between Leads and Drugs? A Historical Perspective , 2001, J. Chem. Inf. Comput. Sci..

[23]  Gareth Jones,et al.  A genetic algorithm for flexible molecular overlay and pharmacophore elucidation , 1995, J. Comput. Aided Mol. Des..

[24]  Thierry Langer,et al.  Comparative Performance Assessment of the Conformational Model Generators Omega and Catalyst: A Large-Scale Survey on the Retrieval of Protein-Bound Ligand Conformations , 2006, J. Chem. Inf. Model..

[25]  P. Charifson,et al.  Conformational analysis of drug-like molecules bound to proteins: an extensive study of ligand reorganization upon binding. , 2004, Journal of medicinal chemistry.

[26]  R. Cramer,et al.  Validation of the general purpose tripos 5.2 force field , 1989 .

[27]  Thomas A. Halgren Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94 , 1996, J. Comput. Chem..

[28]  P. Labute,et al.  Flexible alignment of small molecules. , 2001, Journal of medicinal chemistry.

[29]  Jonas Boström,et al.  Assessing the performance of OMEGA with respect to retrieving bioactive conformations. , 2003, Journal of molecular graphics & modelling.

[30]  Paul N. Mortenson,et al.  Diverse, high-quality test set for the validation of protein-ligand docking performance. , 2007, Journal of medicinal chemistry.

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

[32]  Andrew R. Leach,et al.  Molecular Complexity and Its Impact on the Probability of Finding Leads for Drug Discovery , 2001, J. Chem. Inf. Comput. Sci..

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

[34]  D. Case,et al.  Modification of the Generalized Born Model Suitable for Macromolecules , 2000 .

[35]  C L Brooks,et al.  Do active site conformations of small ligands correspond to low free-energy solution structures? , 1998, Journal of computer-aided molecular design.

[36]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[37]  Omar Haq,et al.  Torsion Angle Preference and Energetics of Small-Molecule Ligands Bound to Proteins , 2007, J. Chem. Inf. Model..

[38]  Evan Bolton,et al.  Assessment of Conformational Ensemble Sizes Necessary for Specific Resolutions of Coverage of Conformational Space , 2007, J. Chem. Inf. Model..

[39]  Christian Senger,et al.  Representation of target-bound drugs by computed conformers: implications for conformational libraries , 2006, BMC Bioinformatics.

[40]  Thierry Langer,et al.  Comparative Analysis of Protein-Bound Ligand Conformations with Respect to Catalyst's Conformational Space Subsampling Algorithms , 2005, J. Chem. Inf. Model..

[41]  Kenneth M. Merz,et al.  Can we separate active from inactive conformations? , 2002, J. Comput. Aided Mol. Des..

[42]  Nicolas Foloppe,et al.  Toward a full characterization of nucleic acid components in aqueous solution: simulations of nucleosides. , 2005, The journal of physical chemistry. B.

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

[44]  Steven L. Teig,et al.  Chemical Function Queries for 3D Database Search , 1994, J. Chem. Inf. Comput. Sci..

[45]  R. Cramer,et al.  Comparative molecular field analysis (CoMFA). 1. Effect of shape on binding of steroids to carrier proteins. , 1988, Journal of the American Chemical Society.

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

[47]  A. Leach A Survey of Methods for Searching the Conformational Space of Small and Medium‐Sized Molecules , 1992 .