From Complex Natural Products to Simple Synthetic Mimetics by Computational De Novo Design

We present the computational de novo design of synthetically accessible chemical entities that mimic the complex sesquiterpene natural product (?)‐Englerin A. We synthesized lead‐like probes from commercially available building blocks and profiled them for activity against a computationally predicted panel of macromolecular targets. Both the design template (?)‐Englerin A and its low‐molecular weight mimetics presented nanomolar binding affinities and antagonized the transient receptor potential calcium channel TRPM8 in a cell‐based assay, without showing target promiscuity or frequent‐hitter properties. This proof‐of‐concept study outlines an expeditious solution to obtaining natural‐product‐inspired chemical matter with desirable properties.     

ALARM NMR: a rapid and robust experimental method to detect reactive false positives in biochemical screens

High-throughput screening (HTS) of large compound collections typically results in numerous small molecule hits that must be carefully evaluated to identify valid drug leads. Although several filtering mechanisms and other tools exist that can assist the chemist in this process, it is often the case that costly synthetic resources are expended in pursuing false positives. We report here a rapid and reliable NMR-based method for identifying reactive false positives including those that oxidize or alkylate a protein target. Importantly, the reactive species need not be the parent compound, as both reactive impurities and breakdown products can be detected. The assay is called ALARM NMR (a La assay to detect reactive molecules by nuclear magnetic resonance) and is based on monitoring DTT-dependent (13)C chemical shift changes of the human La antigen in the presence of a test compound or mixture. Extensive validation has been performed to demonstrate the reliability and utility of using ALARM NMR to assess thiol reactivity. This included comparing ALARM NMR to a glutathione-based fluorescence assay, as well as testing a collection of more than 3500 compounds containing HTS hits from 23 drug targets. The data show that current in silico filtering tools fail to identify more than half of the compounds that can act via reactive mechanisms. Significantly, we show how ALARM NMR data has been critical in identifying reactive compounds that would otherwise have been prioritized for lead optimization. In addition, a new filtering tool has been developed on the basis of the ALARM NMR data that can augment current in silico programs for identifying nuisance compounds and improving the process of hit triage.     

Comparison of luminescence ADP production assay and radiometric scintillation proximity assay for Cdc7 kinase

Several assay technologies have been successfully adapted and used in HTS to screen for protein kinase inhibitors; however, emerging comparative analysis studies report very low hit overlap between the different technologies, which challenges the working assumption that hit identification is not dependent on the assay method of choice. To help address this issue, we performed two screens on the cancer target, Cdc7-Dbf4 heterodimeric protein kinase, using a direct assay detection method measuring [(33)P]-phosphate incorporation into the substrate and an indirect method measuring residual ADP production using luminescence. We conducted the two screens under similar conditions, where in one, we measured [(33)P]-phosphate incorporation using scintillation proximity assay (SPA), and in the other, we detected luminescence signal of the ATP-dependent luciferase after regenerating ATP from residual ADP (LUM). Surprisingly, little or no correlation were observed between the positives identified by the two methods; at a threshold of 30% inhibition, 25 positives were identified in the LUM screen whereas the SPA screen only identified two positives, Tannic acid and Gentian violet, with Tannic acid being common to both. We tested 20 out of the 25 positive compounds in secondary confirmatory study and confirmed 12 compounds including Tannic acid as Cdc7-Dbf4 kinase inhibitors. Gentian violet, which was only positive in the SPA screen, inhibited luminescence detection and categorized as a false positive. This report demonstrates the strong impact in detection format on the success of a screening campaign and the importance of carefully designed confirmatory assays to eliminate those compounds that target the detection part of the assay.     

Structure-based pharmacophore design and virtual screening for novel angiotensin converting enzyme 2 inhibitors

The metallopeptidase Angiotensin Converting Enzyme (ACE) is an important drug target for the treatment of hypertension, heart, kidney, and lung disease. Recently, a close and unique human ACE homologue termed ACE2 has been identified and found to be an interesting new cardiorenal disease target. With the recently resolved inhibitor-bound ACE2 crystal structure available, we have attempted a structure-based approach to identify novel potent and selective inhibitors. Computational approaches focus on pharmacophore-based virtual screening of large compound databases. Selectivity was ensured by initial screening for ACE inhibitors within an internal database and the Derwent World Drug Index, which could be reduced to zero false positives and 0.1% hit rate, respectively. An average hit reduction of 0.44% was achieved with a five feature hypothesis, searching approximately 3.8 million compounds from various commercial databases. Seventeen compounds were selected based on high fit values as well as diverse structure and subjected to experimental validation in a bioassay. We show that all compounds displayed an inhibitory effect on ACE2 activity, the six most promising candidates exhibiting IC50 values in the range of 62-179 microM.     

Construction and validation of a RET TK catalytic domain by homology modeling

RET tyrosine kinase (TK) oncoproteins are potential targets for anticancer therapy. However, the search for novel RET inhibitors has been hampered by the lack of a 3D structure of the receptor. In this study, the "open" and the "closed" structure of the RET TK catalytic domain have been built by homology modeling techniques. The structures were validated by extensive docking studies with practically all the inhibitors reported in the literature and through molecular dynamics simulations of the resulting complexes. All the expected major interactions between the active domain amino acids and the inhibitors have been reproduced in their details. Furthermore, the proposed 3D models are in agreement with the results of available mutation studies. Therefore, these models could be profitably used to filter off from large libraries new potential hit compounds able to target this enzyme.     

Selection of in silico drug screening results for G-protein-coupled receptors by using universal active probes

We developed a new protocol for in silico drug screening for G-protein-coupled receptors (GPCRs) using a set of "universal active probes" (UAPs) with an ensemble docking procedure. UAPs are drug-like compounds, which are actual active compounds of a variety of known proteins. The current targets were nine human GPCRs whose three-dimensional (3D) structures are unknown, plus three GPCRs, namely β(2)-adrenergic receptor (ADRB2), A(2A) adenosine receptor (A(2A)), and dopamine D3 receptor (D(3)), whose 3D structures are known. Homology-based models of the GPCRs were constructed based on the crystal structures with careful sequence inspection. After subsequent molecular dynamics (MD) simulation taking into account the explicit lipid membrane molecules with periodic boundary conditions, we obtained multiple model structures of the GPCRs. For each target structure, docking-screening calculations were carried out via the ensemble docking procedure, using both true active compounds of the target proteins and the UAPs with the multiple target screening (MTS) method. Consequently, the multiple model structures showed various screening results with both poor and high hit ratios, the latter of which could be identified as promising for use in in silico screening to find candidate compounds to interact with the proteins. We found that the hit ratio of true active compounds showed a positive correlation to that of the UAPs. Thus, we could retrieve appropriate target structures from the GPCR models by applying the UAPs, even if no active compound is known for the GPCRs. Namely, the screening result that showed a high hit ratio for the UAPs could be used to identify actual hit compounds for the target GPCRs.     

Enhancement of chemical rules for predicting compound reactivity towards protein thiol groups

Non-specific chemical modification of protein thiol groups continues to be a significant source of false positive hits from high-throughput screening campaigns and can even plague certain protein targets and chemical series well into lead optimization. While experimental tools exist to assess the risk and promiscuity associated with the chemical reactivity of existing compounds, computational tools are desired that can reliably identify substructures that are associated with chemical reactivity to aid in triage of HTS hit lists, external compound purchases, and library design. Here we describe a Bayesian classification model derived from more than 8,800 compounds that have been experimentally assessed for their potential to covalently modify protein targets. The resulting model can be implemented in the large-scale assessment of compound libraries for purchase or design. In addition, the individual substructures identified as highly reactive in the model can be used as look-up tables to guide chemists during hit-to-lead and lead optimization campaigns.     

Improving (Q)SAR predictions by examining bias in the selection of compounds for experimental testing

ABSTRACT

Existing data on structures and biological activities are limited and distributed unevenly across distinct molecular targets and chemical compounds. The question arises if these data represent an unbiased sample of the general population of chemical-biological interactions. To answer this question, we analyzed ChEMBL data for 87,583 molecules tested against 919 protein targets using supervised and unsupervised approaches. Hierarchical clustering of the Murcko frameworks generated using Chemistry Development Toolkit showed that the available data form a big diffuse cloud without apparent structure. In contrast hereto, PASS-based classifiers allowed prediction whether the compound had been tested against the particular molecular target, despite whether it was active or not. Thus, one may conclude that the selection of chemical compounds for testing against specific targets is biased, probably due to the influence of prior knowledge. We assessed the possibility to improve (Q)SAR predictions using this fact: PASS prediction of the interaction with the particular target for compounds predicted as tested against the target has significantly higher accuracy than for those predicted as untested (average ROC AUC are about 0.87 and 0.75, respectively). Thus, considering the existing bias in the data of the training set may increase the performance of virtual screening.     

Evaluation of potential inhibitors of squalene synthase based on virtual screening and in vitro studies

Squalene synthase (SQS) is a potential target for hyperlipidemia treatment. To identify novel chemical scaffolds of SQS inhibitors, we generated 3D-QSAR pharmacophore models using HypoGen. The best quantitative pharmacophore model, Hypo 1, was selected for virtual screening using two chemical databases, Specs and Traditional Chinese Medicine database (TCM). The best-mapped hit compounds were then subjected to filtering by Lipinski^’s rule of five and docking studies to refine the hits. Finally, five compounds were selected from the top-ranked hit compounds for SQS inhibitory assay in vitro. Three of these compounds could inhibit SQS in vitro, and should be further evaluated pre-clinically as a treatment for hyperlipidemia.     

Probing the bioactivity-relevant chemical space of robust reactions and common molecular building blocks

In the search for new bioactive compounds, there is a trend toward increasingly complex compound libraries aiming to target the demanding targets of the future. In contrast, medicinal chemistry and traditional library design rely mainly on a small set of highly established and robust reactions. Here, we probe a set of 58 such reactions for their ability to sample the chemical space of known bioactive molecules, and the potential to create new scaffolds. Combined with ~26,000 common available building blocks, the reactions retrieve around 9% of a scaffold-diverse set of compounds active on human target proteins covering all major pharmaceutical target classes. Almost 80% of generated scaffolds from virtual one-step synthesis products are not present in a large set of known bioactive molecules for human targets, indicating potential for new discoveries. The results suggest that established synthesis resources are well suited to cover the known bioactivity-relevant chemical space and that there are plenty of unexplored regions accessible by these reactions, possibly providing valuable "low-hanging fruit" for hit discovery.     

Tubulin is a molecular target of the Wnt-activating chemical probe

Background

In drug discovery research, cell-based phenotypic screening is an essential method for obtaining potential drug candidates. Revealing the mechanism of action is a key step on the path to drug discovery. However, elucidating the target molecules of hit compounds from phenotypic screening campaigns remains a difficult and troublesome process. Simple and efficient methods for identifying the target molecules are essential.

Results

2-Amino-4-(3,4-(methylenedioxy)benzylamino)-6-(3-methoxyphenyl)pyrimidine (AMBMP) was identified as a senescence inducer from a phenotypic screening campaign. The compound is widely used as a Wnt agonist, although its target molecules remain to be clarified. To identify its target proteins, we compared a series of cellular assay results for the compound with our pathway profiling database. The database comprises the activities of compounds from simple assays of cellular reporter genes and cellular proliferations. In this database, compounds were classified on the basis of statistical analysis of their activities, which corresponded to a mechanism of action by the representative compounds. In addition, the mechanisms of action of the compounds of interest could be predicted using the database. Based on our database analysis, the compound was anticipated to be a tubulin disruptor, which was subsequently confirmed by its inhibitory activity of tubulin polymerization.

Conclusion

These results demonstrate that tubulin is identified for the first time as a target molecule of the Wnt-activating small molecule and that this might have misled the conclusions of some previous studies. Moreover, the present study also emphasizes that our pathway profiling database is a simple and potent tool for revealing the mechanisms of action of hit compounds obtained from phenotypic screenings and off targets of chemical probes.

     

Structural models in the assessment of protein druggability based on HTS data

Insights on the potential of target proteins to bind small molecules with high affinity can be derived from the knowledge of their three-dimensional structural details especially of their binding pockets. The present study uses high-throughput screening (HTS) results on various targets, to obtain mathematical predictive models in which a minimal set of structural parameters significantly contributing to the hit rates or the affinity of the protein binding pockets for small molecular entities, is identified. An emphasis is given to focus on target variation aspect of the data by consideration of commonly tested compounds against the HTS targets. We identify ‘four-parameter’ models with R ², , SEE, and LOO q ² values of 0.70, 0.60, 0.27 and 0.50, respectively, or better. We demonstrate through cross-validation exercises that our regression models apply well on varied data sets. Thus we can use these models to estimate hit rates for HTS campaigns and thereby assign priority to drug targets before they undergo such resource intense experimental screening and follow-up.

Identification of potential inhibitors against nuclear Dam1 complex subunit Ask1 of Candida albicans using virtual screening and MD simulations

Identification of hit compounds against specific target form the starting point for a drug discovery program. A consistent decline of new chemical entities (NCEs) in recent years prompted a challenge to explore newer approaches to discover potential hit compounds that in turn can be converted into leads, and ultimately drug with desired therapeutic efficacy. The vast amount of omics and activity data available in public databases offers an opportunity to identify novel targets and their potential inhibitors. State of the art in silico methods viz., clustering of compounds, virtual screening, molecular docking, MD simulations and MMPBSA calculations were employed in a pipeline to identify potential ‘hits’ against those targets as well whose structures, as of now, could only predict through threading approaches. In the present work, we have started from scratch, amino acid sequence of target and compounds retrieved from PubChem compound database, modeled it in such a way that led to the identification of possible inhibitors of Dam1 complex subunit Ask1 of Candida albicans. We also propose a ligand based binding site determination approach. We have identified potential inhibitors of Ask1 subunit of a Dam1 complex of C. albicans, which is required to prevent precocious spindle elongation in pre-mitotic phases. The proposed scheme may aid to find virtually potential inhibitors of other unique targets against candida.     

Protein flexibility and species specificity in structure-based drug discovery: dihydrofolate reductase as a test system

In structure-based drug discovery, researchers would like to identify all possible scaffolds for a given target. However, techniques that push the boundaries of chemical space could lead to many false positives or inhibitors that lack specificity for the target. Is it possible to broadly identify the appropriate chemical space for the inhibitors and yet maintain target specificity? To address this question, we have turned to dihydrofolate reductase (DHFR), a well-studied metabolic enzyme of pharmacological relevance. We have extended our multiple protein structure (MPS) method for receptor-based pharmacophore models to use multiple X-ray crystallographic structures. Models were created for DHFR from human and Pneumocystis carinii. These models incorporate a fair degree of protein flexibility and are highly selective for known DHFR inhibitors over drug-like non-inhibitors. Despite sharing a highly conserved active site, the pharmacophore models reflect subtle differences between the human and P. carinii forms, which identify species-specific, high-affinity inhibitors. We also use structures of DHFR from Candida albicans as a counter example. The available crystal structures show little flexibility, and the resulting models give poorer performance in identifying species-specific inhibitors. Therapeutic success for this system may depend on achieving species specificity between the related human host and these key fungal targets. The MPS technique is a promising advance for structure-based drug discovery for DHFR and other proteins of biomedical interest.     

Pharmacoinformatics exploration of polyphenol oxidases leading to novel inhibitors by virtual screening and molecular dynamic simulation study

Polyphenol oxidases (PPOs)/tyrosinases are metal-dependent enzymes and known as important targets for melanogenesis. Although considerable attempts have been conducted to control the melanin-associated diseases by using various inhibitors. However, the exploration of the best anti-melanin inhibitor without side effect still remains a challenge in drug discovery. In present study, protein structure prediction, ligand-based pharmacophore modeling, virtual screening, molecular docking and dynamic simulation study were used to screen the strong novel inhibitor to cure melanogenesis. The 3D structures of PPO1 and PPO2 were built through homology modeling, while the 3D crystal structures of PPO3 and PPO4 were retrieved from PDB. Pharmacophore modeling was performed using LigandScout 3.1 software and top five models were selected to screen the libraries (2601 of Aurora and 727, 842 of ZINC). Top 10 hit compounds (C1-10) were short-listed having strong binding affinities for PPO1-4. Drug and synthetic accessibility (SA) scores along with absorption, distribution, metabolism, excretion and toxicity (ADMET) assessment were employed to scrutinize the best lead hit. C4 (name) hit showed the best predicted SA score (5.75), ADMET properties and drug-likeness behavior among the short-listed compounds. Furthermore, docking simulations were performed to check the binding affinity of C1-C10 compounds against target proteins (PPOs). The binding affinity values of complex between C4 and PPOs were higher than those of other complexes (−11.70, −12.1, −9.90 and −11.20 kcal/mol with PPO1, PPO2, PPO3, or PPO4, respectively). From comparative docking energy and binding analyses, PPO2 may be considered as better target for melanogenesis than others. The potential binding modes of C4, C8 and C10 against PPO2 were explored using molecular dynamics simulations. The root mean square deviation and fluctuation (RMSD/RMSF) graphs results depict the significance of C4 over the other compounds. Overall, bioactivity and ligand efficiency profiles suggested that the proposed hit may be more effective inhibitors for melanogenesis.     

Prediction of biological targets for compounds using multiple-category Bayesian models trained on chemogenomics databases

Target identification is a critical step following the discovery of small molecules that elicit a biological phenotype. The present work seeks to provide an in silico correlate of experimental target fishing technologies in order to rapidly fish out potential targets for compounds on the basis of chemical structure alone. A multiple-category Laplacian-modified na?ve Bayesian model was trained on extended-connectivity fingerprints of compounds from 964 target classes in the WOMBAT (World Of Molecular BioAcTivity) chemogenomics database. The model was employed to predict the top three most likely protein targets for all MDDR (MDL Drug Database Report) database compounds. On average, the correct target was found 77% of the time for compounds from 10 MDDR activity classes with known targets. For MDDR compounds annotated with only therapeutic or generic activities such as "antineoplastic", "kinase inhibitor", or "anti-inflammatory", the model was able to systematically deconvolute the generic activities to specific targets associated with the therapeutic effect. Examples of successful deconvolution are given, demonstrating the usefulness of the tool for improving knowledge in chemogenomics databases and for predicting new targets for orphan compounds.     

Identification of novel malarial cysteine protease inhibitors using structure-based virtual screening of a focused cysteine protease inhibitor library

Malaria, in particular that caused by Plasmodium falciparum , is prevalent across the tropics, and its medicinal control is limited by widespread drug resistance. Cysteine proteases of P. falciparum , falcipain-2 (FP-2) and falcipain-3 (FP-3), are major hemoglobinases, validated as potential antimalarial drug targets. Structure-based virtual screening of a focused cysteine protease inhibitor library built with soft rather than hard electrophiles was performed against an X-ray crystal structure of FP-2 using the Glide docking program. An enrichment study was performed to select a suitable scoring function and to retrieve potential candidates against FP-2 from a large chemical database. Biological evaluation of 50 selected compounds identified 21 diverse nonpeptidic inhibitors of FP-2 with a hit rate of 42%. Atomic Fukui indices were used to predict the most electrophilic center and its electrophilicity in the identified hits. Comparison of predicted electrophilicity of electrophiles in identified hits with those in known irreversible inhibitors suggested the soft-nature of electrophiles in the selected target compounds. The present study highlights the importance of focused libraries and enrichment studies in structure-based virtual screening. In addition, few compounds were screened against homologous human cysteine proteases for selectivity analysis. Further evaluation of structure-activity relationships around these nonpeptidic scaffolds could help in the development of selective leads for antimalarial chemotherapy.     

Free resources to assist structure-based virtual ligand screening experiments

In today's research environment, a wealth of experimental/theoretical structural data is available and the number of therapeutically relevant macromolecular structures is growing rapidly. This, coupled with the huge number of small non-peptide potential drug candidates easily available (over 7 million compounds), highlight the need of using computer-aided techniques for the efficient identification and optimization of novel hit compounds. Virtual (or in silico) ligand screening based on the three-dimensional structure of macromolecular targets (SB-VLS) is firmly established as an important approach to identify chemical entities that have a high likelihood of binding to a target molecule to elicit desired biological responses. A myriad of free applications and services facilitating the drug discovery process have been posted on the Web. In this review, we cite over 350 URLs that are useful for SB-VLS projects and essentially free for academic groups. We attempt to provide links for in silico ADME/tox prediction tools, compound collections, some ligand-based methods, characterization/simulation of 3D targets and homology modeling tools, druggable pocket predictions, active site comparisons, analysis of macromolecular interfaces, protein docking tools to help identify binding pockets and protein-ligand docking/scoring methods. As such, we aim at providing both, methods pertaining to the field of Structural Bioinformatics (defined here as tools to study macromolecules) and methods pertaining to the field of Chemoinformatics (defined here as tools to make better decisions faster in the arena of drug/lead identification and optimization). We also report several recent success stories using these free computer methods. This review should help readers finding free computer tools useful for their projects. Overall, we are confident that these tools will facilitate rapid and cost-effective identification of new hit compounds. The URLs presented in this review will be updated regularly at www.vls3d.com in the coming months, "Links" section.