Chemical structural novelty: on-targets and off-targets

Tonantzitlolone B Modulates the Endogenous Opioid System to Promote Antinociception in Mice

Tonantzitlolone B (TZL-B) is a diterpene isolated from the roots of Stillingia loranthacea. Its antinociceptive effects were investigated in male Swiss mice using the following models of pain: formalin test, inflammation induced by Complete Freund's Adjuvant (CFA), tail flick test, and cold plate test. The influence of TZL-B on the opioid system was assessed in vivo, using opioid antagonists; in silico, investigating the chemical similarity among TZL-B and opioid agonists; and ex vivo, measuring preproenkephalin (PENK) gene expression in the spinal cord by RT-qPCR. TZL-B (10-1000 μg/kg) promoted antinociception in the four experimental models without impairing mice's motor function. TZL-B did not alter paw edema during CFA-induced inflammation. The antinociceptive effects of TZL-B in the tail flick and cold plate tests were diminished by the opioid antagonists naloxone (5 mg/kg), NOR-BNI (0.5 mg/kg), naltrindole (3 mg/kg), and CTOP (1 mg/kg), indicating the involvement of κ-, δ-, and μ-opioid receptors. TZL-B showed no significant chemical similarity to opioid agonists, but the treatment with TZL-B (1000 μg/kg) increased PENK gene expression in the spinal cord of mice. These data suggest that TZL-B promotes antinociception by enhancing the transcription of PENK, hence modulating the endogenous opioid system.

Comprehensive assessment of nine target prediction web services: which should we choose for target fishing?

Identification of potential targets for known bioactive compounds and novel synthetic analogs is of considerable significance. In silico target fishing (TF) has become an alternative strategy because of the expensive and laborious wet-lab experiments, explosive growth of bioactivity data and rapid development of high-throughput technologies. However, these TF methods are based on different algorithms, molecular representations and training datasets, which may lead to different results when predicting the same query molecules. This can be confusing for practitioners in practical applications. Therefore, this study systematically evaluated nine popular ligand-based TF methods based on target and ligand-target pair statistical strategies, which will help practitioners make choices among multiple TF methods. The evaluation results showed that SwissTargetPrediction was the best method to produce the most reliable predictions while enriching more targets. High-recall similarity ensemble approach (SEA) was able to find real targets for more compounds compared with other TF methods. Therefore, SwissTargetPrediction and SEA can be considered as primary selection methods in future studies. In addition, the results showed that k = 5 was the optimal number of experimental candidate targets. Finally, a novel ensemble TF method based on consensus voting is proposed to improve the prediction performance. The precision of the ensemble TF method outperforms the individual TF method, indicating that the ensemble TF method can more effectively identify real targets within a given top-k threshold. The results of this study can be used as a reference to guide practitioners in selecting the most effective methods in computational drug discovery.

Random-forest model for drug-target interaction prediction via Kullbeck-Leibler divergence

Virtual screening has significantly improved the success rate of early stage drug discovery. Recent virtual screening methods have improved owing to advances in machine learning and chemical information. Among these advances, the creative extraction of drug features is important for predicting drug–target interaction (DTI), which is a large-scale virtual screening of known drugs. Herein, we report Kullbeck–Leibler divergence (KLD) as a DTI feature and the feature-driven classification model applicable to DTI prediction. For the purpose, E3FP three-dimensional (3D) molecular fingerprints of drugs as a molecular representation allow the computation of 3D similarities between ligands within each target (Q–Q matrix) to identify the uniqueness of pharmacological targets and those between a query and a ligand (Q–L vector) in DTIs. The 3D similarity matrices are transformed into probability density functions via kernel density estimation as a nonparametric estimation. Each density model can exploit the characteristics of each pharmacological target and measure the quasi-distance between the ligands. Furthermore, we developed a random forest model from the KLD feature vectors to successfully predict DTIs for representative 17 targets (mean accuracy: 0.882, out-of-bag score estimate: 0.876, ROC AUC: 0.990). The method is applicable for 2D chemical similarity.

Candidate Therapeutics by Screening for Multitargeting Ligands: Combining the CB2 Receptor With CB1, PPARγ and 5-HT4 Receptors

In recent years, the cannabinoid type 2 receptor (CB2R) has become a major target for treating many disease conditions. The old therapeutic paradigm of “one disease-one target-one drug” is being transformed to “complex disease-many targets-one drug.” Multitargeting, therefore, attracts much attention as a promising approach. We thus focus on designing single multitargeting agents (MTAs), which have many advantages over combined therapies. Using our ligand-based approach, the “Iterative Stochastic Elimination” (ISE) algorithm, we produce activity models of agonists and antagonists for desired therapeutic targets and anti-targets. These models are used for sequential virtual screening and scoring large libraries of molecules in order to pick top-scored candidates for testing in vitro and in vivo. In this study, we built activity models for CB2R and other targets for combinations that could be used for several indications. Those additional targets are the cannabinoid 1 receptor (CB1R), peroxisome proliferator-activated receptor gamma (PPARγ), and 5-Hydroxytryptamine receptor 4 (5-HT4R). All these models have high statistical parameters and are reliable. Many more CB2R/CBIR agonists were found than combined CB2R agonists with CB1R antagonist activity (by 200 fold). CB2R agonism combined with PPARγ or 5-HT4R agonist activity may be used for treating Inflammatory Bowel Disease (IBD). Combining CB2R agonism with 5-HT4R generates more candidates (14,008) than combining CB2R agonism with agonists for the nuclear receptor PPARγ (374 candidates) from an initial set of ∼2.1 million molecules. Improved enrichment of true vs. false positives may be achieved by requiring a better ISE score cutoff or by performing docking. Those candidates can be purchased and tested experimentally to validate their activity. Further, we performed docking to CB2R structures and found lower statistical performance of the docking (“structure-based”) compared to ISE modeling (“ligand-based”). Therefore, ISE modeling may be a better starting point for molecular discovery than docking.

A systems-based method to repurpose marketed therapeutics for antiviral use: a SARS-CoV-2 case study

This study describes two complementary methods that use network-based and sequence similarity tools to identify drug repurposing opportunities predicted to modulate viral proteins. This approach could be rapidly adapted to new and emerging viruses. The first method built and studied a virus-host-physical interaction network; a three-layer multimodal network of drug target proteins, human protein-protein interactions, and viral-host protein-protein interactions. The second method evaluated sequence similarity between viral proteins and other proteins, visualized by constructing a virus-host-similarity interaction network. Methods were validated on the human immunodeficiency virus, hepatitis B, hepatitis C, and human papillomavirus, then deployed on SARS-CoV-2. Comparison of virus-host-physical interaction predictions to known antiviral drugs had AUCs of 0.69, 0.59, 0.78, and 0.67, respectively, reflecting that the scores are predictive of effective drugs. For SARS-CoV-2, 569 candidate drugs were predicted, of which 37 had been included in clinical trials for SARS-CoV-2 (AUC = 0.75, P-value 3.21 × 10-3). As further validation, top-ranked candidate antiviral drugs were analyzed for binding to protein targets in silico; binding scores generated by BindScope indicated a 70% success rate.

5-O-methylcneorumchromone K Exerts Antinociceptive Effects in Mice via Interaction with GABAA Receptors

The proper pharmacological control of pain is a continuous challenge for patients and health care providers. Even the most widely used medications for pain treatment are still ineffective or unsafe for some patients, especially for those who suffer from chronic pain. Substances containing the chromone scaffold have shown a variety of biological activities, including analgesic effects. This work presents for the first time the centrally mediated antinociceptive activity of 5-O-methylcneorumchromone K (5-CK). Cold plate and tail flick tests in mice showed that the 5-CK-induced antinociception was dose-dependent, longer-lasting, and more efficacious than that induced by morphine. The 5-CK-induced antinociception was not reversed by the opioid antagonist naloxone. Topological descriptors (fingerprints) were employed to narrow the antagonist selection to further investigate 5-CK's mechanism of action. Next, based on the results of fingerprints analysis, functional antagonist assays were conducted on nociceptive tests. The effect of 5-CK was completely reversed in both cold plate and tail-flick tests by GABA_A receptor antagonist bicuculline, but not by atropine or glibenclamide. Molecular docking studies suggest that 5-CK binds to the orthosteric binding site, with a similar binding profile to that observed for bicuculline and GABA. These results evidence that 5-CK has a centrally mediated antinociceptive effect, probably involving the activation of GABAergic pathways.

Planomonospora: A Metabolomics Perspective on an Underexplored Actinobacteria Genus

Despite an excellent track record, microbial drug discovery suffers from high rates of rediscovery. Better workflows for the rapid investigation of complex extracts are needed to increase throughput and to allow early prioritization of samples. In addition, systematic characterization of poorly explored strains is seldomly performed. Here, we report a metabolomic study of 72 isolates belonging to the rare actinomycete genus Planomonospora, using a workflow of commonly used open access tools to investigate its secondary metabolites. The results reveal a correlation of chemical diversity and strain phylogeny, with classes of metabolites exclusive to certain phylogroups. We were able to identify previously reported Planomonospora metabolites, including the ureylene-containing oligopeptide antipain, the thiopeptide siomycin including new congeners, and the ribosomally synthesized peptides sphaericin and lantibiotic 97518. In addition, we found that Planomonospora strains can produce the siderophore desferrioxamine or a salinichelin-like peptide. Analysis of the genomes of three newly sequenced strains led to the detection of 59 gene cluster families, of which three were connected to products found by LC-MS/MS profiling. This study demonstrates the value of metabolomic studies to investigate poorly explored taxa and provides a first picture of the biosynthetic capabilities of the genus Planomonospora.

A biaryl-linked tripeptide from Planomonospora reveals a widespread class of minimal RiPP gene clusters

Microbial natural products impress by their bioactivity, structural diversity, and ingenious biosynthesis. While screening the less exploited actinobacterial genus Planomonospora, two cyclopeptides were discovered, featuring an unusual Tyr-His biaryl bridging across a tripeptide scaffold, with the sequences N-acetyl-Tyr-Tyr-His and N-acetyl-Tyr-Phe-His. Planomonospora genomes pointed toward a ribosomal synthesis of the cyclopeptide from a pentapeptide precursor encoded by 18-bp bytA, to our knowledge the smallest coding gene ever reported. Closely linked to bytA is bytO, encoding a cytochrome P450 monooxygenase likely responsible for biaryl installment. In Streptomyces, the bytAO segment was sufficient to direct production of the crosslinked N-acetylated Tyr-Tyr-His tripeptide. Bioinformatic analysis of related cytochrome P450 monooxygenases indicated that they constitute a widespread family of enzymes, and the corresponding genes are closely linked to 5-amino acid coding sequences in approximately 200 (actino)bacterial genomes, all with potential for biaryl linkage between amino acids 1 and 3. We propose the named biarylitides this family of RiPPs.

Inhibition of intracellular Ca2+ mobilization and potassium channels activation are involved in the vasorelaxation induced by 7-hydroxycoumarin

Coumarins exhibit a wide variety of biological effects, including activities in the cardiovascular system and the aim of this study was to evaluate the vascular therapeutic potential of 7-Hydroxicoumarin (7-HC). The vascular effects induced by 7-HC (0.001 μM-300 μM), were investigated by in vitro approaches using isometric tension measurements in rat superior mesenteric arteries and by in silico assays using Ligand-based analysis. Our results suggest that the vasorelaxant effect of 7-HC seems to rely on potassium channels, notably through large conductance Ca²⁺-activated K⁺ (BK_Ca) channels activation. In fact, 7-HC (300 μM) significantly reduced CaCl₂-induced contraction as well as the reduction of intracellular calcium mobilization. However, the relaxation induced by 7-HC was independent of store-operated calcium entry (SOCE). Moreover, in silico analysis suggests that potassium channels have a common binding pocket, where 7-HC may bind and hint that its binding profile is more similar to quinine's than verapamil's. These results are compatible with the inhibition of Ca²⁺ release from intracellular stores, which is prompted by phenylephrine and caffeine. Taken together, these results demonstrate a therapeutic potential of 7-HC on the cardiovascular system, making it a promising lead compound for the development of drugs useful in the treatment of cardiovascular diseases.

Comparing a Query Compound with Drug Target Classes Using 3D-Chemical Similarity

3D similarity is useful in predicting the profiles of unprecedented molecular frameworks that are 2D dissimilar to known compounds. When comparing pairs of compounds, 3D similarity of the pairs depends on conformational sampling, the alignment method, the chosen descriptors, and the similarity coefficients. In addition to these four factors, 3D chemocentric target prediction of an unknown compound requires compound-target associations, which replace compound-to-compound comparisons with compound-to-target comparisons. In this study, quantitative comparison of query compounds to target classes (one-to-group) was achieved via two types of 3D similarity distributions for the respective target class with parameter optimization for the fitting models: (1) maximum likelihood (ML) estimation of queries, and (2) the Gaussian mixture model (GMM) of target classes. While Jaccard-Tanimoto similarity of query-to-ligand pairs with 3D structures (sampled multi-conformers) can be transformed into query distribution using ML estimation, the ligand pair similarity within each target class can be transformed into a representative distribution of a target class through GMM, which is hyperparameterized via the expectation-maximization (EM) algorithm. To quantify the discriminativeness of a query ligand against target classes, the Kullback-Leibler (K-L) divergence of each query was calculated and compared between targets. 3D similarity-based K-L divergence together with the probability and the feasibility index, (F_m), showed discriminative power with regard to some query-class associations. The K-L divergence of 3D similarity distributions can be an additional method for (1) the rank of the 3D similarity score or (2) the p-value of one 3D similarity distribution to predict the target of unprecedented drug scaffolds.

The Road Not Taken with Pyrrole-Imidazole Polyamides: Off-Target Effects and Genomic Binding

The high sequence specificity of minor groove-binding N-methylpyrrole-N-methylimidazole polyamides have made significant advances in cancer and disease biology, yet there have been few comprehensive reports on their off-target effects, most likely as a consequence of the lack of available tools in evaluating genomic binding, an essential aspect that has gone seriously underexplored. Compared to other N-heterocycles, the off-target effects of these polyamides and their specificity for the DNA minor groove and primary base pair recognition require the development of new analytical methods, which are missing in the field today. This review aims to highlight the current progress in deciphering the off-target effects of these N-heterocyclic molecules and suggests new ways that next-generating sequencing can be used in addressing off-target effects.

Electrostatic-field and surface-shape similarity for virtual screening and pose prediction

We introduce a new method for rapid computation of 3D molecular similarity that combines electrostatic field comparison with comparison of molecular surface-shape and directional hydrogen-bonding preferences (called "eSim"). Rather than employing heuristic "colors" or user-defined molecular feature types to represent conformation-dependent molecular electrostatics, eSim calculates the similarity of the electrostatic fields of two molecules (in addition to shape and hydrogen-bonding). We present detailed virtual screening performance data on the standard 102 target DUD-E set. In its moderately fast screening mode, eSim running on a single computing core is capable of processing over 60 molecules per second. In this mode, eSim performed significantly better than all alternate methods for which full DUD-E data were available (mean ROC area of 0.74, p [Formula: see text], by paired t-test, compared with the best performing alternate method). In addition, for 92 targets of the DUD-E set where multiple ligand-bound crystal structures were available, screening performance was assessed using alternate ligands or sets thereof (in their bound poses) as similarity targets. Using the joint alignment of five ligands for each protein target, mean ROC area exceeded 0.82 for the 92 targets. Design-focused application of ligand similarity methods depends on accurate predictions of geometric molecular relationships. We comprehensively assessed pose prediction accuracy by curating nearly 400,000 bound ligand pose pairs across the DUD-E targets. Overall, beginning from agnostic initial poses, we observed an 80% success rate for RMSD [Formula: see text] Å  among the top 20 predicted eSim poses. These examples were split roughly 50/50 into cases with high direct atomic overlap (where a shared scaffold exists between a pair) and low direct atomic overlap (where where a ligand pair has dissimilar scaffolds but largely occupies the same space). Within the high direct atomic overlap subset, the pose prediction success rate was 93%. For the more challenging subset (where dissimilar scaffolds are to be aligned), the success rate was 70%. The eSim approach enables both large-scale screening and rational design of ligands and is rooted in physically meaningful, non-heuristic, molecular comparisons.

Discovery of CNS-Like D3R-Selective Antagonists Using 3D Pharmacophore Guided Virtual Screening

The dopamine D3 receptor is an important CNS target for the treatment of a variety of neurological diseases. Selective dopamine D3 receptor antagonists modulate the improvement of psychostimulant addiction and relapse. In this study, five and six featured pharmacophore models of D3R antagonists were generated and evaluated with the post-hoc score combining two survival scores of active and inactive. Among the Top 10 models, APRRR215 and AHPRRR104 were chosen based on the coefficient of determination (APRRR215: R²training = 0.80; AHPRRR104: R²training = 0.82) and predictability (APRRR215: Q²test = 0.73, R²predictive = 0.82; AHPRRR104: Q²test = 0.86, R²predictive = 0.74) of their 3D-quantitative structure⁻activity relationship models. Pharmacophore-based virtual screening of a large compound library from eMolecules (>3 million compounds) using two optimal models expedited the search process by a 100-fold speed increase compared to the docking-based screening (HTVS scoring function in Glide) and identified a series of hit compounds having promising novel scaffolds. After the screening, docking scores, as an adjuvant predictor, were added to two fitness scores (from the pharmacophore models) and predicted Ki (from PLSs of the QSAR models) to improve accuracy. Final selection of the most promising hit compounds were also evaluated for CNS-like properties as well as expected D3R antagonism.

Quantitative surface field analysis: learning causal models to predict ligand binding affinity and pose

We introduce the QuanSA method for inducing physically meaningful field-based models of ligand binding pockets based on structure-activity data alone. The method is closely related to the QMOD approach, substituting a learned scoring field for a pocket constructed of molecular fragments. The problem of mutual ligand alignment is addressed in a general way, and optimal model parameters and ligand poses are identified through multiple-instance machine learning. We provide algorithmic details along with performance results on sixteen structure-activity data sets covering many pharmaceutically relevant targets. In particular, we show how models initially induced from small data sets can extrapolatively identify potent new ligands with novel underlying scaffolds with very high specificity. Further, we show that combining predictions from QuanSA models with those from physics-based simulation approaches is synergistic. QuanSA predictions yield binding affinities, explicit estimates of ligand strain, associated ligand pose families, and estimates of structural novelty and confidence. The method is applicable for fine-grained lead optimization as well as potent new lead identification.

Most Ligand-Based Classification Benchmarks Reward Memorization Rather than Generalization

Undetected overfitting can occur when there are significant redundancies between training and validation data. We describe AVE, a new measure of training-validation redundancy for ligand-based classification problems, that accounts for the similarity among inactive molecules as well as active ones. We investigated seven widely used benchmarks for virtual screening and classification, and we show that the amount of AVE bias strongly correlates with the performance of ligand-based predictive methods irrespective of the predicted property, chemical fingerprint, similarity measure, or previously applied unbiasing techniques. Therefore, it may be the case that the previously reported performance of most ligand-based methods can be explained by overfitting to benchmarks rather than good prospective accuracy.

90 years of monoamine oxidase: some progress and some confusion

It would not be practical to attempt to deal with all the advances that have informed our understanding of the behavior and functions of this enzyme over the past 90 years. This account concentrates key advances that explain why the monoamine oxidases remain of pharmacological and biochemical interest and on some areas of continuing uncertainty. Some issues that remain to be understood or are in need of further clarification are highlighted.

SDTNBI: an integrated network and chemoinformatics tool for systematic prediction of drug-target interactions and drug repositioning

Computational prediction of drug-target interactions (DTIs) and drug repositioning provides a low-cost and high-efficiency approach for drug discovery and development. The traditional social network-derived methods based on the naïve DTI topology information cannot predict potential targets for new chemical entities or failed drugs in clinical trials. There are currently millions of commercially available molecules with biologically relevant representations in chemical databases. It is urgent to develop novel computational approaches to predict targets for new chemical entities and failed drugs on a large scale. In this study, we developed a useful tool, namely substructure-drug-target network-based inference (SDTNBI), to prioritize potential targets for old drugs, failed drugs and new chemical entities. SDTNBI incorporates network and chemoinformatics to bridge the gap between new chemical entities and known DTI network. High performance was yielded in 10-fold and leave-one-out cross validations using four benchmark data sets, covering G protein-coupled receptors, kinases, ion channels and nuclear receptors. Furthermore, the highest areas under the receiver operating characteristic curve were 0.797 and 0.863 for two external validation sets, respectively. Finally, we identified thousands of new potential DTIs via implementing SDTNBI on a global network. As a proof-of-principle, we showcased the use of SDTNBI to identify novel anticancer indications for nonsteroidal anti-inflammatory drugs by inhibiting AKR1C3, CA9 or CA12. In summary, SDTNBI is a powerful network-based approach that predicts potential targets for new chemical entities on a large scale and will provide a new tool for DTI prediction and drug repositioning. The program and predicted DTIs are available on request.

The comparison of automated clustering algorithms for resampling representative conformer ensembles with RMSD matrix

Background

The accuracy of any 3D-QSAR, Pharmacophore and 3D-similarity based chemometric target fishing models are highly dependent on a reasonable sample of active conformations. Since a number of diverse conformational sampling algorithm exist, which exhaustively generate enough conformers, however model building methods relies on explicit number of common conformers.

Results

In this work, we have attempted to make clustering algorithms, which could find reasonable number of representative conformer ensembles automatically with asymmetric dissimilarity matrix generated from openeye tool kit. RMSD was the important descriptor (variable) of each column of the N × N matrix considered as N variables describing the relationship (network) between the conformer (in a row) and the other N conformers. This approach used to evaluate the performance of the well-known clustering algorithms by comparison in terms of generating representative conformer ensembles and test them over different matrix transformation functions considering the stability. In the network, the representative conformer group could be resampled for four kinds of algorithms with implicit parameters. The directed dissimilarity matrix becomes the only input to the clustering algorithms.

Conclusions

Dunn index, Davies–Bouldin index, Eta-squared values and omega-squared values were used to evaluate the clustering algorithms with respect to the compactness and the explanatory power. The evaluation includes the reduction (abstraction) rate of the data, correlation between the sizes of the population and the samples, the computational complexity and the memory usage as well. Every algorithm could find representative conformers automatically without any user intervention, and they reduced the data to 14–19% of the original values within 1.13 s per sample at the most. The clustering methods are simple and practical as they are fast and do not ask for any explicit parameters. RCDTC presented the maximum Dunn and omega-squared values of the four algorithms in addition to consistent reduction rate between the population size and the sample size. The performance of the clustering algorithms was consistent over different transformation functions. Moreover, the clustering method can also be applied to molecular dynamics sampling simulation results.

Electronic supplementary material

The online version of this article (doi:10.1186/s13321-017-0208-0) contains supplementary material, which is available to authorized users.

The impact of in silico screening in the discovery of novel and safer drug candidates

Drug discovery is a multidisciplinary and multivariate optimization endeavor. As such, in silico screening tools have gained considerable importance to archive, analyze and exploit the vast and ever-increasing amount of experimental data generated throughout the process. The current review will focus on the computer-aided prediction of the numerous properties that need to be controlled during the discovery of a preliminary hit and its promotion to a viable clinical candidate. It does not pretend to the almost impossible task of an exhaustive report but will highlight a few key points that need to be collectively addressed both by chemists and biologists to fuel the drug discovery pipeline with innovative and safe drug candidates.

Metal- and additive-free oxygen-atom transfer reaction: an efficient and chemoselective oxidation of sulfides to sulfoxides with cyclic diacyl peroxides

Metal- and additive-free sulfoxidation was developed using cyclic diacyl peroxides as oxygen sources and a single two-electron transfer mechanism was suggested.

Extrapolative prediction using physically-based QSAR

Surflex-QMOD integrates chemical structure and activity data to produce physically-realistic models for binding affinity prediction . Here, we apply QMOD to a 3D-QSAR benchmark dataset and show broad applicability to a diverse set of targets. Testing new ligands within the QMOD model employs automated flexible molecular alignment, with the model itself defining the optimal pose for each ligand. QMOD performance was compared to that of four approaches that depended on manual alignments (CoMFA, two variations of CoMSIA, and CMF). QMOD showed comparable performance to the other methods on a challenging, but structurally limited, test set. The QMOD models were also applied to test a large and structurally diverse dataset of ligands from ChEMBL, nearly all of which were synthesized years after those used for model construction. Extrapolation across diverse chemical structures was possible because the method addresses the ligand pose problem and provides structural and geometric means to quantitatively identify ligands within a model’s applicability domain. Predictions for such ligands for the four tested targets were highly statistically significant based on rank correlation. Those molecules predicted to be highly active (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {pK}_i \ge 7.5$$\end{document}pKi≥7.5) had a mean experimental \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {pK}_i$$\end{document}pKi of 7.5, with potent and structurally novel ligands being identified by QMOD for each target.

In silico methods to address polypharmacology: current status, applications and future perspectives

Polypharmacology, a new paradigm in drug discovery that focuses on multi-target drugs (MTDs), has potential application for drug repurposing, the process of finding new uses for existing approved drugs, prediction of off-target toxicities and rational design of MTDs. In this scenario, computational strategies have demonstrated great potential in predicting polypharmacology and in facilitating drug repurposing. Here, we provide a comprehensive overview of various computational approaches that enable the prediction and analysis of in vitro and in vivo drug-response phenotypes and outline their potential for drug discovery. We give an outlook on the latest advances in rational design of MTDs and discuss possible future directions of algorithm development in this field.

Relating Essential Proteins to Drug Side-Effects Using Canonical Component Analysis: A Structure-Based Approach

The molecular mechanism of many drug side-effects is unknown and difficult to predict. Previous methods for explaining side-effects have focused on known drug targets and their pathways. However, low affinity binding to proteins that are not usually considered drug targets may also drive side-effects. In order to assess these alternative targets, we used the 3D structures of 563 essential human proteins systematically to predict binding to 216 drugs. We first benchmarked our affinity predictions with available experimental data. We then combined singular value decomposition and canonical component analysis (SVD-CCA) to predict side-effects based on these novel target profiles. Our method predicts side-effects with good accuracy (average AUC: 0.82 for side effects present in <50% of drug labels). We also noted that side-effect frequency is the most important feature for prediction and can confound efforts at elucidating mechanism; our method allows us to remove the contribution of frequency and isolate novel biological signals. In particular, our analysis produces 2768 triplet associations between 50 essential proteins, 99 drugs, and 77 side-effects. Although experimental validation is difficult because many of our essential proteins do not have validated assays, we nevertheless attempted to validate a subset of these associations using experimental assay data. Our focus on essential proteins allows us to find potential associations that would likely be missed if we used recognized drug targets. Our associations provide novel insights about the molecular mechanisms of drug side-effects and highlight the need for expanded experimental efforts to investigate drug binding to proteins more broadly.

PubChem structure-activity relationship (SAR) clusters

BACKGROUND

Developing structure-activity relationships (SARs) of molecules is an important approach in facilitating hit exploration in the early stage of drug discovery. Although information on millions of compounds and their bioactivities is freely available to the public, it is very challenging to infer a meaningful and novel SAR from that information.

RESULTS

Research discussed in the present paper employed a bioactivity-centered clustering approach to group 843,845 non-inactive compounds stored in PubChem according to both structural similarity and bioactivity similarity, with the aim of mining bioactivity data in PubChem for useful SAR information. The compounds were clustered in three bioactivity similarity contexts: (1) non-inactive in a given bioassay, (2) non-inactive against a given protein, and (3) non-inactive against proteins involved in a given pathway. In each context, these small molecules were clustered according to their two-dimensional (2-D) and three-dimensional (3-D) structural similarities. The resulting 18 million clusters, named "PubChem SAR clusters", were delivered in such a way that each cluster contains a group of small molecules similar to each other in both structure and bioactivity.

CONCLUSIONS

The PubChem SAR clusters, pre-computed using publicly available bioactivity information, make it possible to quickly navigate and narrow down the compounds of interest. Each SAR cluster can be a useful resource in developing a meaningful SAR or enable one to design or expand compound libraries from the cluster. It can also help to predict the potential therapeutic effects and pharmacological actions of less-known compounds from those of well-known compounds (i.e., drugs) in the same cluster.

Computational investigations of hERG channel blockers: New insights and current predictive models

Identification of potential human Ether-a-go-go Related-Gene (hERG) potassium channel blockers is an essential part of the drug development and drug safety process in pharmaceutical industries or academic drug discovery centers, as they may lead to drug-induced QT prolongation, arrhythmia and Torsade de Pointes. Recent reports also suggest starting to address such issues at the hit selection stage. In order to prioritize molecules during the early drug discovery phase and to reduce the risk of drug attrition due to cardiotoxicity during pre-clinical and clinical stages, computational approaches have been developed to predict the potential hERG blockage of new drug candidates. In this review, we will describe the current in silico methods developed and applied to predict and to understand the mechanism of actions of hERG blockers, including ligand-based and structure-based approaches. We then discuss ongoing research on other ion channels and hERG polymorphism susceptible to be involved in LQTS and how systemic approaches can help in the drug safety decision.

Knowledge-guided docking: accurate prospective prediction of bound configurations of novel ligands using Surflex-Dock

Prediction of the bound configuration of small-molecule ligands that differ substantially from the cognate ligand of a protein co-crystal structure is much more challenging than re-docking the cognate ligand. Success rates for cross-docking in the range of 20-30 % are common. We present an approach that uses structural information known prior to a particular cutoff-date to make predictions on ligands whose bounds structures were determined later. The knowledge-guided docking protocol was tested on a set of ten protein targets using a total of 949 ligands. The benchmark data set, called PINC ("PINC Is Not Cognate"), is publicly available. Protein pocket similarity was used to choose representative structures for ensemble-docking. The docking protocol made use of known ligand poses prior to the cutoff-date, both to help guide the configurational search and to adjust the rank of predicted poses. Overall, the top-scoring pose family was correct over 60 % of the time, with the top-two pose families approaching a 75 % success rate. Correct poses among all those predicted were identified nearly 90 % of the time. The largest improvements came from the use of molecular similarity to improve ligand pose rankings and the strategy for identifying representative protein structures. With the exception of a single outlier target, the knowledge-guided docking protocol produced results matching the quality of cognate-ligand re-docking, but it did so on a very challenging temporally-segregated cross-docking benchmark.

Prospective performance evaluation of selected common virtual screening tools. Case study: Cyclooxygenase (COX) 1 and 2

Computational methods can be applied in drug development for the identification of novel lead candidates, but also for the prediction of pharmacokinetic properties and potential adverse effects, thereby aiding to prioritize and identify the most promising compounds. In principle, several techniques are available for this purpose, however, which one is the most suitable for a specific research objective still requires further investigation. Within this study, the performance of several programs, representing common virtual screening methods, was compared in a prospective manner. First, we selected top-ranked virtual screening hits from the three methods pharmacophore modeling, shape-based modeling, and docking. For comparison, these hits were then additionally predicted by external pharmacophore- and 2D similarity-based bioactivity profiling tools. Subsequently, the biological activities of the selected hits were assessed in vitro, which allowed for evaluating and comparing the prospective performance of the applied tools. Although all methods performed well, considerable differences were observed concerning hit rates, true positive and true negative hits, and hitlist composition. Our results suggest that a rational selection of the applied method represents a powerful strategy to maximize the success of a research project, tightly linked to its aims. We employed cyclooxygenase as application example, however, the focus of this study lied on highlighting the differences in the virtual screening tool performances and not in the identification of novel COX-inhibitors.

Identification of structural features in chemicals associated with cancer drug response: a systematic data-driven analysis

Motivation: Analysis of relationships of drug structure to biological response is key to understanding off-target and unexpected drug effects, and for developing hypotheses on how to tailor drug therapies. New methods are required for integrated analyses of a large number of chemical features of drugs against the corresponding genome-wide responses of multiple cell models.

Results: In this article, we present the first comprehensive multi-set analysis on how the chemical structure of drugs impacts on genome-wide gene expression across several cancer cell lines [Connectivity Map (CMap) database]. The task is formulated as searching for drug response components across multiple cancers to reveal shared effects of drugs and the chemical features that may be responsible. The components can be computed with an extension of a recent approach called Group Factor Analysis. We identify 11 components that link the structural descriptors of drugs with specific gene expression responses observed in the three cell lines and identify structural groups that may be responsible for the responses. Our method quantitatively outperforms the limited earlier methods on CMap and identifies both the previously reported associations and several interesting novel findings, by taking into account multiple cell lines and advanced 3D structural descriptors. The novel observations include: previously unknown similarities in the effects induced by 15-delta prostaglandin J2 and HSP90 inhibitors, which are linked to the 3D descriptors of the drugs; and the induction by simvastatin of leukemia-specific response, resembling the effects of corticosteroids.

Availability and implementation: Source Code implementing the method is available at: http://research.ics.aalto.fi/mi/software/GFAsparse

Contact:suleiman.khan@aalto.fi or samuel.kaski@aalto.fi

Supplementary Information:Supplementary data are available at Bioinformatics online.

Identifying MurI uncompetitive inhibitors by correlating decomposed binding energies with bioactivity

MurI uncompetitive inhibitors can be virtually identified by a new method that correlates decomposed binding free energies with the bioactivity.

Chemical and protein structural basis for biological crosstalk between PPARα and COX enzymes

We have previously validated a probabilistic framework that combined computational approaches for predicting the biological activities of small molecule drugs. Molecule comparison methods included molecular structural similarity metrics and similarity computed from lexical analysis of text in drug package inserts. Here we present an analysis of novel drug/target predictions, focusing on those that were not obvious based on known pharmacological crosstalk. Considering those cases where the predicted target was an enzyme with known 3D structure allowed incorporation of information from molecular docking and protein binding pocket similarity in addition to ligand-based comparisons. Taken together, the combination of orthogonal information sources led to investigation of a surprising predicted relationship between a transcription factor and an enzyme, specifically, PPARα and the cyclooxygenase enzymes. These predictions were confirmed by direct biochemical experiments which validate the approach and show for the first time that PPARα agonists are cyclooxygenase inhibitors.

A new perspective on sexual mixing among men who have sex with men by body image

Background

“Casual sex” is seldom as non-selective and random as it may sound. During each sexual encounter, people consciously and unconsciously seek their casual sex partners according to different attributes. Influential to a sexual network, research focusing on quantifying the effects of physical appearance on sexual network has been sparse.

Methods

We evaluated the application of Log odds score (LOD) to assess the mixing patterns of 326 men who have sex with men (MSM) in Hong Kong in their networking of casual sex partners by Body Image Type (BIT). This involved an analysis of 1,196 respondents-casual sex partner pairs. Seven BITs were used in the study: Bear, Chubby, Slender, Lean toned, Muscular, Average and Other.

Results

A hierarchical pattern was observed in the preference of MSM for casual sex partners by the latter's BIT. Overall, Muscular men were most preferred, followed by Lean toned while the least preferred was Slender, as illustrated by LOD going down along the hierarchy in the same direction. Marked avoidance was found between men who self-identified as Chubby and men of Other body type (within-group-LOD: 1.25–2.89; between-group-LOD: <−1). None of the respondents reported to have networked a man who self-identified as Average for casual sex.

Conclusions

We have demonstrated the possibility of adopting a mathematical prototype to investigate the influence of BIT in a sexual network of MSM. Construction of matrix based on culture-specific BIT and cross-cultural comparisons would generate new knowledge on the mixing behaviors of MSM.

Comparison of ultra-fast 2D and 3D ligand and target descriptors for side effect prediction and network analysis in polypharmacology

BACKGROUND AND PURPOSE

Some existing computational methods are used to infer protein targets of small molecules and can therefore be used to find new targets for existing drugs, with the goals of re-directing the molecule towards a different therapeutic purpose or explaining off-target effects due to multiple targeting. Inherent limitations, however, arise from the fact that chemical analogy is calculated on the basis of common frameworks or scaffolds and also because target information is neglected. The method we present addresses these issues by taking into account 3D information from both the ligand and the target.

EXPERIMENTAL APPROACH

ElectroShape is an established method for ultra-fast comparison of the shapes and charge distributions of ligands that is validated here for prediction of on-target activities, off-target profiles and adverse effects of drugs and drug-like molecules taken from the DrugBank database.

KEY RESULTS

The method is shown to predict polypharmacology profiles and relate targets from two complementary viewpoints (ligand- and target-based networks).

CONCLUSIONS AND IMPLICATIONS

The open-access web tool presented here (http://ub.cbm.uam.es/chemogenomics/) allows interactive navigation in a unified 'pharmacological space' from the viewpoints of both ligands and targets. It also enables prediction of pharmacological profiles, including likely side effects, for new compounds. We hope this web interface will help many pharmacologists to become aware of this new paradigm (up to now mostly used in the realm of the so-called 'chemical biology') and encourage its use with a view to revealing 'hidden' relationships between new and existing compounds and pharmacologically relevant targets.

Prediction of off-target drug effects through data fusion

We present a probabilistic data fusion framework that combines multiple computational approaches for drawing relationships between drugs and targets. The approach has special relevance to identifying surprising unintended biological targets of drugs. Comparisons between molecules are made based on 2D topological structural considerations, based on 3D surface characteristics, and based on English descriptions of clinical effects. Similarity computations within each modality were transformed into probability scores. Given a new molecule along with a set of molecules sharing some biological effect, a single score based on comparison to the known set is produced, reflecting either 2D similarity, 3D similarity, clinical effects similarity or their combination. The methods were validated within acurated structural pharmacology database (SPDB) and further tested by blind application to data derived from the ChEMBL database. For prediction of off-target effects, 3D-similarity performed best as a single modality, but combining all methods produced performance gains. Striking examples of structurally surprising off-target predictions are presented.

Predicting drug-target interactions using probabilistic matrix factorization

Quantitative analysis of known drug-target interactions emerged in recent years as a useful approach for drug repurposing and assessing side effects. In the present study, we present a method that uses probabilistic matrix factorization (PMF) for this purpose, which is particularly useful for analyzing large interaction networks. DrugBank drugs clustered based on PMF latent variables show phenotypic similarity even in the absence of 3D shape similarity. Benchmarking computations show that the method outperforms those recently introduced provided that the input data set of known interactions is sufficiently large--which is the case for enzymes and ion channels, but not for G-protein coupled receptors (GPCRs) and nuclear receptors. Runs performed on DrugBank after hiding 70% of known interactions show that, on average, 88 of the top 100 predictions hit the hidden interactions. De novo predictions permit us to identify new potential interactions. Drug-target pairs implicated in neurobiological disorders are overrepresented among de novo predictions.

Substrate-driven mapping of the degradome by comparison of sequence logos

Sequence logos are frequently used to illustrate substrate preferences and specificity of proteases. Here, we employed the compiled substrates of the MEROPS database to introduce a novel metric for comparison of protease substrate preferences. The constructed similarity matrix of 62 proteases can be used to intuitively visualize similarities in protease substrate readout via principal component analysis and construction of protease specificity trees. Since our new metric is solely based on substrate data, we can engraft the protease tree including proteolytic enzymes of different evolutionary origin. Thereby, our analyses confirm pronounced overlaps in substrate recognition not only between proteases closely related on sequence basis but also between proteolytic enzymes of different evolutionary origin and catalytic type. To illustrate the applicability of our approach we analyze the distribution of targets of small molecules from the ChEMBL database in our substrate-based protease specificity trees. We observe a striking clustering of annotated targets in tree branches even though these grouped targets do not necessarily share similarity on protein sequence level. This highlights the value and applicability of knowledge acquired from peptide substrates in drug design of small molecules, e.g., for the prediction of off-target effects or drug repurposing. Consequently, our similarity metric allows to map the degradome and its associated drug target network via comparison of known substrate peptides. The substrate-driven view of protein-protein interfaces is not limited to the field of proteases but can be applied to any target class where a sufficient amount of known substrate data is available.

Shaping the interaction landscape of bioactive molecules

MOTIVATION

Most bioactive molecules perform their action by interacting with proteins or other macromolecules. However, for a significant fraction of them, the primary target remains unknown. In addition, the majority of bioactive molecules have more than one target, many of which are poorly characterized. Computational predictions of bioactive molecule targets based on similarity with known ligands are powerful to narrow down the number of potential targets and to rationalize side effects of known molecules.

RESULTS

Using a reference set of 224 412 molecules active on 1700 human proteins, we show that accurate target prediction can be achieved by combining different measures of chemical similarity based on both chemical structure and molecular shape. Our results indicate that the combined approach is especially efficient when no ligand with the same scaffold or from the same chemical series has yet been discovered. We also observe that different combinations of similarity measures are optimal for different molecular properties, such as the number of heavy atoms. This further highlights the importance of considering different classes of similarity measures between new molecules and known ligands to accurately predict their targets.

ChemMapper: a versatile web server for exploring pharmacology and chemical structure association based on molecular 3D similarity method

SUMMARY

ChemMapper is an online platform to predict polypharmacology effect and mode of action for small molecules based on 3D similarity computation. ChemMapper collects >350 000 chemical structures with bioactivities and associated target annotations (as well as >3 000 000 non-annotated compounds for virtual screening). Taking the user-provided chemical structure as the query, the top most similar compounds in terms of 3D similarity are returned with associated pharmacology annotations. ChemMapper is designed to provide versatile services in a variety of chemogenomics, drug repurposing, polypharmacology, novel bioactive compounds identification and scaffold hopping studies.

AVAILABILITY

http://lilab.ecust.edu.cn/chemmapper/.

CONTACT

xfliu@ecust.edu.cn or hlli@ecust.edu.cn

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Bioturbo similarity searching: combining chemical and biological similarity to discover structurally diverse bioactive molecules

Virtual screening using bioactivity profiles has become an integral part of currently applied hit finding methods in pharmaceutical industry. However, a significant drawback of this approach is that it is only applicable to compounds that have been biologically tested in the past and have sufficient activity annotations for meaningful profile comparisons. Although bioactivity data generated in pharmaceutical institutions are growing on an unprecedented scale, the number of biologically annotated compounds still covers only a minuscule fraction of chemical space. For a newly synthesized compound or an isolated natural product to be biologically characterized across multiple assays, it may take a considerable amount of time. Consequently, this chemical matter will not be included in virtual screening campaigns based on bioactivity profiles. To overcome this problem, we herein introduce bioturbo similarity searching that uses chemical similarity to map molecules without biological annotations into bioactivity space and then searches for biologically similar compounds in this reference system. In benchmark calculations on primary screening data, we demonstrate that our approach generally achieves higher hit rates and identifies structurally more diverse compounds than approaches using chemical information only. Furthermore, our method is able to discover hits with novel modes of inhibition that traditional 2D and 3D similarity approaches are unlikely to discover. Test calculations on a set of natural products reveal the practical utility of the approach for identifying novel and synthetically more accessible chemical matter.

3D molecular descriptors important for clinical success

The pharmacokinetic and safety profiles of clinical drug candidates are greatly influenced by their requisite physicochemical properties. In particular, it has been shown that 2D molecular descriptors such as fraction of Sp3 carbon atoms (Fsp3) and number of stereo centers correlate with clinical success. Using the proteomic off-target hit rate of nicotinic ligands, we found that shape-based 3D descriptors such as the radius of gyration and shadow indices discriminate off-target promiscuity better than do Fsp3 and the number of stereo centers. We have deduced the relevant descriptor values required for a ligand to be nonpromiscuous. Investigating the MDL Drug Data Report (MDDR) database as compounds move from the preclinical stage toward the market, we have found that these shape-based 3D descriptors predict clinical success of compounds at preclinical and phase1 stages vs compounds withdrawn from the market better than do Fsp3 and LogD. Further, these computed 3D molecular descriptors correlate well with experimentally observed solubility, which is among well-known physicochemical properties that drive clinical success. We also found that about 84% of launched drugs satisfy either Shadow index or Fsp3 criteria, whereas withdrawn and discontinued compounds fail to meet the same criteria. Our studies suggest that spherical compounds (rather than their elongated counterparts) with a minimal number of aromatic rings may exhibit a high propensity to advance from clinical trials to market.

A targeted library screen reveals a new inhibitor scaffold for protein kinase D

Protein kinase D (PKD) has emerged as a potential therapeutic target in multiple pathological conditions, including cancer and heart diseases. Potent and selective small molecule inhibitors of PKD are valuable for dissecting PKD-mediated cellular signaling pathways and for therapeutic application. In this study, we evaluated a targeted library of 235 small organic kinase inhibitors for PKD1 inhibitory activity at a single concentration. Twenty-eight PKD inhibitory chemotypes were identified and six exhibited excellent PKD1 selectivity. Five of the six lead structures share a common scaffold, with compound 139 being the most potent and selective for PKD vs PKC and CAMK. Compound 139 was an ATP-competitive PKD1 inhibitor with a low double-digit nanomolar potency and was also cell-active. Kinase profiling analysis identified this class of small molecules as pan-PKD inhibitors, confirmed their selectivity again PKC and CAMK, and demonstrated an overall favorable selectivity profile that could be further enhanced through structural modification. Furthermore, using a PKD homology model based on similar protein kinase structures, docking modes for compound 139 were explored and compared to literature examples of PKD inhibition. Modeling of these compounds at the ATP-binding site of PKD was used to rationalize its high potency and provide the foundation for future further optimization. Accordingly, using biochemical screening of a small number of privileged scaffolds and computational modeling, we have identified a new core structure for highly potent PKD inhibition with promising selectivity against closely related kinases. These lead structures represent an excellent starting point for the further optimization and the design of selective and therapeutically effective small molecule inhibitors of PKD.

Does your model weigh the same as a duck?

Computer-aided drug design is a mature field by some measures, and it has produced notable successes that underpin the study of interactions between small molecules and living systems. However, unlike a truly mature field, fallacies of logic lie at the heart of the arguments in support of major lines of research on methodology and validation thereof. Two particularly pernicious ones are cum hoc ergo propter hoc (with this, therefore because of this) and confirmation bias (seeking evidence that is confirmatory of the hypothesis at hand). These fallacies will be discussed in the context of off-target predictive modeling, QSAR, molecular similarity computations, and docking. Examples will be shown that avoid these problems.

Identifying and characterizing promiscuous targets: implications for virtual screening

INTRODUCTION

Ligand-based shape matching approaches have become established as important and popular virtual screening (VS) techniques. However, despite their relative success, the question of how to best choose the initial query compounds and their conformations remains largely unsolved. This issue gains importance when dealing with promiscuous targets, that is, proteins that bind multiple ligand scaffold families in one or more binding site. Conventional shape matching VS approaches assume that there is only one binding mode for a given protein target. This may be true for some targets, but it is certainly not true in all cases. Several recent studies have shown that some protein targets bind to different ligands in different ways.

AREAS COVERED

The authors discuss the concept of promiscuity in the context of virtual drug screening, and present and analyze several examples of promiscuous targets. The article also reports on the impact of the query conformation on the performance of shape-based VS and the potential to improve VS performance by using consensus shape clustering techniques.

EXPERT OPINION

The notion of polypharmacology is becoming highly relevant in drug discovery. Understanding and exploiting promiscuity present challenges and opportunities for drug discovery endeavors. The examples of promiscuity presented here suggest that promiscuous targets and ligands are much more common than previously assumed, and this should be taken into account in practical VS protocols. Although some progress has been made, there is a need to develop more sophisticated computational techniques and protocols that can identify and characterize promiscuous targets on a genomic scale.