Scaffold diversity of exemplified medicinal chemistry space

An Open-Source Implementation of the Scaffold Identification and Naming System (SCINS) and Example Applications

Organizing and partitioning sets of chemical structures is of considerable practical significance, e.g., in compound library analysis and the postprocessing of screening hit lists. Approaches such as unsupervised clustering are computationally demanding and dataset-dependent; on the other hand, rule-based methods, such as those based on Murcko scaffolds, have linear time complexity but are often too fine-grained, leading to a large number of singletons or sparsely populated classes. An alternative rule-based method that seeks to achieve an optimal balance when grouping compounds into sets is the 'Scaffold Identification and Naming System' (SCINS). To facilitate public use of this previously published method, here, we provide an open-source Python implementation of SCINS, dependent only on RDKit. We show that SCINS can be useful in identifying sparsely and densely populated regions in chemical space in large databases, here exemplified with Enamine REAL Diverse and ChEMBL. We find that Enamine REAL Diverse covers a much smaller SCINS space relative to ChEMBL, whereas the opposite is true when Murcko and generic Murcko scaffolds are considered. Additionally, we show that SCINS can result in chemically intuitive grouping of medium-sized sets of bioactive compounds, which can be useful in compound selection from virtual screening campaigns as well as postprocessing of experimental hit lists. Hence, in this work, we provide both an open-source implementation of SCINS and its characterization with relevant use cases.

Chemoinformatic Characterization of NAPROC-13: A Database for Natural Product 13C NMR Dereplication

Natural products (NPs) are secondary metabolites of natural origin with broad applications across various human activities, particularly the discovery of bioactive compounds. Structural elucidation of new NPs entails significant cost and effort. On the other hand, the dereplication of known compounds is crucial for the early exclusion of irrelevant compounds in contemporary pharmaceutical research. NAPROC-13 stands out as a publicly accessible database, providing structural and 13C NMR spectroscopic information for over 25 000 compounds, rendering it a pivotal resource in natural product (NP) research, favoring open science. This study seeks to quantitatively analyze the chemical content, structural diversity, and chemical space coverage of NPs within NAPROC-13, compared to FDA-approved drugs and a very diverse subset of NPs, UNPD-A. Findings indicated that NPs in NAPROC-13 exhibit properties comparable to those in UNPD-A, albeit showcasing a notably diverse array of structural content, scaffolds, ring systems of pharmaceutical interest, and molecular fragments. NAPROC-13 covers a specific region of the chemical multiverse (a generalization of the chemical space from different chemical representations) regarding physicochemical properties and a region as broad as UNPD-A in terms of the structural features represented by fingerprints.

Prediction of Drug-Induced Liver Injury: From Molecular Physicochemical Properties and Scaffold Architectures to Machine Learning Approaches

The process of developing new drugs is widely acknowledged as being time-intensive and requiring substantial financial investment. Despite ongoing efforts to reduce time and expenses in drug development, ensuring medication safety remains an urgent problem. One of the major problems involved in drug development is hepatotoxicity, specifically known as drug-induced liver injury (DILI). The popularity of new drugs often poses a significant barrier during development and frequently leads to their recall after launch. In silico methods have many advantages compared with traditional in vivo and in vitro assays. To establish a more precise and reliable prediction model, it is necessary to utilize an extensive and high-quality database consisting of information on drug molecule properties and structural patterns. In addition, we should also carefully select appropriate molecular descriptors that can be used to accurately depict compound characteristics. The aim of this study was to conduct a comprehensive investigation into the prediction of DILI. First, we conducted a comparative analysis of the physicochemical properties of extensively well-prepared DILI-positive and DILI-negative compounds. Then, we used classic substructure dissection methods to identify structural pattern differences between these two different types of chemical molecules. These findings indicate that it is not feasible to establish property or substructure-based rules for distinguishing between DILI-positive and DILI-negative compounds. Finally, we developed quantitative classification models for predicting DILI using the naïve Bayes classifier (NBC) and recursive partitioning (RP) machine learning techniques. The optimal DILI prediction model was obtained using NBC, which combines 21 physicochemical properties, the VolSurf descriptors and the LCFP_10 fingerprint set. This model achieved a global accuracy (GA) of 0.855 and an area under the curve (AUC) of 0.704 for the training set, while the corresponding values were 0.619 and 0.674 for the test set, respectively. Moreover, indicative substructural fragments favorable or unfavorable for DILI were identified from the best naïve Bayesian classification model. These findings may help prioritize lead compounds in the early stage of drug development pipelines.

Occurrence of "Natural Selection" in Successful Small Molecule Drug Discovery

Published compounds from ChEMBL version 32 are used to seek evidence for the occurrence of "natural selection" in drug discovery. Three measures of natural product (NP) character were applied, to compare time- and target-matched compounds reaching the clinic (clinical compounds in phase 1-3 development and approved drugs) with background compounds (reference compounds). Pseudo-NPs (PNPs), containing NP fragments combined in ways inaccessible by nature, are increasing over time, reaching 67% of clinical compounds first disclosed since 2010. PNPs are 54% more likely to be found in post-2008 clinical versus reference compounds. The majority of target classes show increased clinical compound NP character versus their reference compounds. Only 176 NP fragments appear in >1000 clinical compounds published since 2008, yet these make up on average 63% of the clinical compound's core scaffolds. There is untapped potential awaiting exploitation, by applying nature's building blocks─"natural intelligence"─to drug design.

The Goldilocks paradigm: comparing classical machine learning, large language models, and few-shot learning for drug discovery applications

Recent advances in machine learning (ML) have led to newer model architectures including transformers (large language models, LLMs) showing state of the art results in text generation and image analysis as well as few-shot learning (FSLC) models which offer predictive power with extremely small datasets. These new architectures may offer promise, yet the 'no-free lunch' theorem suggests that no single model algorithm can outperform at all possible tasks. Here, we explore the capabilities of classical (SVR), FSLC, and transformer models (MolBART) over a range of dataset tasks and show a 'goldilocks zone' for each model type, in which dataset size and feature distribution (i.e. dataset "diversity") determines the optimal algorithm strategy. When datasets are small ( < 50 molecules), FSLC tend to outperform both classical ML and transformers. When datasets are small-to-medium sized (50-240 molecules) and diverse, transformers outperform both classical models and few-shot learning. Finally, when datasets are of larger and of sufficient size, classical models then perform the best, suggesting that the optimal model to choose likely depends on the dataset available, its size and diversity. These findings may help to answer the perennial question of which ML algorithm is to be used when faced with a new dataset.

Exploring SureChEMBL from a drug discovery perspective

In the pharmaceutical industry, the patent protection of drugs and medicines is accorded importance because of the high costs involved in the development of novel drugs. Over the years, researchers have analyzed patent documents to identify freedom-to-operate spaces for novel drug candidates. To assist this, several well-established public patent document data repositories have enabled automated methodologies for extracting information on therapeutic agents. In this study, we delve into one such publicly available patent database, SureChEMBL, which catalogues patent documents related to life sciences. Our exploration begins by identifying patent compounds across public chemical data resources, followed by pinpointing sections in patent documents where the chemical annotations were found. Next, we exhibit the potential of compounds to serve as drug candidates by evaluating their conformity to drug-likeness criteria. Lastly, we examine the drug development stage reported for these compounds to understand their clinical success. In summary, our investigation aims at providing a comprehensive overview of the patent compounds catalogued in SureChEMBL, assessing their relevance to pharmaceutical drug discovery.

Expanding Complex Morpholines Using Systematic Chemical Diversity

The morpholine heterocycle is a structural unit found in many bioactive compounds and FDA-approved drugs, but the generation of more complex C-functionalized morpholine derivatives remains considerably underexplored. Using systematic chemical diversity (SCD), a concept that guides the expansion of saturated drug-like scaffolds through regiochemical and stereochemical variation, we describe the synthesis of a collection of methyl-substituted morpholine acetic acid esters starting from enantiomerically pure amino acids and amino alcohols. In total, 24 diverse substituted morpholines were produced that vary systematically in regiochemistry and stereochemistry (relative and absolute). These diverse C-substituted morpholines can be directly applied in fragment screening or incorporated as building blocks in medicinal chemistry and library synthesis.

Chemical Multiverse and Diversity of Food Chemicals

Food chemicals have a fundamental role in our lives, with an extended impact on nutrition, disease prevention, and marked economic implications in the food industry. The number of food chemical compounds in public databases has substantially increased in the past few years, which can be characterized using chemoinformatics approaches. We and other groups explored public food chemical libraries containing up to 26,500 compounds. This study aimed to analyze the chemical contents, diversity, and coverage in the chemical space of food chemicals and additives and, from here on, food components. The approach to food components addressed in this study is a public database with more than 70,000 compounds, including those predicted via omics techniques. It was concluded that food components have distinctive physicochemical properties and constitutional descriptors despite sharing many chemical structures with natural products. Food components, on average, have large molecular weights and several apolar structures with saturated hydrocarbons. Compared to reference databases, food component structures have low scaffold and fingerprint-based diversity and high structural complexity, as measured by the fraction of sp3 carbons. These structural features are associated with a large fraction of macronutrients as lipids. Lipids in food components were decompiled by an analysis of the maximum common substructures. The chemical multiverse representation of food chemicals showed a larger coverage of chemical space than natural products and FDA-approved drugs by using different sets of representations.

Techniques and Strategies in Drug Design and Discovery

The process of drug discovery constitutes a highly intricate and formidable undertaking, encompassing the identification and advancement of novel therapeutic entities [...].

Exploring Scoring Function Space: Developing Computational Models for Drug Discovery

BACKGROUND

The idea of scoring function space established a systems-level approach to address the development of models to predict the affinity of drug molecules by those interested in drug discovery.

OBJECTIVE

Our goal here is to review the concept of scoring function space and how to explore it to develop machine learning models to address protein-ligand binding affinity.

METHODS

We searched the articles available in PubMed related to the scoring function space. We also utilized crystallographic structures found in the protein data bank (PDB) to represent the protein space.

RESULTS

The application of systems-level approaches to address receptor-drug interactions allows us to have a holistic view of the process of drug discovery. The scoring function space adds flexibility to the process since it makes it possible to see drug discovery as a relationship involving mathematical spaces.

CONCLUSION

The application of the concept of scoring function space has provided us with an integrated view of drug discovery methods. This concept is useful during drug discovery, where we see the process as a computational search of the scoring function space to find an adequate model to predict receptor-drug binding affinity.

Identification of Novel β-Tubulin Inhibitors Using a Combined In Silico/In Vitro Approach

Due to their potential as leads for various therapeutic applications, including as antimitotic and antiparasitic agents, the development of tubulin inhibitors offers promise for drug discovery. In this study, an in silico pharmacophore-based virtual screening approach targeting the colchicine binding site of β-tubulin was employed. Several structure- and ligand-based models for known tubulin inhibitors were generated. Compound databases were virtually screened against the models, and prioritized hits from the SPECS compound library were tested in an in vitro tubulin polymerization inhibition assay for their experimental validation. Out of the 41 SPECS compounds tested, 11 were active tubulin polymerization inhibitors, leading to a prospective true positive hit rate of 26.8%. Two novel inhibitors displayed IC50 values in the range of colchicine. The most potent of which was a novel acetamide-bridged benzodiazepine/benzimidazole derivative with an IC50 = 2.9 μM. The screening workflow led to the identification of diverse inhibitors active at the tubulin colchicine binding site. Thus, the pharmacophore models show promise as valuable tools for the discovery of compounds and as potential leads for the development of cancer therapeutic agents.

Identification of Potential Aldose Reductase Inhibitors Using Convolutional Neural Network-Based in Silico Screening

Aldose reductase (ALR2) is a notable enzyme of the polyol pathway responsible for aggravating diabetic neuropathy complications. The first step begins when it catalyzes the reduction of glucose to sorbitol with NADPH as a coenzyme. Elevated concentrations of sorbitol damage the tissues, leading to complications like neuropathy. Though considerable effort has been pushed toward the successful discovery of potent inhibitors, its discovery still remains an elusive task. To this end, we present a 3D convolutional neural network (3D-CNN) based ALR2 inhibitor classification technique by dealing with snapshots of images captured from 3D chemical structures with multiple rotations as input data. The CNN-based architecture was trained on the 360 sets of image data along each axis and further prediction on the Maybridge library by each of the models. Subjecting the retrieved hits to molecular docking leads to the identification of the top 10 molecules with high binding affinity. The hits displayed a better blood-brain barrier penetration (BBB) score (90% with more than four scores) as compared to standard inhibitors (38%), reflecting the superior BBB penetrating efficiency of the hits. Followed by molecular docking, the biological evaluation spotlighted five compounds as promising ALR2 inhibitors and can be considered as a likely prospect for further structural optimization with medicinal chemistry efforts to improve their inhibition efficacy and consolidate them as new ALR2 antagonists in the future. In addition, the study also demonstrated the usefulness of scaffold analysis of the molecules as a method for investigating the significance of structurally diverse compounds in data-driven studies. For reproducibility and accessibility purposes, all of the source codes used in our study are publicly available.

Systematic Evaluation of Local and Global Machine Learning Models for the Prediction of ADME Properties

Machine learning (ML) has become an indispensable tool to predict absorption, distribution, metabolism, and excretion (ADME) properties in pharmaceutical research. ML algorithms are trained on molecular structures and corresponding ADME assay data to develop quantitative structure-property relationship (QSPR) models. Traditional QSPR models were trained on compound sets of limited size. With the advent of more complex ML algorithms and data availability, training sets have become larger and more diverse. Most common training approaches consist in either training a model with a small set of similar compounds, namely, compounds designed for the same drug discovery project or chemical series (local model approach) or with a larger set of diverse compounds (global model approach). Global models are built with all experimental data available for an assay, combining compound data from different projects and disease areas. Despite the ML progress made so far, the choice of the appropriate data composition for building ML models is still unclear. Herein, a systematic evaluation of local and global ML models was performed for 10 different experimental assays and 112 drug discovery projects. Results show a consistent superior performance of global models for ADME property predictions. Diagnostic analyses were also carried out to investigate the influence of training set size, structural diversity, and data shift in the relative performance of local and global ML models. Training set and structural diversity did not have an impact in the relative performance on the methods. Instead, data shift helped to identify the projects with larger performance differences between local and global models. Results presented in this work can be leveraged to improve ML-based ADME properties predictions and thus decision-making in drug discovery projects.

Scaffold and Structural Diversity of the Secondary Metabolite Space of Medicinal Fungi

Medicinal fungi, including mushrooms, have well-documented therapeutic uses. In this study, we perform a cheminformatics-based investigation of the scaffold and structural diversity of the secondary metabolite space of medicinal fungi and, moreover, perform a detailed comparison with approved drugs, other natural product libraries, and semi-synthetic libraries. We find that the secondary metabolite space of medicinal fungi has similar or higher scaffold diversity in comparison to other natural product libraries analyzed here. Notably, 94% of the scaffolds in the secondary metabolite space of medicinal fungi are not present in the approved drugs. Further, we find that the secondary metabolites, on the one hand, are structurally far from the approved drugs, while, on the other hand, they are close in terms of molecular properties to the approved drugs. Lastly, chemical space visualization using dimensionality reduction methods showed that the secondary metabolite space has minimal overlap with the approved drug space. In a nutshell, our results underscore that the secondary metabolite space of medicinal fungi is a valuable resource for identifying potential lead molecules for natural product-based drug discovery.

Comparative analysis of an anthraquinone and chalcone derivatives-based virtual combinatorial library. A cheminformatics "proof-of-concept" study

A Laplacian scoring algorithm for gene selection and the Gini coefficient to identify the genes whose expression varied least across a large set of samples were the state-of-the-art methods used here. These methods have not been trialed for their feasibility in cheminformatics. This was a maiden attempt to investigate a complete comparative analysis of an anthraquinone and chalcone derivatives-based virtual combinatorial library. This computational "proof-of-concept" study illustrated the combinatorial approach used to explain how the structure of the selected natural products (NPs) undergoes molecular diversity analysis. A virtual combinatorial library (1.6 M) based on 20 anthraquinones and 24 chalcones was enumerated. The resulting compounds were optimized to the near drug-likeness properties, and the physicochemical descriptors were calculated for all datasets including FDA, Non-FDA, and NPs from ZINC 15. UMAP and PCA were applied to compare and represent the chemical space coverage of each dataset. Subsequently, the Laplacian score and Gini coefficient were applied to delineate feature selection and selectivity among properties, respectively. Finally, we demonstrated the diversity between the datasets by employing Murcko's and the central scaffolds systems, calculating three fingerprint descriptors and analyzing their diversity by PCA and SOM. The optimized enumeration resulted in 1,610,268 compounds with NP-Likeness, and synthetic feasibility mean scores close to FDA, Non-FDA, and NPs datasets. The overlap between the chemical space of the 1.6 M database was more prominent than with the NPs dataset. A Laplacian score prioritized NP-likeness and hydrogen bond acceptor properties (1.0 and 0.923), respectively, while the Gini coefficient showed that all properties have selective effects on datasets (0.81-0.93). Scaffold and fingerprint diversity indicated that the descending order for the tested datasets was FDA, Non-FDA, NPs and 1.6 M. Virtual combinatorial libraries based on NPs can be considered as a source of the combinatorial compound with NP-likeness properties. Furthermore, measuring molecular diversity is supposed to be performed by different methods to allow for comparison and better judgment.

In Silico Drug Repurposing Framework Predicts Repaglinide, Agomelatine and Protokylol as TRPV1 Modulators with Analgesic Activity

Pain is one of the most common symptoms experienced by patients. The use of current analgesics is limited by low efficacy and important side effects. Transient receptor potential vanilloid-1 (TRPV1) is a non-selective cation channel, activated by capsaicin, heat, low pH or pro-inflammatory agents. Since TRPV1 is a potential target for the development of novel analgesics due to its distribution and function, we aimed to develop an in silico drug repositioning framework to predict potential TRPV1 ligands among approved drugs as candidates for treating various types of pain. Structures of known TRPV1 agonists and antagonists were retrieved from ChEMBL databases and three datasets were established: agonists, antagonists and inactive molecules (pIC50 or pEC50 < 5 M). Structures of candidates for repurposing were retrieved from the DrugBank database. The curated active/inactive datasets were used to build and validate ligand-based predictive models using Bemis−Murcko structural scaffolds, plain ring systems, flexophore similarities and molecular descriptors. Further, molecular docking studies were performed on both active and inactive conformations of the TRPV1 channel to predict the binding affinities of repurposing candidates. Variables obtained from calculated scaffold-based activity scores, molecular descriptors criteria and molecular docking were used to build a multi-class neural network as an integrated machine learning algorithm to predict TRPV1 antagonists and agonists. The proposed predictive model had a higher accuracy for classifying TRPV1 agonists than antagonists, the ROC AUC values being 0.980 for predicting agonists, 0.972 for antagonists and 0.952 for inactive molecules. After screening the approved drugs with the validated algorithm, repaglinide (antidiabetic) and agomelatine (antidepressant) emerged as potential TRPV1 antagonists, and protokylol (bronchodilator) as an agonist. Further studies are required to confirm the predicted activity on TRPV1 and to assess the candidates’ efficacy in alleviating pain.

Novel 5-Nitrofuran-Tagged Imidazo-Fused Azines and Azoles Amenable by the Groebke–Blackburn–Bienaymé Multicomponent Reaction: Activity Profile against ESKAPE Pathogens and Mycobacteria

A chemically diverse set of 13 5-nitrofuran-tagged heterocyclic compounds has been prepared via the Groebke–Blackburn–Bienaymé multicomponent reaction. The testing of these compounds against the so-called ESKAPE panel of pathogens identified an apparent lead compound—N-cyclohexyl-2-(5-nitrofuran-2-yl)imidazo[1,2-a]pyridine-3-amine (4a)—which showed an excellent profile against Enterobacter cloacae, Staphylococcus aureus, Klebsiella pneumoniae, and Enterococcus faecalis (MIC 0.25, 0.06, 0.25 and 0.25 µg/mL, respectively). Its antibacterial profile and practically convenient synthesis warrant further pre-clinical development. Certain structure-activity relationships were established in the course of this study which were rationalized by the flexible docking experiments in silico. The assessment of antitubercular potential of the compounds synthesized against drug sensitive H37v strain of Mycobacterium tuberculosis revealed little potential of the imidazo-fused products of the Groebke–Blackburn–Bienaymé multicomponent reaction as chemotherapeutic agents against this pathogen.

Synergy between machine learning and natural products cheminformatics: Application to the lead discovery of anthraquinone derivatives

Cheminformatics utilizing machine learning (ML) techniques have opened up a new horizon in drug discovery. This is owing to vast chemical space expansion with rocketing numbers of expected hits and lead compounds that match druggable macromolecular targets, in particular from natural compounds. Due to the natural products' (NP) structural complexity, uniqueness, and diversity, they could occupy a bigger space in pharmaceuticals, allowing the industry to pursue more selective leads in the nanomolar range of binding affinity. ML is an essential part of each step of the drug design pipeline, such as target prediction, compound library preparation, and lead optimization. Notably, molecular mechanic and dynamic simulations, induced docking, and free energy perturbations are essential in predicting best binding poses, binding free energy values, and molecular mechanics force fields. Those applications have leveraged from artificial intelligence (AI), which decreases the computational costs required for such costly simulations. This review aimed to describe chemical space and compound libraries related to NPs. High-throughput screening utilized for fractionating NPs and high-throughput virtual screening and their strategies, and significance, are reviewed. Particular emphasis was given to AI approaches, ML tools, algorithms, and techniques, especially in drug discovery of macrocyclic compounds and approaches in computer-aided and ML-based drug discovery. Anthraquinone derivatives were discussed as a source of new lead compounds that can be developed using ML tools for diverse medicinal uses such as cancer, infectious diseases, and metabolic disorders. Furthermore, the power of principal component analysis in understanding relevant protein conformations, and molecular modeling of protein-ligand interaction were also presented. Apart from being a concise reference for cheminformatics, this review is a useful text to understand the application of ML-based algorithms to molecular dynamics simulation and in silico absorption, distribution, metabolism, excretion, and toxicity prediction.

A unified "top-down" approach for the synthesis of diverse lead-like molecular scaffolds

A "top-down" synthetic approach enabled the step-efficient synthesis of 21 diverse novel molecular scaffolds. The scaffolds were derived from four complex intermediates that had been prepared using cycloaddition chemistry. Scaffold-hopping of these intermediates was achieved through attachment of an additional ring, ring cleavage, ring expansion and/or ring fusion. It was shown that the resulting scaffolds could be decorated to yield diverse lead-like screening compounds.

Large-Scale Evaluation of Collision Cross Sections to Investigate Blood-Brain Barrier Permeation of Drugs

Successful drug administration to the central nervous system requires accurate adjustment of the drugs’ molecular properties. Therefore, structure-derived descriptors of potential brain therapeutic agents are essential for an early evaluation of pharmacokinetics during drug development. The collision cross section (CCS) of molecules was recently introduced as a novel measurable parameter to describe blood-brain barrier (BBB) permeation. This descriptor combines molecular information about mass, structure, volume, branching and flexibility. As these chemical properties are known to influence cerebral pharmacokinetics, CCS determination of new drug candidates may provide important additional spatial information to support existing models of BBB penetration of drugs. Besides measuring CCS, calculation is also possible; but however, the reliability of computed CCS values for an evaluation of BBB permeation has not yet been fully investigated. In this work, prediction tools based on machine learning were used to compute CCS values of a large number of compounds listed in drug libraries as negative or positive with respect to brain penetration (BBB⁺ and BBB⁻ compounds). Statistical evaluation of computed CCS and several other descriptors could prove the high value of CCS. Further, CCS-deduced maximum molecular size of BBB⁺ drugs matched the dimensions of BBB pores. A threshold for transcellular penetration and possible permeation through pore-like openings of cellular tight-junctions is suggested. In sum, CCS evaluation with modern in silico tools shows high potential for its use in the drug development process.

SynthI: A New Open-Source Tool for Synthon-Based Library Design

Most of the existing computational tools for de novo library design are focused on the generation, rational selection, and combination of promising structural motifs to form members of the new library. However, the absence of a direct link between the chemical space of the retrosynthetically generated fragments and the pool of available reagents makes such approaches appear as rather theoretical and reality-disconnected. In this context, here we present Synthons Interpreter (SynthI), a new open-source toolkit for de novo library design that allows merging those two chemical spaces into a single synthons space. Here synthons are defined as actual fragments with valid valences and special labels, specifying the position and the nature of reactive centers. They can be issued from either the "breakup" of reference compounds according to 38 retrosynthetic rules or real reagents, after leaving group withdrawal or transformation. Such an approach not only enables the design of synthetically accessible libraries and analog generation but also facilitates reagents (building blocks) analysis in the medicinal chemistry context. SynthI code is publicly available at https://github.com/Laboratoire-de-Chemoinformatique/SynthI.

Application of Networking Approaches to Assess the Chemical Diversity, Biogeography, and Pharmaceutical Potential of Verongiida Natural Products

This study provides a review of all isolated natural products (NPs) reported for sponges within the order Verongiida (1960 to May 2020) and includes a comprehensive compilation of their geographic and physico-chemical parameters. Physico-chemical parameters were used in this study to infer pharmacokinetic properties as well as the potential pharmaceutical potential of NPs from this order of marine sponge. In addition, a network analysis for the NPs produced by the Verongiida sponges was applied to systematically explore the chemical space relationships between taxonomy, secondary metabolite and drug score variables, allowing for the identification of differences and correlations within a dataset. The use of scaffold networks as well as bipartite relationship networks provided a platform to explore chemical diversity as well as the use of chemical similarity networks to link pharmacokinetic properties with structural similarity. This study paves the way for future applications of network analysis procedures in the field of natural products for any order or family.

Magic Rings: Navigation in the Ring Chemical Space Guided by the Bioactive Rings

The large majority of bioactive molecules contain a more or less complex ring system as a central structural element. This central core determines the basic molecule shape, keeps substituents in their proper positions, and often also contributes to the biological activity itself. In this study the ring systems extracted from one billion molecules are processed and differences between rings from bioactive molecules and common synthetic molecules are analyzed. The bioactive rings seem to be distributed throughout the large portion of chemical space, but not uniformly; one can see several more dense regions, where the bioactive rings often appear in small clusters, as well as empty areas. A web tool offering an interactive navigation in the ring chemical space and supporting identification of bioisosteric ring analogs available at https://bit.ly/magicrings is also described.

Exploring the Effect of Structure-Based Scaffold Hopping on the Inhibition of Coxsackievirus A24v Transduction by Pentavalent N-Acetylneuraminic Acid Conjugates

Coxsackievirus A24 variant (CVA24v) is the primary causative agent of the highly contagious eye infection designated acute hemorrhagic conjunctivitis (AHC). It is solely responsible for two pandemics and several recurring outbreaks of the disease over the last decades, thus affecting millions of individuals throughout the world. To date, no antiviral agents or vaccines are available for combating this disease, and treatment is mainly supportive. CVA24v utilizes Neu5Ac-containing glycans as attachment receptors facilitating entry into host cells. We have previously reported that pentavalent Neu5Ac conjugates based on a glucose-scaffold inhibit CVA24v infection of human corneal epithelial cells. In this study, we report on the design and synthesis of scaffold-replaced pentavalent Neu5Ac conjugates and their effect on CVA24v cell transduction and the use of cryogenic electron microscopy (cryo-EM) to study the binding of these multivalent conjugates to CVA24v. The results presented here provide insights into the development of Neu5Ac-based inhibitors of CVA24v and, most significantly, the first application of cryo-EM to study the binding of a multivalent ligand to a lectin.

Pseudo Natural Products-Chemical Evolution of Natural Product Structure

Pseudo-natural products (PNPs) combine natural product (NP) fragments in novel arrangements not accessible by current biosynthesis pathways. As such they can be regarded as non-biogenic fusions of NP-derived fragments. They inherit key biological characteristics of the guiding natural product, such as chemical and physiological properties, yet define small molecule chemotypes with unprecedented or unexpected bioactivity. We iterate the design principles underpinning PNP scaffolds and highlight their syntheses and biological investigations. We provide a cheminformatic analysis of PNP collections assessing their molecular properties and shape diversity. We propose and discuss how the iterative analysis of NP structure, design, synthesis, and biological evaluation of PNPs can be regarded as a human-driven branch of the evolution of natural products, that is, a chemical evolution of natural product structure.

Design of gp120 HIV-1 entry inhibitors by scaffold hopping via isosteric replacements

We present the development of alternative scaffolds and validation of their synthetic pathways as a tool for the exploration of new HIV gp120 inhibitors based on the recently discovered inhibitor of this class, NBD-14136. The new synthetic routes were based on isosteric replacements of the amine and acid precursors required for the synthesis of NBD-14136, guided by molecular modeling and chemical feasibility analysis. To ensure that these synthetic tools and new scaffolds had the potential for further exploration, we eventually tested few representative compounds from each newly designed scaffold against the gp120 inhibition assay and cell viability assays.

Synthesis of Three-Dimensional (Di)Azatricyclododecene Scaffold and Its Application to Peptidomimetics

A novel sp³ carbon-rich tricyclic 3D scaffold-based peptide mimetic compound library was constructed to target protein-protein interactions. Tricyclic framework 7 was synthesized from 9-azabicyclo[3,3,1]nonan-3-one (11) via a gold(I)-catalyzed Conia-ene reaction. The electron-donating group on the pendant alkyne of cyclization precursor 12 b-e was the key to forming 6-endo-dig cyclized product 7 with complete regioselectivity. Using the synthetic strategy for regioselective construction of bridged tricyclic framework 7, a diazatricyclododecene 3D-scaffold 8 a, which enables the introduction of substituents into the scaffold to mimic amino acid side chains, was designed and synthesized. The peptide mimetics 21 a-u were synthesized via step-by-step installation of three substituents on diazatricyclododecene scaffold 8 a. Compounds 21 a-h were synthesized as α-helix peptide mimics of hydrophobic ZZxxZ and ZxxZZ sequences (Z=Leu or Phe) and subjected to cell-based assays: antiproliferative activity, HIF-1 transcriptional activity which is considered to affect cancer malignancy, and antiviral activity against rabies virus. Compound 21 a showed the strongest inhibitory activity of HIF-1 transcriptional activity (IC₅₀ =4.1±0.8 μM), whereas compounds 21 a-g showed antiviral activity with IC₅₀ values of 4.2-12.4 μM, suggesting that the 3D-scaffold 8 a has potential as a versatile peptide mimic skeleton.

Chemical Space Exploration of DprE1 Inhibitors Using Chemoinformatics and Artificial Intelligence

Tuberculosis (TB), entrained by Mycobacterium tuberculosis, continues to be an enfeebling disease, killing nearly 1.5 million people in 2019, with 2 billion people worldwide affected by latent TB. The multidrug-resistant and totally drug-resistant emerging strains further exacerbate the TB infection. The cell wall of bacteria provides critical virulence components such as cell surface proteins, regulators, signal transduction proteins, and toxins. The cell wall biosynthesis pathway of Mycobacterium tuberculosis is exhaustively studied to discover novel drug targets. Decaprenylphosphoryl-β-d-ribose-2'-epimerase (DprE1) is an important enzyme involved in the arabinogalactan biosynthetic pathway of Mycobacterium tuberculosis cell wall and is essential for both latent and persistent bacterial infection. We analyzed all known ∼1300 DprE1 inhibitors to gain deep insights into the chemogenomic space of DprE1-ligand complexes. Physicochemical descriptors of the DprE1 inhibitors showed a marked lipophilic character forming a cluster distinct from the existing TB drugs, as revealed by the principal component analysis. Similarity analysis using Murcko scaffolds and rubber band scaling revealed scarce representation of the chemical space. Further, Murcko scaffold analysis uncovered favorable and unfavorable scaffolds, where benzo and pyridine-based core scaffolds exhibit the highest biological activity, as evidenced by their MIC and IC₅₀ values. Automatic SAR and R-group decomposition analysis resulted in the identification of substructures responsible for the inhibitory activity of the DprE1 enzyme. Further, with activity cliff analysis, we observed prominent discontinuity in the SAR of DprE1 inhibitors, where even simple structural modification in the chemical scaffold resulted in significant potency difference, presumably due to the binding orientation and interaction in the active site. Thiophene, 6-membered aromatic rings, and unsubstituted benzene ring-based toxicophores were identified in the DprE1 chemical space using an artificial intelligence approach based on inductive logic programming. This paper, hence, ushers in new insights for the design and development of potent covalent and non-covalent DprE1 inhibitors and guides hit and lead optimization for the development of non-hazardous small molecule therapeutics for Mycobacterium tuberculosis.

Natural product-informed exploration of chemical space to enable bioactive molecular discovery

The search for new bioactive molecules remains an open challenge limiting our ability to discover new drugs to treat disease and chemical probes to comprehensively study biological processes. The vastness of chemical space renders its exploration unfeasible by synthesis alone. Historically, chemists have tended to explore chemical space unevenly without committing to systematic frameworks for navigation. This minireview covers a range of approaches that take inspiration from the structure or origin of natural products, and help focus molecular discovery on biologically-relevant regions of chemical space. All these approaches have enabled the discovery of distinctive and novel bioactive small molecules such as useful chemical probes of biological mechanisms. This minireview comments on how such approaches may be developed into more general frameworks for the systematic identification of currently unexplored regions of biologically-relevant chemical space, a challenge that is central to both chemical biology and medicinal chemistry.

Derivatization Design of Synthetically Accessible Space for Optimization: In Silico Synthesis vs Deep Generative Design

Molecular design is of utmost importance in lead optimization programs ultimately determining the fate of the project and the speed to reach preclinical stage. Newly designed lead analogues or new chemotypes must successfully address the challenges in the multidimensional optimization process throughout several optimization cycles. The speed, quality, and creativity of the designs can have a major impact on the cycle time, the number of required cycles, and the number of compounds needed to be synthesized and evaluated that in combination affect the overall timeline and cost of the lead optimization phase. Recently, a new concept, generative design with deep learning, has become popular for de novo design of project relevant analogue sets. We have developed a de novo design technology called "derivatization design" that applies artificial-intelligence-assisted forward in silico synthesis for the generation of near neighbor lead analogues as well as scaffold variations. The several attractive features of the methodology include synthetic feasibility, reagent availability and cost data associated with each new molecule; thus, detailed synthetic assessment is automatically generated during the design. As a result, these practically important data types can become an early part of the ranking and selection process for cycle time reduction. The power of derivatization design is demonstrated in a simple design study of DDR1 inhibitors and comparison of the produced molecules to a recently published data set obtained with deep generative design.

Efficient unified synthesis of diverse bridged polycyclic scaffolds using a complexity-generating 'stitching' annulation approach

Regioselective and stereospecific directed C-H arylation of simple amine substrates, and cyclisation, delivered 30 diverse, three-dimensional scaffolds. The unified approach significantly expanded the range of bridged ring systems that contain both a nitrogen atom and an aromatic ring.

Exploiting cheminformatic and machine learning to navigate the available chemical space of potential small molecule inhibitors of SARS-CoV-2

The current life-threatening and tenacious pandemic eruption of coronavirus disease in 2019 (COVID-19) has posed a significant global hazard concerning high mortality rate, economic meltdown, and everyday life distress. The rapid spread of COVID-19 demands countermeasures to combat this deadly virus. Currently, there are no drugs approved by the FDA to treat COVID-19. Therefore, discovering small molecule therapeutics for treating COVID-19 infection is essential. So far, only a few small molecule inhibitors are reported for coronaviruses. There is a need to expand the small chemical space of coronaviruses inhibitors by adding potent and selective scaffolds with anti-COVID activity. In this context, the huge antiviral chemical space already available can be analysed using cheminformatic and machine learning to unearth new scaffolds. We created three specific datasets called "antiviral dataset" (N = 38,428) "drug-like antiviral dataset" (N = 20,963) and "anticorona dataset" (N = 433) for this purpose. We analyzed the 433 molecules of "anticorona dataset" for their scaffold diversity, physicochemical distributions, principal component analysis, activity cliffs, R-group decomposition, and scaffold mapping. The scaffold diversity of the "anticorona dataset" in terms of Murcko scaffold analysis demonstrates a thorough representation of diverse chemical scaffolds. However, physicochemical descriptor analysis and principal component analysis demonstrated negligible drug-like features for the "anticorona dataset" molecules. The "antiviral dataset" and "drug-like antiviral dataset" showed low scaffold diversity as measured by the Gini coefficient. The hierarchical clustering of the "antiviral dataset" against the "anticorona dataset" demonstrated little molecular similarity. We generated a library of frequent fragments and polypharmacological ligands targeting various essential viral proteins such as main protease, helicase, papain-like protease, and replicase polyprotein 1ab. Further structural and chemical features of the "anticorona dataset" were compared with SARS-CoV-2 repurposed drugs, FDA-approved drugs, natural products, and drugs currently in clinical trials. Using machine learning tool DCA (DMax Chemistry Assistant), we converted the "anticorona dataset" into an elegant hypothesis with significant functional biological relevance. Machine learning analysis uncovered that FDA approved drugs, Tizanidine HCl, Cefazolin, Raltegravir, Azilsartan, Acalabrutinib, Luliconazole, Sitagliptin, Meloxicam (Mobic), Succinyl sulfathiazole, Fluconazole, and Pranlukast could be repurposed as effective drugs for COVID-19. Fragment-based scaffold analysis and R-group decomposition uncovered pyrrolidine and the indole molecular scaffolds as the potent fragments for designing and synthesizing the novel drug-like molecules for targeting SARS-CoV-2. This comprehensive and systematic assessment of small-molecule viral therapeutics' entire chemical space realised critical insights to potentially privileged scaffolds that could aid in enrichment and rapid discovery of efficacious antiviral drugs for COVID-19.

Uncertainty quantification in drug design

Machine learning and artificial intelligence are increasingly being applied to the drug-design process as a result of the development of novel algorithms, growing access, the falling cost of computation and the development of novel technologies for generating chemically and biologically relevant data. There has been recent progress in fields such as molecular de novo generation, synthetic route prediction and, to some extent, property predictions. Despite this, most research in these fields has focused on improving the accuracy of the technologies, rather than on quantifying the uncertainty in the predictions. Uncertainty quantification will become a key component in autonomous decision making and will be crucial for integrating machine learning and chemistry automation to create an autonomous design-make-test-analyse cycle. This review covers the empirical, frequentist and Bayesian approaches to uncertainty quantification, and outlines how they can be used for drug design. We also outline the impact of uncertainty quantification on decision making.

Cheminformatics to Characterize Pharmacologically Active Natural Products

Natural products have a significant role in drug discovery. Natural products have distinctive chemical structures that have contributed to identifying and developing drugs for different therapeutic areas. Moreover, natural products are significant sources of inspiration or starting points to develop new therapeutic agents. Natural products such as peptides and macrocycles, and other compounds with unique features represent attractive sources to address complex diseases. Computational approaches that use chemoinformatics and molecular modeling methods contribute to speed up natural product-based drug discovery. Several research groups have recently used computational methodologies to organize data, interpret results, generate and test hypotheses, filter large chemical databases before the experimental screening, and design experiments. This review discusses a broad range of chemoinformatics applications to support natural product-based drug discovery. We emphasize profiling natural product data sets in terms of diversity; complexity; acid/base; absorption, distribution, metabolism, excretion, and toxicity (ADME/Tox) properties; and fragment analysis. Novel techniques for the visual representation of the chemical space are also discussed.

Autonomous Discovery in the Chemical Sciences Part II: Outlook

This two-part Review examines how automation has contributed to different aspects of discovery in the chemical sciences. In this second part, we reflect on a selection of exemplary studies. It is increasingly important to articulate what the role of automation and computation has been in the scientific process and how that has or has not accelerated discovery. One can argue that even the best automated systems have yet to "discover" despite being incredibly useful as laboratory assistants. We must carefully consider how they have been and can be applied to future problems of chemical discovery in order to effectively design and interact with future autonomous platforms. The majority of this Review defines a large set of open research directions, including improving our ability to work with complex data, build empirical models, automate both physical and computational experiments for validation, select experiments, and evaluate whether we are making progress towards the ultimate goal of autonomous discovery. Addressing these practical and methodological challenges will greatly advance the extent to which autonomous systems can make meaningful discoveries.

In Silico Design and Analysis of a Kinase-Focused Combinatorial Library Considering Diversity and Quality

A structurally diverse, high-quality, and kinase-focused database plays a critical role in finding hits or leads in kinase drug discovery. Here, we propose a workflow for designing a virtual kinase-focused combinatorial library using existing structures. Based on the analysis of known protein kinase inhibitors (PKIs), detailed fragment optimization, fragment selection, fragment linking, and a molecular filtering scheme were defined. Quick recognition of core fragments that can possibly form dual hydrogen bonds with the hinge region of the ATP-pocket was proposed. Furthermore, three diversity and four quality metrics were chosen for compound library analysis, which can be applied to databases with over 30 million structures. Compared with 13 commercial libraries, our protocol demonstrates a special advantage in terms of good skeleton diversity, acceptable fingerprint diversity, balanced scaffold distribution, and high quality, which can work well not only on existing PKIs, but also on four chosen commercial libraries. Overall, the strategy can greatly facilitate the expansion of a desirable chemical space for kinase drug discovery.

A general approach for retrosynthetic molecular core analysis

Scaffold analysis of compound data sets has reemerged as a chemically interpretable alternative to machine learning for chemical space and structure–activity relationships analysis. In this context, analog series-based scaffolds (ASBS) are synthetically relevant core structures that represent individual series of analogs. As an extension to ASBS, we herein introduce the development of a general conceptual framework that considers all putative cores of molecules in a compound data set, thus softening the often applied “single molecule–single scaffold” correspondence. A putative core is here defined as any substructure of a molecule complying with two basic rules: (a) the size of the core is a significant proportion of the whole molecule size and (b) the substructure can be reached from the original molecule through a succession of retrosynthesis rules. Thereafter, a bipartite network consisting of molecules and cores can be constructed for a database of chemical structures. Compounds linked to the same cores are considered analogs. We present case studies illustrating the potential of the general framework. The applications range from inter- and intra-core diversity analysis of compound data sets, structure–property relationships, and identification of analog series and ASBS. The molecule–core network herein presented is a general methodology with multiple applications in scaffold analysis. New statistical methods are envisioned that will be able to draw quantitative conclusions from these data. The code to use the method presented in this work is freely available as an additional file. Follow-up applications include analog searching and core structure–property relationships analyses.

Computational Drug Repurposing Algorithm Targeting TRPA1 Calcium Channel as a Potential Therapeutic Solution for Multiple Sclerosis

Multiple sclerosis (MS) is a chronic autoimmune disease affecting the central nervous system (CNS) through neurodegeneration and demyelination, leading to physical/cognitive disability and neurological defects. A viable target for treating MS appears to be the Transient Receptor Potential Ankyrin 1 (TRPA1) calcium channel, whose inhibition has been shown to have beneficial effects on neuroglial cells and protect against demyelination. Using computational drug discovery and data mining methods, we performed an in silico screening study combining chemical graph mining, quantitative structure-activity relationship (QSAR) modeling, and molecular docking techniques in a global prediction model in order to identify repurposable drugs as potent TRPA1 antagonists that may serve as potential treatments for MS patients. After screening the DrugBank database with the combined generated algorithm, 903 repurposable structures were selected, with 97 displaying satisfactory inhibition probabilities and pharmacokinetics. Among the top 10 most probable inhibitors of TRPA1 with good blood brain barrier (BBB) permeability, desvenlafaxine, paliperidone, and febuxostat emerged as the most promising repurposable agents for treating MS. Molecular docking studies indicated that desvenlafaxine, paliperidone, and febuxostat are likely to induce allosteric TRPA1 channel inhibition. Future in vitro and in vivo studies are needed to confirm the biological activity of the selected hit molecules.

Expansion of the structure-activity relationships of BACE1 inhibitors by harnessing diverse building blocks prepared using a unified synthetic approach

The structural diversity of β-site amyloid precursor protein cleaving enzyme 1 (BACE1) inhibitors was expanded by harnessing diverse building blocks that had been prepared via a unified lead-oriented synthetic approach. It was shown that the lipophilic cyclohexylmethyl group within a known series of BACE1 inhibitors could be productively replaced with a range of alternative ring systems.

Screening Library Design

Thanks to technological advances and a greater understanding of the biological and chemical natures of targets and related diseases, high-throughput screening (HTS) has been allowed to be faster, cheaper, and more accessible. Yet, despite these increased technologies and understandings, the frequency of novel and drugs are being approved each year has not being increasing over the years. 2017 was considered a "bumper" year with a total of 46 approved drugs, over double that of the previous year. However, it is thought that part of the problem that HTS has not lived up to expectations is because of the contents of current chemical libraries. Therefore, new methods to design screening libraries are of great interest.