The purchasable chemical space: a detailed picture

Protein Structure-Based Organic Chemistry-Driven Ligand Design from Ultralarge Chemical Spaces

Ultralarge chemical spaces describing several billion compounds are revolutionizing hit identification in early drug discovery. Because of their size, such chemical spaces cannot be fully enumerated and require ad-hoc computational tools to navigate them and pick potentially interesting hits. We here propose a structure-based approach to ultralarge chemical space screening in which commercial chemical reagents are first docked to the target of interest and then directly connected according to organic chemistry and topological rules, to enumerate drug-like compounds under three-dimensional constraints of the target. When applied to bespoke chemical spaces of different sizes and chemical complexity targeting two receptors of pharmaceutical interest (estrogen β receptor, dopamine D3 receptor), the computational method was able to quickly enumerate hits that were either known ligands (or very close analogs) of targeted receptors as well as chemically novel candidates that could be experimentally confirmed by in vitro binding assays. The proposed approach is generic, can be applied to any docking algorithm, and requires few computational resources to prioritize easily synthesizable hits from billion-sized chemical spaces.

Chemical Library Space: Definition and DNA-Encoded Library Comparison Study Case

The development of DNA-encoded library (DEL) technology introduced new challenges for the analysis of chemical libraries. It is often useful to consider a chemical library as a stand-alone chemoinformatic object─represented both as a collection of independent molecules, and yet an individual entity─in particular, when they are inseparable mixtures, like DELs. Herein, we introduce the concept of chemical library space (CLS), in which resident items are individual chemical libraries. We define and compare four vectorial library representations obtained using generative topographic mapping. These allow for an effective comparison of libraries, with the ability to tune and chemically interpret the similarity relationships. In particular, property-tuned CLS encodings enable to simultaneously compare libraries with respect to both property and chemotype distributions. We apply the various CLS encodings for the selection problem of DELs that optimally "match" a reference collection (here ChEMBL28), showing how the choice of the CLS descriptors may help to fine-tune the "matching" (overlap) criteria. Hence, the proposed CLS may represent a new efficient way for polyvalent analysis of thousands of chemical libraries. Selection of an easily accessible compound collection for drug discovery, as a substitute for a difficult to produce reference library, can be tuned for either primary or target-focused screening, also considering property distributions of compounds. Alternatively, selection of libraries covering novel regions of the chemical space with respect to a reference compound subspace may serve for library portfolio enrichment.

Accelerators for Classical Molecular Dynamics Simulations of Biomolecules

Atomistic Molecular Dynamics (MD) simulations provide researchers the ability to model biomolecular structures such as proteins and their interactions with drug-like small molecules with greater spatiotemporal resolution than is otherwise possible using experimental methods. MD simulations are notoriously expensive computational endeavors that have traditionally required massive investment in specialized hardware to access biologically relevant spatiotemporal scales. Our goal is to summarize the fundamental algorithms that are employed in the literature to then highlight the challenges that have affected accelerator implementations in practice. We consider three broad categories of accelerators: Graphics Processing Units (GPUs), Field-Programmable Gate Arrays (FPGAs), and Application Specific Integrated Circuits (ASICs). These categories are comparatively studied to facilitate discussion of their relative trade-offs and to gain context for the current state of the art. We conclude by providing insights into the potential of emerging hardware platforms and algorithms for MD.

A Close-up Look at the Chemical Space of Commercially Available Building Blocks for Medicinal Chemistry

The ability to efficiently synthesize desired compounds can be a limiting factor for chemical space exploration in drug discovery. This ability is conditioned not only by the existence of well-studied synthetic protocols but also by the availability of corresponding reagents, so-called building blocks (BBs). In this work, we present a detailed analysis of the chemical space of 400 000 purchasable BBs. The chemical space was defined by corresponding synthons─fragments contributed to the final molecules upon reaction. They allow an analysis of BB physicochemical properties and diversity, unbiased by the leaving and protective groups in actual reagents. The main classes of BBs were analyzed in terms of their availability, rule-of-two-defined quality, and diversity. Available BBs were eventually compared to a reference set of biologically relevant synthons derived from ChEMBL fragmentation, in order to illustrate how well they cover the actual medicinal chemistry needs. This was performed on a newly constructed universal generative topographic map of synthon chemical space that enables visualization of both libraries and analysis of their overlapped and library-specific regions.

NP Navigator: A New Look at the Natural Product Chemical Space

Natural products (NPs), being evolutionary selected over millions of years to bind to biological macromolecules, remained an important source of inspiration for medicinal chemists even after the advent of efficient drug discovery technologies such as combinatorial chemistry and high-throughput screening. Thus, there is a strong demand for efficient and user-friendly computational tools that allow to analyze large libraries of NPs. In this context, we introduce NP Navigator - a freely available intuitive online tool for visualization and navigation through the chemical space of NPs and NP-like molecules. It is based on the hierarchical ensemble of generative topographic maps, featuring NPs from the COlleCtion of Open NatUral producTs (COCONUT), bioactive compounds from ChEMBL and commercially available molecules from ZINC. NP Navigator allows to efficiently analyze different aspects of NPs - chemotype distribution, physicochemical properties, biological activity and commercial availability of NPs. The latter concerns not only purchasable NPs but also their close analogs that can be considered as synthetic mimetics of NPs or pseudo-NPs.

Targeting Xanthine Oxidase by Natural Products as a Therapeutic Approach for Mental Disorders

Mental disorders comprise diverse human pathologies, including depression, bipolar affective disorder, schizophrenia, and dementia that affect millions of people around the world. The causes of mental disorders are unclear, but growing evidence suggests that oxidative stress and the purine/adenosine system play a key role in their development and progression. Xanthine oxidase (XO) is a flavoprotein enzyme essential for the catalysis of the oxidative hydroxylation of purines -hypoxanthine and xanthine- to generate uric acid. As a consequence of the oxidative reaction of XO, reactive oxygen species (ROS) such as superoxide and hydrogen peroxide are produced and, further, contribute to the pathogenesis of mental disorders. Altered XO activity has been associated with free radical-mediated neurotoxicity inducing cell damage and inflammation. Diverse studies reported a direct association between an increased activity of XO and diverse mental diseases including depression or schizophrenia. Small-molecule inhibitors, such as the well-known allopurinol, and dietary flavonoids, can modulate the XO activity and subsequent ROS production. In the present work, we review the available literature on XO inhibition by small molecules and their potential therapeutic application in mental disorders. In addition, we discuss the chemistry and molecular mechanism of XO inhibitors, as well as the use of structure-based and computational methods to design specific inhibitors with the capability of modulating XO activity.

FragNet, a Contrastive Learning-Based Transformer Model for Clustering, Interpreting, Visualizing, and Navigating Chemical Space

The question of molecular similarity is core in cheminformatics and is usually assessed via a pairwise comparison based on vectors of properties or molecular fingerprints. We recently exploited variational autoencoders to embed 6M molecules in a chemical space, such that their (Euclidean) distance within the latent space so formed could be assessed within the framework of the entire molecular set. However, the standard objective function used did not seek to manipulate the latent space so as to cluster the molecules based on any perceived similarity. Using a set of some 160,000 molecules of biological relevance, we here bring together three modern elements of deep learning to create a novel and disentangled latent space, viz transformers, contrastive learning, and an embedded autoencoder. The effective dimensionality of the latent space was varied such that clear separation of individual types of molecules could be observed within individual dimensions of the latent space. The capacity of the network was such that many dimensions were not populated at all. As before, we assessed the utility of the representation by comparing clozapine with its near neighbors, and we also did the same for various antibiotics related to flucloxacillin. Transformers, especially when as here coupled with contrastive learning, effectively provide one-shot learning and lead to a successful and disentangled representation of molecular latent spaces that at once uses the entire training set in their construction while allowing "similar" molecules to cluster together in an effective and interpretable way.

Chemography: Searching for Hidden Treasures

The days when medicinal chemistry was limited to a few series of compounds of therapeutic interest are long gone. Nowadays, no human may succeed to acquire a complete overview of more than a billion existing or feasible compounds within which the potential "blockbuster drugs" are well hidden and yet only a few mouse clicks away. To reach these "hidden treasures", we adapted the generative topographic mapping method to enable efficient navigation through the chemical space, from a global overview to a structural pattern detection, covering, for the first time, the complete ZINC library of purchasable compounds, relative to 1.6 million biologically relevant ChEMBL molecules. About 40 000 hierarchical maps of the chemical space were constructed. Structural motifs inherent to only one library were identified. Roughly 20 000 off-market ChEMBL compound families represent incentives to enrich commercial catalogs. Alternatively, 125 000 ZINC-specific compound classes, absent in structure-activity bases, are novel paths to explore in medicinal chemistry. The complete list of these chemotypes can be downloaded using the link https://forms.gle/B6bUJj82t9EfmttV6.

Cheminformatics in Natural Product-based Drug Discovery

This review seeks to provide a timely survey of the scope and limitations of cheminformatics methods in natural product-based drug discovery. Following an overview of data resources of chemical, biological and structural information on natural products, we discuss, among other aspects, in silico methods for (i) data curation and natural products dereplication, (ii) analysis, visualization, navigation and comparison of the chemical space, (iii) quantification of natural product-likeness, (iv) prediction of the bioactivities (virtual screening, target prediction), ADME and safety profiles (toxicity) of natural products, (v) natural products-inspired de novo design and (vi) prediction of natural products prone to cause interference with biological assays. Among the many methods discussed are rule-based, similarity-based, shape-based, pharmacophore-based and network-based approaches, docking and machine learning methods.

The Symbiotic Relationship Between Drug Discovery and Organic Chemistry

All pharmaceutical products contain organic molecules; the source may be a natural product or a fully synthetic molecule, or a combination of both. Thus, it follows that organic chemistry underpins both existing and upcoming pharmaceutical products. The reverse relationship has also affected organic synthesis, changing its landscape towards increasingly complex targets. This Review article sets out to give a concise appraisal of this symbiotic relationship between organic chemistry and drug discovery, along with a discussion of the design concepts and highlighting key milestones along the journey. In particular, criteria for a high-quality compound library design enabling efficient virtual navigation of chemical space, as well as rise and fall of concepts for its synthetic exploration (such as combinatorial chemistry; diversity-, biology-, lead-, or fragment-oriented syntheses; and DNA-encoded libraries) are critically surveyed.

Fragment-based drug design of nature-inspired compounds

Natural product (NP)-derived drugs can be extracts, biological macromolecules, or purified small molecule substances. Small molecule drugs can be originally purified from NPs, can represent semisynthetic molecules, natural fragments containing small molecules, or are fully synthetic molecules that mimic natural compounds. New semisynthetic NP-like drugs are entering the pharmaceutical market almost every year and reveal growing interests in the application of fragment-based approaches for NPs. Thus, several NP databases were constructed to be implemented in the fragment-based drug design (FBDD) workflows. FBDD has been established previously as an approach for hit identification and lead generation. Several biophysical and computational methods are used for fragment screening to identify potential hits. Once the fragments within the binding pocket of the protein are identified, they can be grown, linked, or merged to design more active compounds. This work discusses applications of NPs and NP scaffolds to FBDD. Moreover, it briefly reviews NP databases containing fragments and reports on case studies where the approach has been successfully applied for the design of antimalarial and anticancer drug candidates.

SAR by Space: Enriching Hit Sets from the Chemical Space

We introduce SAR by Space, a concept to drastically accelerate structure-activity relationship (SAR) elucidation by synthesizing neighboring compounds that originate from vast chemical spaces. The space navigation is accomplished within minutes on affordable standard computer hardware using a tree-based molecule descriptor and dynamic programming. Maximizing the synthetic accessibility of the results from the computer is achieved by applying a careful selection of building blocks in combination with suitably chosen reactions; a decade of in-house quality control shows that this is a crucial part in the process. The REAL Space is the largest chemical space of commercially available compounds, counting 11 billion molecules as of today. It was used to mine actives against bromodomain 4 (BRD4). Before synthesis, compounds were docked into the binding site using a scoring function, which incorporates intrinsic desolvation terms, thus avoiding time-consuming simulations. Five micromolar hits have been identified and verified within less than six weeks, including the measurement of IC₅₀ values. We conclude that this procedure is a substantial time-saver, accelerating both ligand and structure-based approaches in hit generation and lead optimization stages.

A primer on natural product-based virtual screening

Databases play an important role in various computational techniques, including virtual screening (VS) and molecular modeling in general. These collections of molecules can contain a large amount of information, making them suitable for several drug discovery applications. For example, vendor, bioactivity data or target type can be found when searching a database. The introduction of these data resources and their characteristics is used for the design of an experiment. The description of the construction of a database can also be a good advisor for the creation of a new one. There are free available databases and commercial virtual libraries of molecules. Furthermore, a computational chemist can find databases for a general purpose or a specific subset such as natural products (NPs). In this chapter, NP database resources are presented, along with some guidelines when preparing an NP database for drug discovery purposes.

The next level in chemical space navigation: going far beyond enumerable compound libraries

Recent innovations have brought pharmacophore-driven methods for navigating virtual chemical spaces, the size of which can reach into the billions of molecules, to the fingertips of every chemist. There has been a paradigm shift in the underlying computational chemistry that drives chemical space search applications, incorporating intelligent reaction knowledge into their core so that they can readily deliver commercially available molecules as nearest neighbor hits from within giant virtual spaces. These vast resources enable medicinal chemists to execute rapid scaffold-hopping experiments, rapid hit expansion, and structure-activity relationship (SAR) exploitation in largely intellectual property (IP)-free territory and at unparalleled low cost.

Resources for Chemical, Biological, and Structural Data on Natural Products

Natural products from plants, marine life, animals, fungi, bacteria, and other organisms remain the most productive source of inspiration for small-molecule drug discovery. Today, a wealth of information on natural products that is particularly valuable to applications in cheminformatics is at our disposal. In this contribution, we provide a timely overview of relevant resources for measured chemical, biological, and structural data on natural products. In particular, we comment on the accessibility, scope, chemical space, and limitations of the individual data sources. The bottleneck of natural products remains the limited availability of material for testing. In this context, we analyze the number of natural products readily obtainable from commercial and other sources.

Chemical Space Expansion of Bromodomain Ligands Guided by in Silico Virtual Couplings (AutoCouple)

Expanding the chemical space and simultaneously ensuring synthetic accessibility is of upmost importance, not only for the discovery of effective binders for novel protein classes but, more importantly, for the development of compounds against hard-to-drug proteins. Here, we present AutoCouple, a de novo approach to computational ligand design focused on the diversity-oriented generation of chemical entities via virtual couplings. In a benchmark application, chemically diverse compounds with low-nanomolar potency for the CBP bromodomain and high selectivity against the BRD4(1) bromodomain were achieved by the synthesis of about 50 derivatives of the original fragment. The binding mode was confirmed by X-ray crystallography, target engagement in cells was demonstrated, and antiproliferative activity was showcased in three cancer cell lines. These results reveal AutoCouple as a useful in silico coupling method to expand the chemical space in hit optimization campaigns resulting in potent, selective, and cell permeable bromodomain ligands.

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

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

Predicting novel substrates for enzymes with minimal experimental effort with active learning

Enzymatic substrate promiscuity is more ubiquitous than previously thought, with significant consequences for understanding metabolism and its application to biocatalysis. This realization has given rise to the need for efficient characterization of enzyme promiscuity. Enzyme promiscuity is currently characterized with a limited number of human-selected compounds that may not be representative of the enzyme's versatility. While testing large numbers of compounds may be impractical, computational approaches can exploit existing data to determine the most informative substrates to test next, thereby more thoroughly exploring an enzyme's versatility. To demonstrate this, we used existing studies and tested compounds for four different enzymes, developed support vector machine (SVM) models using these datasets, and selected additional compounds for experiments using an active learning approach. SVMs trained on a chemically diverse set of compounds were discovered to achieve maximum accuracies of ~80% using ~33% fewer compounds than datasets based on all compounds tested in existing studies. Active learning-selected compounds for testing resolved apparent conflicts in the existing training data, while adding diversity to the dataset. The application of these algorithms to wide arrays of metabolic enzymes would result in a library of SVMs that can predict high-probability promiscuous enzymatic reactions and could prove a valuable resource for the design of novel metabolic pathways.

Harnessing the Properties of Natural Products

Natural products (NPs) have been used as traditional medicines since antiquity. With more than 10⁶⁰ estimated compounds with molecular weights less than 500 Da representing chemical space, NPs occupy a very small percentage; however, they are significantly overrepresented in biologically relevant chemical space. The classical approach concentrates on identifying one or more NPs with biological activity from a source organism. There is much more to be learned from NPs than we can discover this narrow view. In this review, we discuss ways to harness the global properties of NPs.

In silico prediction of chemical mechanism of action via an improved network-based inference method

BACKGROUND AND PURPOSE

Deciphering chemical mechanism of action (MoA) enables the development of novel therapeutics (e.g. drug repositioning) and evaluation of drug side effects. Development of novel computational methods for chemical MoA assessment under a systems pharmacology framework would accelerate drug discovery and development with greater efficiency and low cost.

EXPERIMENTAL APPROACH

In this study, we proposed an improved network-based inference method, balanced substructure-drug-target network-based inference (bSDTNBI), to predict MoA for old drugs, clinically failed drugs and new chemical entities. Specifically, three parameters were introduced into network-based resource diffusion processes to adjust the initial resource allocation of different node types, the weighted values of different edge types and the influence of hub nodes. The performance of the method was systematically validated by benchmark datasets and bioassays.

KEY RESULTS

High performance was yielded for bSDTNBI in both 10-fold and leave-one-out cross validations. A global drug-target network was built to explore MoA of anticancer drugs and repurpose old drugs for 15 cancer types/subtypes. In a case study, 27 predicted candidates among 56 commercially available compounds were experimentally validated to have binding affinities on oestrogen receptor α with IC₅₀ or EC₅₀ values ≤10 μM. Furthermore, two dual ligands with both agonistic and antagonistic activities ≤1 μM would provide potential lead compounds for the development of novel targeted therapy in breast cancer or osteoporosis.

CONCLUSION AND IMPLICATIONS

In summary, bSDTNBI would provide a powerful tool for the MoA assessment on both old drugs and novel compounds in drug discovery and development.

Chemoinformatic expedition of the chemical space of fungal products

AIM

Fungi are valuable resources for bioactive secondary metabolites. However, the chemical space of fungal secondary metabolites has been studied only on a limited basis. Herein, we report a comprehensive chemoinformatic analysis of a unique set of 207 fungal metabolites isolated and characterized in a USA National Cancer Institute funded drug discovery project.

RESULTS

Comparison of the molecular complexity of the 207 fungal metabolites with approved anticancer and nonanticancer drugs, compounds in clinical studies, general screening compounds and molecules Generally Recognized as Safe revealed that fungal metabolites have high degree of complexity. Molecular fingerprints showed that fungal metabolites are as structurally diverse as other natural products and have, in general, drug-like physicochemical properties.

CONCLUSION

Fungal products represent promising candidates to expand the medicinally relevant chemical space. This work is a significant expansion of an analysis reported years ago for a smaller set of compounds (less than half of the ones included in the present work) from filamentous fungi using different structural properties.

The Glycome Analytics Platform: an integrative framework for glycobioinformatics

MOTIVATION

Complex carbohydrates play a central role in cellular communication and in disease development. O- and N-glycans, which are post-translationally attached to proteins and lipids, are sugar chains that are rooted, tree structures. Independent efforts to develop computational tools for analyzing complex carbohydrate structures have been designed to exploit specific databases requiring unique formatting and limited transferability. Attempts have been made at integrating these resources, yet it remains difficult to communicate and share data across several online resources. A disadvantage of the lack of coordination between development efforts is the inability of the user community to create reproducible analyses (workflows). The latter results in the more serious unreliability of glycomics metadata.

RESULTS

In this paper, we realize the significance of connecting multiple online glycan resources that can be used to design reproducible experiments for obtaining, generating and analyzing cell glycomes. To address this, a suite of tools and utilities, have been integrated into the analytic functionality of the Galaxy bioinformatics platform to provide a Glycome Analytics Platform (GAP).Using this platform, users can design in silico workflows to manipulate various formats of glycan sequences and analyze glycomes through access to web data and services. We illustrate the central functionality and features of the GAP by way of example; we analyze and compare the features of the N-glycan glycome of monocytic cells sourced from two separate data depositions.This paper highlights the use of reproducible research methods for glycomics analysis and the GAP presents an opportunity for integrating tools in glycobioinformatics.

AVAILABILITY AND IMPLEMENTATION

This software is open-source and available online at https://bitbucket.org/scientificomputing/glycome-analytics-platform

CONTACTS

chris.barnett@uct.ac.za or kevin.naidoo@uct.ac.za

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

4-Acyl Pyrrole Derivatives Yield Novel Vectors for Designing Inhibitors of the Acetyl-Lysine Recognition Site of BRD4(1)

Several human diseases, including cancer, show altered signaling pathways resulting from changes in the activity levels of epigenetic modulators. In the past few years, small-molecule inhibitors against specific modulators, including the bromodomain and extra-terminal (BET) bromodomain family of acetylation readers, have shown early promise in the treatment of the genetically defined midline carcinoma and hematopoietic malignancies. We have recently developed a novel potent inhibitor of BET proteins, 1 (XD14[ Angew. Chem., Int. Ed. 2013, 52, 14055]), which exerts a strong inhibitory potential on the proliferation of specific leukemia cell lines. In the study presented here, we designed analogues of 1 to study the potential of substitutions on the 4-acyl pyrrole backbone to occupy additional sites within the substrate recognition site of BRD4(1). The compounds were profiled using ITC, DSF, and X-ray crystallography. We could introduce several substitutions that address previously untargeted areas of the substrate recognition site. This work may substantially contribute to the development of therapeutics with increased target specificity against BRD4-related malignancies.

A chemical space odyssey of inhibitors of histone deacetylases and bromodomains

The interest in epigenetic drug and probe discovery is growing as reflected in the large amount of structure-epigenetic activity information available.

StreptomeDB 2.0--an extended resource of natural products produced by streptomycetes

Over the last decades, the genus Streptomyces has stirred huge interest in the scientific community as a source of bioactive compounds. The majority of all known antibiotics is isolated from these bacterial strains, as well as a variety of other drugs such as antitumor agents, immunosuppressants and antifungals. To the best of our knowledge, StreptomeDB was the first database focusing on compounds produced by streptomycetes. The new version presented herein represents a major step forward: its content has been increased to over 4000 compounds and more than 2500 host organisms. In addition, we have extended the background information and included hundreds of new manually curated references to literature. The latest update features a unique scaffold-based navigation system, which enables the exploration of the chemical diversity of StreptomeDB on a structural basis. We have included a phylogenetic tree, based on 16S rRNA sequences, which comprises more than two-thirds of the included host organisms. It enables visualizing the frequency, appearance, and persistence of compounds and scaffolds in an evolutionary context. Additionally, we have included predicted MS- and NMR-spectra of thousands of compounds for assignment of experimental data. The database is freely accessible via http://www.pharmaceutical-bioinformatics.org/streptomedb.