Training based on ligand efficiency improves prediction of bioactivities of ligands and drug target proteins in a machine learning approach

Antifungal peptides produced by actinomycetes and their biological activities against plant diseases

Antibacterial peptides are a class of naturally occurring peptides produced by eukaryotes and prokaryotes. Some of them exhibit broad-spectrum antifungal activity. Antifungal peptides (AFPs) can be developed as antibiotic to control fungal infections in agriculture due to their different antifungal mechanisms. As actinomycetes are still one of the most important sources of novel antibiotics, in this review, the mechanisms of action of AFPs are explained. Characterization of several AFPs produced by actinomycetes and their biological activities against plant diseases are summarized. Furthermore, the pathway for total synthesis of naturally occurring cyclodepsipeptide, valinomycin, is proposed. Finally, the pathway for biosynthesis of kutzneride 2 is proposed and the structure-activity relationship of kutznerides is discussed.

An Application of Fit Quality to Screen MDM2/p53 Protein-Protein Interaction Inhibitors

The judicious application of ligand or binding efficiency (LE) metrics, which quantify the molecular properties required to obtain binding affinity for a drug target, is gaining traction in the selection and optimization of fragments, hits and leads. Here we report for the first time the use of LE based metric, fit quality (FQ), in virtual screening (VS) of MDM2/p53 protein-protein interaction inhibitors (PPIIs). Firstly, a Receptor-Ligand pharmacophore model was constructed on multiple MDM2/ligand complex structures to screen the library. The enrichment factor (EF) for screening was calculated based on a decoy set to define the screening threshold. Finally, 1% of the library, 335 compounds, were screened and re-filtered with the FQ metric. According to the statistical results of FQ vs. activity of 156 MDM2/p53 PPIIs extracted from literatures, the cut-off was defined as FQ = 0.8. After the second round of VS, six compounds with the FQ > 0.8 were picked out for assessing their antitumor activity. At the cellular level, the six hits exhibited a good selectivity (larger than 3) against HepG2 (wt-p53) vs. Hep3B (p53 null) cell lines. On the further study, the six hits exhibited an acceptable affinity (range of K_i from 10² to 10³ nM) to MDM2 when comparing to Nutlin-3a. Based on our work, FQ based VS strategy could be applied to discover other PPIIs.

Predicted Biological Activity of Purchasable Chemical Space

Whereas 400 million distinct compounds are now purchasable within the span of a few weeks, the biological activities of most are unknown. To facilitate access to new chemistry for biology, we have combined the Similarity Ensemble Approach (SEA) with the maximum Tanimoto similarity to the nearest bioactive to predict activity for every commercially available molecule in ZINC. This method, which we label SEA+TC, outperforms both SEA and a naïve-Bayesian classifier via predictive performance on a 5-fold cross-validation of ChEMBL’s bioactivity data set (version 21). Using this method, predictions for over 40% of compounds (>160 million) have either high significance (pSEA ≥ 40), high similarity (ECFP4MaxTc ≥ 0.4), or both, for one or more of 1382 targets well described by ligands in the literature. Using a further 1347 less-well-described targets, we predict activities for an additional 11 million compounds. To gauge whether these predictions are sensible, we investigate 75 predictions for 50 drugs lacking a binding affinity annotation in ChEMBL. The 535 million predictions for over 171 million compounds at 2629 targets are linked to purchasing information and evidence to support each prediction and are freely available via https://zinc15.docking.org and https://files.docking.org.

Consensus scoring model for the molecular docking study of mTOR kinase inhibitor

The discovery of mammalian target of rapamycin (mTOR) kinase inhibitors has always been a research hotspot of antitumor drugs. Consensus scoring used in the docking study of mTOR kinase inhibitors usually improves hit rate of virtual screening. Herein, we attempt to build a series of consensus scoring models based on a set of the common scoring functions. In this paper, twenty-five kinds of mTOR inhibitors (16 clinical candidate compounds and 9 promising preclinical compounds) are carefully collected, and selected for the molecular docking study used by the Glide docking programs within the standard precise (SP) mode. The predicted poses of these ligands are saved, and revaluated by twenty-six available scoring functions, respectively. Subsequently, consensus scoring models are trained based on the obtained rescoring results by the partial least squares (PLS) method, and validated by Leave-one-out (LOO) method. In addition, three kinds of ligand efficiency indices (BEI, SEI, and LLE) instead of pIC₅₀ as the activity could greatly improve the statistical quality of build models. Two best calculated models 10 and 22 using the same BEI indice have following statistical parameters, respectively: for model 10, training set R²=0.767, Q²=0.647, RMSE=0.024, and for test set R²=0.932, RMSE=0.026; for model 22, raining set R²=0.790, Q²=0.627, RMSE=0.023, and for test set R²=0.955, RMSE=0.020. These two consensus scoring model would be used for the docking virtual screening of novel mTOR inhibitors.

Large-Scale Prediction of Drug-Target Interaction: a Data-Centric Review

The prediction of drug-target interactions (DTIs) is of extraordinary significance to modern drug discovery in terms of suggesting new drug candidates and repositioning old drugs. Despite technological advances, large-scale experimental determination of DTIs is still expensive and laborious. Effective and low-cost computational alternatives remain in strong need. Meanwhile, open-access resources have been rapidly growing with massive amount of bioactivity data becoming available, creating unprecedented opportunities for the development of novel in silico models for large-scale DTI prediction. In this work, we review the state-of-the-art computational approaches for identifying DTIs from a data-centric perspective: what the underlying data are and how they are utilized in each study. We also summarize popular public data resources and online tools for DTI prediction. It is found that various types of data were employed including properties of chemical structures, drug therapeutic effects and side effects, drug-target binding, drug-drug interactions, bioactivity data of drug molecules across multiple biological targets, and drug-induced gene expressions. More often, the heterogeneous data were integrated to offer better performance. However, challenges remain such as handling data imbalance, incorporating negative samples and quantitative bioactivity data, as well as maintaining cross-links among different data sources, which are essential for large-scale and automated information integration.

Ligand efficiency metrics in drug discovery: the pros and cons from a practical perspective

INTRODUCTION

Ligand efficiency metrics are almost universally accepted as a valuable indicator of compound quality and an aid to reduce attrition. Areas covered: In this review, the authors describe ligand efficiency metrics giving a balanced overview on their merits and points of weakness in order to enable the readers to gain an informed opinion. Relevant theoretical breakthroughs and drug-like properties are also illustrated. Several recent exemplary case studies are discussed in order to illustrate the main fields of application of ligand efficiency metrics. Expert opinion: As a medicinal chemist guide, ligand efficiency metrics perform in a context- and chemotype-dependent manner; thus, they should not be used as a magic box. Since the 'big bang' of efficiency metrics occurred more or less ten years ago and the average time to develop a new drug is over the same period, the next few years will give a clearer outlook on the increased rate of success, if any, gained by means of these new intriguing tools.

Between Descriptors and Properties: Understanding the Ligand Efficiency Trends for G Protein-Coupled Receptor and Kinase Structure-Activity Data Sets

The chemical meaning of the ligand efficiency (LE) metrics is explained in this paper using a large G protein-coupled receptor (GPCR) and kinase structure-activity (IC₅₀, K_i) data set. Although there is a controversy in the literature regarding both the mathematical validity and the performance of LE, it is in common use as an early estimator for drug optimization. Apparently, the numerous con arguments are not convincing enough. We show here for the first time that the main misunderstanding of the chemical meaning of LE is its interpretation as a molecular descriptor connected with a single molecule. Instead, LE should be interpreted as a statistical property. We show that the LE, which is designed as a regression of a binding property on the heavy atom count (HAC), is correlated to the reciprocal of the molecular weight because of Avogadro statistics. This indicates that the hyperbolic model of LE is basically a consequence of a nonbinding effect, an increase in the number of ligands that are available to a receptor for smaller molecules, and not a real increase in the binding potency for a single HAC as interpreted in the literature. Accordingly, we need to revisit and carefully reevaluate LE-based molecular comparisons.

Combination therapeutics in complex diseases

The biological redundancies in molecular networks of complex diseases limit the efficacy of many single drug therapies. Combination therapeutics, as a common therapeutic method, involve pharmacological intervention using several drugs that interact with multiple targets in the molecular networks of diseases and may achieve better efficacy and/or less toxicity than monotherapy in practice. The development of combination therapeutics is complicated by several critical issues, including identifying multiple targets, targeting strategies and the drug combination. This review summarizes the current achievements in combination therapeutics, with a particular emphasis on the efforts to develop combination therapeutics for complex diseases.

Accurate and efficient target prediction using a potency-sensitive influence-relevance voter

BACKGROUND

A number of algorithms have been proposed to predict the biological targets of diverse molecules. Some are structure-based, but the most common are ligand-based and use chemical fingerprints and the notion of chemical similarity. These methods tend to be computationally faster than others, making them particularly attractive tools as the amount of available data grows.

RESULTS

Using a ChEMBL-derived database covering 490,760 molecule-protein interactions and 3236 protein targets, we conduct a large-scale assessment of the performance of several target-prediction algorithms at predicting drug-target activity. We assess algorithm performance using three validation procedures: standard tenfold cross-validation, tenfold cross-validation in a simulated screen that includes random inactive molecules, and validation on an external test set composed of molecules not present in our database.

CONCLUSIONS

We present two improvements over current practice. First, using a modified version of the influence-relevance voter (IRV), we show that using molecule potency data can improve target prediction. Second, we demonstrate that random inactive molecules added during training can boost the accuracy of several algorithms in realistic target-prediction experiments. Our potency-sensitive version of the IRV (PS-IRV) obtains the best results on large test sets in most of the experiments. Models and software are publicly accessible through the chemoinformatics portal at http://chemdb.ics.uci.edu/.

Identification of potential inhibitors based on compound proposal contest: Tyrosine-protein kinase Yes as a target

A search of broader range of chemical space is important for drug discovery. Different methods of computer-aided drug discovery (CADD) are known to propose compounds in different chemical spaces as hit molecules for the same target protein. This study aimed at using multiple CADD methods through open innovation to achieve a level of hit molecule diversity that is not achievable with any particular single method. We held a compound proposal contest, in which multiple research groups participated and predicted inhibitors of tyrosine-protein kinase Yes. This showed whether collective knowledge based on individual approaches helped to obtain hit compounds from a broad range of chemical space and whether the contest-based approach was effective.

Ligand efficiency-based support vector regression models for predicting bioactivities of ligands to drug target proteins

The concept of ligand efficiency (LE) indices is widely accepted throughout the drug design community and is frequently used in a retrospective manner in the process of drug development. For example, LE indices are used to investigate LE optimization processes of already-approved drugs and to re-evaluate hit compounds obtained from structure-based virtual screening methods and/or high-throughput experimental assays. However, LE indices could also be applied in a prospective manner to explore drug candidates. Here, we describe the construction of machine learning-based regression models in which LE indices are adopted as an end point and show that LE-based regression models can outperform regression models based on pIC50 values. In addition to pIC50 values traditionally used in machine learning studies based on chemogenomics data, three representative LE indices (ligand lipophilicity efficiency (LLE), binding efficiency index (BEI), and surface efficiency index (SEI)) were adopted, then used to create four types of training data. We constructed regression models by applying a support vector regression (SVR) method to the training data. In cross-validation tests of the SVR models, the LE-based SVR models showed higher correlations between the observed and predicted values than the pIC50-based models. Application tests to new data displayed that, generally, the predictive performance of SVR models follows the order SEI > BEI > LLE > pIC50. Close examination of the distributions of the activity values (pIC50, LLE, BEI, and SEI) in the training and validation data implied that the performance order of the SVR models may be ascribed to the much higher diversity of the LE-based training and validation data. In the application tests, the LE-based SVR models can offer better predictive performance of compound-protein pairs with a wider range of ligand potencies than the pIC50-based models. This finding strongly suggests that LE-based SVR models are better than pIC50-based models at predicting bioactivities of compounds that could exhibit a much higher (or lower) potency.