SHAFTS: a hybrid approach for 3D molecular similarity calculation. 1. Method and assessment of virtual screening

We developed a novel approach called SHAFTS (SHApe-FeaTure Similarity) for 3D molecular similarity calculation and ligand-based virtual screening. SHAFTS adopts a hybrid similarity metric combined with molecular shape and colored (labeled) chemistry groups annotated by pharmacophore features for 3D similarity calculation and ranking, which is designed to integrate the strength of pharmacophore matching and volumetric overlay approaches. A feature triplet hashing method is used for fast molecular alignment poses enumeration, and the optimal superposition between the target and the query molecules can be prioritized by calculating corresponding "hybrid similarities". SHAFTS is suitable for large-scale virtual screening with single or multiple bioactive compounds as the query "templates" regardless of whether corresponding experimentally determined conformations are available. Two public test sets (DUD and Jain's sets) including active and decoy molecules from a panel of useful drug targets were adopted to evaluate the virtual screening performance. SHAFTS outperformed several other widely used virtual screening methods in terms of enrichment of known active compounds as well as novel chemotypes, thereby indicating its robustness in hit compounds identification and potential of scaffold hopping in virtual screening.     

Combining ethnopharmacology and virtual screening for lead structure discovery: COX-inhibitors as application example

Virtual screening of large libraries of organic compounds combined with pharmacological high throughput screening is widely used for drug discovery in the pharmaceutical industry. Our aim was to explore the efficiency of using a biased 3D database comprising secondary metabolites from antiinflammatory medicinal plants as a source for the virtual screening. For this study pharmacophore models of cyclooxygenase I and II (COX-1, COX-2), key enzymes in the inflammation process, were generated with structure-based as well as common feature based modeling, resulting in three COX hypotheses. Four different multiconfomational 3D databases limited in molecular weight between 300 and 700 Da were applied to the screening in order to compare and analyze the obtained hit rates. Two of them were created in-house (DIOS, NPD). The database DIOS consists of 2752 compounds from phytochemical reports of antiinflammatory medicinal plants described by the ethnopharmacological source 'De material medica' of Pedanius Dioscorides, whereas NPD contains almost 80,000 compounds gathered arbitrarily from natural sources. In addition, two available multiconformational 3D libraries comprising marketed and development drug substances (DWI and NCI), mainly originating from synthesis, were used for comparison. As a test of the pharmacophore models' capability in natural sources, the models were used to search for known COX inhibitory natural products. This was achieved with some exceptions, which are discussed in the paper. Depending on the hypothesis used, DWI and NCI library searches produced hit rates in the range of 6.6% to 13.7%. A slight increase of the number of molecules assessed for binding was achieved with the database of natural products (NPD). Using the biased 3D database DIOS, however, the average increase of efficiency reached 77% to 133% compared to the hit rates resulting from WDI and NCI. The statistical benefit of a combination of an ethnopharmacological approach with the potential of computer aided drug discovery by in silico screening was demonstrated exemplified on the applied targets COX-1 and COX-2.     

A pseudo-ligand approach to virtual screening

A virtual screening method is presented that is grounded on a receptor-derived pharmacophore model termed "virtual ligand" or "pseudo-ligand". The model represents an idealized constellation of potential ligand sites that interact with residues of the binding pocket. For rapid virtual screening of compound libraries the potential pharmacophore points of the virtual ligand are encoded as an alignment-free correlation vector, avoiding spatial alignment of pharmacophore features between the pharmacophore query (i.e., the virtual ligand) and the candidate molecule. The method was successfully applied to retrieving factor Xa inhibitors from a Ugi three-component combinatorial library, and yielded high enrichment of actives in a retrospective search for cyclooxygenase-2 (COX-2) inhibitors. The approach provides a concept for "de-orphanizing" potential drug targets and identifying ligands for hitherto unexplored or allosteric binding pockets.     

The discovery of Kv1.5 blockers as a case study for the application of virtual screening approaches

Different virtual screening techniques are available as alternatives to high throughput screening. These different techniques have been rarely used together on the same target. We had the opportunity to do so in order to discover novel blockers of the voltage-dependent potassium channel Kv1.5, a potential target for the treatment of atrial fibrillation. Our corporate database was searched, using a protein-based pharmacophore, derived from a homology model, as query. As a result, 244 molecules were screened in vitro, 19 of them (7.8%) were found to be active. Five of them, belonging to five different chemical classes, exhibited IC50 values under 10 microM. The performance of this structure-based virtual screening protocol has been compared with those of similarity and ligand-based pharmacophore searches. The analysis of the results supports the conventional wisdom of using as many virtual screening techniques as possible in order to maximize the chance of finding as many chemotypes as possible.     

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.     

Automated pharmacophore query optimization with genetic algorithms - a case study using the MC4R system

Due to the recent availability of high quality small molecule databases, such as ZINC and PubChem,1,2 virtual screening is playing an even more important role in identifying biologically relevant molecules in drug discovery campaigns. The success of pharmacophore-based virtual screening (PBVS) relies largely on the accuracy and specificity of the pharmacophore query employed. Deriving a pharmacophore query from a single structure inevitably introduces uncertainty, and the derived query is unlikely to be optimal against every collection of input compounds, especially when it is desired to discriminate among compounds with similar chemical structures. In this study, we present an optimization approach empowered by genetic algorithms (GA) to enhance the accuracy and specificity of a primary pharmacophore query. The example utilized is the human melanocortin type 4 receptor (hMC4R), for which the pharmacophore query was built on the basis of the structure of a rigid cyclic peptide agonist.(3) The optimized query is shown to be capable of identifying 37 positive hMC4R agonists with no false positives from a training set containing 55 agonists and 51 nonagonists. This represents a significant improvement from the initial query which exhibited a 37/32 hit rate. The final, optimized query is challenged with a testing set comprising of 55 hMC4R agonists and 50 nonagonists and achieves a hit rate of 33/8, that improved from 40/31. The impact of GA controlling parameters, including mutation rate, crossover rate, fitness function, population size, and convergence criterion, on performance of optimization are examined and discussed.     

Novel approach to structure-based pharmacophore search using computational geometry and shape matching techniques

Computationally efficient structure-based virtual screening methods have recently been reported that seek to find effective means to utilize experimental structure information without employing detailed molecular docking calculations. These tools can be coupled with efficient experimental screening technologies to improve the probability of identifying hits and leads for drug discovery research. Commercial software ROCS (rapid overlay of chemical structures) from Open Eye Scientific is such an example, which is a shape-based virtual screening method using the 3D structure of a ligand, typically from a bound X-ray costructure, as the query. We report here the development of a new structure-based pharmacophore search method (called Shape4) for virtual screening. This method adopts a variant of the ROCS shape technology and expands its use to work with an empty crystal structure. It employs a rigorous computational geometry method and a deterministic geometric casting algorithm to derive the negative image (i.e., pseudoligand) of a target binding site. Once the negative image (or pseudoligand) is generated, an efficient shape comparison algorithm in the commercial OE SHAPE Toolkit is adopted to compare and match small organic molecules with the shape of the pseudoligand. We report the detailed computational protocol and its computational validation using known biologically active compounds extracted from the WOMBAT database. Models derived for five selected targets were used to perform the virtual screening experiments to obtain the enrichment data for various virtual screening methods. It was found that our approach afforded similar or better enrichment ratios than other related methods, often with better diversity among the top ranking computational hits.     

血必净抗炎作用药效物质基础和多靶点作用效应

研究了血必净抗炎作用的药效物质基础, 并在分子层面上对复方多靶点作用进行阐释. 借助于计算机辅助药物设计技术, 构建血必净化学成分数据库, 综合应用同源模建、分子对接、药效团模型、数据库搜索等方法, 探讨其与炎症靶点5-脂氧合酶(5-LOX)、环氧合酶-2(COX-2)、IKK-2受体的相互作用关系. 血必净化学成分中与靶点5-LOX、COX-2、IKK-2结合效应较好的分别有30、36、8个分子; 有16个分子与2个或3个靶点存在作用, 其中15个分子对靶点5-LOX和COX-2有抑制作用, 迷迭香酸对3个靶点均有作用. 从分子层次上阐释了复方血必净抗炎作用的药效物质基础和多靶点作用效应, 为血必净复方的临床应用提供了科学依据; 同时, 也为寻找新型抗炎药物提供一定的参考和借鉴.     

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.     

Comparative analysis of pharmacophore screening tools

The pharmacophore concept is of central importance in computer-aided drug design (CADD) mainly because of its successful application in medicinal chemistry and, in particular, high-throughput virtual screening (HTVS). The simplicity of the pharmacophore definition enables the complexity of molecular interactions between ligand and receptor to be reduced to a handful set of features. With many pharmacophore screening softwares available, it is of the utmost interest to explore the behavior of these tools when applied to different biological systems. In this work, we present a comparative analysis of eight pharmacophore screening algorithms (Catalyst, Unity, LigandScout, Phase, Pharao, MOE, Pharmer, and POT) for their use in typical HTVS campaigns against four different biological targets by using default settings. The results herein presented show how the performance of each pharmacophore screening tool might be specifically related to factors such as the characteristics of the binding pocket, the use of specific pharmacophore features, and the use of these techniques in specific steps/contexts of the drug discovery pipeline. Algorithms with rmsd-based scoring functions are able to predict more compound poses correctly as overlay-based scoring functions. However, the ratio of correctly predicted compound poses versus incorrectly predicted poses is better for overlay-based scoring functions that also ensure better performances in compound library enrichments. While the ensemble of these observations can be used to choose the most appropriate class of algorithm for specific virtual screening projects, we remarked that pharmacophore algorithms are often equally good, and in this respect, we also analyzed how pharmacophore algorithms can be combined together in order to increase the success of hit compound identification. This study provides a valuable benchmark set for further developments in the field of pharmacophore search algorithms, e.g., by using pose predictions and compound library enrichment criteria.     

Using halogen bonds to address the protein backbone: a systematic evaluation

Halogen bonds are specific embodiments of the sigma hole bonding paradigm. They represent directional interactions between the halogens chlorine, bromine, or iodine and an electron donor as binding partner. Using quantum chemical calculations at the MP2 level, we systematically explore how they can be used in molecular design to address the omnipresent carbonyls of the protein backbone. We characterize energetics and directionality and elucidate their spatial variability in sub-optimal geometries that are expected to occur in protein-ligand complexes featuring a multitude of concomitant interactions. By deriving simple rules, we aid medicinal chemists and chemical biologists in easily exploiting them for scaffold decoration and design. Our work shows that carbonyl-halogen bonds may be used to expand the patentable medicinal chemistry space, redefining halogens as key features. Furthermore, this data will be useful for implementing halogen bonds into pharmacophore models or scoring functions making the QM information available for automatic molecular recognition in virtual high throughput screening.     

Virtual screening for compounds that mimic protein-protein interface epitopes

Modulation of protein-protein interactions (PPI) has emerged as a new concept in rational drug design. Here, we present a computational protocol for identifying potential PPI inhibitors. Relevant regions of interfaces (epitopes) are predicted for three-dimensional protein models and serve as queries for virtual compound screening. We present a computational screening protocol that incorporates two different pharmacophore models. One model is based on the mathematical concept of autocorrelation vectors and the other utilizes fuzzy labeled graphs. In a proof-of-concept study, we were able to identify serine protease inhibitors using a predicted trypsin epitope as query. Our virtual screening framework may be suited for rapid identification of PPI inhibitors and suggesting bioactive tool compounds.