Calculating the knowledge-based similarity of functional groups using crystallographic data

A knowledge-based method for calculating the similarity of functional groups is described and validated. The method is based on experimental information derived from small molecule crystal structures. These data are used in the form of scatterplots that show the likelihood of a non-bonded interaction being formed between functional group A (the `central group') and functional group B (the `contact group' or `probe'). The scatterplots are converted into three-dimensional maps that show the propensity of the probe at different positions around the central group. Here we describe how to calculate the similarity of a pair of central groups based on these maps. The similarity method is validated using bioisosteric functional group pairs identified in the Bioster database and Relibase. The Bioster database is a critical compilation of thousands of bioisosteric molecule pairs, including drugs, enzyme inhibitors and agrochemicals. Relibase is an object-oriented database containing structural data about protein-ligand interactions. The distributions of the similarities of the bioisosteric functional group pairs are compared with similarities for all the possible pairs in IsoStar, and are found to be significantly different. Enrichment factors are also calculated showing the similarity method is statistically significantly better than random in predicting bioisosteric functional group pairs.     

Energetic analysis of fragment docking and application to structure-based pharmacophore hypothesis generation

We have developed a method that uses energetic analysis of structure-based fragment docking to elucidate key features for molecular recognition. This hybrid ligand- and structure-based methodology uses an atomic breakdown of the energy terms from the Glide XP scoring function to locate key pharmacophoric features from the docked fragments. First, we show that Glide accurately docks fragments, producing a root mean squared deviation (RMSD) of <1.0 Å for the top scoring pose to the native crystal structure. We then describe fragment-specific docking settings developed to generate poses that explore every pocket of a binding site while maintaining the docking accuracy of the top scoring pose. Next, we describe how the energy terms from the Glide XP scoring function are mapped onto pharmacophore sites from the docked fragments in order to rank their importance for binding. Using this energetic analysis we show that the most energetically favorable pharmacophore sites are consistent with features from known tight binding compounds. Finally, we describe a method to use the energetically selected sites from fragment docking to develop a pharmacophore hypothesis that can be used in virtual database screening to retrieve diverse compounds. We find that this method produces viable hypotheses that are consistent with known active compounds. In addition to retrieving diverse compounds that are not biased by the co-crystallized ligand, the method is able to recover known active compounds from a database screen, with an average enrichment of 8.1 in the top 1% of the database.     

Automated site-directed drug design: An assessment of the transferability of atomic residual charges (CNDO) for molecular fragments

Summary In this paper a database of atomic residual charges has been constructed for all the molecular fragments defined previously in a combinatorial search of the Cambridge Structural Database. The charges generated for the atoms in each fragment are compared with charges calculated for whole molecules containing those fragments. The fragment atomic charges lie within 1 S.D. of the mean for 68%, and within 2 S.D. for 91%, of the atoms whose charges were computed for whole molecules. The actual charges on any atom are strongly influenced by the adjacent connected atoms. There is a large spread of atomic residual charge within the fragments database.     

Developing an antituberculosis compounds database and data mining in the search of a motif responsible for the activity of a diverse class of antituberculosis agents

A novel data mining procedure to look for new antitubercular agents and targets as well as to find a minimum common bioactive substructure (MCBS), has been reported here. The methodology extracts MCBS, both across the diverse chemical classes and within the particular chemical class, known to be present in the various marketed drugs alongside antimycobacterial compounds with known MICs. For this purpose a small in-house database of compounds has been created, for which MICs against Mycobacterium are known. The compounds have been collected from literature available on the synthetic compounds, having known MICs against Mycobacterium tuberculosis. An elaborate HQSAR (Hologram QSAR) study has been attempted to extract active fragment from a diverse class of compounds, in combination with the clustering technique to select a homogeneous group of compounds having good a profile toward the activity. The 2D pharmacophore (the 2D fragments extracted from HQSAR) has been validated searching the database. It has been found further that this validated 2D pharmacophore could be used for searching the orphan target in Mycobacterium effectively.     

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.     

Modeling flexible pharmacophores with distance geometry, scoring, and bound stretching

The study of pharmacophores, i.e., of common features between different ligands, is important for the quantitative identification of "compatible" enzymes and binding species. A pharmacophore-based technique is developed that combines multiple conformations with a distance geometry method to create flexible pharmacophore representations. It uses a set of low-energy conformations combined with a new process we call bound stretching to create sets of distance bounds, which contain all or most of the low-energy conformations. The bounds can be obtained using the exact distances between pairs of atoms from the different low-energy conformations. To avoid missing conformations, we can take advantage of the triangle distance inequality between sets of three points to logically expand a set of upper and lower distance bounds (bound stretching). The flexible pharmacophore can be found using a 3-D maximal common subgraph method, which uses the overlap of distance bounds to determine the overlapping structure. A scoring routine is implemented to select the substructures with the largest overlap because there will typically be many overlaps with the maximum number of overlapping bounds. A case study is presented in which 3-D flexible pharmacophores are generated and used to eliminate potential binding species identified by a 2-D pharmacophore method. A second case study creates flexible pharmacophores from a set of thrombin ligands. These are used to compare the new method with existing pharmacophore identification software.     

Development and NMR validation of minimal pharmacophore hypotheses for the generation of fragment libraries enriched in heparanase inhibitors

A combined strategy based on the development of pharmacophore hypotheses and NMR approaches is reported for the identification of novel inhibitors of heparanase, a key enzyme involved in tumor metastasis through the remodeling of the subepithelial and subendothelial basement membranes, resulting in the dissemination of metastatic cancer cells. Several pharmacophore hypotheses were initially developed from the most active heparanase inhibitors known to date and, after their application to a pool of 27 known heparanase inhibitors and a database of 1,120 compounds approved by the FDA, a four-point pharmacophore model was selected as the most predictive. This model was subsequently applied to a database of 686 chemical fragments, and a subset of 100 fragments accomplishing completely or partially the four-point model was selected to perform nuclear magnetic resonance experiments to validate the hypothesis. The experimental studies confirmed the reliability of our pharmacophore model, its applicability to in silico databases in order to reduce the number of compounds to be experimentally screened, and the possibility of generating fragment libraries enriched in heparanase inhibitors. Rafael Gozalbes and Silvia Mosulén contributed equally to this work.     

CONFIRM: connecting fragments found in receptor molecules

A novel algorithm for the connecting of fragment molecules is presented and validated for a number of test systems. Within the CONFIRM (Connecting Fragments Found in Receptor Molecules) approach a pre-prepared library of bridges is searched to extract those which match a search criterion derived from known experimental or computational binding information about fragment molecules within a target binding site. The resulting bridge 'hits' are then connected, in an automated fashion, to the fragments and docked into the target receptor. Docking poses are assessed in terms of root-mean-squared deviation from the known positions of the fragment molecules, as well as docking score should known inhibitors be available. The creation of the bridge library, the full details and novelty of the CONFIRM algorithm, and the general applicability of this approach within the field of fragment-based de novo drug design are discussed.     

New leads for selective GSK-3 inhibition: pharmacophore mapping and virtual screening studies

Summary Glycogen Synthase Kinase-3 is a regulatory serine/threonine kinase, which is being targeted for the treatment of a number of human diseases including type-2 diabetes mellitus, neurodegenerative diseases, cancer and chronic inflammation. Selective GSK-3 inhibition is an important requirement owing to the possibility of side effects arising from other kinases. A pharmacophore mapping strategy is employed in this work to identify new leads for selective GSK-3 inhibition. Ligands known to show selective GSK-3 inhibition were employed in generating a pharmacophore map using distance comparison method (DISCO). The derived pharmacophore map was validated using (i) important interactions involved in selective GSK-3 inhibitions, and (ii) an in-house database containing different classes of GSK-3 selective, non-selective and inactive molecules. New Lead identification was carried out by performing virtual screening using validated pharmacophoric query and three chemical databases namely NCI, Maybridge and Leadquest. Further data reduction was carried out by employing virtual filters based on (i) Lipinski’s rule of 5 (ii) van der Waals bumps and (iii) restricting the number of rotatable bonds to seven. Final screening was carried out using FlexX based molecular docking study.     

Reaction site mapping of xenobiotic biotransformations

Predictive metabolism methods can be used in drug discovery projects to enhance the understanding of structure-metabolism relationships. The present study uses data mining methods to exploit biotransformation data that have been recorded in the MDL Metabolite database. Reacting center fingerprints were derived from a comparison of substrates and their corresponding products listed in the database. This process yields two fingerprint databases: all atoms in all substrates and all reacting centers. The metabolic reaction data are then mined by submitting a new molecule and searching for fingerprint matches to every atom in the new molecule in both databases. An "occurrence ratio" is derived from the fingerprint matches between the submitted compound and the reacting center and substrate fingerprint databases. Normalization of the occurrence ratio within each submitted molecule enables the results of the search to be rank-ordered as a measure of the relative frequency of a reaction occurring at a specific site within the submitted molecule. Predictive performance that would allow this method to be used by drug discovery teams to generate useful hypotheses regarding structure metabolism relationships was observed.     

FLAME: a program to flexibly align molecules

Herein, we describe a method to flexibly align molecules (FLAME = FLexibly Align MolEcules). FLAME aligns two molecules by first finding maximum common pharmacophores between them using a genetic algorithm. The resulting alignments are then subjected to simultaneous optimizations of their internal energies and an alignment score. The utility of the method in pairwise alignment, multiple molecule flexible alignment, and database searching was examined. For pairwise alignment, two carboxypeptidase ligands (Protein Data Bank codes and ), two estrogen receptor ligands ( and ), and two thrombin ligands ( and ) were used as test sets. Alignments generated by FLAME starting from CONCORD structures compared very well to the X-ray structures (average root-mean-square deviation = 0.36 A) even without further minimization in the presence of the protein. For multiple flexible alignments, five structurally diverse D3 receptor ligands were used as a test set. The FLAME alignment automatically identified three common pharmacophores: a base, a hydrogen-bond acceptor, and a hydrophobe/aromatic ring. The best alignment was then used to search the MDDR database. The search results were compared to the results using atom pair and Daylight fingerprint similarity. A similar database search comparison was also performed using estrogen receptor modulators. In both cases, hits identified by FLAME were structurally more diverse compared to those from the atom pair and Daylight fingerprint methods.     

Introduction of the conditional correlated Bernoulli model of similarity value distributions and its application to the prospective prediction of fingerprint search performance

A statistical approach named the conditional correlated Bernoulli model is introduced for modeling of similarity scores and predicting the potential of fingerprint search calculations to identify active compounds. Fingerprint features are rationalized as dependent Bernoulli variables and conditional distributions of Tanimoto similarity values of database compounds given a reference molecule are assessed. The conditional correlated Bernoulli model is utilized in the context of virtual screening to estimate the position of a compound obtaining a certain similarity value in a database ranking. Through the generation of receiver operating characteristic curves from cumulative distribution functions of conditional similarity values for known active and random database compounds, one can predict how successful a fingerprint search might be. The comparison of curves for different fingerprints makes it possible to identify fingerprints that are most likely to identify new active molecules in a database search given a set of known reference molecules.     

A self-organizing algorithm for molecular alignment and pharmacophore development

We present a method for simultaneous three-dimensional (3D) structure generation and pharmacophore-based alignment using a self-organizing algorithm called Stochastic Proximity Embedding (SPE). Current flexible molecular alignment methods either start from a single low-energy structure for each molecule and tweak bonds or torsion angles, or choose from multiple conformations of each molecule. Methods that generate structures and align them iteratively (e.g., genetic algorithms) are often slow. In earlier work, we used SPE to generate good-quality 3D conformations by iteratively adjusting pairwise distances between atoms based on a set of geometric rules, and showed that it samples conformational space better and runs faster than earlier programs. In this work, we run SPE on the entire ensemble of molecules to be aligned. Additional information about which atoms or groups of atoms in each molecule correspond to points in the pharmacophore can come from an automatically generated hypothesis or be specified manually. We add distance terms to SPE to bring pharmacophore points from different molecules closer in space, and also to line up normal/direction vectors associated with these points. We also permit pharmacophore points to be constrained to lie near external coordinates from a binding site. The aligned 3D molecular structures are nearly correct if the pharmacophore hypothesis is chemically feasible; postprocessing by minimization of suitable distance and energy functions further improves the structures and weeds out infeasible hypotheses. The method can be used to test 3D pharmacophores for a diverse set of active ligands, starting from only a hypothesis about corresponding atoms or groups.     

Protein structure prediction from inaccurate and sparse NMR data using an enhanced genetic algorithm

Nuclear Magnetic Resonance Spectroscopy (most commonly known as NMR Spectroscopy) is used to generate approximate and partial distances between pairs of atoms of the native structure of a protein. To predict protein structure from these partial distances by solving the Euclidean distance geometry problem from the partial distances obtained from NMR Spectroscopy, we can predict three-dimensional (3D) structure of a protein. In this paper, a new genetic algorithm is proposed to efficiently address the Euclidean distance geometry problem towards building 3D structure of a given protein applying NMR's sparse data. Our genetic algorithm uses (i) a greedy mutation and crossover operator to intensify the search; (ii) a twin removal technique for diversification in the population; (iii) a random restart method to recover from stagnation; and (iv) a compaction factor to reduce the search space. Reducing the search space drastically, our approach improves the quality of the search. We tested our algorithms on a set of standard benchmarks. Experimentally, we show that our enhanced genetic algorithms significantly outperforms the traditional genetic algorithms and a previously proposed state-of-the-art method. Our method is capable of producing structures that are very close to the native structures and hence, the experimental biologists could adopt it to determine more accurate protein structures from NMR data.     

Structure-based fragment hopping for lead optimization using predocked fragment database

In this work, we describe a structure-based de novo optimization process, called "LeadOp" (short for lead optimization), that decomposes a compound into fragments of different molecular components either by chemical or user-defined rules. Each fragment is evaluated through a predocked fragment database that ranks fragments according to specific fragment-receptor binding interactions, replacing fragments that contribution the least to binding and finally reassembling the fragments to form a new ligand. The fundamental idea is to replace "bad" fragments of a ligand with "good" fragments while leaving the core of the ligand intact, thus improving the compound's activity. The molecular fragments were selected from a collection of 27,417 conformers that are the fragments of compounds in the DrugBank database. The collection of molecular fragments are docked to the target's binding site and evaluated using group efficiency (calculated binding affinity divided by the number of heavy atoms), and the "strongest" binder is selected. The LeadOp method was tested with two biomolecular systems: mutant B-Raf kinase and human 5-lipoxygenase. The LeadOp methodology was able to optimize the query molecules and systematically developed improved analogs for each of our example systems.