A genetic algorithm for flexible molecular overlay and pharmacophore elucidation

SummaryA genetic algorithm (GA) has been developed for the superimposition of sets of flexible molecules. Molecules are represented by a chromosome that encodes angles of rotation about flexible bonds and mappings between hydrogen-bond donor proton, acceptor lone pair and ring centre features in pairs of molecules. The molecule with the smallest number of features in the data set is used as a template, onto which the remaining molecules are fitted with the objective of maximising structural equivalences. The fitness function of the GA is a weighted combination of: (i) the number and the similarity of the features that have been overlaid in this way; (ii) the volume integral of the overlay; and (iii) the van der Waals energy of the molecular conformations defined by the torsion angles encoded in the chromosomes. The algorithm has been applied to a number of pharmacophore elucidation problems, i.e., angiotensin II receptor antagonists, Leu-enkephalin and a hybrid morphine molecule, 5-HT1D agonists, benzodiazepine receptor ligands, 5-HT3 antagonists, dopamine D2 antagonists, dopamine reuptake blockers and FKBP12 ligands. The resulting pharmacophores are generated rapidly and are in good agreement with those derived from alternative means.

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

[2]  P. W. Codding,et al.  Molecular structure of Ro15-1788 and a model for the binding of benzodiazepine receptor ligands. Structural identification of common features in antagonists. , 1985, Molecular pharmacology.

[3]  D. E. Goldberg,et al.  Genetic Algorithms in Search , 1989 .

[4]  S Forrest,et al.  Genetic algorithms , 1996, CSUR.

[5]  Gareth Jones,et al.  Matching two-dimensional chemical graphs using genetic algorithms , 1994, J. Chem. Inf. Comput. Sci..

[6]  David A. Frenkel,et al.  CoMFA validation of the superposition of six classes of compounds which block GABA receptors non-competitively , 1993, J. Comput. Aided Mol. Des..

[7]  P. G. N. Digby,et al.  Multivariate Analysis of Ecological Communities , 1987 .

[8]  J C White,et al.  Novel antagonists of the 5-HT3 receptor. Synthesis and structure-activity relationships of (2-alkoxybenzoyl)ureas. , 1992, Journal of medicinal chemistry.

[9]  Tommy Liljefors,et al.  A textbook of drug design and development , 1996 .

[10]  R M Eglen,et al.  2-(Quinuclidin-3-yl)pyrido[4,3-b]indol-1-ones and isoquinolin-1-ones. Potent conformationally restricted 5-HT3 receptor antagonists. , 1993, Journal of medicinal chemistry.

[11]  Ralph E. Christoffersen,et al.  Computer-assisted drug design , 1979 .

[12]  Mark Froimowitz.  Conformational analysis of cocaine, the potent analog 2β‐carbomethoxy‐3β‐(4‐fluorophenyl)tropane (CFT), and other dopamine reuptake blockers , 1993, J. Comput. Chem..

[13]  C. F. Curtiss,et al.  Molecular Theory Of Gases And Liquids , 1954 .

[14]  Jon Clardy,et al.  DESIGN, SYNTHESIS, AND KINETIC EVALUATION OF HIGH-AFFINITY FKBP LIGANDS AND THE X-RAY CRYSTAL-STRUCTURES OF THEIR COMPLEXES WITH FKBP12. , 1994 .

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

[16]  J J Baldwin,et al.  Application of the three-dimensional structures of protein target molecules in structure-based drug design. , 1994, Journal of medicinal chemistry.

[17]  Garland R. Marshall,et al.  The Conformational Parameter in Drug Design: The Active Analog Approach , 1979 .

[18]  B. Masek,et al.  Molecular shape comparison of angiotensin II receptor antagonists. , 1993, Journal of medicinal chemistry.

[19]  David E. Goldberg,et al.  Genetic Algorithms in Search Optimization and Machine Learning , 1988 .

[20]  R F Standaert,et al.  Atomic structures of the human immunophilin FKBP-12 complexes with FK506 and rapamycin. , 1993, Journal of molecular biology.

[21]  Philip M. Dean,et al.  An exploration of a novel strategy for superposing several flexible molecules , 1993, J. Comput. Aided Mol. Des..

[22]  Lawrence. Davis,et al.  Handbook Of Genetic Algorithms , 1990 .

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

[24]  G J Williams,et al.  The Protein Data Bank: a computer-based archival file for macromolecular structures. , 1977, Journal of molecular biology.

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

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