Scaffold-Hopping Potential of Fragment-Based De Novo Design: The Chances and Limits of Variation

Comparative study of the mechanism of natural compounds with similar structures using docking and transcriptome data for improving in silico herbal medicine experimentations

Natural products have successfully treated several diseases using a multi-component, multi-target mechanism. However, a precise mechanism of action (MOA) has not been identified. Systems pharmacology methods have been used to overcome these challenges. However, there is a limitation as those similar mechanisms of similar components cannot be identified. In this study, comparisons of physicochemical descriptors, molecular docking analysis and RNA-seq analysis were performed to compare the MOA of similar compounds and to confirm the changes observed when similar compounds were mixed and used. Various analyses have confirmed that compounds with similar structures share similar MOA. We propose an advanced method for in silico experiments in herbal medicine research based on the results. Our study has three novel findings. First, an advanced network pharmacology research method was suggested by partially presenting a solution to the difficulty in identifying multi-component mechanisms. Second, a new natural product analysis method was proposed using large-scale molecular docking analysis. Finally, various biological data and analysis methods were used, such as in silico system pharmacology, docking analysis and drug response RNA-seq. The results of this study are meaningful in that they suggest an analysis strategy that can improve existing systems pharmacology research analysis methods by showing that natural product-derived compounds with the same scaffold have the same mechanism.

Discovery of Novel and Highly Selective Cyclopropane ALK Inhibitors through a Fragment-Assisted, Structure-Based Drug Design

Fragment screening is frequently used for hit identification. However, there was no report starting from a small fragment for the development of an anaplastic lymphoma kinase (ALK) inhibitor, despite the number of ALK inhibitors reported. We began our research with the fragment hit F-1 and our subsequent linker design, and its docking analysis yielded novel cis-1,2,2-trisubstituted cyclopropane 1. The fragment information was integrated with a structure-based approach to improve upon the selectivity over tropomyosin receptor kinase A, leading to the potent and highly selective ALK inhibitor, 4-trifluoromethylphenoxy-cis-1,2,2-trisubstituted cyclopropane 12. This work shows that fragments become a powerful tool for both lead generation and optimization, such as the improvement of selectivity, by combining them with a structure-based drug design approach, resulting in the fast and efficient development of a novel, potent, and highly selective compound.

Unraveling Plant Natural Chemical Diversity for Drug Discovery Purposes

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

Analysis of drug-endogenous human metabolite similarities in terms of their maximum common substructures

In previous work, we have assessed the structural similarities between marketed drugs (‘drugs’) and endogenous natural human metabolites (‘metabolites’ or ‘endogenites’), using ‘fingerprint’ methods in common use, and the Tanimoto and Tversky similarity metrics, finding that the fingerprint encoding used had a dramatic effect on the apparent similarities observed. By contrast, the maximal common substructure (MCS), when the means of determining it is fixed, is a means of determining similarities that is largely independent of the fingerprints, and also has a clear chemical meaning. We here explored the utility of the MCS and metrics derived therefrom. In many cases, a shared scaffold helps cluster drugs and endogenites, and gives insight into enzymes (in particular transporters) that they both share. Tanimoto and Tversky similarities based on the MCS tend to be smaller than those based on the MACCS fingerprint-type encoding, though the converse is also true for a significant fraction of the comparisons. While no single molecular descriptor can account for these differences, a machine learning-based analysis of the nature of the differences (MACCS_Tanimoto vs MCS_Tversky) shows that they are indeed deterministic, although the features that are used in the model to account for this vary greatly with each individual drug. The extent of its utility and interpretability vary with the drug of interest, implying that while MCS is neither ‘better’ nor ‘worse’ for every drug–endogenite comparison, it is sufficiently different to be of value. The overall conclusion is thus that the use of the MCS provides an additional and valuable strategy for understanding the structural basis for similarities between synthetic, marketed drugs and natural intermediary metabolites.

Protocol for fragment hopping

Fragment hopping is a fragment-based approach to designing biologically active small molecules. The key of this approach is the determination of the minimal pharmacophoric elements in the three-dimensional space. Based on the derived minimal pharmacophoric elements, new fragments with different chemotypes can be generated and positioned to the active site of the target protein. Herein, we detail a protocol for performing fragment hopping. This approach can not only explore a wide chemical space to produce new ligands with novel scaffolds but also characterize and utilize the delicate differences in the active sites between isofunctional proteins to produce new ligands with high target selectivity/specificity.

De novo design - hop(p)ing against hope

Current trends in computational de novo design provide a fresh approach to 'scaffold-hopping' in drug discovery. The methodological repertoire is no longer limited to receptor-based methods, but specifically ligand-based techniques that consider multiple properties in parallel, including the synthetic feasibility of the computer-generated molecules and their polypharmacology, provide innovative ideas for the discovery of new chemical entities. The concept of fragment-based and virtual reaction-driven design enables rapid compound optimization from scratch with a manageable complexity of the search. Starting from known drugs as a reference, such algorithms suggest drug-like molecules with motivated scaffold variations, and advanced mathematical models of structure-activity landscapes and multi-objective design techniques have created new opportunities for hit and lead finding.

Using molecular features of xenobiotics to predict hepatic gene expression response

Despite recent advances in molecular medicine and rational drug design, many drugs still fail because toxic effects arise at the cellular and tissue level. In order to better understand these effects, cellular assays can generate high-throughput measurements of gene expression changes induced by small molecules. However, our understanding of how the chemical features of small molecules influence gene expression is very limited. Therefore, we investigated the extent to which chemical features of small molecules can reliably be associated with significant changes in gene expression. Specifically, we analyzed the gene expression response of rat liver cells to 170 different drugs and searched for genes whose expression could be related to chemical features alone. Surprisingly, we can predict the up-regulation of 87 genes (increased expression of at least 1.5 times compared to controls). We show an average cross-validation predictive area under the receiver operating characteristic curve (AUROC) of 0.7 or greater for each of these 87 genes. We applied our method to an external data set of rat liver gene expression response to a novel drug and achieved an AUROC of 0.7. We also validated our approach by predicting up-regulation of Cytochrome P450 1A2 (CYP1A2) in three drugs known to induce CYP1A2 that were not in our data set. Finally, a detailed analysis of the CYP1A2 predictor allowed us to identify which fragments made significant contributions to the predictive scores.

Enrichment of virtual hits by progressive shape-matching and docking

The main applications of virtual chemical screening include the selection of a minimal receptor-relevant subset of a chemical library with a maximal chemical diversity. We have previously reported that the combination of ligand-centric and receptor-centric virtual screening methods may provide a compromise between computational time and accuracy during the hit enrichment process. In the present work, we propose a "progressive distributed docking" method that improves the virtual screening process using an iterative combination of shape-matching and docking steps. Known ligands with low docking scores were used as initial 3D templates for the shape comparisons with the chemical library. Next, new compounds with good template shape matches and low receptor docking scores were selected for the next round of shape searching and docking. The present iterative virtual screening process was tested for enriching peroxisome proliferator-activated receptor and phosphoinositide 3-kinase relevant compounds from a selected subset of the chemical libraries. It was demonstrated that the iterative combination improved the lead-hopping practice by improving the chemical diversity in the selected list of virtual hits.

Structural diversity of biologically interesting datasets: a scaffold analysis approach

Background

The recent public availability of the human metabolome and natural product datasets has revitalized "metabolite-likeness" and "natural product-likeness" as a drug design concept to design lead libraries targeting specific pathways. Many reports have analyzed the physicochemical property space of biologically important datasets, with only a few comprehensively characterizing the scaffold diversity in public datasets of biological interest. With large collections of high quality public data currently available, we carried out a comparative analysis of current day leads with other biologically relevant datasets.

Results

In this study, we note a two-fold enrichment of metabolite scaffolds in drug dataset (42%) as compared to currently used lead libraries (23%). We also note that only a small percentage (5%) of natural product scaffolds space is shared by the lead dataset. We have identified specific scaffolds that are present in metabolites and natural products, with close counterparts in the drugs, but are missing in the lead dataset. To determine the distribution of compounds in physicochemical property space we analyzed the molecular polar surface area, the molecular solubility, the number of rings and the number of rotatable bonds in addition to four well-known Lipinski properties. Here, we note that, with only few exceptions, most of the drugs follow Lipinski's rule. The average values of the molecular polar surface area and the molecular solubility in metabolites is the highest while the number of rings is the lowest. In addition, we note that natural products contain the maximum number of rings and the rotatable bonds than any other dataset under consideration.

Conclusions

Currently used lead libraries make little use of the metabolites and natural products scaffold space. We believe that metabolites and natural products are recognized by at least one protein in the biosphere therefore, sampling the fragment and scaffold space of these compounds, along with the knowledge of distribution in physicochemical property space, can result in better lead libraries. Hence, we recommend the greater use of metabolites and natural products while designing lead libraries. Nevertheless, metabolites have a limited distribution in chemical space that limits the usage of metabolites in library design.

The role of computational methods in the identification of bioactive compounds

Computational methods play an ever increasing role in lead finding. A vast repertoire of molecular design and virtual screening methods emerged in the past two decades and are today routinely used. There is increasing awareness that there is no single best computational protocol and correspondingly there is a shift recommending the combination of complementary methods. A promising trend for the application of computational methods in lead finding is to take advantage of the vast amounts of HTS (High Throughput Screening) data to allow lead assessment by detailed systems-based data analysis, especially for phenotypic screens where the identification of compound-target pairs is the primary goal. Herein, we review trends and provide examples of successful applications of computational methods in lead finding.

De novo design: balancing novelty and confined chemical space

IMPORTANCE OF THE FIELD

De novo drug design serves as a tool for the discovery of new ligands for macromolecular targets as well as optimization of known ligands. Recently developed tools aim to address the multi-objective nature of drug design in an unprecedented manner.

AREAS COVERED IN THIS REVIEW

This article discusses recent advances in de novo drug design programs and accessory programs used to evaluate compounds post-generation.

WHAT THE READER WILL GAIN

The reader is introduced to the challenges inherent in de novo drug design and will become familiar with current trends in de novo design. Furthermore, the reader will be better prepared to assess the value of a tool, and be equipped to design more elegant tools in the future.

TAKE HOME MESSAGE

De novo drug design can assist in the efficient discovery of new compounds with a high affinity for a given target. The inclusion of existing chemoinformatic methods with current structure-based de novo design tools provides a means of enhancing the therapeutic value of these generated compounds.

Bioisosteric Replacement and Scaffold Hopping in Lead Generation and Optimization

Bioisosteric replacement and scaffold hopping are twin methods used in drug design to improve the synthetic accessibility, potency and drug like properties of a compound and to move into novel chemical space. Bioisosteric replacement involves swapping functional groups of a molecule with other functional groups that have similar biological properties. Scaffold hopping is the replacement of the core framework of a molecule with another scaffold that will improve the properties of the molecule or to find similar potent compounds that exist in novel chemical space. This review outlines the key concepts, importance and challenges of both methods using examples and comparisons of techniques available for finding bioisosteric replacements and scaffold hops. There are many methods available for bioisosteric replacement and scaffold hopping, all with their own advantages and disadvantages. Drug design projects would benefit from a combination of these methods to retrieve diverse and complimentary results. Continuing progress in these fields will allow further validation of both methods as well as the accumulation of knowledge on bioisosteres and possible scaffold replacements.