State-of-the-art in ligand-based virtual screening

HybridDock: A Hybrid Protein-Ligand Docking Protocol Integrating Protein- and Ligand-Based Approaches

Structure-based molecular docking and ligand-based similarity search are two commonly used computational methods in computer-aided drug design. Structure-based docking tries to utilize the structural information on a drug target like protein, and ligand-based screening takes advantage of the information on known ligands for a target. Given their different advantages, it would be desirable to use both protein- and ligand-based approaches in drug discovery when information for both the protein and known ligands is available. Here, we have presented a general hybrid docking protocol, referred to as HybridDock, to utilize both the protein structures and known ligands by combining the molecular docking program MDock and the ligand-based similarity search method SHAFTS, and evaluated our hybrid docking protocol on the CSAR 2013 and 2014 exercises. The results showed that overall our hybrid docking protocol significantly improved the performance in both binding affinity and binding mode predictions, compared to the sole MDock program. The efficacy of the hybrid docking protocol was further confirmed using the combination of DOCK and SHAFTS, suggesting an alternative docking approach for modern drug design/discovery.

State-of-the-art strategies for targeting protein-protein interactions by small-molecule inhibitors

Targeting protein-protein interactions (PPIs) has emerged as a viable approach in modern drug discovery. However, the identification of small molecules enabling us to effectively interrupt their interactions presents significant challenges. In the recent past, significant advances have been made in the development of new biological and chemical strategies to facilitate the discovery process of small-molecule PPI inhibitors. This review aims to highlight the state-of-the-art technologies and the achievements made recently in this field. The "hot spots" of PPIs have been proved to be critical for small molecules to bind. Three strategies including screening, designing, and synthetic approaches have been explored for discovering PPI inhibitors by targeting the "hot spots". Although the classic high throughput screening approach can be used, fragment screening, fragment-based drug design and newly improved virtual screening are demonstrated to be more effective in the discovery of PPI inhibitors. In addition to screening approaches, design strategies including anchor-based and small molecule mimetics of secondary structures involved in PPIs have become powerful tools as well. Finally, constructing new chemically spaced libraries with high diversity and complexity is becoming an important area of interest for PPI inhibitors. The successful cases from the recent five year studies are used to illustrate how these approaches are implemented to uncover and optimize small molecule PPI inhibitors and notably some of them have become promising therapeutics.

Molecular docking and structure-based drug design strategies

Pharmaceutical research has successfully incorporated a wealth of molecular modeling methods, within a variety of drug discovery programs, to study complex biological and chemical systems. The integration of computational and experimental strategies has been of great value in the identification and development of novel promising compounds. Broadly used in modern drug design, molecular docking methods explore the ligand conformations adopted within the binding sites of macromolecular targets. This approach also estimates the ligand-receptor binding free energy by evaluating critical phenomena involved in the intermolecular recognition process. Today, as a variety of docking algorithms are available, an understanding of the advantages and limitations of each method is of fundamental importance in the development of effective strategies and the generation of relevant results. The purpose of this review is to examine current molecular docking strategies used in drug discovery and medicinal chemistry, exploring the advances in the field and the role played by the integration of structure- and ligand-based methods.

Virtual Screening for HIV Protease Inhibitors Using a Novel Database Filtering Procedure

A virtual screening to find novel inhibitors for HIV protease was performed on the ZINC database.1 A critical part in virtual screening and associated techniques is preliminary database filtering and size reduction and for that purpose a novel feature matrix matching procedure was used. The reduction of ∼14 million available ligands to a subset of 14299 ligands was achieved with a structure based approach where the analysis of the 3D structure of the active site of the protease produced a graph with hydrogen bond donor, hydrogen bond acceptor and hydrophobic subsites represented as graph nodes. A similar treatment was also applied to the compound database content and the comparison of binding site and ligand graphs was used to preselect potentially active ligands. The resulting set was further subjected to docking. The algorithm used was able to find several novel as well as previously known and experimentally tested ligands, demonstrating the validity of the approach.

[Chemical libraries dedicated to protein-protein interactions]

The identification of complete networks of protein-protein interactions (PPI) within a cell has contributed to major breakthroughs in understanding biological pathways, host-pathogen interactions and cancer development. As a consequence, PPI have emerged as a new class of promising therapeutic targets. However, they are still considered as a challenging class of targets for drug discovery programs. Recent successes have allowed the characterization of structural and physicochemical properties of protein-protein interfaces leading to a better understanding of how they can be disrupted with small molecule compounds. In addition, characterization of the profiles of PPI inhibitors has allowed the development of PPI-focused libraries. In this review, we present the current efforts at developing chemical libraries dedicated to these innovative targets.

Discovery of TRAF-2 and NCK-interacting kinase (TNIK) inhibitors by ligand-based virtual screening methods

TRAF-2 and NCK-interacting kinase (TNIK) is a serine–threonine kinase with a proposed role in Wnt/β-catenin and JNK pathways.

Practical Considerations in Virtual Screening and Molecular Docking

Molecular docking has become an important common component of the drug discovery toolbox, and its relative low-cost implications and perceived simplicity of use has stimulated an everincreasing popularity within academic communities. The inherent “garbage-in-garbage-out” defect of molecular docking, however, leads a lot of researchers to dedicate countless hours to the identification of hit compounds that later prove to be inactive. Several considerations that can greatly improve the success and enrichment of true bioactive hit compounds are commonly overlooked at the initial stages of a molecular docking study. This chapter will cover several of these considerations, including protonation states, active site waters, separating actives from decoys, consensus docking and molecular mechanics generalized-Born/surface area (MM-GBSA) rescoring, and incorporation of pharmacophoric constraints, in an attempt to clarify what is, in fact, very complicated and inherent difficulties of a structure-based drug design study.

Benchmarking methods and data sets for ligand enrichment assessment in virtual screening

Retrospective small-scale virtual screening (VS) based on benchmarking data sets has been widely used to estimate ligand enrichments of VS approaches in the prospective (i.e. real-world) efforts. However, the intrinsic differences of benchmarking sets to the real screening chemical libraries can cause biased assessment. Herein, we summarize the history of benchmarking methods as well as data sets and highlight three main types of biases found in benchmarking sets, i.e. "analogue bias", "artificial enrichment" and "false negative". In addition, we introduce our recent algorithm to build maximum-unbiased benchmarking sets applicable to both ligand-based and structure-based VS approaches, and its implementations to three important human histone deacetylases (HDACs) isoforms, i.e. HDAC1, HDAC6 and HDAC8. The leave-one-out cross-validation (LOO CV) demonstrates that the benchmarking sets built by our algorithm are maximum-unbiased as measured by property matching, ROC curves and AUCs.

Scaffold hopping with virtual screening from IP3 to a drug-like partial agonist of the inositol trisphosphate receptor

Inositol 1,4,5-trisphosphate (IP3 ) is a universal signalling molecule that releases calcium from stores within cells by activating the IP3 receptor. Although chemical tools that modulate the IP3 receptor exist, none is ideal due to trade offs between potency, selectivity and cell permeability, and their chemical properties make them challenging starting points for optimisation. Therefore, to find new leads, we used virtual screening to scaffold hop from IP3 by using the program ROCS to perform a 3D ligand-based screen of the ZINC database of purchasable compounds. We then used the program FRED to dock the top-ranking hits into the IP3 binding pocket of the receptor. We tested the 12 highest-scoring hits in a calcium-release bioassay and identified SI-9 as a partial agonist. SI-9 competed with [(3) H]IP3 binding, and reduced histamine-induced calcium signalling in HeLa cells. SI-9 has a novel 2D scaffold that represents a tractable lead for designing improved IP3 receptor modulators.

Pharmacophore Modeling, Ensemble Docking, Virtual Screening, and Biological Evaluation on Glycogen Synthase Kinase-3β

Glycogen synthase kinase-3 (GSK-3) is a multifunctional serine/threonine protein kinase which is engaged in a variety of signaling pathways, regulating a wide range of cellular processes. GSK-3β, also known as tau protein kinase I (TPK-I), is one of the most important kinases implicated in the hyperphosphorylation of tau that leads to neurodegenerative diseases. Hence, GSK-3β has emerged as an important therapeutic target. To identify compounds that are structurally novel and diverse compared to previously reported ATP-competitive GSK-3β inhibitors, we performed virtual screening by implementing a mixed ligand/structure-based approach, which included pharmacophore modeling, diversity analysis, and ensemble docking. The sensitivities of different docking protocols to induced-fit effects were explored. An enrichment study was employed to verify the robustness of ensemble docking, using 13 X-ray structures of GSK-3β, compared to individual docking in terms of retrieving active compounds from a decoy dataset. A total of 24 structurally diverse compounds obtained from the virtual screening underwent biological validation. The bioassay results showed that 15 out of the 24 hit compounds are indeed GSK-3β inhibitors, and among them, one compound exhibiting sub-micromolar inhibitory activity is a reasonable starting point for further optimization.

Ligand efficiency based approach for efficient virtual screening of compound libraries

Here we report for the first time the use of fit quality (FQ), a ligand efficiency (LE) based measure for virtual screening (VS) of compound libraries. The LE based VS protocol was used to screen an in-house database of 125,000 compounds to identify aurora kinase A inhibitors. First, 20 known aurora kinase inhibitors were docked to aurora kinase A crystal structure (PDB ID: 2W1C); and the conformations of docked ligand were used to create a pharmacophore (PH) model. The PH model was used to screen the database compounds, and rank (PH rank) them based on the predicted IC50 values. Next, LE_Scale, a weight-dependant LE function, was derived from 294 known aurora kinase inhibitors. Using the fit quality (FQ = LE/LE_Scale) score derived from the LE_Scale function, the database compounds were reranked (PH_FQ rank) and the top 151 (0.12% of database) compounds were assessed for aurora kinase A inhibition biochemically. This VS protocol led to the identification of 7 novel hits, with compound 5 showing aurora kinase A IC50 = 1.29 μM. Furthermore, testing of 5 against a panel of 31 kinase reveals that it is selective toward aurora kinase A & B, with <50% inhibition for other kinases at 10 μM concentrations and is a suitable candidate for further development. Incorporation of FQ score in the VS protocol not only helped identify a novel aurora kinase inhibitor, 5, but also increased the hit rate of the VS protocol by improving the enrichment factor (EF) for FQ based screening (EF = 828), compared to PH based screening (EF = 237) alone. The LE based VS protocol disclosed here could be applied to other targets for hit identification in an efficient manner.

Novel colchicine-site binders with a cyclohexanedione scaffold identified through a ligand-based virtual screening approach

Vascular disrupting agents (VDAs) constitute an innovative anticancer therapy that targets the tumor endothelium, leading to tumor necrosis. Our approach for the identification of new VDAs has relied on a ligand 3-D shape similarity virtual screening (VS) approach using the ROCS program as the VS tool and as query colchicine and TN-16, which both bind the α,β-tubulin dimer. One of the hits identified, using TN-16 as query, has been explored by the synthesis of its structural analogues, leading to 2-(1-((2-methoxyphenyl)amino)ethylidene)-5-phenylcyclohexane-1,3-dione (compound 16c) with an IC50 = 0.09 ± 0.01 μM in HMEC-1 and BAEC, being 100-fold more potent than the initial hit. Compound 16c caused cell cycle arrest in the G2/M phase and interacted with the colchicine-binding site in tubulin, as confirmed by a competition assay with N,N'-ethylenebis(iodoacetamide) and by fluorescence spectroscopy. Moreover, 16c destroyed an established endothelial tubular network at 1 μM and inhibited the migration and invasion of human breast carcinoma cells at 0.4 μM. In conclusion, our approach has led to a new chemotype of promising antiproliferative compounds with antimitotic and potential VDA properties.

Condorcet and borda count fusion method for ligand-based virtual screening

Background

It is known that any individual similarity measure will not always give the best recall of active molecule structure for all types of activity classes. Recently, the effectiveness of ligand-based virtual screening approaches can be enhanced by using data fusion. Data fusion can be implemented using two different approaches: group fusion and similarity fusion. Similarity fusion involves searching using multiple similarity measures. The similarity scores, or ranking, for each similarity measure are combined to obtain the final ranking of the compounds in the database.

Results

The Condorcet fusion method was examined. This approach combines the outputs of similarity searches from eleven association and distance similarity coefficients, and then the winner measure for each class of molecules, based on Condorcet fusion, was chosen to be the best method of searching. The recall of retrieved active molecules at top 5% and significant test are used to evaluate our proposed method. The MDL drug data report (MDDR), maximum unbiased validation (MUV) and Directory of Useful Decoys (DUD) data sets were used for experiments and were represented by 2D fingerprints.

Conclusions

Simulated virtual screening experiments with the standard two data sets show that the use of Condorcet fusion provides a very simple way of improving the ligand-based virtual screening, especially when the active molecules being sought have a lowest degree of structural heterogeneity. However, the effectiveness of the Condorcet fusion was increased slightly when structural sets of high diversity activities were being sought.

Focused chemical libraries--design and enrichment: an example of protein-protein interaction chemical space

One of the many obstacles in the development of new drugs lies in the limited number of therapeutic targets and in the quality of screening collections of compounds. In this review, we present general strategies for building target-focused chemical libraries with a particular emphasis on protein-protein interactions (PPIs). We describe the chemical spaces spanned by nine commercially available PPI-focused libraries and compare them to our 2P2I3D academic library, dedicated to orthosteric PPI modulators. We show that although PPI-focused libraries have been designed using different strategies, they share common subspaces. PPI inhibitors are larger and more hydrophobic than standard drugs; however, an effort has been made to improve the drug-likeness of focused chemical libraries dedicated to this challenging class of targets.

The use of 2D fingerprint methods to support the assessment of structural similarity in orphan drug legislation

BACKGROUND

In the European Union, medicines are authorised for some rare disease only if they are judged to be dissimilar to authorised orphan drugs for that disease. This paper describes the use of 2D fingerprints to show the extent of the relationship between computed levels of structural similarity for pairs of molecules and expert judgments of the similarities of those pairs. The resulting relationship can be used to provide input to the assessment of new active compounds for which orphan drug authorisation is being sought.

RESULTS

143 experts provided judgments of the similarity or dissimilarity of 100 pairs of drug-like molecules from the DrugBank 3.0 database. The similarities of these pairs were also computed using BCI, Daylight, ECFC4, ECFP4, MDL and Unity 2D fingerprints. Logistic regression analyses demonstrated a strong relationship between the human and computed similarity assessments, with the resulting regression models having significant predictive power in experiments using data from submissions of orphan drug medicines to the European Medicines Agency. The BCI fingerprints performed best overall on the DrugBank dataset while the BCI, Daylight, ECFP4 and Unity fingerprints performed comparably on the European Medicines Agency dataset.

CONCLUSIONS

Measures of structural similarity based on 2D fingerprints can provide a useful source of information for the assessment of orphan drug status by regulatory authorities.

SCISSORS: practical considerations

Molecular similarity has been effectively applied to many problems in cheminformatics and computational drug discovery, but modern methods can be prohibitively expensive for large-scale applications. The SCISSORS method rapidly approximates measures of pairwise molecular similarity such as ROCS and LINGO Tanimotos, acting as a filter to quickly reduce the size of a problem. We report an in-depth analysis of SCISSORS performance, including a mapping of the SCISSORS error distribution, benchmarking, and investigation of several algorithmic modifications. We show that SCISSORS can accurately predict multiconformer similarity and suggest a method for estimating optimal SCISSORS parameters in a data set-specific manner. These results are a useful resource for researchers seeking to incorporate SCISSORS into molecular similarity applications.

Receptor-ligand molecular docking

Docking methodology aims to predict the experimental binding modes and affinities of small molecules within the binding site of particular receptor targets and is currently used as a standard computational tool in drug design for lead compound optimisation and in virtual screening studies to find novel biologically active molecules. The basic tools of a docking methodology include a search algorithm and an energy scoring function for generating and evaluating ligand poses. In this review, we present the search algorithms and scoring functions most commonly used in current molecular docking methods that focus on protein-ligand applications. We summarise the main topics and recent computational and methodological advances in protein-ligand docking. Protein flexibility, multiple ligand binding modes and the free-energy landscape profile for binding affinity prediction are important and interconnected challenges to be overcome by further methodological developments in the docking field.

Molecular crime scene investigation - dusting for fingerprints

In chemoinformatics and drug design, fingerprints (FPs) are defined as string representations of molecular structure and properties and are popular descriptors for similarity searching. FPs are generally characterized by the simplicity of their design and ease of use. Despite a long history in chemoinformatics, the potential and limitations of FP searching are often not well under- stood. Standard FPs can also be subjected to engineering techniques to tune them for specific search applications.

Integrating virtual and biochemical screening for protein tyrosine phosphatase inhibitor discovery

Protein tyrosine phosphatases (PTPs) represent an important class of enzymes that mediate signal transduction and control diverse aspects of cell behavior. The importance of their activity is exemplified by their significant contribution to disease etiology with over half of all human PTP genes implicated in at least one disease. Small molecule inhibitors targeting individual PTPs are important biological tools, and are needed to fully characterize the function of these enzymes. Moreover, potent and selective PTP inhibitors hold the promise to transform the treatment of many diseases. While numerous methods exist to develop PTP-directed small molecules, we have found that complimentary use of both virtual (in silico) and biochemical (in vitro) screening approaches expedite compound identification and drug development. Here, we summarize methods pertinent to our work and others. Focusing on specific challenges and successes we have experienced, we discuss the considerable caution that must be taken to avoid enrichment of inhibitors that function by non-selective oxidation. We also discuss the utility of using "open" PTP structures to identify active-site directed compounds, a rather unconventional choice for virtual screening. When integrated closely, virtual and biochemical screening can be used in a productive workflow to identify small molecules targeting PTPs.

Hit identification and optimization in virtual screening: practical recommendations based on a critical literature analysis

A critical analysis of virtual screening results published between 2007 and 2011 was performed. The activity of reported hit compounds from over 400 studies was compared to their hit identification criteria. Hit rates and ligand efficiencies were calculated to assist in these analyses, and the results were compared with factors such as the size of the virtual library and the number of compounds tested. A series of promiscuity, druglike, and ADMET filters were applied to the reported hits to assess the quality of compounds reported, and a careful analysis of a subset of the studies that presented hit optimization was performed. These data allowed us to make several practical recommendations with respect to selection of compounds for experimental testing, definition of hit identification criteria, and general virtual screening hit criteria to allow for realistic hit optimization. A key recommendation is the use of size-targeted ligand efficiency values as hit identification criteria.

The future of virtual compound screening

We provide a future perspective of the virtual screening field. A number of challenges will be highlighted that virtual screening will likely face when compound data will further grow at or beyond current rates and when much more target information will become available. These challenges go beyond computational efficiency issues (that will of course also play a critical role). For example, for structure-based approaches, the accuracy of scoring functions and energy calculations will need to be improved. For ligand-based approaches, the compound class-dependence of similarity methods needs to be further explored and relationships between molecular similarity and activity similarity need to be established. We also comment on the current and future value of virtual screening. Opportunities for further development in a postgenome era are also discussed. It is hoped that some of the views and hypotheses we articulate might stimulate further discussion about the virtual screening field going forward.

Shifting from the single to the multitarget paradigm in drug discovery

Increasing evidence that several drug compounds exert their effects through interactions with multiple targets is boosting the development of research fields that challenge the data reductionism approach. In this article, we review and discuss the concepts of drug repurposing, polypharmacology, chemogenomics, phenotypic screening and high-throughput in vivo testing of mixture-based libraries in an integrated manner. These research fields offer alternatives to the current paradigm of drug discovery, from a one target-one drug model to a multiple-target approach. Furthermore, the goals of lead identification are being expanded accordingly to identify not only 'key' compounds that fit with a single-target 'lock', but also 'master key' compounds that favorably interact with multiple targets (i.e. operate a set of desired locks to gain access to the expected clinical effects).

Fingerprint design and engineering strategies: rationalizing and improving similarity search performance

Fingerprints (FPs) are bit or integer string representations of molecular structure and properties, and are popular descriptors for chemical similarity searching. A major goal of similarity searching is the identification of novel active compounds on the basis of known reference molecules. In this review recent FP design and engineering strategies are discussed. New types of FPs continue to be replaced, often applying different design principles. FP engineering techniques have recently been introduced to further improve search performance and computational efficiency and elucidate mechanisms by which FPs recognize active compounds. In addition, through feature selection and hybridization techniques, standard FPs have been transformed into compound class-specific versions with further increased search performance. Moreover, scaffold hopping mechanisms have been explored. FPs will continue to play an important role in the search for novel active compounds.

Fusing similarity rankings in ligand-based virtual screening

Data fusion is the name given to a range of methods for combining multiple sources of evidence. This mini-review summarizes the use of one such class of methods for combining the rankings obtained when similarity searching is used for ligand-based virtual screening. Two main approaches are described: similarity fusion involves combining rankings from single searches based on multiple similarity measures; and group fusion involves combining rankings from multiple searches based on a single similarity measure. The review then focuses on the rules that are available for combining similarity rankings, and on the evidence that exists for the superiority of fusion-based methods over conventional similarity searching.

The holistic integration of virtual screening in drug discovery

During the past decade, virtual screening (VS) has come of age. In this review, we document the evolution and maturation of VS from a rather exotic, stand-alone method toward a versatile hit and lead identification technology. VS campaigns have become fully integrated into drug discovery campaigns, evenly matched and complementary to high-throughput screening (HTS) methods. Here, we propose a novel classification of VS applications to help to monitor the advances in VS and to support future improvement of computational hit and lead identification methods. Several relevant VS studies from recent publications, in both academic and industrial settings, were selected to demonstrate the progress in this area. Furthermore, we identify challenges that lie ahead for the development of integrated VS campaigns.

Comparison of chemical similarity measures using different numbers of query structures

Many different similarity measures have been described for searching chemical databases. Drawing on previous work in textual information retrieval, this paper investigates the numbers of queries that are required to make robust statements as to the relative retrieval effectiveness of different similarity measures. Experiments with the MDL Drug Data Report database suggest that much larger numbers of queries are ideally required for this purpose than has been the case in previous comparative studies in chemoinformatics.

Virtual screening in drug design

Virtual screening has become a standard tool in drug discovery to identify novel lead compounds that target a biomolecule of interest. I present several concepts in ligand-based and structure-based virtual screening and discuss some of the current shortcomings and new developments. I also highlight approaches that combine concepts from structure- and ligand-based design.

Molecular simulations of drug–receptor complexes in anticancer research

Molecular modeling and computer simulation techniques have matured significantly in recent years and proved their value in the study of drug–DNA, drug–DNA–protein, drug–protein and protein–protein interactions. Evolution in this area has gone hand-in-hand with an increased availability of structural data on biological macromolecules, major advances in molecular mechanics force fields and considerable improvements in computer technologies, most significantly processing speeds, multiprocessor programming and data-storage capacity. The information derived from molecular simulations of drug–receptor complexes can be used to extract structural and energetic information that is usually beyond current experimental possibilities, provide independent accounts of experimentally observed behavior, help in the interpretation of biochemical or pharmacological results, and open new avenues for research by posing novel relevant questions that can guide the design of new experiments. As drug-screening tools, ligand- and fragment-docking platforms stand out as powerful techniques that can provide candidate molecules for hit and lead development. This review provides an overall perspective of the main methods and focuses on some selected applications to both classical and novel anticancer targets.

Shape-based reprofiling of FDA-approved drugs for the H₁ histamine receptor

Reprofiling of existing drugs to treat conditions not originally targeted is an attractive means of addressing the problem of a decreasing stream of approved drugs. To determine if 3D shape similarity can be used to rationalize an otherwise serendipitous process, we employed 3D shape-based virtual screening to reprofile existing FDA-approved drugs. The study was conducted in two phases. First, multiple histamine H(1) receptor antagonists were identified to be used as query molecules, and these were compared to a database of approved drugs. Second, the hits were ranked according to 3D similarity and the top drugs evaluated in a cell-based assay. The virtual screening methodology proved highly successful, as 13 of 23 top drugs tested selectively inhibited histamine-induced calcium release with the best being chlorprothixene (IC(50) 1 nM). Finally, we confirmed that the drugs identified using the cell-based assay were all acting at the receptor level by conducting a radioligand-binding assay using rat membrane.

Analysis of structure-based virtual screening studies and characterization of identified active compounds

Structure-based virtual screening makes explicit or implicit use of 3D target structure information to detect novel active compounds. Results of nearly 300 currently available original applications have been analyzed to characterize the state-of-the-art in this field. Compound selection from docking calculations is much influenced by subjective criteria. Although submicromolar compounds are identified, the majority of docking hits are only weakly potent. However, only a small percentage of docking hits can be reproduced by ligand-based methods. When docking calculations identify potent hits, they often originate from specialized compound sources (e.g., pharmaceutical compound decks or target-focused libraries) and also display a notable bias towards kinase targets. Structure-based virtual screening is the dominant approach to computational hit identification. Docking calculations frequently identify active compounds. Limited accuracy of compound scoring and ranking currently presents a major caveat of the approach that is often compensated for by chemical intuition and knowledge.

Virtual compound screening in drug discovery

Virtual screening (VS) methods are applied in both academia and drug discovery, and can be divided into ligand- and target structure-based approaches. The VS field is still evolving and is characterized by scientific heterogeneity. The value of virtual compound screening for drug discovery is often debated, in particular, given the large investments made in experimental high-throughput screening technologies. The current state-of-the-art in the VS field is discussed. Despite its limitations, VS applications have often succeeded in identifying novel hits including first-in-class active compounds and novel chemotypes. VS has its place in pharmaceutical research, but there is still much room for further improvements including method evaluation and drug discovery applications. The potential of VS is currently underutilized because its complementarity to high-throughput screening is not sufficiently exploited. Building close interfaces between computational and experimental screening would further streamline the hit identification process.

Haptic-driven applications to molecular modeling: state-of-the-art and perspectives

Drug design is a creative process that combines different scientific expertise. With the development of increasingly powerful computers, disciplines such as molecular modeling and, in particular, drug design, are becoming an important component of drug discovery. However, modern software often limits the user interaction with the computer calculation, reducing the potential for researchers to use their knowledge in the design process. For this reason, interactive methodologies have been investigated in recent years. In particular, haptic-driven simulators offer the possibility for users to drive and control the modeling simulations, efficiently combining human knowledge and computational power. In this article, we will discuss the state-of-the-art and future perspectives of such methodologies.

Established and emerging trends in computational drug discovery in the structural genomics era

Bioinformatics and chemoinformatics approaches contribute to hit discovery, hit-to-lead optimization, safety profiling, and target identification and enhance our overall understanding of the health and disease states. A vast repertoire of computational methods has been reported and increasingly combined in order to address more and more challenging targets or complex molecular mechanisms in the context of large-scale integration of structure and bioactivity data produced by private and public drug research. This review explores some key computational methods directly linked to drug discovery and chemical biology with a special emphasis on compound collection preparation, virtual screening, protein docking, and systems pharmacology. A list of generally freely available software packages and online resources is provided, and examples of successful applications are briefly commented upon.

Survey of public domain software for docking simulations and virtual screening

Progress in functional genomics and structural studies on biological macromolecules are generating a growing number of potential targets for therapeutics, adding to the importance of computational approaches for small molecule docking and virtual screening of candidate compounds. In this review, recent improvements in several public domain packages that are widely used in the context of drug development, including DOCK, AutoDock, AutoDock Vina and Screening for Ligands by Induced-fit Docking Efficiently (SLIDE) are surveyed. The authors also survey methods for the analysis and visualisation of docking simulations, as an important step in the overall assessment of the results. In order to illustrate the performance and limitations of current docking programs, the authors used the National Center for Toxicological Research (NCTR) oestrogen receptor benchmark set of 232 oestrogenic compounds with experimentally measured strength of binding to oestrogen receptor alpha. The methods tested here yielded a correlation coefficient of up to 0.6 between the predicted and observed binding affinities for active compounds in this benchmark.

On setting up and assessing docking simulations for virtual screening

Small molecule docking and virtual screening of candidate compounds have become an integral part of drug discovery pipelines, complementing and streamlining experimental efforts in that regard. In this chapter, we describe specific software packages and protocols that can be used to efficiently set up a computational screening using a library of compounds and a docking program. We also discuss consensus- and clustering-based approaches that can be used to assess the results, and potentially re-rank the hits. While docking programs share many common features, they may require tailored implementation of virtual screening pipelines for specific computing platforms. Here, we primarily focus on solutions for several public domain packages that are widely used in the context of drug development.