GAPE: An Improved Genetic Algorithm for Pharmacophore Elucidation

Prediction of the binding mode for a series of active compounds, in the absence of known protein structure, is a problem of paramount importance in rational drug design. GAPE (genetic algorithm for pharmacophore elucidation) is an automated multicompound overlay creation program, based on the original GASP program, that uses a genetic algorithm to fully explore the conformational space of the input structures and their alignments, so as to elucidate a pharmacophore. The software was evaluated on 13 test systems from nine protein targets using overlaid ligands extracted from the PDB. Using objective rmsd criteria and starting from 2D structures, in the absence of any protein information, GAPE was observed in eight systems to approximate the crystallographically observed binding mode. In the predicted alignments for each of those eight systems, at least half the input structures were within 2 Å rmsd of the crystal structure coordinates. Further analysis, using stricter subjective criteria, showed considerable success in five systems. For example, the prediction for a set of 12 ligands targeting P38 had 11 ligands with a 1.8 Å rmsd to crystal structure coordinates. Finally, the algorithm was favorably compared with the current GASP and Galahad programs.

[1]  Andrew Smellie,et al.  Identification of Common Functional Configurations Among Molecules , 1996, J. Chem. Inf. Comput. Sci..

[2]  Antonio Zamora,et al.  An Algorithm for Finding the Smallest Set of Smallest Rings , 1976, J. Chem. Inf. Comput. Sci..

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

[4]  R C Glen,et al.  Molecular recognition using a binary genetic search algorithm. , 1993, Journal of molecular graphics.

[5]  Christian Lemmen,et al.  Computational methods for the structural alignment of molecules , 2000, J. Comput. Aided Mol. Des..

[6]  Gareth Jones,et al.  Further Development of a Genetic Algorithm for Ligand Docking and Its Application to Screening Combinatorial Libraries , 1999 .

[7]  Yvonne C. Martin,et al.  A fast new approach to pharmacophore mapping and its application to dopaminergic and benzodiazepine agonists , 1993, J. Comput. Aided Mol. Des..

[8]  Thierry Langer,et al.  Comparative Analysis of Protein‐Bound Ligand Conformations with Respect to Catalyst′s Conformational Space Subsampling Algorithms. , 2005 .

[9]  Mark A. Murcko,et al.  Virtual screening : an overview , 1998 .

[10]  Sung Jin Cho,et al.  FLAME: A Program to Flexibly Align Molecules. , 2006 .

[11]  Philip M. Dean,et al.  Three-dimensional hydrogen-bond geometry and probability information from a crystal survey , 1996, J. Comput. Aided Mol. Des..

[12]  J. A. Grant,et al.  A fast method of molecular shape comparison: A simple application of a Gaussian description of molecular shape , 1996, J. Comput. Chem..

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

[14]  Peter Willett,et al.  GALAHAD: 1. Pharmacophore identification by hypermolecular alignment of ligands in 3D , 2006, J. Comput. Aided Mol. Des..

[15]  Valerie J. Gillet,et al.  Incorporating partial matches within multiobjective pharmacophore identification , 2006, J. Comput. Aided Mol. Des..

[16]  P Willett,et al.  Development and validation of a genetic algorithm for flexible docking. , 1997, Journal of molecular biology.

[17]  Robert D. Clark,et al.  SYBYL Line Notation (SLN): A Single Notation To Represent Chemical Structures, Queries, Reactions, and Virtual Libraries , 2008, J. Chem. Inf. Model..

[18]  Michal Vieth,et al.  Geometric Accuracy of Three-Dimensional Molecular Overlays , 2006, J. Chem. Inf. Model..

[19]  Ralph Kühne,et al.  Prediction of the Intrinsic Hydrogen Bond Acceptor Strength of Organic Compounds by Local Molecular Parameters , 2009, J. Chem. Inf. Model..

[20]  Dagmar Ringe,et al.  Analogous inhibitors of elastase do not always bind analogously , 1994, Nature Structural Biology.

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

[22]  Sun Choi,et al.  The Effects of Biasing Torsional Mutations in a Conformational GA , 2006, J. Chem. Inf. Model..

[23]  Valerie J. Gillet,et al.  Generation of multiple pharmacophore hypotheses using multiobjective optimisation techniques , 2004, J. Comput. Aided Mol. Des..

[24]  Jun Feng,et al.  PharmID: Pharmacophore Identification Using Gibbs Sampling , 2006, J. Chem. Inf. Model..

[25]  Andrew R. Leach,et al.  A comparison of the pharmacophore identification programs: Catalyst, DISCO and GASP , 2002, J. Comput. Aided Mol. Des..

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

[27]  L. Hardy,et al.  The impact of structure-guided drug design on clinical agents , 2003 .

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

[29]  Magdalena Bacilieri,et al.  Ligand-based drug design methodologies in drug discovery process: an overview. , 2006, Current drug discovery technologies.

[30]  Ajay N. Jain,et al.  Ligand-based structural hypotheses for virtual screening. , 2004, Journal of medicinal chemistry.

[31]  R C Glen,et al.  Characterisation of the solution conformation of a cyclic RGD peptide analogue by NMR spectroscopy allied with a genetic algorithm approach and constrained molecular dynamics. , 2009, International journal of peptide and protein research.

[32]  R. Glen,et al.  Molecular recognition of receptor sites using a genetic algorithm with a description of desolvation. , 1995, Journal of molecular biology.

[33]  Arthur Dalby,et al.  Description of several chemical structure file formats used by computer programs developed at Molecular Design Limited , 1992, J. Chem. Inf. Comput. Sci..

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

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

[36]  Valerie J. Gillet,et al.  Multiobjective Optimization of Pharmacophore Hypotheses: Bias Toward Low-Energy Conformations , 2009, J. Chem. Inf. Model..

[37]  J. A. Grant,et al.  A Gaussian Description of Molecular Shape , 1995 .

[38]  Anil K. Jain,et al.  Data clustering: a review , 1999, CSUR.

[39]  David E. Shaw,et al.  PHASE: a new engine for pharmacophore perception, 3D QSAR model development, and 3D database screening: 1. Methodology and preliminary results , 2006, J. Comput. Aided Mol. Des..

[40]  Gareth Jones,et al.  Elucidating Molecular Overlays from Pairwise Alignments Using a Genetic Algorithm , 2009, J. Chem. Inf. Model..