Searching for pharmacophoric patterns in databases of three-dimensional chemical structures

This paper provides an overview of the research that has been carried out in Sheffield over the last decade into searching techniques for databases of three-dimensional (3D) chemical structures. A 3D structure or query pattern is represented by a labelled graph, in which the nodes and the edges of the graph are used to represent atoms and the associated inter-atomic distances, respectively. The presence of a pharmacophore in each of the structures in a database can then be tested by means of a subgraph isomorphism algorithm, the computational requirements of which are minimized by the use of an initial screening procedure that eliminates the majority of the structures from the subgraph-isomorphism search. Analogous graph-based representation and searching methods can also be used with flexible 3D structures: in this case, the edges of the graphs represent inter-atomic distance ranges and a final conformational search needs to be carried out for those molecules that match the query pharmacophore in the subgraph-isomorphism search. The paper also reviews related work on the automatic identification of pharmacophoric patterns and on 3D similarity searching.     

Pharmacophoric pattern matching in files of 3-D chemical structures: selection of interatomic distance screens

This paper describes an algorithmic selection procedure for the identification of interatomic distance screens. The screens are used to increase the search efficiency of an experimental pharmacophoric pattern matching system for 3-D structures taken from the Cambridge Crystallographic Data Bank. The selection procedure takes account of the variant frequencies of occurrence of the atomic types present in the file, and involves the identification of distance ranges so that the resulting screens occur with approximately equal frequencies.     

Pharmacophoric pattern matching in files of three-dimensional chemical structures: use of bounded distance matrices for the representation and searching of conformationally flexible molecules

This paper discusses the use of bounded distance matrices for the representation of conformationally flexible threedimensional (3D) molecules. It is shown that pharmacophoric pattern searches of databases of flexible 3D molecules represented in this way can be carried out using screen and geometric searching algorithms that are analogous to those used for searching databases of rigid 3D structures. Molecules matching a query pattern after the geometric search must then undergo a final conformational search to determine whether they can, in fact, adopt a conformation that matches the query. An analysis of this three-stage searching procedure shows that searching databases of flexible 3D molecules is extremely demanding of computational resources.     

Construction of a three-dimensional pharmacophore for Bcl-2 inhibitors by flexible docking and the multiple copy simultaneous search method

B-Cell lymphoma-2 (Bcl-2) protein is a new promising target for anticancer drugs. A number of anticancer Bcl-2 inhibitors with diverse chemical structures have been discovered in recent years. In this paper, the flexible docking was performed to determine the binding modes of the representative inhibitors from different structural types. Subsequently, the binding modes of inhibitor were used to construct a primary three- dimensional (3D) pharmacophore model. It proved that this model can effectively disrupt the binding of the BH3 domain of proapoptotic Bcl-2 family members to Bcl-2, and match the structural requirement of a new type of Bcl-2 inhibitors. However, these distances between pharmacophoric points are not optimal due to the fact that not all of individual functional groups are located in the ideal position when inhibitors bind to its receptor. In this paper, we proposed a new idea to improve the quality of the pharmacophore model: the multiple copy simultaneous search (MCSS) method was performed to determine the energetically favorable distribution of functional groups with similar features to these pharmacophoric points in the active site of Bcl-2 first. Then their most energetically favorable minima in the positions near the pharmacophoric points were used to optimize the distances between pharmacophoric points. By examining the binding modes of several inhibitors from the same structural type, it was found that the more potent the inhibitor was, the closer it was to the optimized distances between pharmacophoric points. The optimized 3D pharmacophore model obtained in this paper may provide a good starting point for further rational design of Bcl-2 inhibitors.     

CYP isoform specificity toward drug metabolism: analysis using common feature hypothesis

Three dimensional pharmacophoric maps were generated for each isoforms of CYP2C9, CYP2D6 and CYP3A4 separately using independent training sets consist of highly potent substrates (seven substrates for each isoform). HipHop module of CATALYST software was used in the generation of pharmacophore models. The best pharmacophore model was chosen out of the several models on the basis of (i) highest ranking score, (ii) better fit value among training set, (iii) capability to screen substrates from data set and (iv) efficiency to identify the isoform specificity. The individual pharmacophore models (CYP2C9-hypo1, CYP2D6-hypo1 and CYP3A4-hypo1) are characterized by the pharmacophoric features XZDH, RPZH and XYZHH for the CYP2C9, CYP2D6 and CYP3A4 respectively. Each of the chosen models was validated by using data sets of CYP substrates. This comparative study of CYP substrates demonstrates the importance of acidic character along with HBD and HBAl features for CYP2C9, basic character with ring aromatic features for CYP2D6 and hydrophobic features for CYP3A4. Acidity, basicity and hydrophobicity features arising from the functional groups of the substrates are also responsible for demonstrating CYP isoform specificity. Hence, these chemical features are incorporated in the decision tree along with pharmacophore maps. Finally, a decision tree based on chemical features and pharmacophore features was generated to identify the isoform specificity of novel query molecule toward the three isoforms.     

Pharmacophoric studies of in vitro inhibition of Plasmodium falciparum growth

Malaria continues to be a scourge in India and the situation has been compounded by the emergence of resistant strains of Plasmodium falciparum which is the primary cause of fatality in this disease. Therefore, there is an urgent need to develop newer drugs. Molecular modeling and pharmacophoric determination have become predominant methods today in the design and synthesis of newer and more effective drugs. Many Plasmodium specific enzymes and proteins involved in crucial biochemical pathways have been identified and their structures have been determined by X-ray crystallography. These enzymes and proteins are excellent targets for newer antimalarial agents. Bisphosphonates have shown potent inhibitory activity against Plasmodium farneysl diphosphate synthase (FPPS) enzyme, which is vital to the protein prenylation pathway of the organism. In this study, a set of 26 bisphosphonate inhibitors, synthesized by Oldfield et al [J Med Chem (2008) 51, 7827-7833] were subjected to rigorous 3D-QSAR studies using the PHASE computational package. In vitro Plasmodium growth inhibition rather than direct enzyme inhibition was considered in the study for a more realistic approach. Good statistical correlations were obtained for the pharmacophoric model as revealed by the regression values, indicating good stability of the model. Three hydrogen bond acceptors and a hydrogen bond donor defined the pharmacophore from the present study. This pharmacophore, AAAD (A = Hydrogen bond acceptor and D = Hydrogen bond donor) was put through a search-run for matching structures from the SPECS database yielding four matching structures, which could function as starting points for more novel and potent antimalarials.     

Use of techniques derived from graph theory to compare secondary structure motifs in proteins

A substructure matching algorithm is described that can be used for the automatic identification of secondary structural motifs in three-dimensional protein structures from the Protein Data Bank. The proteins and motifs are stored for searching as labelled graphs, with the nodes of a graph corresponding to linear representations of helices and strands and the edges to the inter-line angles and distances. A modification of Ullman's subgraph isomorphism algorithm is described that can be used to search these graph representations. Tests with patterns from the protein structure literature demonstrate both the efficiency and the effectiveness of the search procedure, which has been implemented in FORTRAN 77 on a MicroVAX-II system, coupled to the molecular fitting program FRODO on an Evans and Sutherland PS300 graphics system.     

Pharmacophore mapping and in silico screening to identify new potent leads for A(2A) adenosine receptor as antagonists

A(2A) adenosine receptor (AR) antagonists play an important role in neurodegenerative diseases like Parkinson's disease. A 3D-QSAR study of A(2A) AR antagonists, was taken up to design best pharmacophore model. The pharmacophoric features (ADHRR) containing a hydrogen bond acceptor (A), a hydrogen bond donor (D), a hydrophobic group (H) and two aromatic rings (R), is projected as the best predictive pharmacophore model. The QSAR model was further treated as a template for in silico search of databases to identify new scaffolds. The binding patterns of the leads with A(2A) AR are analysed using docking studies and novel potent ligands of A(2A) AR are projected.     

CO: A chemical ontology for identification of functional groups and semantic comparison of small molecules

A novel chemical ontology based on chemical functional groups automatically, objectively assigned by a computer program, was developed to categorize small molecules. It has been applied to PubChem and the small molecule interaction database to demonstrate its utility as a basic pharmacophore search system. Molecules can be compared using a semantic similarity score based on functional group assignments rather than 3D shape, which succeeds in identifying small molecules known to bind a common binding site. This ontology will serve as a powerful tool for searching chemical databases and identifying key functional groups responsible for biological activities.     

Pharmacophore mapping and in silico screening to identify new potent leads for A2A adenosine receptor as antagonists

A_2A adenosine receptor (AR) antagonists play an important role in neurodegenerative diseases like Parkinson’s disease. A 3D-QSAR study of A_2A AR antagonists, was taken up to design best pharmacophore model. The pharmacophoric features (ADHRR) containing a hydrogen bond acceptor (A), a hydrogen bond donor (D), a hydrophobic group (H) and two aromatic rings (R), is projected as the best predictive pharmacophore model. The QSAR model was further treated as a template for in silico search of databases to identify new scaffolds. The binding patterns of the leads with A_2A AR are analysed using docking studies and novel potent ligands of A_2A AR are projected.     

Pharmacophore and Virtual Screening of JAK3 inhibitors

Janus kinase 3 (JAK3) is a non-receptor tyrosine kinases family of protein which is comprised of JAK1, JAK2, JAK3 and TYK2. It plays an important role in immune function and lymphoid development and it only resides in the hematopoietic system. Therefore, selective targeting JAK3 is a rational approach in developing new therapeutic molecule. In this study, about 116 JAK3 inhibitors were collected from the literature and were used to build four-point pharmacophore model using Phase (Schrodinger module). The statistically significant pharmacophore hypothesis of AAHR.92 with r2 value of 0.942 was used as 3D query to search against 3D database namely Zincpharmer. A total of 2, 27,483 compounds obtained as hit were subjected to high throughput virtual screening (HTVS module of Schrodinger). Among the hits, ten compounds with good G-score ranging from -12.96 to -11.18 with good binding energy to JAK3 were identified.     

Discovery of novel fatty acid synthase (FAS) inhibitors based on the structure of ketoaceyl synthase (KS) domain

Development of fatty acid synthase (FAS) inhibitors has increasingly attracted much attention in recent years due to their potential therapeutic use in obesity and cancers. In this investigation, pharmacophore modeling based on the first crystal structure of human KS domain of FAS was carried out. The established pharmacophore model was taken as a 3D query for retrieving potent FAS inhibitors from the chemical database Specs. Docking study was further carried out to refine the obtained hit compounds. Finally, a total of 28 compounds were selected based on the ranking order and visual examination, which were first evaluated by a cell line-based assay. Seven compounds that have good inhibition activity against two FAS overexpressing cancer cell lines were further evaluated by an enzyme-based assay. One compound with a new chemical scaffold was found to have low micromolar inhibition potency against FAS, which has been subjected to further chemical structural modification.     

Pharmacophore modeling study based on known Spleen tyrosine kinase inhibitors together with virtual screening for identifying novel inhibitors

In this investigation, chemical features based 3D pharmacophore models were developed based on the known inhibitors of Spleen tyrosine kinase (Syk) with the aid of hiphop and hyporefine modules within catalyst. The best quantitative pharmacophore model, Hypo1, was used as a 3D structural query for retrieving potential inhibitors from chemical databases including Specs, NCI, MayBridge, and Chinese Nature Product Database (CNPD). The hit compounds were subsequently subjected to filtering by Lipinski’s rule of five and docking studies to refine the retrieved hits. Finally 30 compounds were selected from the top ranked hit compounds and conducted an in vitro kinase inhibitory assay. Six compounds showed a good inhibitory potency against Syk, which have been selected for further investigation.     

Query Chem: a Google-powered web search combining text and chemical structures

Query Chem (www.QueryChem.com) is a Web program that integrates chemical structure and text-based searching using publicly available chemical databases and Google's Web Application Program Interface (API). Query Chem makes it possible to search the Web for information about chemical structures without knowing their common names or identifiers. Furthermore, a structure can be combined with textual query terms to further restrict searches. Query Chem's search results can retrieve many interesting structure-property relationships of biomolecules on the Web.     

Rapid 3D protein structure database searching using information retrieval techniques

MOTIVATION: As the sizes of three-dimensional (3D) protein structure databases are growing rapidly nowadays, exhaustive database searching, in which a 3D query structure is compared to each and every structure in the database, becomes inefficient. We propose a rapid 3D protein structure retrieval system named 'ProtDex2', in which we adopt the techniques used in information retrieval systems in order to perform rapid database searching without having access to every 3D structure in the database. The retrieval process is based on the inverted-file index constructed on the feature vectors of the relationships between the secondary structure elements (SSEs) of all the 3D protein structures in the database. ProtDex2 is a significant improvement, both in terms of speed and accuracy, upon its predecessor system, ProtDex. RESULTS: The experimental results show that ProtDex2 is very much faster than two well-known protein structure comparison methods, DALI and CE, yet not sacrificing on the accuracy of the comparison. When comparing with a similar SSE-based method, namely TopScan, ProtDex2 is much faster with comparable degree of accuracy. AVAILABILITY: The software is available at: http://xena1.ddns.comp.nus.edu.sg/~genesis/PD2.htm     

The SuMo server: 3D search for protein functional sites

SUMMARY: We provide the scientific community with a web server which gives access to SuMo, a bioinformatic system for finding similarities in arbitrary 3D structures or substructures of proteins. SuMo is based on a unique representation of macromolecules using selected triplets of chemical groups having their own geometry and symmetry, regardless of the restrictive notions of main chain and lateral chains of amino acids. The heuristic for extracting similar sites was used to drive two major large-scale approaches. First, searching for ligand binding sites onto a query structure has been made possible by comparing the structure against each of the ligand binding sites found in the Protein Data Bank (PDB). Second, the reciprocal process, i.e. searching for a given 3D site of interest among the structures of the PDB is also possible and helps detect cross-reacting targets in drug design projects. AVAILABILITY: The web server is freely accessible to academia through http://sumo-pbil.ibcp.fr and full support is available from MEDIT (http://www.medit.fr). CONTACT: mjambon@burnham.org.