Drug-likeness analysis of traditional Chinese medicines: 2. Characterization of scaffold architectures for drug-like compounds, non-drug-like compounds, and natural compounds from traditional Chinese medicines

Recent advances in the development of P2Y14R inhibitors: a patent and literature review (2018-present)

INTRODUCTION

The P2Y14 receptor (P2Y14R), a member of the G protein-coupled receptor family, is activated by extracellular nucleotides. Due to its involvement in inflammatory, immunological and other associated processes, P2Y14R has emerged as a promising therapeutic target. Despite lacking a determined three-dimensional crystal structure, the homology modeling technique based on closely related P2Y receptors' crystallography has been extensively utilized for developing active compounds targeting P2Y14R. Recent discoveries have unveiled numerous highly effective and subtype-specific P2Y14R inhibitors. This study presents an overview of the latest advancements in P2Y14R inhibitors.

AREAS COVERED

This review presents an overview of the advancements in P2Y14R inhibitor research over the past five years, encompassing new patents, journal articles, and highlighting the therapeutic prospects inherent in these compounds.

EXPERT OPINION

The recent revelation of the vast potential of P2Y14R inhibitors has led to the development of novel compounds that exhibit promising capabilities for the treatment of sterile inflammation of the kidney, potentially diabetes, and asthma. Despite being a relatively nascent class of compounds, certain members have already exhibited their capacity to surmount specific challenges posed by conventional P2Y14R inhibitors. Targeting P2Y14R through small molecules may present a promising therapeutic strategy for effectively managing diverse inflammatory diseases.

Scopy: an integrated negative design python library for desirable HTS/VS database design

BACKGROUND

High-throughput screening (HTS) and virtual screening (VS) have been widely used to identify potential hits from large chemical libraries. However, the frequent occurrence of 'noisy compounds' in the screened libraries, such as compounds with poor drug-likeness, poor selectivity or potential toxicity, has greatly weakened the enrichment capability of HTS and VS campaigns. Therefore, the development of comprehensive and credible tools to detect noisy compounds from chemical libraries is urgently needed in early stages of drug discovery.

RESULTS

In this study, we developed a freely available integrated python library for negative design, called Scopy, which supports the functions of data preparation, calculation of descriptors, scaffolds and screening filters, and data visualization. The current version of Scopy can calculate 39 basic molecular properties, 3 comprehensive molecular evaluation scores, 2 types of molecular scaffolds, 6 types of substructure descriptors and 2 types of fingerprints. A number of important screening rules are also provided by Scopy, including 15 drug-likeness rules (13 drug-likeness rules and 2 building block rules), 8 frequent hitter rules (four assay interference substructure filters and four promiscuous compound substructure filters), and 11 toxicophore filters (five human-related toxicity substructure filters, three environment-related toxicity substructure filters and three comprehensive toxicity substructure filters). Moreover, this library supports four different visualization functions to help users to gain a better understanding of the screened data, including basic feature radar chart, feature-feature-related scatter diagram, functional group marker gram and cloud gram.

CONCLUSION

Scopy provides a comprehensive Python package to filter out compounds with undesirable properties or substructures, which will benefit the design of high-quality chemical libraries for drug design and discovery. It is freely available at https://github.com/kotori-y/Scopy.

Plant Metabolite Databases: From Herbal Medicines to Modern Drug Discovery

Traditional herbal medicine has been an inseparable part of the traditional medical science in many countries throughout history. Nowadays, the popularity of using herbal medicines in daily life, as well as clinical practices, has gradually expanded to numerous Western countries with positive impacts and acceptance. The continuous growth of the herbal consumption market has promoted standardization and modernization of herbal-derived products with present pharmacological criteria. To store and extensively share this knowledge with the community and serve scientific research, various herbal metabolite databases have been developed with diverse focuses under the support of modern advances. The advent of these databases has contributed to accelerating research on pharmaceuticals of natural origins. In the scope of this study, we critically review 30 herbal metabolite databases, discuss different related perspectives, and provide a comparative analysis of 18 accessible noncommercial ones. We hope to provide you with fundamental information and multidimensional perspectives from herbal medicines to modern drug discovery.

Multi-Target Drugs Against Metabolic Disorders

BACKGROUND

Metabolic disorders are a major cause of illness and death worldwide. Metabolism is the process by which the body makes energy from proteins, carbohydrates, and fats; chemically breaking these down in the digestive system towards sugars and acids which constitute the human body's fuel for immediate use, or to store in body tissues, such as the liver, muscles, and body fat.

OBJECTIVE

The efficiency of treatments for multifactor diseases has not been proved. It is accepted that to manage multifactor diseases, simultaneous modulation of multiple targets is required leading to the development of new strategies for discovery and development of drugs against metabolic disorders.

METHODS

In silico studies are increasingly being applied by researchers due to reductions in time and costs for new prototype synthesis; obtaining substances that present better therapeutic profiles.

DISCUSSION

In the present work, in addition to discussing multi-target drug discovery and the contributions of in silico studies to rational bioactive planning against metabolic disorders such as diabetes and obesity, we review various in silico study contributions to the fight against human metabolic pathologies.

CONCLUSION

In this review, we have presented various studies involved in the treatment of metabolic disorders; attempting to obtain hybrid molecules with pharmacological activity against various targets and expanding biological activity by using different mechanisms of action to treat a single pathology.

Cheminformatic Insight into the Differences between Terrestrial and Marine Originated Natural Products

This is a new golden age for drug discovery based on natural products derived from both marine and terrestrial sources. Herein, a straightforward but important question is "what are the major structural differences between marine natural products (MNPs) and terrestrial natural products (TNPs)?" To answer this question, we analyzed the important physicochemical properties, structural features, and drug-likeness of the two types of natural products and discussed their differences from the perspective of evolution. In general, MNPs have lower solubility and are often larger than TNPs. On average, particularly from the perspective of unique fragments and scaffolds, MNPs usually possess more long chains and large rings, especially 8- to 10-membered rings. MNPs also have more nitrogen atoms and halogens, notably bromines, and fewer oxygen atoms, suggesting that MNPs may be synthesized by more diverse biosynthetic pathways than TNPs. Analysis of the frequently occurring Murcko frameworks in MNPs and TNPS also reveals a striking difference between MNPs and TNPs. The scaffolds of the former tend to be longer and often contain ester bonds connected to 10-membered rings, while the scaffolds of the latter tend to be shorter and often bear more stable ring systems and bond types. Besides, the prediction from the naïve Bayesian drug-likeness classification model suggests that most compounds in MNPs and TNPs are drug-like, although MNPs are slightly more drug-like than TNPs. We believe that MNPs and TNPs with novel drug-like scaffolds have great potential to be drug leads or drug candidates in drug discovery campaigns.

HawkRank: a new scoring function for protein-protein docking based on weighted energy terms

Deciphering the structural determinants of protein–protein interactions (PPIs) is essential to gain a deep understanding of many important biological functions in the living cells. Computational approaches for the structural modeling of PPIs, such as protein–protein docking, are quite needed to complement existing experimental techniques. The reliability of a protein–protein docking method is dependent on the ability of the scoring function to accurately distinguish the near-native binding structures from a huge number of decoys. In this study, we developed HawkRank, a novel scoring function designed for the sampling stage of protein–protein docking by summing the contributions from several energy terms, including van der Waals potentials, electrostatic potentials and desolvation potentials. First, based on the solvation free energies predicted by the Generalized Born model for ~ 800 proteins, a SASA (solvent accessible surface area)-based solvation model was developed, which can give the aqueous solvation free energies for proteins by summing the contributions of 21 atom types. Then, the van der Waals potentials and electrostatic potentials based on the Amber ff14SB force field were computed. Finally, the HawkRank scoring function was derived by determining the most optimal weights for five energy terms based on the training set. Here, MSR (modified success rate), a novel protein–protein scoring quality index, was used to assess the performance of HawkRank and three other popular protein–protein scoring functions, including ZRANK, FireDock and dDFIRE. The results show that HawkRank outperformed the other three scoring functions according to the total number of hits and MSR. HawkRank is available at http://cadd.zju.edu.cn/programs/hawkrank.

ADMET Evaluation in Drug Discovery. Part 17: Development of Quantitative and Qualitative Prediction Models for Chemical-Induced Respiratory Toxicity

As a dangerous end point, respiratory toxicity can cause serious adverse health effects and even death. Meanwhile, it is a common and traditional issue in occupational and environmental protection. Pharmaceutical and chemical industries have a strong urge to develop precise and convenient computational tools to evaluate the respiratory toxicity of compounds as early as possible. Most of the reported theoretical models were developed based on the respiratory toxicity data sets with one single symptom, such as respiratory sensitization, and therefore these models may not afford reliable predictions for toxic compounds with other respiratory symptoms, such as pneumonia or rhinitis. Here, based on a diverse data set of mouse intraperitoneal respiratory toxicity characterized by multiple symptoms, a number of quantitative and qualitative predictions models with high reliability were developed by machine learning approaches. First, a four-tier dimension reduction strategy was employed to find an optimal set of 20 molecular descriptors for model building. Then, six machine learning approaches were used to develop the prediction models, including relevance vector machine (RVM), support vector machine (SVM), regularized random forest (RRF), extreme gradient boosting (XGBoost), naïve Bayes (NB), and linear discriminant analysis (LDA). Among all of the models, the SVM regression model shows the most accurate quantitative predictions for the test set (q²_ext = 0.707), and the XGBoost classification model achieves the most accurate qualitative predictions for the test set (MCC of 0.644, AUC of 0.893, and global accuracy of 82.62%). The application domains were analyzed, and all of the tested compounds fall within the application domain coverage. We also examined the structural features of the compounds and important fragments with large prediction errors. In conclusion, the SVM regression model and the XGBoost classification model can be employed as accurate prediction tools for respiratory toxicity.

Chinese Herbal Medicine Meets Biological Networks of Complex Diseases: A Computational Perspective

With the rapid development of cheminformatics, computational biology, and systems biology, great progress has been made recently in the computational research of Chinese herbal medicine with in-depth understanding towards pharmacognosy. This paper summarized these studies in the aspects of computational methods, traditional Chinese medicine (TCM) compound databases, and TCM network pharmacology. Furthermore, we chose arachidonic acid metabolic network as a case study to demonstrate the regulatory function of herbal medicine in the treatment of inflammation at network level. Finally, a computational workflow for the network-based TCM study, derived from our previous successful applications, was proposed.

Comparative analyses of structural features and scaffold diversity for purchasable compound libraries

Large purchasable screening libraries of small molecules afforded by commercial vendors are indispensable sources for virtual screening (VS). Selecting an optimal screening library for a specific VS campaign is quite important to improve the success rates and avoid wasting resources in later experimental phases. Analysis of the structural features and molecular diversity for different screening libraries can provide valuable information to the decision making process when selecting screening libraries for VS. In this study, the structural features and scaffold diversity of eleven purchasable screening libraries and Traditional Chinese Medicine Compound Database (TCMCD) were analyzed and compared. Their scaffold diversity represented by the Murcko frameworks and Level 1 scaffolds was characterized by the scaffold counts and cumulative scaffold frequency plots, and visualized by Tree Maps and SAR Maps. The analysis demonstrates that, based on the standardized subsets with similar molecular weight distributions, Chembridge, ChemicalBlock, Mucle, TCMCD and VitasM are more structurally diverse than the others. Compared with all purchasable screening libraries, TCMCD has the highest structural complexity indeed but more conservative molecular scaffolds. Moreover, we found that some representative scaffolds were important components of drug candidates against different drug targets, such as kinases and guanosine-binding protein coupled receptors, and therefore the molecules containing pharmacologically important scaffolds found in screening libraries might be potential inhibitors against the relevant targets. This study may provide valuable perspective on which purchasable compound libraries are better for you to screen.Graphical abstract

Drug-likeness prediction of chemical constituents isolated from Chinese materia medica Ciwujia

ETHNOPHARMACOLOGICAL RELEVANCE

Ciwujia (CWJ), one of the most commonly used Chinese materia medicas (CMMs), is derived from the roots, rhizomes, and stems of Acanthopanax senticosus harms (AS). CWJ has been used for the treatment of various central nervous system (CNS) and peripheral system diseases. Drug-likeness prediction can help to analyze the absorption, distribution, metabolism, and excretion (ADME) processes of the compounds in CWJ, as well as their potential therapeutic and toxic effects, which is of significance in the confirmation of the active material bases of CWJ.

MATERIALS AND METHODS

The ADME properties of the compounds were calculated through web based PreADMET program and ACD/I-Lab 2.0. The potential therapeutic and toxicity targets of these compounds were screened by the ChemQuery tool in DrugBank and T3DB.

RESULTS

14/39 compounds had moderate or good oral bioavailability (OB). 29/39 compounds bound weakly to the plasma proteins. 18/39 compounds might pass across the blood-brain barrier (BBB). Most of these compounds showed low renal excretion ability. 25/39 compounds had 99 structurally similar drugs and 158 potential therapeutic targets. Additionally, 17/39 compounds had 53 structurally similar toxins and 126 potential toxicity targets.

CONCLUSION

Our study suggests that these compounds have a certain drug-likeness potentials, which are also likely to be the material bases of CWJ. These results may provide a reference for the safe use of CWJ and the expansion of its application scope.

ADMET evaluation in drug discovery: 15. Accurate prediction of rat oral acute toxicity using relevance vector machine and consensus modeling

Background

Determination of acute toxicity, expressed as median lethal dose (LD₅₀), is one of the most important steps in drug discovery pipeline. Because in vivo assays for oral acute toxicity in mammals are time-consuming and costly, there is thus an urgent need to develop in silico prediction models of oral acute toxicity.

Results

In this study, based on a comprehensive data set containing 7314 diverse chemicals with rat oral LD₅₀ values, relevance vector machine (RVM) technique was employed to build the regression models for the prediction of oral acute toxicity in rate, which were compared with those built using other six machine learning approaches, including k-nearest-neighbor regression, random forest (RF), support vector machine, local approximate Gaussian process, multilayer perceptron ensemble, and eXtreme gradient boosting. A subset of the original molecular descriptors and structural fingerprints (PubChem or SubFP) was chosen by the Chi squared statistics. The prediction capabilities of individual QSAR models, measured by q_ext² for the test set containing 2376 molecules, ranged from 0.572 to 0.659.

Conclusion

Considering the overall prediction accuracy for the test set, RVM with Laplacian kernel and RF were recommended to build in silico models with better predictivity for rat oral acute toxicity. By combining the predictions from individual models, four consensus models were developed, yielding better prediction capabilities for the test set (q_ext² = 0.669–0.689). Finally, some essential descriptors and substructures relevant to oral acute toxicity were identified and analyzed, and they may be served as property or substructure alerts to avoid toxicity. We believe that the best consensus model with high prediction accuracy can be used as a reliable virtual screening tool to filter out compounds with high rat oral acute toxicity. Graphical abstract

Workflow of combinatorial QSAR modelling to predict rat oral acute toxicity

Collapse

Exploration of Scaffolds from Natural Products with Antiplasmodial Activities, Currently Registered Antimalarial Drugs and Public Malarial Screen Data

In light of current resistance to antimalarial drugs, there is a need to discover new classes of antimalarial agents with unique mechanisms of action. Identification of unique scaffolds from natural products with in vitro antiplasmodial activities may be the starting point for such new classes of antimalarial agents. We therefore conducted scaffold diversity and comparison analysis of natural products with in vitro antiplasmodial activities (NAA), currently registered antimalarial drugs (CRAD) and malaria screen data from Medicine for Malaria Ventures (MMV). The scaffold diversity analyses on the three datasets were performed using scaffold counts and cumulative scaffold frequency plots. Scaffolds from the NAA were compared to those from CRAD and MMV. A Scaffold Tree was also generated for each of the datasets and the scaffold diversity of NAA was found to be higher than that of MMV. Among the NAA compounds, we identified unique scaffolds that were not contained in any of the other compound datasets. These scaffolds from NAA also possess desirable drug-like properties making them ideal starting points for antimalarial drug design considerations. The Scaffold Tree showed the preponderance of ring systems in NAA and identified virtual scaffolds, which may be potential bioactive compounds.

Structure-Activity Relationships and Anti-inflammatory Activities of N-Carbamothioylformamide Analogues as MIF Tautomerase Inhibitors

Macrophage migration inhibitory factor (MIF), a proinflammatory cytokine, is an attractive therapeutic target for the treatment of inflammatory diseases. In our previous study, 3-[(biphenyl-4-ylcarbonyl)carbamothioyl]amino benzoic acid (compound 1) was discovered as a potent inhibitor of MIF by docking-based virtual screening and bioassays. Here, a series of analogues of compound 1 derived from similarity search and chemical synthesis were evaluated for their MIF tautomerase activities, and their structure-activity relationships were then analyzed. The most potent inhibitor (compound 5) with an IC50 of 370 nM strongly suppressed lipopolysaccharide (LPS)-induced production of TNF-α and IL-6 in a dose-dependent manner and significantly enhanced the survival rate of mice with LPS-induced endotoxic shock from 0 to 35% at 0.5 mg/kg and to 45% at 1 mg/kg, highlighting the therapeutic potential of the MIF tautomerase inhibition in inflammatory diseases.

The application of in silico drug-likeness predictions in pharmaceutical research

The concept of drug-likeness, established from the analyses of the physiochemical properties or/and structural features of existing small organic drugs or/and drug candidates, has been widely used to filter out compounds with undesirable properties, especially poor ADMET (absorption, distribution, metabolism, excretion, and toxicity) profiles. Here, we summarize various approaches for drug-likeness evaluations, including simple rules/filters based on molecular properties/structures and quantitative prediction models based on sophisticated machine learning methods, and provide a comprehensive review of recent advances in this field. Moreover, the strengths and weaknesses of these approaches are briefly outlined. Finally, the drug-likeness analyses of natural products and traditional Chinese medicines (TCM) are discussed.

Capturing nature's diversity

Natural products are universally recognized to contribute valuable chemical diversity to the design of molecular screening libraries. The analysis undertaken in this work, provides a foundation for the generation of fragment screening libraries that capture the diverse range of molecular recognition building blocks embedded within natural products. Physicochemical properties were used to select fragment-sized natural products from a database of known natural products (Dictionary of Natural Products). PCA analysis was used to illustrate the positioning of the fragment subset within the property space of the non-fragment sized natural products in the dataset. Structural diversity was analysed by three distinct methods: atom function analysis, using pharmacophore fingerprints, atom type analysis, using radial fingerprints, and scaffold analysis. Small pharmacophore triplets, representing the range of chemical features present in natural products that are capable of engaging in molecular interactions with small, contiguous areas of protein binding surfaces, were analysed. We demonstrate that fragment-sized natural products capture more than half of the small pharmacophore triplet diversity observed in non fragment-sized natural product datasets. Atom type analysis using radial fingerprints was represented by a self-organizing map. We examined the structural diversity of non-flat fragment-sized natural product scaffolds, rich in sp3 configured centres. From these results we demonstrate that 2-ring fragment-sized natural products effectively balance the opposing characteristics of minimal complexity and broad structural diversity when compared to the larger, more complex fragment-like natural products. These naturally-derived fragments could be used as the starting point for the generation of a highly diverse library with the scope for further medicinal chemistry elaboration due to their minimal structural complexity. This study highlights the possibility to capture a high proportion of the individual molecular interaction motifs embedded within natural products using a fragment screening library spanning 422 structural clusters and comprised of approximately 2800 natural products.

Drug/nondrug classification using Support Vector Machines with various feature selection strategies

In conjunction with the advance in computer technology, virtual screening of small molecules has been started to use in drug discovery. Since there are thousands of compounds in early-phase of drug discovery, a fast classification method, which can distinguish between active and inactive molecules, can be used for screening large compound collections. In this study, we used Support Vector Machines (SVM) for this type of classification task. SVM is a powerful classification tool that is becoming increasingly popular in various machine-learning applications. The data sets consist of 631 compounds for training set and 216 compounds for a separate test set. In data pre-processing step, the Pearson's correlation coefficient used as a filter to eliminate redundant features. After application of the correlation filter, a single SVM has been applied to this reduced data set. Moreover, we have investigated the performance of SVM with different feature selection strategies, including SVM-Recursive Feature Elimination, Wrapper Method and Subset Selection. All feature selection methods generally represent better performance than a single SVM while Subset Selection outperforms other feature selection methods. We have tested SVM as a classification tool in a real-life drug discovery problem and our results revealed that it could be a useful method for classification task in early-phase of drug discovery.

Virtual screening of traditional Chinese medicine (TCM) database: identification of fragment-like lead molecules for filariasis target asparaginyl-tRNA synthetase

Lymphatic filariasis (LF) is a vector borne infectious disease caused by the nematode Wuchereria bancrofti, Brugia malayi, and Brugia timori. Over 120 million people are affected by LF in the world, of which two-thirds are in Asia. The infection restricts the normal flow of lymph from the infected area resulting in swelling of the extremities and causing permanent disability. As the available drugs for the treatment of LF are becoming ineffective due to the development of resistance, there is an urgent need to find new leads for drug development. In this study, asparaginyl-tRNA synthetase (AsnRS; PDB ID: 2XGT) essential for the protein bio-synthesis in the filarial nematode was used to carry out virtual screening (VS) of plant constituents from traditional Chinese medicine (TCM) database. Docking as well as E-pharmacophore based VS were carried out to identify the hits. The top scoring hits, Agri 1 (1,3,8-trihydroxy-4,5-dimethoxyxanthen-9-one-3-O-beta-D-glucopyranoside) and Agri 2 (5,7-dihydroxy-2-propylchromone 7-O-beta-D-glucopyranoside), constituents of Agrimonia pilosa, were selected for molecular dynamics (MD) simulation study for 10 ns. MD simulation showed that both the glycosides Agri 1 and Agri 2 were forming stable interactions with the target protein. Moreover, docking and MD simulation of the lead A (1,3,8-trihydroxy-4,5-dimethoxyxanthen-9-one; Mol. Wt.: 304.25; CLogP: 3.07) and lead B (5,7-dihydroxy-2-propylchromone; Mol. Wt.: 220.22; CLogP: 3.02), the aglycones of Agri 1 and Agri 2, respectively, were carried out with the target AsnRS. The in silico investigations of the aglycones suggest that the lead B could be a suitable fragment-like lead molecule for anti-filarial drug discovery.

A Drug-Target Network-Based Approach to Evaluate the Efficacy of Medicinal Plants for Type II Diabetes Mellitus

The use of plants as natural medicines in the treatment of type II diabetes mellitus (T2DM) has long been of special interest. In this work, we developed a docking score-weighted prediction model based on drug-target network to evaluate the efficacy of medicinal plants for T2DM. High throughput virtual screening from chemical library of natural products was adopted to calculate the binding affinity between natural products contained in medicinal plants and 33 T2DM-related proteins. The drug-target network was constructed according to the strength of the binding affinity if the molecular docking score satisfied the threshold. By linking the medicinal plant with T2DM through drug-target network, the model can predict the efficacy of natural products and medicinal plant for T2DM. Eighteen thousand nine hundred ninety-nine natural products and 1669 medicinal plants were predicted to be potentially bioactive.