Virtual Screening and Scaffold Hopping Based on GRID Molecular Interaction Fields

Molecular medicinal insights into scaffold hopping-based drug discovery success

In both academia and the pharmaceutical industry, innovative hypotheses, methodologies and technologies that can shorten the drug research and development, leading to higher success rates, are vital. In this review, we demonstrate how innovative variations of the scaffold-hopping strategy have been used to create new druggable molecular spaces, drugs, clinical candidates, preclinical candidates, and bioactive agents. We also analyze molecular modulations that enabled improvements of the pharmacodynamic (PD), physiochemical, and pharmacokinetic (PK) properties (P3 properties) of the drugs resulting from these scaffold-hopping strategies.

Computational Drug Repurposing Resources and Approaches for Discovering Novel Antifungal Drugs against Candida albicans N-Myristoyl Transferase

Candida albicans is a yeast that is an opportunistic fungal pathogen and also identified as ubiquitous polymorphic species that is mainly linked with major fungal infections in humans, particularly in the immunocompromised patients including transplant recipients, chemotherapy patients, HIV-infected patients as well as in low-birth-weight infants. Systemic Candida infections have a high mortality rate of around 29 to 76%. For reducing its infection, limited drugs are existing such as caspofungin, fluconazole, terbinafine, and amphotericin B, etc. which contain unlikable side effects and also toxic. This review intends to utilize advanced bioinformatics technologies such as Molecular docking, Scaffold hopping, Virtual screening, Pharmacophore modeling, Molecular dynamics (MD) simulation for the development of potentially new drug candidates with a drug-repurpose approach against Candida albicans within a limited time frame and also cost reductive.

Screening of Bacterial Quorum Sensing Inhibitors in a Vibrio fischeri LuxR-Based Synthetic Fluorescent E. coli Biosensor

A library of 23 pure compounds of varying structural and chemical characteristics was screened for their quorum sensing (QS) inhibition activity using a synthetic fluorescent Escherichia coli biosensor that incorporates a modified version of lux regulon of Vibrio fischeri. Four such compounds exhibited QS inhibition activity without compromising bacterial growth, namely, phenazine carboxylic acid (PCA), 2-heptyl-3-hydroxy-4-quinolone (PQS), 1H-2-methyl-4-quinolone (MOQ) and genipin. When applied at 50 µM, these compounds reduced the QS response of the biosensor to 33.7% ± 2.6%, 43.1% ± 2.7%, 62.2% ± 6.3% and 43.3% ± 1.2%, respectively. A series of compounds only showed activity when tested at higher concentrations. This was the case of caffeine, which, when applied at 1 mM, reduced the QS to 47% ± 4.2%. In turn, capsaicin, caffeic acid phenethyl ester (CAPE), furanone and polygodial exhibited antibacterial activity when applied at 1mM, and reduced the bacterial growth by 12.8% ± 10.1%, 24.4% ± 7.0%, 91.4% ± 7.4% and 97.5% ± 3.8%, respectively. Similarly, we confirmed that trans-cinnamaldehyde and vanillin, when tested at 1 mM, reduced the QS response to 68.3% ± 4.9% and 27.1% ± 7.4%, respectively, though at the expense of concomitantly reducing cell growth by 18.6% ± 2.5% and 16% ± 2.2%, respectively. Two QS natural compounds of Pseudomonas aeruginosa, namely PQS and PCA, and the related, synthetic compounds MOQ, 1H-3-hydroxyl-4-quinolone (HOQ) and 1H-2-methyl-3-hydroxyl-4-quinolone (MHOQ) were used in molecular docking studies with the binding domain of the QS receptor TraR as a target. We offer here a general interpretation of structure-function relationships in this class of compounds that underpins their potential application as alternatives to antibiotics in controlling bacterial virulence.

The Role of APOSTART in Switching between Sexuality and Apomixis in Poa pratensis

The production of seeds without sex is considered the holy grail of plant biology. The transfer of apomixis to various crop species has the potential to transform plant breeding, since it will allow new varieties to retain valuable traits thorough asexual reproduction. Therefore, a greater molecular understanding of apomixis is fundamental. In a previous work we identified a gene, namely APOSTART, that seemed to be involved in this asexual mode of reproduction, which is very common in Poa pratensis L., and here we present a detailed work aimed at clarifying its role in apomixis. In situ hybridization showed that PpAPOSTART is expressed in reproductive tissues from pre-meiosis to embryo development. Interestingly, it is expressed early in few nucellar cells of apomictic individuals possibly switching from a somatic to a reproductive cell as in aposporic apomixis. Moreover, out of 13 APOSTART members, we identified one, APOSTART_6, as specifically expressed in flower tissue. APOSTART_6 also exhibited delayed expression in apomictic genotypes when compared with sexual types. Most importantly, the SCAR (Sequence Characterized Amplified Region) derived from the APOSTART_6 sequence completely co-segregated with apomixis.

Characterization of PD-L1 binding sites by a combined FMO/GRID-DRY approach

The programmed cell death protein 1 (PD-1) and its ligand, PD-L1, constitute an important co-inhibitory immune checkpoint leading to downregulation of immune system. Tumor cells developed a strategy to trigger PD-1/PD-L1 pathway reducing the T cell anticancer activity. Anti-PD-L1 small drugs, generally with improved pharmacokinetic and technological profiles than monoclonal antibodies, became an attractive research topic. Nevertheless, still few works have been published on the chemical features of possible binding sites. In this work, we applied a novel computational protocol based on the combination of the ab initio Fragment Molecular Orbital (FMO) method and a newly developed GRID-DRY approach in order to characterize the PD-L1 binding sites, starting from PD-1/PD-L1 and PD-L1/BMS-ligands (Bristol-Mayers Squibb ligands) complexes. The FMO method allows the calculation of the pair-residues as well as the ligand-residues interactions with ab initio accuracy, whereas the GRID-DRY approach is an effective tool to investigate hydrophobic interactions, not easily detectable by ab initio methods. The present GRID-DRY protocol is able to determine the energy contributions of each ligand atoms to each hydrophobic interaction, both qualitatively and quantitatively. We were also able to identify the three specific hot regions involved in PD-1/PD-L1 protein-protein interaction and in PD-L1/BMS-ligand interactions, in agreement with preceding theoretical/experimental results, and to suggest a specific pharmacophore for PD-L1 inhibitors.

WebMetabase: cleavage sites analysis tool for natural and unnatural substrates from diverse data source

SUMMARY

More than 150 peptide therapeutics are globally in clinical development. Many enzymatic barriers should be crossed by a successful drug to be prosperous in such a process. Therefore, the new peptide drugs must be designed preventing the potential protease cleavage to make the compound less susceptible to protease reaction. We present a new data analysis tool developed in WebMetabase, an approach that stores the information from liquid chromatography mass spectrometry-based experimental data or from external sources such as the MEROPS database. The tool is a chemically aware system where each peptide substrate is presented as a sequence of structural blocks (SBs) connected by amide bonds and not being limited to the natural amino acids. Each SB is characterized by its pharmacophoric and physicochemical properties including a similarity score that describes likelihood between a SB and each one of the other SBs in the database. This methodology can be used to perform a frequency analysis to discover the most frequent cleavage sites for similar amide bonds, defined based on the similarity of the SB that participate in such a bond within the experimentally derived and/or public database.

AVAILABILITY AND IMPLEMENTATION

http://webmetabase.com:8182/WebMetabaseBioinformatics/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

All in One: Cavity Detection, Druggability Estimate, Cavity-Based Pharmacophore Perception, and Virtual Screening

Discovering the very first ligands of pharmacologically important targets in a fast and cost-efficient manner is an important issue in drug discovery. In the absence of structural information on either endogenous or synthetic ligands, computational chemists classically identify the very first hits by docking compound libraries to a binding site of interest, with well-known biases arising from the usage of scoring functions. We herewith propose a novel computational method tailored to ligand-free protein structures and consisting in the generation of simple cavity-based pharmacophores to which potential ligands could be aligned by the use of a smooth Gaussian function. The method, embedded in the IChem toolkit, automatically detects ligand-binding cavities, then predicts their structural druggability, and last creates a structure-based pharmacophore for predicted druggable binding sites. A companion tool (Shaper2) was designed to align ligands to cavity-derived pharmacophoric features. The proposed method is as efficient as state-of-the-art virtual screening methods (ROCS, Surflex-Dock) in both posing and virtual screening challenges. Interestingly, IChem-Shaper2 is clearly orthogonal to these latter methods in retrieving unique chemotypes from high-throughput virtual screening data.

Design, synthesis and pharmacological evaluation of N4,N6-disubstituted pyrimidine-4,6-diamine derivatives as potent EGFR inhibitors in non-small cell lung cancer

A novel series of 4, 6-disubstituted pyrimidines derivatives were designed, synthesized, and evaluated as epidermal growth factor receptor (EGFR) inhibitors for non-small cell lung cancer(NSCLC). 4, 6-disubstituted pyrimidines as core structure was utilized to substitute the lead structure AZD3759 of the quinazoline basic skeleton via an approach involving scaffold hopping. It was found that compound Yfq07 exhibited the best inhibitory effect compared with AZD3759 in vitro and in vivo: Yfq07 exhibited a competitive ATP inhibitory effect, multiple target effects, and further featured a stronger activity against H3255, A431, HCC827, PC-9 and H1975 compared to AZD3759. Moreover, a stronger pro-apoptotic effect, inhibition of cell G2/M phase on A431, H3255, HCC827 and H1975 could also be observed. In this study, the ultimate goal was changing the core structure to improve other epidermal growth factor receptor tyrosine kinase inhibitors (EGFR-TKIs) properties while retaining the overall potency. Yfq07 was further explored as an effective 4, 6-pyrimidine anticancer agent for the treatment of human NSCLC.

Structure-based design of small-molecule protein–protein interaction modulators: the story so far

As the pivotal role of protein–protein interactions in cell growth, transcriptional activity, intracellular trafficking, signal transduction and pathological conditions has been assessed, experimental and in silico strategies have been developed to design protein–protein interaction modulators. State-of-the-art structure-based design methods, mainly pharmacophore modeling and docking, which have succeeded in the identification of enzyme inhibitors, receptor agonists and antagonists, and new tools specifically conceived to target protein–protein interfaces (e.g., hot-spot and druggable pocket prediction methods) have been applied in the search for small-molecule protein–protein interaction modulators. Many successful applications of structure-based design approaches that, despite the challenge of targeting protein–protein interfaces with small molecules, have led to the identification of micromolar and submicromolar hits are reviewed here.

Probing the S1' site for the identification of non-zinc-binding MMP-2 inhibitors

Matrix metalloproteinases (MMPs) are zinc-dependent enzymes involved in several pathological states. Among them, MMP-2 is a relevant therapeutic target because of its role in cancer development and progression. Many MMP inhibitors (MMPIs) have been discovered over the last 30 years, and the majority of them contain a functional group that binds the zinc ion (zinc-binding group; ZBG). Unfortunately, no MMPIs have reached the market yet, owing to toxic effects due to unselective interactions of the ZBG. The new generation of MMPIs that do not bind the zinc ion could overcome problems of selectivity and toxicity, but have so far been developed only for MMP-8, -12, and -13. In this work, a virtual screening protocol was established by combining ligand- and structure-based methods to identify non-zinc-binding MMP-2 inhibitors using a new-generation MMP-8 inhibitor as a probe to find unexplored interactions in the MMP-2 S1' site. The screening allowed the identification of micromolar MMP-2 inhibitors that putatively avoid binding the zinc ion, as demonstrated by docking calculations. The LIA model, built to correlate predicted and experimental binding energies of the identified non-zinc-binding MMP-2 hits, underpins the reliability of the predicted docking poses.

From the protein's perspective: the benefits and challenges of protein structure-based pharmacophore modeling

Protein structure-based pharmacophore (SBP) models derive the molecular features a ligand must contain to be biologically active by conversion of protein properties to reciprocal ligand space. SBPs improve molecular understanding of ligand–protein interactions and can be used as valuable tools for hit and lead optimization, compound library design, and target hopping.

A novel approach for predicting P-glycoprotein (ABCB1) inhibition using molecular interaction fields

P-glycoprotein (Pgp or ABCB1) is an ABC transporter protein involved in intestinal absorption, drug metabolism, and brain penetration, and its inhibition can seriously alter a drug's bioavailability and safety. In addition, inhibitors of Pgp can be used to overcome multidrug resistance. Given this dual purpose, reliable in silico procedures to predict Pgp inhibition are of great interest. A large and accurate literature collection yielded more than 1200 structures; a model was then constructed using various molecular interaction field-based technologies, considering pharmacophoric features and those physicochemical properties related to membrane partitioning. High accuracy was demonstrated internally with two different validation sets and, moreover, using a number of molecules, for which Pgp inhibition was not experimentally available but was evaluated in-house. All of the validations confirmed the robustness of the model and its suitability to help medicinal chemists in drug discovery. The information derived from the model was rationalized as a pharmacophore for competitive Pgp inhibition.

FLAP: GRID molecular interaction fields in virtual screening. validation using the DUD data set

The performance of FLAP (Fingerprints for Ligands and Proteins) in virtual screening is assessed using a subset of the DUD (Directory of Useful Decoys) benchmarking data set containing 13 targets each with more than 15 different chemotype classes. A variety of ligand and receptor-based virtual screening approaches are examined, using combinations of individual templates 2D structures of known actives, a cocrystallized ligand, a receptor structure, or a cocrystallized ligand-biased receptor structure. We examine several data fusion approaches to combine the results of the individual virtual screens. In doing so, we show that excellent chemotype enrichment is achieved in both single target ligand-based and receptor-based approaches, of approximately 17-fold over random on average at a false positive rate of 1%. We also show that using as much starting knowledge as possible improves chemotype enrichment, and that data fusion using Pareto ranking is an effective method to do this giving up to 50% improvement in enrichment over the single methods. Finally we show that if inactivity or decoy data is incorporated, automatically training the scoring function in FLAP improves recovery still further, with almost 2-fold improvement over the enrichments shown by the single methods. The results clearly demonstrate the utility of FLAP for virtual screening when either a limited or wide range of prior knowledge is available.

Discovery of novel inhibitors of the NorA multidrug transporter of Staphylococcus aureus

Four novel inhibitors of the NorA efflux pump of Staphylococcus aureus, discovered through a virtual screening process, are reported. The four compounds belong to different chemical classes and were tested for their in vitro ability to block the efflux of a well-known NorA substrate, as well as for their ability to potentiate the effect of ciprofloxacin (CPX) on several strains of S. aureus, including a NorA overexpressing strain. Additionally, the MIC values of each of the compounds individually are reported. A structure-activity relationship study was also performed on these novel chemotypes, revealing three new compounds that are also potent NorA inhibitors. The virtual screening procedure employed FLAP, a new methodology based on GRID force field descriptors.

Modeling approaches for ligand-based 3D similarity

3D ligand-based similarity approaches are widely used in the early phases of drug discovery for tasks such as hit finding by virtual screening or compound design with quantitative structure–activity relationships. Here in we review widely used software for performing such tasks. Some techniques are based on relatively mature technology, shape-based similarity for instance. Typically, these methods remained in the realm of the expert user, the experienced modeler. However, advances in implementation and speed have improved usability and allow these methods to be applied to databases comprising millions of compounds. There are now many reports of such methods impacting drug-discovery projects. As such, the medicinal chemistry community has become the intended market for some of these new tools, yet they may consider the wide array and choice of approaches somewhat disconcerting. Each method has subtle differences and is better suited to certain tasks than others. In this article we review some of the widely used computational methods via application, provide straightforward background on the underlying theory and provide examples for the interested reader to pursue in more detail. In the new era of preclinical drug discovery there will be ever more pressure to move faster and more efficiently, and computational approaches based on 3D ligand similarity will play an increasing role in in this process.

Lead generation and optimization based on protein-ligand complementarity

This work proposes a computational procedure for structure-based lead generation and optimization, which relies on the complementarity of the protein-ligand interactions. This procedure takes as input the known structure of a protein-ligand complex. Retaining the positions of the ligand heavy atoms in the protein binding site it designs structurally similar compounds considering all possible combinations of atomic species (N, C, O, CH_3, NH, etc). Compounds are ranked based on a score which incorporates energetic contributions evaluated using molecular mechanics force fields. This procedure was used to design new inhibitor molecules for three serine/threonine protein kinases (p38 MAP kinase, p42 MAP kinase (ERK2), and c-Jun N-terminal kinase 3 (JNK3)). For each enzyme, the calculations produce a set of potential inhibitors whose scores are in agreement with IC50 data and Ki values. Furthermore, the native ligands for each protein target, scored within the five top-ranking compounds predicted by our method, one of the top-ranking compounds predicted to inhibit JNK3 was synthesized and his inhibitory activity confirmed against ATP hydrolysis. Our computational procedure is therefore deemed to be a useful tool for generating chemically diverse molecules active against known target proteins.

A similarity search using molecular topological graphs

A molecular similarity measure has been developed using molecular topological graphs and atomic partial charges. Two kinds of topological graphs were used. One is the ordinary adjacency matrix and the other is a matrix which represents the minimum path length between two atoms of the molecule. The ordinary adjacency matrix is suitable to compare the local structures of molecules such as functional groups, and the other matrix is suitable to compare the global structures of molecules. The combination of these two matrices gave a similarity measure. This method was applied to in silico drug screening, and the results showed that it was effective as a similarity measure.

Optimal assignment methods for ligand-based virtual screening

BACKGROUND

Ligand-based virtual screening experiments are an important task in the early drug discovery stage. An ambitious aim in each experiment is to disclose active structures based on new scaffolds. To perform these "scaffold-hoppings" for individual problems and targets, a plethora of different similarity methods based on diverse techniques were published in the last years. The optimal assignment approach on molecular graphs, a successful method in the field of quantitative structure-activity relationships, has not been tested as a ligand-based virtual screening method so far.

RESULTS

We evaluated two already published and two new optimal assignment methods on various data sets. To emphasize the "scaffold-hopping" ability, we used the information of chemotype clustering analyses in our evaluation metrics. Comparisons with literature results show an improved early recognition performance and comparable results over the complete data set. A new method based on two different assignment steps shows an increased "scaffold-hopping" behavior together with a good early recognition performance.

CONCLUSION

The presented methods show a good combination of chemotype discovery and enrichment of active structures. Additionally, the optimal assignment on molecular graphs has the advantage to investigate and interpret the mappings, allowing precise modifications of internal parameters of the similarity measure for specific targets. All methods have low computation times which make them applicable to screen large data sets.

Hot-spots-guided receptor-based pharmacophores (HS-Pharm): a knowledge-based approach to identify ligand-anchoring atoms in protein cavities and prioritize structure-based pharmacophores

The design of biologically active compounds from ligand-free protein structures using a structure-based approach is still a major challenge. In this paper, we present a fast knowledge-based approach (HS-Pharm) that allows the prioritization of cavity atoms that should be targeted for ligand binding, by training machine learning algorithms with atom-based fingerprints of known ligand-binding pockets. The knowledge of hot spots for ligand binding is here used for focusing structure-based pharmacophore models. Three targets of pharmacological interest (neuraminidase, beta2 adrenergic receptor, and cyclooxygenase-2) were used to test the evaluated methodology, and the derived structure-based pharmacophores were used in retrospective virtual screening studies. The current study shows that structure-based pharmacophore screening is a powerful technique for the fast identification of potential hits in a chemical library, and that it is a valid alternative to virtual screening by molecular docking.

Computational chemistry approaches to drug discovery in signal transduction

The advent of therapeutic strategies aimed at targeting specific macromolecular components of deregulated signaling pathways associated with particular disease states has given rise to the idea that it should be possible to design ligands as drug candidates to these targets from first principles. This concept has been beckoning for a long time but structure-based ligand design only became feasible once it was possible to determine the 3-D structures of molecular targets at atomic resolution. However, structure-based design turned out to be difficult, chiefly because under physiological conditions both receptors and ligands are not static but they behave dynamically. While it is possible to design ligands with high steric and electronic complementarity to a receptor site, it is always uncertain how biologically relevant the assumed conformations of both ligand and receptor actually are. The fact that it remains beyond our current abilities to predict with sufficient accuracy the affinity between hypothetical ligand and receptor poses is in part connected with this problem and continues to confound the reliable prediction of drug-like ligands for therapeutic targets. Nevertheless, significant progress has been made and so-called virtual screening methods that use computational methods to dock candidate ligands into receptor sites and to score the resulting complexes are now used routinely as one of the components in drug discovery screening campaigns. Here an overview is given of the underlying principles, implementations, and applications of structure-guided computational design technologies. Although the emphasis is on receptor-based strategies, mention will also be made of some of the more established ligand-based approaches, such as similarity analyses and quantitative structure-activity relationship methods.

ParaFrag--an approach for surface-based similarity comparison of molecular fragments

A frequent task in computer-aided drug design is to identify novel chemotypes similar in activity but structurally different to a given reference structure. Here we report the development of a novel method for atom-independent similarity comparison of molecular fragments (substructures of drug-like molecules). The fragments are characterized by their local surface properties coded in the form of 3D pharmacophores. As surface properties, we used the electrostatic potential (MEP), the local ionization energy (IE(L)), local electron affinity (EA(L)) and local polarizability (POL) calculated on isodensity surfaces. A molecular fragment can then be represented by a minimal set of extremes for each surface property. We defined a tolerance sphere for each of these extremes, thus allowing us to assess the similarity of fragments in an analogous manner to classical pharmacophore comparison. As a first application of this method we focused on comparing rigid fragments suitable for scaffold hopping. A retrospective analysis of successful scaffold hopping reported for Factor Xa inhibitors [Wood MR et al (2006) J Med Chem 49:1231] showed that our method performs well where atom-based similarity metrics fail.

New directions in library design and analysis

The high costs associated with high-throughput screening (HTS) coupled with the limited coverage and bias of current screening collections is such that diversity analysis continues to be an important criterion in lead generation. Whereas early approaches to diversity analysis were based on traditional descriptors such as two-dimensional fingerprints a recent emphasis has been on assessing scaffold coverage to ensure that a variety of different chemotypes are represented. Moreover, whether designing diverse or focused libraries, it is widely recognised that designs should aim to achieve a balance in a number of different properties and multiobjective optimisation provides an effective way of achieving such designs.

In-silico drug screening method based on the protein–compound affinity matrix using the factor selection technique

We have developed a new in-silico drug screening method, a modified version of a docking score index (DSI) method, based on a protein-compound docking affinity matrix. By using this method, the docking scores are converted to the docking score indexes by the principal component analysis (PCA) method and each compound is projected into a PCA space. In this study, we propose a method to select a set of suitable principal component axes and evaluate the database enrichment for 12 target proteins. This method selects the new active compounds or hits, which are close to the known active compounds, thereby enhancing the database enrichment.

SHOP: scaffold HOPping by GRID-based similarity searches

A new GRID-based method for scaffold hopping (SHOP) is presented. In a fully automatic manner, scaffolds were identified in a database based on three types of 3D-descriptors. SHOP's ability to recover scaffolds was assessed and validated by searching a database spiked with fragments of known ligands of three different protein targets relevant for drug discovery using a rational approach based on statistical experimental design. Five out of eight and seven out of eight thrombin scaffolds and all seven HIV protease scaffolds were recovered within the top 10 and 31 out of 31 neuraminidase scaffolds were in the 31 top-ranked scaffolds. SHOP also identified new scaffolds with substantially different chemotypes from the queries. Docking analysis indicated that the new scaffolds would have similar binding modes to those of the respective query scaffolds observed in X-ray structures. The databases contained scaffolds from published combinatorial libraries to ensure that identified scaffolds could be feasibly synthesized.

Comparison of shape-matching and docking as virtual screening tools

Ligand docking is a widely used approach in virtual screening. In recent years a large number of publications have appeared in which docking tools are compared and evaluated for their effectiveness in virtual screening against a wide variety of protein targets. These studies have shown that the effectiveness of docking in virtual screening is highly variable due to a large number of possible confounding factors. Another class of method that has shown promise in virtual screening is the shape-based, ligand-centric approach. Several direct comparisons of docking with the shape-based tool ROCS have been conducted using data sets from some of these recent docking publications. The results show that a shape-based, ligand-centric approach is more consistent than, and often superior to, the protein-centric approach taken by docking.

Chapter 9 Molecular Similarity: Advances in Methods, Applications and Validations in Virtual Screening and QSAR

This chapter discusses recent developments in some of the areas that exploit the molecular similarity principle, novel approaches to capture molecular properties by the use of novel descriptors, focuses on a crucial aspect of computational models-their validity, and discusses additional ways to examine data available, such as those from high-throughput screening (HTS) campaigns and to gain more knowledge from this data. The chapter also presents some of the recent applications of methods discussed focusing on the successes of virtual screening applications, database clustering and comparisons (such as drug- and in-house-likeness), and the recent large-scale validations of docking and scoring programs. While a great number of descriptors and modeling methods has been proposed until today, the recent trend toward proper model validation is very much appreciated. Although some of their limitations are surely because of underlying principles and limitations of fundamental concepts, others will certainly be eliminated in the future.

Unsupervised 3D Ring Template Searching as an Ideas Generator for Scaffold Hopping: Use of the LAMDA, RigFit, and Field-Based Similarity Search (FBSS) Methods

Crystal structures taken from the Cambridge Structural Database were used to build a ring scaffold database containing 19 050 3D structures, with each such scaffold then being used to generate a centroid connecting path (CCP) representation. The CCP is a novel object that connects ring centroids, ring linker atoms, and other important points on the connection path between ring centroids. Unsupervised searching in the scaffold and CCP data sets was carried out using the atom-based LAMDA and RigFit search methods and the field-based similarity search method. The performance of these methods was tested with three different ring scaffold queries. These searches demonstrated that unsupervised 3D scaffold searching methods can find not only the types of ring systems that might be retrieved in carefully defined pharmacophore searches (supervised approach) but also additional, structurally diverse ring systems that could form the starting point for lead discovery programs or other scaffold-hopping applications. Not only are the methods effective but some are sufficiently rapid to permit scaffold searching in large chemical databases on a routine basis.