The state of the art in secondary pharmacology and its impact on the safety of new medicines

Secondary pharmacology screening of investigational small-molecule drugs for potentially adverse off-target activities has become standard practice in pharmaceutical research and development, and regulatory agencies are increasingly requesting data on activity against targets with recognized adverse effect relationships. However, the screening strategies and target panels used by pharmaceutical companies may vary substantially. To help identify commonalities and differences, as well as to highlight opportunities for further optimization of secondary pharmacology assessment, we conducted a broad-ranging survey across 18 companies under the auspices of the DruSafe leadership group of the International Consortium for Innovation and Quality in Pharmaceutical Development. Based on our analysis of this survey and discussions and additional research within the group, we present here an overview of the current state of the art in secondary pharmacology screening. We discuss best practices, including additional safety-associated targets not covered by most current screening panels, and present approaches for interpreting and reporting off-target activities. We also provide an assessment of the safety impact of secondary pharmacology screening, and a perspective on opportunities and challenges in this rapidly developing field.

TERT activation targets DNA methylation and multiple aging hallmarks

Insufficient telomerase activity, stemming from low telomerase reverse transcriptase (TERT) gene transcription, contributes to telomere dysfunction and aging pathologies. Besides its traditional function in telomere synthesis, TERT acts as a transcriptional co-regulator of genes pivotal in aging and age-associated diseases. Here, we report the identification of a TERT activator compound (TAC) that upregulates TERT transcription via the MEK/ERK/AP-1 cascade. In primary human cells and naturally aged mice, TAC-induced elevation of TERT levels promotes telomere synthesis, blunts tissue aging hallmarks with reduced cellular senescence and inflammatory cytokines, and silences p16INK4a expression via upregulation of DNMT3B-mediated promoter hypermethylation. In the brain, TAC alleviates neuroinflammation, increases neurotrophic factors, stimulates adult neurogenesis, and preserves cognitive function without evident toxicity, including cancer risk. Together, these findings underscore TERT's critical role in aging processes and provide preclinical proof of concept for physiological TERT activation as a strategy to mitigate multiple aging hallmarks and associated pathologies.

Computer-aided pattern scoring - A multitarget dataset-driven workflow to predict ligands of orphan targets

The identification of lead molecules and the exploration of novel pharmacological drug targets are major challenges of medical life sciences today. Genome-wide association studies, multi-omics, and systems pharmacology steadily reveal new protein networks, extending the known and relevant disease-modifying proteome. Unfortunately, the vast majority of the disease-modifying proteome consists of 'orphan targets' of which intrinsic ligands/substrates, (patho)physiological roles, and/or modulators are unknown. Undruggability is a major challenge in drug development today, and medicinal chemistry efforts cannot keep up with hit identification and hit-to-lead optimization studies. New 'thinking-outside-the-box' approaches are necessary to identify structurally novel and functionally distinctive ligands for orphan targets. Here we present a unique dataset that includes critical information on the orphan target ABCA1, from which a novel cheminformatic workflow - computer-aided pattern scoring (C@PS) - for the identification of novel ligands was developed. Providing a hit rate of 95.5% and molecules with high potency and molecular-structural diversity, this dataset represents a suitable template for general deorphanization studies.

Medicinal Polypharmacology in the Clinic - Translating the Polypharmacolome into Therapeutic Benefit

Drugs with multiple targets, often annotated as 'unselective', 'promiscuous', 'multitarget', or 'polypharmacological', are widely considered in both academic and industrial research as a high risk due to the likelihood of adverse effects. However, retrospective analyses have shown that particularly approved drugs bear rich polypharmacological profiles. This raises the question whether our perception of the specificity paradigm ('one drug-one target concept') is correct - and if specifically multitarget drugs should be developed instead of being rejected. These questions provoke a paradigm shift - regarding the development of polypharmacological drugs not as a 'waste of investment', but acknowledging the existence of a 'lack of investment'. This perspective provides an insight into modern drug development highlighting latest drug candidates that have not been assessed in a broader polypharmacology-based context elsewhere embedded in a historic framework of classical and modern approved multitarget drugs. The article shall be an inspiration to the scientific community to re-consider current standards, and more, to evolve to a better understanding of polypharmacology from a challenge to an opportunity.

Medicinal polypharmacology: Exploration and exploitation of the polypharmacolome in modern drug development

At the core of complex and multifactorial human diseases, such as cancer, metabolic syndrome, or neurodegeneration, are multiple players that cross-talk in robust biological networks which are intrinsically resilient to alterations. These multifactorial diseases are characterized by sophisticated feedback mechanisms which manifest cellular imbalance and resistance to drug therapy. By adhering to the specificity paradigm ("one target-one drug concept"), research focused for many years on drugs with very narrow mechanisms of action. This narrow focus promoted therapy ineffectiveness and resistance. However, modern drug discovery has evolved over the last years, increasingly emphasizing integral strategies for the development of clinically effective drugs. These integral strategies include the controlled engagement of multiple targets to overcome therapy resistance. Apart from the additive or even synergistic effects in therapy, multitarget drugs harbor molecular-structural attributes to explore orphan targets of which intrinsic substrates/physiological role(s) and/or modulators are unknown for future therapy purposes. We designated this multidisciplinary and translational research field between medicinal chemistry, chemical biology, and molecular pharmacology as 'medicinal polypharmacology'. Medicinal polypharmacology emerged as alternative approach to common single-targeted pharmacology stretching from basic drug and target identification processes to clinical evaluation of multitarget drugs, and the exploration and exploitation of the 'polypharmacolome' is at the forefront of modern drug development research.

Multi-level advances in databases related to systems pharmacology in traditional Chinese medicine: a 60-year review

The therapeutic effects of traditional Chinese medicine (TCM) involve intricate interactions among multiple components and targets. Currently, computational approaches play a pivotal role in simulating various pharmacological processes of TCM. The application of network analysis in TCM research has provided an effective means to explain the pharmacological mechanisms underlying the actions of herbs or formulas through the lens of biological network analysis. Along with the advances of network analysis, computational science has coalesced around the core chain of TCM research: formula-herb-component-target-phenotype-ZHENG, facilitating the accumulation and organization of the extensive TCM-related data and the establishment of relevant databases. Nonetheless, recent years have witnessed a tendency toward homogeneity in the development and application of these databases. Advancements in computational technologies, including deep learning and foundation model, have propelled the exploration and modeling of intricate systems into a new phase, potentially heralding a new era. This review aims to delves into the progress made in databases related to six key entities: formula, herb, component, target, phenotype, and ZHENG. Systematically discussions on the commonalities and disparities among various database types were presented. In addition, the review raised the issue of research bottleneck in TCM computational pharmacology and envisions the forthcoming directions of computational research within the realm of TCM.

Promiscuity in drug discovery on the verge of the structural revolution: recent advances and future chances

INTRODUCTION

Promiscuity denotes the ability of ligands and targets to specifically interact with multiple binding partners. Despite negative aspects like side effects, promiscuity is receiving increasing attention in drug discovery as it can enhance drug efficacy and provides a molecular basis for drug repositioning. The three-dimensional structure of ligand-target complexes delivers exclusive insights into the molecular mechanisms of promiscuity and structure-based methods enable the identification of promiscuous interactions. With the recent breakthrough in protein structure prediction, novel possibilities open up to reveal unknown connections in ligand-target interaction networks.

AREAS COVERED

This review highlights the significance of structure in the identification and characterization of promiscuity and evaluates the potential of protein structure prediction to advance our knowledge of drug-target interaction networks. It discusses the definition and relevance of promiscuity in drug discovery and explores different approaches to detecting promiscuous ligands and targets.

EXPERT OPINION

Examination of structural data is essential for understanding and quantifying promiscuity. The recent advancements in structure prediction have resulted in an abundance of targets that are well-suited for structure-based methods like docking. In silico approaches may eventually completely transform our understanding of drug-target networks by complementing the millions of predicted protein structures with billions of predicted drug-target interactions.

Prospects for the application of traditional Chinese medicine network pharmacology in food science research

There has always been a particular difficulty with in-depth research on the mechanisms of food nutrition and bioactivity. The main function of food is to meet the nutritional needs of the human body, rather than to exert a therapeutic effect. Its relatively modest biological activity makes it difficult to study from the perspective of general pharmacological models. With the popularity of functional foods and the concept of dietary therapy, and the development of information and multi-omics technology in food research, research into these mechanisms is moving towards a more microscopic future. Network pharmacology has accumulated nearly 20 years of research experience in traditional Chinese medicine (TCM), and there has been no shortage of work from this perspective on the medicinal functions of food. Given the similarity between the concept of 'multi-component-multi-target' properties of food and TCM, we think that network pharmacology is applicable to the study of the complex mechanisms of food. Here we review the development of network pharmacology, summarize its application to 'medicine and food homology', and propose a methodology based on food characteristics for the first time, demonstrating its feasibility for food research. © 2023 Society of Chemical Industry.

Targeting RNA with small molecules-A safety perspective

RNA is a major player in cellular function, and consequently can drive a number of disease pathologies. Over the past several years, small molecule-RNA targeting (smRNA targeting) has developed into a promising drug discovery approach. Numerous techniques, tools, and assays have been developed to support this field, and significant investments have been made by pharmaceutical and biotechnology companies. To date, the focus has been on identifying disease validated primary targets for smRNA drug development, yet RNA as a secondary (off) target for all small molecule drug programs largely has been unexplored. In this perspective, we discuss structure, target, and mechanism-driven safety aspects of smRNAs and highlight how these parameters can be evaluated in drug discovery programs to produce potentially safer drugs.

Communalistic use of psychoactive plants as a bridge between traditional healing practices and Western medicine: A new path for the Global Mental Health movement

The Global Mental Health (GMH) movement aims to provide urgently needed treatment to those with mental illness, especially in low- and middle-income countries. Due to the complexity of providing mental health services to people from various cultures, there is much debate among GMH advocates regarding the best way to proceed. While biomedical interventions offer some degree of help, complementary approaches should focus on the social/community aspects. Many cultures conduct traditional rituals involving the communal use of psychoactive plants. We propose that these practices should be respected, protected, and promoted as valuable tools with regard to mental health care at the community level. The traditional use of psychoactive plants promotes community engagement and participation, and they are relatively affordable. Furthermore, the worldviews and meaning-making systems of local population are respected. The medical systems surrounding the use of psychoactive plants can be explained in biomedical terms, and many recently published clinical trials have demonstrated their therapeutic potential. Psychoactive plants and associated rituals offer potential benefits as complementary aspects of mental health services. They should be considered as such by international practitioners and advocates of the GMH movement.

The Use of Cannabis sativa L. for Pest Control: From the Ethnobotanical Knowledge to a Systematic Review of Experimental Studies

Background: Despite the benefits that synthetic pesticides have provided in terms of pest and disease control, they cause serious long-term consequences for both the environment and living organisms. Interest in eco-friendly products has subsequently increased in recent years. Methods: This article briefly analyzes the available ethnobotanical evidence regarding the use of Cannabis sativa as a pesticide and offers a systematic review of experimental studies. Results: Our findings indicate that both ethnobotanical and experimental procedures support the use of C. sativa as a pesticide, as remarkable toxicity has been observed against pest organisms. The results included in the systematic review of experimental studies (n=30) show a high degree of heterogeneity, but certain conclusions can be extracted to guide further research. For instance, promising pesticide properties were reported for most of the groups of species tested, especially Arachnida and Insecta; the efficacy of C. sativa as a pesticide can be derived from a wide variety of compounds that it contains and possible synergistic effects; it is crucial to standardize the phytochemical profile of C. sativa plants used as well as to obtain easily reproducible results; appropriate extraction methods should be explored; and upper inflorescences of the plant may be preferred for the production of the essential oil, but further studies should explore better other parts of the plant. Conclusion: In the coming years, as new findings are produced, the promising potential of C. sativa as a pesticide will be elucidated, and reviews such as the present one constitute useful basic tools to make these processes easier.

Comprehensive Assessment of Indian Variations in the Druggable Kinome Landscape Highlights Distinct Insights at the Sequence, Structure and Pharmacogenomic Stratum

India confines more than 17% of the world’s population and has a diverse genetic makeup with several clinically relevant rare mutations belonging to many sub-group which are undervalued in global sequencing datasets like the 1000 Genome data (1KG) containing limited samples for Indian ethnicity. Such databases are critical for the pharmaceutical and drug development industry where diversity plays a crucial role in identifying genetic disposition towards adverse drug reactions. A qualitative and comparative sequence and structural study utilizing variant information present in the recently published, largest curated Indian genome database (IndiGen) and the 1000 Genome data was performed for variants belonging to the kinase coding genes, the second most targeted group of drug targets. The sequence-level analysis identified similarities and differences among different populations based on the nsSNVs and amino acid exchange frequencies whereas a comparative structural analysis of IndiGen variants was performed with pathogenic variants reported in UniProtKB Humsavar data. The influence of these variations on structural features of the protein, such as structural stability, solvent accessibility, hydrophobicity, and the hydrogen-bond network was investigated. In-silico screening of the known drugs to these Indian variation-containing proteins reveals critical differences imparted in the strength of binding due to the variations present in the Indian population. In conclusion, this study constitutes a comprehensive investigation into the understanding of common variations present in the second largest population in the world and investigating its implications in the sequence, structural and pharmacogenomic landscape. The preliminary investigation reported in this paper, supporting the screening and detection of ADRs specific to the Indian population could aid in the development of techniques for pre-clinical and post-market screening of drug-related adverse events in the Indian population.

The Tunicate Metabolite 2-(3,5-Diiodo-4-methoxyphenyl)ethan-1-amine Targets Ion Channels of Vertebrate Sensory Neurons

Marine tunicates produce defensive amino-acid-derived metabolites, including 2-(3,5-diiodo-4-methoxyphenyl)ethan-1-amine (DIMTA), but their mechanisms of action are rarely known. Using an assay-guided approach, we found that out of the many different sensory cells in the mouse dorsal root ganglion (DRG), DIMTA selectively affected low-threshold cold thermosensors. Whole-cell electrophysiology experiments using DRG cells, channels expressed in Xenopus oocytes, and human cell lines revealed that DIMTA blocks several potassium channels, reducing the magnitude of the afterhyperpolarization and increasing the baseline intracellular calcium concentration [Ca²⁺]_i of low-threshold cold thermosensors. When injected into mice, DIMTA increased the threshold of cold sensation by >3 °C. DIMTA may thus serve as a lead in the further design of compounds that inhibit problems in the cold-sensory system, such as cold allodynia and other neuropathic pain conditions.

GPCR_LigandClassify.py; a rigorous machine learning classifier for GPCR targeting compounds

The current study describes the construction of various ligand-based machine learning models to be used for drug-repurposing against the family of G-Protein Coupled Receptors (GPCRs). In building these models, we collected > 500,000 data points, encompassing experimentally measured molecular association data of > 160,000 unique ligands against > 250 GPCRs. These data points were retrieved from the GPCR-Ligand Association (GLASS) database. We have used diverse molecular featurization methods to describe the input molecules. Multiple supervised ML algorithms were developed, tested and compared for their accuracy, F scores, as well as for their Matthews' correlation coefficient scores (MCC). Our data suggest that combined with molecular fingerprinting, ensemble decision trees and gradient boosted trees ML algorithms are on the accuracy border of the rather sophisticated deep neural nets (DNNs)-based algorithms. On a test dataset, these models displayed an excellent performance, reaching a ~ 90% classification accuracy. Additionally, we showcase a few examples where our models were able to identify interesting connections between known drugs from the Drug-Bank database and members of the GPCR family of receptors. Our findings are in excellent agreement with previously reported experimental observations in the literature. We hope the models presented in this paper synergize with the currently ongoing interest of applying machine learning modeling in the field of drug repurposing and computational drug discovery in general.

Melatonin controversies, an update

Melatonin was discovered more than 60 years ago. Since then, several seminal discoveries have allowed us to define its function as a neuroendocrine hormone and its molecular targets in mammals and many other species. However, many fundamental issues have not yet been solved such as the subcellular localization of melatonin synthesis and the full spectrum of its molecular targets. In addition, a considerable number of controversies persist in the field, mainly concerning how many functions melatonin has. Altogether, this illustrates how "immature" the field still is. The intention of this opinion article is to note the controversies and limitations in the field, to initiate a discussion and to make proposals/guidelines to overcome them and move the field forward.

Therapeutic Potential of Natural Psychoactive Drugs for Central Nervous System Disorders: A Perspective from Polypharmacology

In the drug development, the formation of highly selective ligands has been unsuccessful in the treatment of central nervous system disorders. Multi-target ligands, from the polypharmacology paradigm, are being proposed as treatments for these complex disorders, since they offer enhanced efficacy and a strong safety profile. Natural products are the best examples of multi-target compounds, so they are of high interest within this paradigm. Additionally, recent research on psychoactive drugs of natural origin, such as ayahuasca and cannabis, has demonstrated the promising therapeutic potential for the treatment of some psychiatric and neurological disorders. In this text, we describe how research on psychoactive drugs can be effectively combined with the polypharmacology paradigm, providing ayahuasca and cannabis research as examples. The advantages and disadvantages are also discussed.

Finding New Molecular Targets of Familiar Natural Products Using In Silico Target Prediction

Natural products comprise a rich reservoir for innovative drug leads and are a constant source of bioactive compounds. To find pharmacological targets for new or already known natural products using modern computer-aided methods is a current endeavor in drug discovery. Nature's treasures, however, could be used more effectively. Yet, reliable pipelines for the large-scale target prediction of natural products are still rare. We developed an in silico workflow consisting of four independent, stand-alone target prediction tools and evaluated its performance on dihydrochalcones (DHCs)-a well-known class of natural products. Thereby, we revealed four previously unreported protein targets for DHCs, namely 5-lipoxygenase, cyclooxygenase-1, 17β-hydroxysteroid dehydrogenase 3, and aldo-keto reductase 1C3. Moreover, we provide a thorough strategy on how to perform computational target predictions and guidance on using the respective tools.

Drugs, host proteins and viral proteins: how their promiscuities shape antiviral design

The reciprocal nature of drug specificity and target specificity implies that the same is true for their respective promiscuities. Protein promiscuity has two broadly different types of footprint in drug design. The first is relaxed specificity of binding sites for substrates, inhibitors, effectors or cofactors. The second involves protein-protein interactions of regulatory processes such as signal transduction and transcription, and here protein intrinsic disorder plays an important role. Both viruses and host cells exploit intrinsic disorder for their survival, as do the design and discovery programs for antivirals. Drug action, strictly speaking, always relies upon promiscuous activity, with drug promiscuity enlarging its scope. Drug repurposing searches for additional promiscuity on the part of both the drug and the target in the host. Understanding the subtle nuances of these promiscuities is critical in the design of novel and more effective antivirals.

Daurichromenic Acid from the Chinese Traditional Medicinal Plant Rhododendron dauricum Inhibits Sphingomyelin Synthase and Aβ Aggregation

Species of the genus Rhododendron have been used in traditional Chinese medicine, with the medicinal herb "Manshanfong" used as an expectorant and for the treatment of acute bronchitis. Daurichromenic acid (DCA), a constituent of Rhododendron dauricum, is a meroterpenoid with antibacterial, anti-HIV, and anti-inflammatory activities. However, the mechanisms underlying these pharmacologic activities are poorly understood. To develop new drugs based on DCA, more information is required regarding its interactions with biomolecules. The present study showed that DCA inhibits the activity of the enzyme sphingomyelin synthase, with an IC₅₀ of 4 µM. The structure-activity relationships between DCA and sphingomyelin synthase were evaluated using derivatives and cyclized hongoquercin A. In addition, DCA was found to inhibit amyloid β aggregation. These results may help in the design of effective drugs based on DCA.

X-ray Structure-Based Chemoinformatic Analysis Identifies Promiscuous Ligands Binding to Proteins from Different Classes with Varying Shapes

(1) Background: Compounds with multitarget activity are of interest in basic research to explore molecular foundations of promiscuous binding and in drug discovery as agents eliciting polypharmacological effects. Our study has aimed to systematically identify compounds that form complexes with proteins from distinct classes and compare their bioactive conformations and molecular properties. (2) Methods: A large-scale computational investigation was carried out that combined the analysis of complex X-ray structures, ligand binding modes, compound activity data, and various molecular properties. (3) Results: A total of 515 ligands with multitarget activity were identified that included 70 organic compounds binding to proteins from different classes. These multiclass ligands (MCLs) were often flexible and surprisingly hydrophilic. Moreover, they displayed a wide spectrum of binding modes. In different target structure environments, binding shapes of MCLs were often similar, but also distinct. (4) Conclusions: Combined structural and activity data analysis identified compounds with activity against proteins with distinct structures and functions. MCLs were found to have greatly varying shape similarity when binding to different protein classes. Hence, there were no apparent canonical binding shapes indicating multitarget activity. Rather, conformational versatility characterized MCL binding.

Unraveling Plant Natural Chemical Diversity for Drug Discovery Purposes

The screening and testing of extracts against a variety of pharmacological targets in order to benefit from the immense natural chemical diversity is a concern in many laboratories worldwide. And several successes have been recorded in finding new actives in natural products, some of which have become new drugs or new sources of inspiration for drugs. But in view of the vast amount of research on the subject, it is surprising that not more drug candidates were found. In our view, it is fundamental to reflect upon the approaches of such drug discovery programs and the technical processes that are used, along with their inherent difficulties and biases. Based on an extensive survey of recent publications, we discuss the origin and the variety of natural chemical diversity as well as the strategies to having the potential to embrace this diversity. It seemed to us that some of the difficulties of the area could be related with the technical approaches that are used, so the present review begins with synthetizing some of the more used discovery strategies, exemplifying some key points, in order to address some of their limitations. It appears that one of the challenges of natural product-based drug discovery programs should be an easier access to renewable sources of plant-derived products. Maximizing the use of the data together with the exploration of chemical diversity while working on reasonable supply of natural product-based entities could be a way to answer this challenge. We suggested alternative ways to access and explore part of this chemical diversity with in vitro cultures. We also reinforced how important it was organizing and making available this worldwide knowledge in an "inventory" of natural products and their sources. And finally, we focused on strategies based on synthetic biology and syntheses that allow reaching industrial scale supply. Approaches based on the opportunities lying in untapped natural plant chemical diversity are also considered.

Dose-Response Metabolomics To Understand Biochemical Mechanisms and Off-Target Drug Effects with the TOXcms Software

Small-molecule drugs and toxicants commonly interact with more than a single protein target, each of which may have unique effects on cellular phenotype. Although untargeted metabolomics is often applied to understand the mode of action of these chemicals, simple pairwise comparisons of treated and untreated samples are insufficient to resolve the effects of disrupting two or more independent protein targets. Here, we introduce a workflow for dose-response metabolomics to evaluate chemicals that potentially affect multiple proteins with different potencies. Our approach relies on treating samples with various concentrations of compound prior to analysis with mass spectrometry-based metabolomics. Data are then processed with software we developed called TOXcms, which statistically evaluates dose-response trends for each metabolomic signal according to user-defined tolerances and subsequently groups those that follow the same pattern. Although TOXcms was built upon the XCMS framework, it is compatible with any metabolomic data-processing software. Additionally, to enable correlation of dose responses beyond those that can be measured by metabolomics, TOXcms also accepts data from respirometry, cell death assays, other omic platforms, etc. In this work, we primarily focus on applying dose-response metabolomics to find off-target effects of drugs. Using metformin and etomoxir as examples, we demonstrate that each group of dose-response patterns identified by TOXcms signifies a metabolic response to a different protein target with a unique drug binding affinity. TOXcms is freely available on our laboratory website at http://pattilab.wustl.edu/software/toxcms .

Identification of Target Associations for Polypharmacology from Analysis of Crystallographic Ligands of the Protein Data Bank

The design of a chemical entity that potently and selectively binds to a biological target of therapeutic relevance has dominated the scene of drug discovery so far. However, recent findings suggest that multitarget ligands may be endowed with superior efficacy and be less prone to drug resistance. The Protein Data Bank (PDB) provides experimentally validated structural information about targets and bound ligands. Therefore, it represents a valuable source of information to help identifying active sites, understanding pharmacophore requirements, designing novel ligands, and inferring structure-activity relationships. In this study, we performed a large-scale analysis of the PDB by integrating different ligand-based and structure-based approaches, with the aim of identifying promising target associations for polypharmacology based on reported crystal structure information. First, the 2D and 3D similarity profiles of the crystallographic ligands were evaluated using different ligand-based methods. Then, activity data of pairs of similar ligands binding to different targets were inspected by comparing structural information with bioactivity annotations reported in the ChEMBL, BindingDB, BindingMOAD, and PDBbind databases. Afterward, extensive docking screenings of ligands in the identified cross-targets were made in order to validate and refine the ligand-based results. Finally, the therapeutic relevance of the identified target combinations for polypharmacology was evaluated from comparison with information on therapeutic targets reported in the Therapeutic Target Database (TTD). The results led to the identification of several target associations with high therapeutic potential for polypharmacology.

Promiscuity of in Vitro Secondary Pharmacology Assays and Implications for Lead Optimization Strategies

We conducted an analysis on screening data generated from 1445 compounds against a panel of 130 enzymes, ion channels, and receptors to assess secondary pharmacological risks. Hit rates of these targets as well as physicochemical properties for those hits were evaluated. A majority of targets yielded hits with higher clogP, molecular weight, and more basic character than inactive compounds. Although most targets favored lipophilic hits, the average clogP of hits at a given target did not correlate with its hit rate. Furthermore, a matched pair analysis was completed to determine structural changes that impacted off-target activities. A correlation of binding assays used in this analysis illustrated that some pharmacologically related binding assays are highly correlative and may be substituted for a smaller set of surrogate assays.

Molecular Inhibitor of QSOX1 Suppresses Tumor Growth In Vivo

Quiescin sulfhydryl oxidase 1 (QSOX1) is an enzyme overexpressed by many different tumor types. QSOX1 catalyzes the formation of disulfide bonds in proteins. Because short hairpin knockdowns (KD) of QSOX1 have been shown to suppress tumor growth and invasion in vitro and in vivo, we hypothesized that chemical compounds inhibiting QSOX1 enzymatic activity would also suppress tumor growth, invasion, and metastasis. High throughput screening using a QSOX1-based enzymatic assay revealed multiple potential QSOX1 inhibitors. One of the inhibitors, known as "SBI-183," suppresses tumor cell growth in a Matrigel-based spheroid assay and inhibits invasion in a modified Boyden chamber, but does not affect viability of nonmalignant cells. Oral administration of SBI-183 inhibits tumor growth in 2 independent human xenograft mouse models of renal cell carcinoma. We conclude that SBI-183 warrants further exploration as a useful tool for understanding QSOX1 biology and as a potential novel anticancer agent in tumors that overexpress QSOX1.

The human endogenous metabolome as a pharmacology baseline for drug discovery

We have limited understanding of the variation in in vitro affinities of drugs for their targets. An analysis of a highly curated set of 815 interactions between 566 drugs and 129 primary targets reveals that 71% of drug-target affinities have values above that of the corresponding endogenous ligand, 96% of them fitting within a range of two orders of magnitude. Our findings suggest that the evolutionary optimised affinity of endogenous ligands for their native proteins can serve as a baseline for the primary pharmacology of drugs. We show that the degree of off-target selectivity and safety risks of drugs derived from their secondary pharmacology depend very much on that baseline. Thus, we propose a new approach for estimating safety margins.

SAR Matrix Method for Large-Scale Analysis of Compound Structure-Activity Relationships and Exploration of Multitarget Activity Spaces

As the number of compounds and the volume of bioactivity data rapidly grow, advanced computational methods are required to study structure-activity relationships (SARs) on a large scale. Herein, the SAR matrix (SARM) methodology is described that was designed to systematically extract structural relationships between bioactive compounds from large databases, explore structure-activity relationships, and navigate multitarget activity spaces, which is one of the core tasks in chemogenomics. In addition, the SARM approach was designed to visualize structural and structure-activity relationships, which is often of critical importance for making this information available in an intuitive form for practical applications.

Inhibition of Oncogenic Kinases: An In Vitro Validated Computational Approach Identified Potential Multi-Target Anticancer Compounds

Tumorigenesis in humans is a multistep progression that imitates genetic changes leading to cell transformation and malignancy. Oncogenic kinases play a central role in cancer progression, rendering them putative targets for the design of anti-cancer drugs. The presented work aims to identify the potential multi-target inhibitors of oncogenic receptor tyrosine kinases (RTKs) and serine/threonine kinases (STKs). For this, chemoinformatics and structure-based virtual screening approaches were combined with an in vitro validation of lead hits on both cancerous and non-cancerous cell lines. A total of 16 different kinase structures were screened against ~739,000 prefiltered compounds using diversity selection, after which the top hits were filtered for promising pharmacokinetic properties. This led to the identification of 12 and 9 compounds against RTKs and STKs, respectively. Molecular dynamics (MD) simulations were carried out to better comprehend the stability of the predicted hit kinase-compound complexes. Two top-ranked compounds against each kinase class were tested in vitro for cytotoxicity, with compound F34 showing the most promising inhibitory activity in HeLa, HepG2, and Vero cell lines with IC₅₀ values of 145.46 μM, 175.48 μM, and 130.52 μM, respectively. Additional docking of F34 against various RTKs was carried out to support potential multi-target inhibition. Together with reliable MD simulations, these results suggest the promising potential of identified multi-target STK and RTK scaffolds for further kinase-specific anti-cancer drug development toward combinatorial therapies.

In Silico Target Prediction for Small Molecules

Drugs modulate disease states through their actions on targets in the body. Determining these targets aids the focused development of new treatments, and helps to better characterize those already employed. One means of accomplishing this is through the deployment of in silico methodologies, harnessing computational analytical and predictive power to produce educated hypotheses for experimental verification. Here, we provide an overview of the current state of the art, describe some of the well-established methods in detail, and reflect on how they, and emerging technologies promoting the incorporation of complex and heterogeneous data-sets, can be employed to improve our understanding of (poly)pharmacology.

Synthesis, cytotoxic evaluation, and molecular docking studies of the semi-synthetic "triterpenoid-steroid" hybrids

Synthetic transformations of steroids for drug discovery and improvement of drug effectiveness have been an important part of modern medicinal chemistry and pharmaceutical sciences. Pentacyclic triterpenoids, being represented in the nature by various structures and biogenetically related to steroids, can largely expand the spectrum of biologically active steroidal agents via synthesis of the so-called "triterpenoid-steroid" hybrids. In the presented work, the nitrile anion cyclizations of 3,4-secolupane and 3,4-seco-oleanane nitriles and follow-up synthetic transformations of the cyclized products with formation of the gemm-dimethyl-free A ring "triterpenoid-steroid" hybrids were studied. Furthermore, the resulting cyclic compounds were modified at C3, C4, and/or C5 positions of ring A, as well as at C20, C28, and C30 positions of the isopropylidene moiety in the case of lupane triterpenoids. The cytotoxic effect of the synthesized compounds against seven cancer cell lines HEp-2, HCT 116, MS, RD TE32, A549, MCF7, and PC3 was evaluated. The in silico identification of potential anticancer protein targets with regard to the compounds, which were active at micromolar concentrations against tested cell lines, was carried out. The molecular docking studies showed that compound 19, which demonstrated most pronounced cytotoxicity (IC₅₀ 0.64-3.17 μM) against all tested cell lines, fits well the active sites of CDK6 and HER2/neu.

Dual Inhibitors of PARPs and ROCKs

Recent network and system biology analyses suggest that most complex diseases are regulated by robust and highly interconnected pathways that could be better modulated by small molecules binding to multiple biological targets. These pieces of evidence recently led to devote efforts on identifying single chemical entities that bind to two different disease-relevant targets. Here, we first predicted in silico and later confirmed in vitro that UPF 1069, a known bioactive poly(ADP-ribose) polymerase-1/2 (PARP1/2) molecule, and hydroxyfasudil, a known bioactive Rho-associated protein kinase-1/2 (ROCK1/2) molecule, have low-micromolar cross-affinity for ROCK1/2 and PARP1/2, respectively. These molecules can now be regarded as chemical seeds from which pharmacological tools could be generated to study the impact of dual inhibition of PARPs and ROCKs in preclinical models of a variety of complex diseases where both targets are involved.

Drugs as habitable planets in the space of dark chemical matter

A recent study demonstrated antifungal activity of dark chemical matter (DCM) compounds that were otherwise inactive in more than 100 HTS assays. These compounds were proposed to possess unique activity and 'clean' safety profiles. Here, we present an outlook of the promiscuity and safety of these compounds by retrospectively comparing their chemical and biological spaces with those of drugs. Significant amounts of marketed drugs (16%), withdrawn drugs (16.5%) and natural compounds (3.5%) share structural identity with DCM. Compound promiscuity assessment indicates that dark matter compounds could potentially interact with multiple biological targets. Further, thousands of DCM compounds showed presence of frequent-hitting pan-assay interference compound (PAINS) substructures. In light of these observations, filtering these compounds from screening libraries can be an irrevocable loss.

Network-based Approaches in Pharmacology

In drug discovery, network-based approaches are expected to spotlight our understanding of drug action across multiple layers of information. On one hand, network pharmacology considers the drug response in the context of a cellular or phenotypic network. On the other hand, a chemical-based network is a promising alternative for characterizing the chemical space. Both can provide complementary support for the development of rational drug design and better knowledge of the mechanisms underlying the multiple actions of drugs. Recent progress in both concepts is discussed here. In addition, a network-based approach using drug-target-therapy data is introduced as an example.

Entering the 'big data' era in medicinal chemistry: molecular promiscuity analysis revisited

The 'big data' concept plays an increasingly important role in many scientific fields. Big data involves more than unprecedentedly large volumes of data that become available. Different criteria characterizing big data must be carefully considered in computational data mining, as we discuss herein focusing on medicinal chemistry. This is a scientific discipline where big data is beginning to emerge and provide new opportunities. For example, the ability of many drugs to specifically interact with multiple targets, termed promiscuity, forms the molecular basis of polypharmacology, a hot topic in drug discovery. Compound promiscuity analysis is an area that is much influenced by big data phenomena. Different results are obtained depending on chosen data selection and confidence criteria, as we also demonstrate.

The scoring bias in reverse docking and the score normalization strategy to improve success rate of target fishing

Target fishing often relies on the use of reverse docking to identify potential target proteins of ligands from protein database. The limitation of reverse docking is the accuracy of current scoring funtions used to distinguish true target from non-target proteins. Many contemporary scoring functions are designed for the virtual screening of small molecules without special optimization for reverse docking, which would be easily influenced by the properties of protein pockets, resulting in scoring bias to the proteins with certain properties. This bias would cause lots of false positives in reverse docking, interferring the identification of true targets. In this paper, we have conducted a large-scale reverse docking (5000 molecules to 100 proteins) to study the scoring bias in reverse docking by DOCK, Glide, and AutoDock Vina. And we found that there were actually some frequency hits, namely interference proteins in all three docking procedures. After analyzing the differences of pocket properties between these interference proteins and the others, we speculated that the interference proteins have larger contact area (related to the size and shape of protein pockets) with ligands (for all three docking programs) or higher hydrophobicity (for Glide), which could be the causes of scoring bias. Then we applied the score normalization method to eliminate this scoring bias, which was effective to make docking score more balanced between different proteins in the reverse docking of benchmark dataset. Later, the Astex Diver Set was utilized to validate the effect of score normalization on actual cases of reverse docking, showing that the accuracy of target prediction significantly increased by 21.5% in the reverse docking by Glide after score normalization, though there was no obvious change in the reverse docking by DOCK and AutoDock Vina. Our results demonstrate the effectiveness of score normalization to eliminate the scoring bias and improve the accuracy of target prediction in reverse docking. Moreover, the properties of protein pockets causing scoring bias to certain proteins we found here can provide the theory basis to further optimize the scoring functions of docking programs for future research.

Large-scale bioactivity analysis of the small-molecule assayed proteome

This study presents an analysis of the small molecule bioactivity profiles across large quantities of diverse protein families represented in PubChem BioAssay. We compared the bioactivity profiles of FDA approved drugs to non-FDA approved compounds, and report several distinct patterns characteristic of the approved drugs. We found that a large fraction of the previously reported higher target promiscuity among FDA approved compounds, compared to non-FDA approved bioactives, was frequently due to cross-reactivity within rather than across protein families. We identified 804 potentially novel protein target candidates for FDA approved drugs, as well as 901 potentially novel target candidates with active non-FDA approved compounds, but no FDA approved drugs with activity against these targets. We also identified 486348 potentially novel compounds active against the same targets as FDA approved drugs, as well as 153402 potentially novel compounds active against targets without active FDA approved drugs. By quantifying the agreement among replicated screens, we estimated that more than half of these novel outcomes are reproducible. Using biclustering, we identified many dense clusters of FDA approved drugs with enriched activity against a common set of protein targets. We also report the distribution of compound promiscuity using a Bayesian statistical model, and report the sensitivity and specificity of two common methods for identifying promiscuous compounds. Aggregator assays exhibited greater accuracy in identifying highly promiscuous compounds, while PAINS substructures were able to identify a much larger set of "middle range" promiscuous compounds. Additionally, we report a large number of promiscuous compounds not identified as aggregators or PAINS. In summary, the results of this study represent a rich reference for selecting novel drug and target protein candidates, as well as for eliminating candidate compounds with unselective activities.

Discovery of a Chemical Probe Bisamide (CCT251236): An Orally Bioavailable Efficacious Pirin Ligand from a Heat Shock Transcription Factor 1 (HSF1) Phenotypic Screen

Phenotypic screens, which focus on measuring and quantifying discrete cellular changes rather than affinity for individual recombinant proteins, have recently attracted renewed interest as an efficient strategy for drug discovery. In this article, we describe the discovery of a new chemical probe, bisamide (CCT251236), identified using an unbiased phenotypic screen to detect inhibitors of the HSF1 stress pathway. The chemical probe is orally bioavailable and displays efficacy in a human ovarian carcinoma xenograft model. By developing cell-based SAR and using chemical proteomics, we identified pirin as a high affinity molecular target, which was confirmed by SPR and crystallography.

Highly Promiscuous Small Molecules from Biological Screening Assays Include Many Pan-Assay Interference Compounds but Also Candidates for Polypharmacology

In PubChem screening assays, 466 highly promiscuous compounds were identified that were examined for known pan-assay interference compounds (PAINS) and aggregators using publicly available filters. These filters detected 210 PAINS and 67 aggregators. Compounds passing the filters included additional PAINS that were not detected, mostly due to tautomerism, and a variety of other potentially reactive compounds currently not encoded as PAINS. For a subset of compounds passing the filters, there was no evidence of potential artifacts. These compounds are considered candidates for further exploring multitarget activities and the molecular basis of polypharmacology.

Comparing atom-based with residue-based descriptors in predicting binding site similarity: do backbone atoms matter?

AIM

We question the level of detail required in protein 3D-representation to detect site similarity which is relevant for polypharmacology prediction.

RESULTS

We modified the in-house program SiteAlign to replace generic pharmacophoric descriptors of cavity-lining amino acids by descriptors accounting for solvent exposure. Benchmarking the novel, atom-based, method (SiteAlign2) revealed no global improvement of performance. However, in the rare cases of no sequence or global structure similarities between the compared proteins, SiteAlign2 was more successful if backbone atoms are key determinants of ligand binding.

CONCLUSION

SiteAlign suits the comparison of binding sites for close or distant homologs. SiteAlign2 provides a better insight into the physical model of site similarity between nonhomologs, but at the expense of an increased sensitivity to atomic coordinates.

Computational Method for the Systematic Identification of Analog Series and Key Compounds Representing Series and Their Biological Activity Profiles

A computational methodology is introduced for detecting all unique series of analogs in large compound data sets, regardless of chemical relationships between analogs. No prior knowledge of core structures or R-groups is required, which are automatically determined. The approach is based upon the generation of retrosynthetic matched molecular pairs and analog networks from which distinct series are isolated. The methodology was applied to systematically extract more than 17 000 distinct series from the ChEMBL database. For comparison, analog series were also isolated from screening compounds and drugs. Known biological activities were mapped to series from ChEMBL, and in more than 13 000 of these series, key compounds were identified that represented substitution sites of all analogs within a series and its complete activity profile. The analog series, key compounds, and activity profiles are made freely available as a resource for medicinal chemistry applications.

Phenotypic Assessment and the Discovery of Topiramate

The role of phenotypic assessment in drug discovery is discussed, along with the discovery and development of TOPAMAX (topiramate), a billion-dollar molecule for the treatment of epilepsy and migraine.

Geomfinder: a multi-feature identifier of similar three-dimensional protein patterns: a ligand-independent approach

Background

Since the structure of proteins is more conserved than the sequence, the identification of conserved three-dimensional (3D) patterns among a set of proteins, can be important for protein function prediction, protein clustering, drug discovery and the establishment of evolutionary relationships. Thus, several computational applications to identify, describe and compare 3D patterns (or motifs) have been developed. Often, these tools consider a 3D pattern as that described by the residues surrounding co-crystallized/docked ligands available from X-ray crystal structures or homology models. Nevertheless, many of the protein structures stored in public databases do not provide information about the location and characteristics of ligand binding sites and/or other important 3D patterns such as allosteric sites, enzyme-cofactor interaction motifs, etc. This makes necessary the development of new ligand-independent methods to search and compare 3D patterns in all available protein structures.

Results

Here we introduce Geomfinder, an intuitive, flexible, alignment-free and ligand-independent web server for detailed estimation of similarities between all pairs of 3D patterns detected in any two given protein structures. We used around 1100 protein structures to form pairs of proteins which were assessed with Geomfinder. In these analyses each protein was considered in only one pair (e.g. in a subset of 100 different proteins, 50 pairs of proteins can be defined). Thus: (a) Geomfinder detected identical pairs of 3D patterns in a series of monoamine oxidase-B structures, which corresponded to the effectively similar ligand binding sites at these proteins; (b) we identified structural similarities among pairs of protein structures which are targets of compounds such as acarbose, benzamidine, adenosine triphosphate and pyridoxal phosphate; these similar 3D patterns are not detected using sequence-based methods; (c) the detailed evaluation of three specific cases showed the versatility of Geomfinder, which was able to discriminate between similar and different 3D patterns related to binding sites of common substrates in a range of diverse proteins.

Conclusions

Geomfinder allows detecting similar 3D patterns between any two pair of protein structures, regardless of the divergency among their amino acids sequences. Although the software is not intended for simultaneous multiple comparisons in a large number of proteins, it can be particularly useful in cases such as the structure-based design of multitarget drugs, where a detailed analysis of 3D patterns similarities between a few selected protein targets is essential.

Electronic supplementary material

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

Determining the Degree of Promiscuity of Extensively Assayed Compounds

In the context of polypharmacology, an emerging concept in drug discovery, promiscuity is rationalized as the ability of compounds to specifically interact with multiple targets. Promiscuity of drugs and bioactive compounds has thus far been analyzed computationally on the basis of activity annotations, without taking assay frequencies or inactivity records into account. Most recent estimates have indicated that bioactive compounds interact on average with only one to two targets, whereas drugs interact with six or more. In this study, we have further extended promiscuity analysis by identifying the most extensively assayed public domain compounds and systematically determining their promiscuity. These compounds were tested in hundreds of assays against hundreds of targets. In our analysis, assay promiscuity was distinguished from target promiscuity and separately analyzed for primary and confirmatory assays. Differences between the degree of assay and target promiscuity were surprisingly small and average and median degrees of target promiscuity of 2.6 to 3.4 and 2.0 were determined, respectively. Thus, target promiscuity remained at a low level even for most extensively tested active compounds. These findings provide further evidence that bioactive compounds are less promiscuous than drugs and have implications for pharmaceutical research. In addition to a possible explanation that drugs are more extensively tested for additional targets, the results would also support a “promiscuity enrichment model” according to which promiscuous compounds might be preferentially selected for therapeutic efficacy during clinical evaluation to ultimately become drugs.

Polypharmacology: in silico methods of ligand design and development

How to design a ligand to bind multiple targets, rather than to a single target, is the focus of this review. Rational polypharmacology draws on knowledge that is both broad ranging and hierarchical. Computer-aided multitarget ligand design methods are described according to their nested knowledge level. Ligand-only and then receptor-ligand strategies are first described; followed by the metabolic network viewpoint. Subsequently strategies that view infectious diseases as multigenomic targets are discussed, and finally the disease level interpretation of medicinal therapy is considered. As yet there is no consensus on how best to proceed in designing a multitarget ligand. The current methodologies are bought together in an attempt to give a practical overview of how polypharmacology design might be best initiated.

The In Vitro Pharmacological Profile of Drugs as a Proxy Indicator of Potential In Vivo Organ Toxicities

The potential of a drug to cause certain organ toxicities is somehow implicitly contained in its full pharmacological profile, provided the drug reaches and accumulates at the various organs where the different interacting proteins in its profile, both targets and off-targets, are expressed. Under this assumption, a computational approach was implemented to obtain a projected anatomical profile of a drug from its in vitro pharmacological profile linked to protein expression data across 47 organs. It was observed that the anatomical profiles obtained when using only the known primary targets of the drugs reflected roughly the intended organ targets. However, when both known and predicted secondary pharmacology was considered, the projected anatomical profiles of the drugs were able to clearly highlight potential organ off-targets. Accordingly, when applied to sets of drugs known to cause cardiotoxicity and hepatotoxicity, the approach is able to identify heart and liver, respectively, as the organs where the proteins in the pharmacological profile of the corresponding drugs are specifically expressed. When applied to a set of drugs linked to a risk of Torsades de Pointes, heart is again the organ clearly standing out from the rest and a potential protein profile hazard is proposed. The approach can be used as a proxy indicator of potential in vivo organ toxicities.

Assessing the Growth of Bioactive Compounds and Scaffolds over Time: Implications for Lead Discovery and Scaffold Hopping

The increase in compounds with activity against five major therapeutic target families has been quantified on a time scale and investigated employing a compound-scaffold-cyclic skeleton (CSK) hierarchy. The analysis was designed to better understand possible reasons for target-dependent growth of bioactive compounds. There was strong correlation between compound and scaffold growth across all target families. Active compounds becoming available over time were mostly represented by new scaffolds. On the basis of scaffold-to-compound ratios, new active compounds were structurally diverse and, on the basis of CSK-to-scaffold ratios, often had previously unobserved topologies. In addition, novel targets emerged that complemented major families. The analysis revealed that compound growth is associated with increasing chemical diversity and that current pharmaceutical targets are capable of recognizing many structurally different compounds, which provides a rationale for the rapid increase in the number of bioactive compounds over the past decade. In light of these findings, it is likely that new chemical entities will be discovered for many small molecule targets including relatively unexplored ones as well as for popular and well-studied therapeutic targets. Moreover, given the wealth of new "active scaffolds" that have been increasingly identified for many targets over time, computational scaffold-hopping exercises should generally have a high likelihood of success.