Pharmacophore Fingerprint-Based Approach to Binding Site Subpocket Similarity and Its Application to Bioisostere Replacement

Bioisosteres have been defined as structurally different molecules or substructures that can form comparable intermolecular interactions, and therefore, fragments that bind to similar protein structures exhibit a degree of bioisosterism. We present KRIPO (Key Representation of Interaction in POckets): a new method for quantifying the similarities of binding site subpockets based on pharmacophore fingerprints. The binding site fingerprints have been optimized to improve their performance for both intra- and interprotein family comparisons. A range of attributes of the fingerprints was considered in the optimization, including the placement of pharmacophore features, whether or not the fingerprints are fuzzified, and the resolution and complexity of the pharmacophore fingerprints (2-, 3-, and 4-point fingerprints). Fuzzy 3-point pharmacophore fingerprints were found to represent the optimal balance between computational resource requirements and the identification of potential replacements. The complete PDB was converted into a database comprising almost 300,000 optimized fingerprints of local binding sites together with their associated ligand fragments. The value of the approach is demonstrated by application to two crystal structures from the Protein Data Bank: (1) a MAP kinase P38 structure in complex with a pyridinylimidazole inhibitor (1A9U) and (2) a complex of thrombin with melagatran (1K22). Potentially valuable bioisosteric replacements for all subpockets of the two studied protein are identified.

[1]  G Klebe,et al.  Factorising ligand affinity: a combined thermodynamic and crystallographic study of trypsin and thrombin inhibition. , 2001, Journal of molecular biology.

[2]  Kalidas Yeturu,et al.  PocketMatch: A new algorithm to compare binding sites in protein structures , 2008, BMC Bioinformatics.

[3]  Y. Martin,et al.  A practitioner's perspective of the role of quantitative structure-activity analysis in medicinal chemistry. , 1981, Journal of medicinal chemistry.

[4]  Gerhard Klebe,et al.  From Structure to Function: A New Approach to Detect Functional Similarity among Proteins Independent from Sequence and Fold Homology. , 2001, Angewandte Chemie.

[5]  Izhar Wallach,et al.  Prediction of sub-cavity binding preferences using an adaptive physicochemical structure representation , 2009, Bioinform..

[6]  G. Klebe,et al.  From the Similarity Analysis of Protein Cavities to the Functional Classification of Protein Families Using Cavbase , 2006, Journal of Molecular Biology.

[7]  Gabriele Cruciani,et al.  A Common Reference Framework for Analyzing/Comparing Proteins and Ligands. Fingerprints for Ligands And Proteins (FLAP): Theory and Application , 2007, J. Chem. Inf. Model..

[8]  G. Klebe,et al.  A new method to detect related function among proteins independent of sequence and fold homology. , 2002, Journal of molecular biology.

[9]  Joseph S. Verducci,et al.  A Modification of the Jaccard–Tanimoto Similarity Index for Diverse Selection of Chemical Compounds Using Binary Strings , 2002, Technometrics.

[10]  Artem Vorotyntsev,et al.  A Computational Fragment Approach by Mining the Protein Data Bank: Library Design and Bioisosterism , 2011 .

[11]  Eyke Hüllermeier,et al.  Functional Classification of Protein Kinase Binding Sites Using Cavbase , 2007, ChemMedChem.

[12]  E. LaVoie,et al.  Bioisosterism: A Rational Approach in Drug Design. , 1996, Chemical reviews.

[13]  P. Lam,et al.  Structure-based design of novel guanidine/benzamidine mimics: potent and orally bioavailable factor Xa inhibitors as novel anticoagulants. , 2003, Journal of medicinal chemistry.

[14]  Dragos Horvath,et al.  Neighborhood Behavior of in Silico Structural Spaces with Respect to In Vitro Activity Spaces-A Benchmark for Neighborhood Behavior Assessment of Different in Silico Similarity Metrics , 2003, J. Chem. Inf. Comput. Sci..

[15]  S. Henikoff,et al.  Amino acid substitution matrices from protein blocks. , 1992, Proceedings of the National Academy of Sciences of the United States of America.

[16]  G. Moreau,et al.  Rational approach to drug design , 2008 .

[17]  Gert Vriend,et al.  Increasing the precision of comparative models with YASARA NOVA—a self‐parameterizing force field , 2002, Proteins.

[18]  A. Konagurthu,et al.  MUSTANG: A multiple structural alignment algorithm , 2006, Proteins.

[19]  Jean-Philippe Vert,et al.  A new protein binding pocket similarity measure based on comparison of clouds of atoms in 3D: application to ligand prediction , 2010, BMC Bioinformatics.

[20]  Robert D. Clark,et al.  Managing bias in ROC curves , 2008, J. Comput. Aided Mol. Des..

[21]  W L Jorgensen,et al.  Estimation of binding affinities for selective thrombin inhibitors via Monte Carlo simulations. , 2001, Journal of medicinal chemistry.

[22]  Gerhard Klebe,et al.  KNOBLE: a knowledge-based approach for the design and synthesis of readily accessible small-molecule chemical probes to test protein binding. , 2007, Angewandte Chemie.

[23]  James E. J. Mills,et al.  High-Throughput Virtual Screening of Proteins Using GRID Molecular Interaction Fields , 2010, J. Chem. Inf. Model..

[24]  Christopher A. Lipinski,et al.  Chapter 27. Bioisosterism in Drug Design , 1986 .

[25]  Markus Wagener,et al.  The Quest for Bioisosteric Replacements , 2006, J. Chem. Inf. Model..

[26]  C. Sander,et al.  Detection of common three‐dimensional substructures in proteins , 1991, Proteins.

[27]  Christophe Combet,et al.  The SuMo server: 3D search for protein functional sites , 2005, Bioinform..

[28]  Naomie Salim,et al.  Analysis and Display of the Size Dependence of Chemical Similarity Coefficients , 2003, J. Chem. Inf. Comput. Sci..

[29]  M. Jambon,et al.  A new bioinformatic approach to detect common 3D sites in protein structures , 2003, Proteins.

[30]  E. Myers,et al.  Basic local alignment search tool. , 1990, Journal of molecular biology.

[31]  V. Ramensky,et al.  A novel approach to local similarity of protein binding sites substantially improves computational drug design results , 2007, Proteins.

[32]  David Weininger,et al.  SMILES, a chemical language and information system. 1. Introduction to methodology and encoding rules , 1988, J. Chem. Inf. Comput. Sci..

[33]  J. Mason,et al.  New 4-point pharmacophore method for molecular similarity and diversity applications: overview of the method and applications, including a novel approach to the design of combinatorial libraries containing privileged substructures. , 1999, Journal of medicinal chemistry.

[34]  S Wold,et al.  Statistical molecular design, parallel synthesis, and biological evaluation of a library of thrombin inhibitors. , 2001, Journal of medicinal chemistry.

[35]  Nicola D. Gold,et al.  SitesBase: a database for structure-based protein–ligand binding site comparisons , 2005, Nucleic Acids Res..

[36]  Ksenia Oguievetskaia,et al.  Computational Fragment-Based Approach at PDB Scale by Protein Local Similarity , 2009, J. Chem. Inf. Model..

[37]  J. Kraus,et al.  Isosterism and molecular modification in drug design: tetrazole analogue of GABA: effects on enzymes of the gamma-aminobutyrate system. , 1983, Pharmacological research communications.

[38]  István Ujváry,et al.  Extended Summary: BIOSTER—a database of structurally analogous compounds , 1997 .

[39]  C. Thornber,et al.  Isosterism and molecular modification in drug design , 1979 .

[40]  P. Olesen,et al.  The use of bioisosteric groups in lead optimization. , 2001, Current opinion in drug discovery & development.

[41]  Jacob D. Durrant,et al.  CrystalDock: A Novel Approach to Fragment-Based Drug Design , 2011, J. Chem. Inf. Model..

[42]  E. Goldsmith,et al.  Structural basis of inhibitor selectivity in MAP kinases. , 1998, Structure.

[43]  A Burger,et al.  Isosterism and bioisosterism in drug design. , 1991, Progress in drug research. Fortschritte der Arzneimittelforschung. Progres des recherches pharmaceutiques.

[44]  Nathanael Weill,et al.  Alignment-Free Ultra-High-Throughput Comparison of Druggable Protein-Ligand Binding Sites , 2010, J. Chem. Inf. Model..

[45]  Alex M Aronov,et al.  Toward a pharmacophore for kinase frequent hitters. , 2004, Journal of medicinal chemistry.

[46]  Jonathan S. Mason,et al.  Ligand-Receptor 3-D Similarity Studies Using Multiple 4-Point Pharmacophores , 1998, Pacific Symposium on Biocomputing.

[47]  L. M. Lima,et al.  Bioisosterism: a useful strategy for molecular modification and drug design. , 2005, Current medicinal chemistry.

[48]  Peter Willett,et al.  Identification of target-specific bioisosteric fragments from ligand–protein crystallographic data , 2006, J. Comput. Aided Mol. Des..

[49]  Lawrence C Kuo,et al.  Unexpected enhancement of thrombin inhibitor potency with o-aminoalkylbenzylamides in the P1 position. , 2003, Bioorganic & medicinal chemistry letters.