What Can We Learn from Bioactivity Data? Chemoinformatics Tools and Applications in Chemical Biology Research

QSARtuna: An Automated QSAR Modeling Platform for Molecular Property Prediction in Drug Design

Machine-learning (ML) and deep-learning (DL) approaches to predict the molecular properties of small molecules are increasingly deployed within the design-make-test-analyze (DMTA) drug design cycle to predict molecular properties of interest. Despite this uptake, there are only a few automated packages to aid their development and deployment that also support uncertainty estimation, model explainability, and other key aspects of model usage. This represents a key unmet need within the field, and the large number of molecular representations and algorithms (and associated parameters) means it is nontrivial to robustly optimize, evaluate, reproduce, and deploy models. Here, we present QSARtuna, a molecule property prediction modeling pipeline, written in Python and utilizing the Optuna, Scikit-learn, RDKit, and ChemProp packages, which enables the efficient and automated comparison between molecular representations and machine learning models. The platform was developed by considering the increasingly important aspect of model uncertainty quantification and explainability by design. We provide details for our framework and provide illustrative examples to demonstrate the capability of the software when applied to simple molecular property, reaction/reactivity prediction, and DNA encoded library enrichment classification. We hope that the release of QSARtuna will further spur innovation in automatic ML modeling and provide a platform for education of best practices in molecular property modeling. The code for the QSARtuna framework is made freely available via GitHub.

Acute behavioral and Neurochemical Effects of Novel N-Benzyl-2-Phenylethylamine Derivatives in Adult Zebrafish

Hallucinogenic drugs potently affect brain and behavior and have also recently emerged as potentially promising agents in pharmacotherapy. Complementing laboratory rodents, the zebrafish (Danio rerio) is a powerful animal model organism for screening neuroactive drugs, including hallucinogens. Here, we test a battery of ten novel N-benzyl-2-phenylethylamine (NBPEA) derivatives with the 2,4- and 3,4-dimethoxy substitutions in the phenethylamine moiety and the -OCH₃, -OCF₃, -F, -Cl, and -Br substitutions in the ortho position of the phenyl ring of the N-benzyl moiety, assessing their acute behavioral and neurochemical effects in the adult zebrafish. Overall, substitutions in the Overall, substitutions in the N-benzyl moiety modulate locomotion, and substitutions in the phenethylamine moiety alter zebrafish anxiety-like behavior, also affecting the brain serotonin and/or dopamine turnover. The 24H-NBOMe(F) and 34H-NBOMe(F) treatment also reduced zebrafish despair-like behavior. Computational analyses of zebrafish behavioral data by artificial intelligence identified several distinct clusters for these agents, including anxiogenic/hypolocomotor (24H-NBF, 24H-NBOMe, and 34H-NBF), behaviorally inert (34H-NBBr, 34H-NBCl, and 34H-NBOMe), anxiogenic/hallucinogenic-like (24H-NBBr, 24H-NBCl, and 24H-NBOMe(F)), and anxiolytic/hallucinogenic-like (34H-NBOMe(F)) drugs. Our computational analyses also revealed phenotypic similarity of the behavioral activity of some NBPEAs to that of selected conventional serotonergic and antiglutamatergic hallucinogens. In silico functional molecular activity modeling further supported the overlap of the drug targets for NBPEAs tested here and the conventional serotonergic and antiglutamatergic hallucinogens. Overall, these findings suggest potent neuroactive properties of several novel synthetic NBPEAs, detected in a sensitive in vivo vertebrate model system, the zebrafish, raising the possibility of their potential clinical use and abuse.

Discovery of an Unexpected Similarity in Ligand Binding between BRD4 and PPARγ

Knowledge about interrelationships between different proteins is crucial in fundamental research for the elucidation of protein networks and pathways. Furthermore, it is especially critical in chemical biology to identify further key regulators of a disease and to take advantage of polypharmacology effects. Here, we present a new concept that combines a scaffold-based analysis of bioactivity data with a subsequent screening to identify novel inhibitors for a protein target of interest. The initial scaffold-based analysis revealed a flavone-like scaffold that can be found in ligands of different unrelated proteins indicating a similarity in ligand binding. This similarity was further investigated by testing compounds on bromodomain-containing protein 4 (BRD4) that were similar to known ligands of the other identified protein targets. Several new BRD4 inhibitors were identified and proven to be validated hits based on orthogonal assays and X-ray crystallography. The most important discovery was an unexpected relationship between BRD4 and peroxisome-proliferator activated receptor gamma (PPARγ). Both proteins share binding site similarities near a common hydrophobic subpocket which should allow the design of a polypharmacology-based ligand targeting both proteins. Such dual-BRD4-PPARγ modulators open up new therapeutic opportunities, because both are important drug targets for cancer therapy and many more important diseases. Thereon, a complex structure of sulfasalazine was obtained that involves two bromodomains and could be a potential starting point for the design of a bivalent BRD4 inhibitor.

Schistosomiasis Drug Discovery in the Era of Automation and Artificial Intelligence

Schistosomiasis is a parasitic disease caused by trematode worms of the genus Schistosoma and affects over 200 million people worldwide. The control and treatment of this neglected tropical disease is based on a single drug, praziquantel, which raises concerns about the development of drug resistance. This, and the lack of efficacy of praziquantel against juvenile worms, highlights the urgency for new antischistosomal therapies. In this review we focus on innovative approaches to the identification of antischistosomal drug candidates, including the use of automated assays, fragment-based screening, computer-aided and artificial intelligence-based computational methods. We highlight the current developments that may contribute to optimizing research outputs and lead to more effective drugs for this highly prevalent disease, in a more cost-effective drug discovery endeavor.

C-SPADE: a web-tool for interactive analysis and visualization of drug screening experiments through compound-specific bioactivity dendrograms

The advent of polypharmacology paradigm in drug discovery calls for novel chemoinformatic tools for analyzing compounds’ multi-targeting activities. Such tools should provide an intuitive representation of the chemical space through capturing and visualizing underlying patterns of compound similarities linked to their polypharmacological effects. Most of the existing compound-centric chemoinformatics tools lack interactive options and user interfaces that are critical for the real-time needs of chemical biologists carrying out compound screening experiments. Toward that end, we introduce C-SPADE, an open-source exploratory web-tool for interactive analysis and visualization of drug profiling assays (biochemical, cell-based or cell-free) using compound-centric similarity clustering. C-SPADE allows the users to visually map the chemical diversity of a screening panel, explore investigational compounds in terms of their similarity to the screening panel, perform polypharmacological analyses and guide drug-target interaction predictions. C-SPADE requires only the raw drug profiling data as input, and it automatically retrieves the structural information and constructs the compound clusters in real-time, thereby reducing the time required for manual analysis in drug development or repurposing applications. The web-tool provides a customizable visual workspace that can either be downloaded as figure or Newick tree file or shared as a hyperlink with other users. C-SPADE is freely available at http://cspade.fimm.fi/.

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.

Exploring G Protein-Coupled Receptors (GPCRs) Ligand Space via Cheminformatics Approaches: Impact on Rational Drug Design

The primary goal of rational drug discovery is the identification of selective ligands which act on single or multiple drug targets to achieve the desired clinical outcome through the exploration of total chemical space. To identify such desired compounds, computational approaches are necessary in predicting their drug-like properties. G Protein-Coupled Receptors (GPCRs) represent one of the largest and most important integral membrane protein families. These receptors serve as increasingly attractive drug targets due to their relevance in the treatment of various diseases, such as inflammatory disorders, metabolic imbalances, cardiac disorders, cancer, monogenic disorders, etc. In the last decade, multitudes of three-dimensional (3D) structures were solved for diverse GPCRs, thus referring to this period as the "golden age for GPCR structural biology." Moreover, accumulation of data about the chemical properties of GPCR ligands has garnered much interest toward the exploration of GPCR chemical space. Due to the steady increase in the structural, ligand, and functional data of GPCRs, several cheminformatics approaches have been implemented in its drug discovery pipeline. In this review, we mainly focus on the cheminformatics-based paradigms in GPCR drug discovery. We provide a comprehensive view on the ligand- and structure-based cheminformatics approaches which are best illustrated via GPCR case studies. Furthermore, an appropriate combination of ligand-based knowledge with structure-based ones, i.e., integrated approach, which is emerging as a promising strategy for cheminformatics-based GPCR drug design is also discussed.

Spectrophores as one-dimensional descriptors calculated from three-dimensional atomic properties: applications ranging from scaffold hopping to multi-target virtual screening

Spectrophores are novel descriptors that are calculated from the three-dimensional atomic properties of molecules. In our current implementation, the atomic properties that were used to calculate spectrophores include atomic partial charges, atomic lipophilicity indices, atomic shape deviations and atomic softness properties. This approach can easily be widened to also include additional atomic properties. Our novel methodology finds its roots in the experimental affinity fingerprinting technology developed in the 1990’s by Terrapin Technologies. Here we have translated it into a purely virtual approach using artificial affinity cages and a simplified metric to calculate the interaction between these cages and the atomic properties. A typical spectrophore consists of a vector of 48 real numbers. This makes it highly suitable for the calculation of a wide range of similarity measures for use in virtual screening and for the investigation of quantitative structure–activity relationships in combination with advanced statistical approaches such as self-organizing maps, support vector machines and neural networks. In our present report we demonstrate the applicability of our novel methodology for scaffold hopping as well as virtual screening.

CHIPMUNK: A Virtual Synthesizable Small-Molecule Library for Medicinal Chemistry, Exploitable for Protein-Protein Interaction Modulators

A common issue during drug design and development is the discovery of novel scaffolds for protein targets. On the one hand the chemical space of purchasable compounds is rather limited; on the other hand artificially generated molecules suffer from a grave lack of accessibility in practice. Therefore, we generated a novel virtual library of small molecules which are synthesizable from purchasable educts, called CHIPMUNK (CHemically feasible In silico Public Molecular UNiverse Knowledge base). Altogether, CHIPMUNK covers over 95 million compounds and encompasses regions of the chemical space that are not covered by existing databases. The coverage of CHIPMUNK exceeds the chemical space spanned by the Lipinski rule of five to foster the exploration of novel and difficult target classes. The analysis of the generated property space reveals that CHIPMUNK is well suited for the design of protein-protein interaction inhibitors (PPIIs). Furthermore, a recently developed structural clustering algorithm (StruClus) for big data was used to partition the sub-libraries into meaningful subsets and assist scientists to process the large amount of data. These clustered subsets also contain the target space based on ChEMBL data which was included during clustering.

Polypharmacology of conformationally locked methanocarba nucleosides

A single molecular scaffold can be adapted to interact with diverse targets, either separately or simultaneously. Nucleosides and nucleotides in which ribose is substituted with bicyclo[3.1.0]hexane are an example of a versatile drug-like scaffold for increasing selectivity at their classical targets: kinases, polymerases, adenosine and P2 receptors. Also, by applying structure-based functional group manipulations, rigidified adenosine derivatives can be repurposed to satisfy pharmacophoric requirements of various GPCRs, ion channels, enzymes and transporters, initially detected as off-target activities. Recent examples include 5HT2B serotonin receptor antagonists and novel dopamine transporter allosteric modulators. This directable target diversity establishes rigid nucleosides as privileged scaffolds.

MouseTox: An online toxicity assessment tool for small molecules through Enalos Cloud platform

Advances in the drug discovery research substantially depend on in silico methods and techniques that capitalize on experimental data to enable the accurate property/activity assessment by employing a variety of computational techniques. These in silico tools can significantly reduce expensive and time consuming experimental procedures required and are strongly recommended to avoid animal testing, especially as far as toxicity evaluation and risk assessment is concerned. In this context, in the present work we aim to develop a predictive model for the cytotoxic effects of a wide range of compounds based solely on calculated molecular descriptors that account for their topological, geometric and structural characteristics. The developed model was fully validated and was released online via Enalos Cloud platform accessible through http://enalos.insilicotox.com/MouseTox/. This ready-to-use web service offers, through a user-friendly interface, free access to the model results and therefore can act as a toxicity prediction tool for the risk assessment of novel compounds, without any special requirements or prior programming skills.

Screening and identification of novel biologically active natural compounds

With the advent of very rapid and cheap genome analyses and the linkage of these plus microbial metabolomics to potential compound structures came the realization that there was an immense sea of novel agents to be mined and tested. In addition, it is now recognized that there is significant microbial involvement in many natural products isolated from “nominally non-microbial sources”.

This short review covers the current screening methods that have evolved and one might even be tempted to say “devolved” in light of the realization that target-based screens had problems when the products entered clinical testing, with off-target effects being the major ones. Modern systems include, but are not limited to, screening in cell lines utilizing very modern techniques (a high content screen) that are designed to show interactions within cells when treated with an “agent”. The underlying principle(s) used in such systems dated back to unpublished attempts in the very early 1980s by the pharmaceutical industry to show toxic interactions within animal cells by using automated light microscopy. Though somewhat successful, the technology was not adequate for any significant commercialization. Somewhat later, mammalian cell lines that were “genetically modified” to alter signal transduction cascades, either up or down, and frequently linked to luciferase readouts, were then employed in a 96-well format. In the case of microbes, specific resistance parameters were induced in isogenic cell lines from approximately the mid-1970s. In the latter two cases, comparisons against parent and sibling cell lines were used in order that a rapid determination of potential natural product “hits” could be made. Obviously, all of these assay systems could also be, and were, used for synthetic molecules.

These methods and their results have led to a change in what the term “screening for bioactivity” means. In practice, versions of phenotypic screening are returning, but in a dramatically different scientific environment from the 1970s, as I hope to demonstrate in the short article that follows.

Scaffold Hunter: a comprehensive visual analytics framework for drug discovery

The era of big data is influencing the way how rational drug discovery and the development of bioactive molecules is performed and versatile tools are needed to assist in molecular design workflows. Scaffold Hunter is a flexible visual analytics framework for the analysis of chemical compound data and combines techniques from several fields such as data mining and information visualization. The framework allows analyzing high-dimensional chemical compound data in an interactive fashion, combining intuitive visualizations with automated analysis methods including versatile clustering methods. Originally designed to analyze the scaffold tree, Scaffold Hunter is continuously revised and extended. We describe recent extensions that significantly increase the applicability for a variety of tasks.