Next Article in Journal
Bufotalin Suppresses Proliferation and Metastasis of Triple-Negative Breast Cancer Cells by Promoting Apoptosis and Inhibiting the STAT3/EMT Axis
Previous Article in Journal
Enzyme-Assisted Extraction of Ulvan from the Green Macroalgae Ulva fenestrata
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Discovering the Active Ingredients of Medicine and Food Homologous Substances for Inhibiting the Cyclooxygenase-2 Metabolic Pathway by Machine Learning Algorithms

1
State Key Laboratory of Chemical Resource Engineering, Department of Pharmaceutical Engineering, Beijing University of Chemical Technology, 15 Beisanhuan East Road, Beijing 100029, China
2
Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, China
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(19), 6782; https://doi.org/10.3390/molecules28196782
Submission received: 10 July 2023 / Revised: 20 September 2023 / Accepted: 21 September 2023 / Published: 23 September 2023
(This article belongs to the Section Computational and Theoretical Chemistry)

Abstract

:
Cyclooxygenase-2 (COX-2) and microsomal prostaglandin E2 synthase (mPGES-1) are two key targets in anti-inflammatory therapy. Medicine and food homology (MFH) substances have both edible and medicinal properties, providing a valuable resource for the development of novel, safe, and efficient COX-2 and mPGES-1 inhibitors. In this study, we collected active ingredients from 503 MFH substances and constructed the first comprehensive MFH database containing 27,319 molecules. Subsequently, we performed Murcko scaffold analysis and K-means clustering to deeply analyze the composition of the constructed database and evaluate its structural diversity. Furthermore, we employed four supervised machine learning algorithms, including support vector machine (SVM), random forest (RF), deep neural networks (DNNs), and eXtreme Gradient Boosting (XGBoost), as well as ensemble learning, to establish 640 classification models and 160 regression models for COX-2 and mPGES-1 inhibitors. Among them, ModelA_ensemble_RF_1 emerged as the optimal classification model for COX-2 inhibitors, achieving predicted Matthews correlation coefficient (MCC) values of 0.802 and 0.603 on the test set and external validation set, respectively. ModelC_RDKIT_SVM_2 was identified as the best regression model based on COX-2 inhibitors, with root mean squared error (RMSE) values of 0.419 and 0.513 on the test set and external validation set, respectively. ModelD_ECFP_SVM_4 stood out as the top classification model for mPGES-1 inhibitors, attaining MCC values of 0.832 and 0.584 on the test set and external validation set, respectively. The optimal regression model for mPGES-1 inhibitors, ModelF_3D_SVM_1, exhibited predictive RMSE values of 0.253 and 0.35 on the test set and external validation set, respectively. Finally, we proposed a ligand-based cascade virtual screening strategy, which integrated the well-performing supervised machine learning models with unsupervised learning: the self-organized map (SOM) and molecular scaffold analysis. Using this virtual screening workflow, we discovered 10 potential COX-2 inhibitors and 15 potential mPGES-1 inhibitors from the MFH database. We further verified candidates by molecular docking, investigated the interaction of the candidate molecules upon binding to COX-2 or mPGES-1. The constructed comprehensive MFH database has laid a solid foundation for the further research and utilization of the MFH substances. The series of well-performing machine learning models can be employed to predict the COX-2 and mPGES-1 inhibitory capabilities of unknown compounds, thereby aiding in the discovery of anti-inflammatory medications. The COX-2 and mPGES-1 potential inhibitor molecules identified through the cascade virtual screening approach provide insights and references for the design of highly effective and safe novel anti-inflammatory drugs.

1. Introduction

Eicosanoids, derived from arachidonic acid (AA) and related polyunsaturated fatty acids (PUFAs), play a crucial role in regulating a wide range of homeostatic and inflammatory processes associated with many diseases, such as atherosclerosis, Alzheimer’s disease, and cancer [1]. Studies of eicosanoids have primarily focused on prostaglandin E2 (PGE2), due to the fact that PEG2 primarily contributes to the fundamental symptoms of inflammation, such as fever, swelling, redness, pain, and loss of function [2]. In the production of PGE2, cyclooxygenase-2 (COX-2) catalyzes the conversion of arachidonic acid (AA) into prostaglandin H2 (PGH2). Subsequently, microsomal prostaglandin E2 synthase (mPGES-1) converts PGH2 into PGE2. Therefore, suppressing the overexpression of COX-2 and mPGES-1 is a reasonable strategy in the treatment of inflammation.
Traditional non-steroidal anti-inflammatory drugs (NSAIDs), such as Aspirin, inhibit both COX-1 and COX-2. This reduces the production of the pro-inflammatory factor PEG2 along with the cytoprotective prostaglandin I2, leading to an increased risk of severe gastrointestinal and cardiovascular disease. Although highly selective COX2 inhibitors, such as Celecoxib and Rofecoxib, can alleviate some of the side effects associated with NSAIDs, their long-term use does not completely avoid the risk of cardiovascular disease [3]. Consequently, current research is focused on discovering potential COX-2 inhibitors that exhibit both high selectivity and low toxicity. In recent years, mPGES-1 inhibitors have been suggested as attractive anti-inflammation therapeutics since they can selectively reduce the PEG2 without affecting the cytoprotective PGs that regulate homeostasis. However, the wide variation in mPGES-1 inhibitor efficacy between humans and mice poses a challenge to fully exploit rodent disease models, and the lack of effective inhibitors with cross-species activity further hinders the development of safe and effective mPGES-1 inhibitors [4].
Medicine and food homology (MFH) refers to substances that are both food and herbal medicine. It reflects the traditional Chinese idea of health care, and includes the contents of food therapy, health care, and medicinal food in Chinese medicine [5]. MFHs have passed the food safety risk assessment of the Chinese National Health Commission and have been proven to be both for long-term use as daily food and for disease treatment and health benefits [6]. Natural products have served as a valuable source of drugs and lead compounds due to their complexity, structural diversity, and extensive pharmacological effects. Growing studies have revealed that many active ingredients of plant-derived and marine-derived natural products are involved in the regulation of inflammatory responses in humans [7], especially with COX-2 and mPGES-1 inhibitory activity. A series of secondary metabolites isolated from marine products have shown potency in inhibiting COX-2, their chemical structures are shown in sub-figure A of Figure 1. Vaccinal A, Botryoisocoumarin A [8], Axinelline A, Capnellene, Stachybogrisephenone B, and Actinoquinoline A inhibit COX-2 in vitro with IC50 values of 1.8 μM, 6.5 μM, 2.8 μM, 6.2 μM, 8.9 μM, and 2.1 μM, respectively [9]. Some compounds from medicinal plants have been verified to exhibit the prominent mPGES-1 inhibitory effect (seen in sub-figure B of Figure 1), Curcumin, Hyperforin, Garcinol, Arzanol, Myrtucommulone, Carnosic acid, and Emeblin inhibit mPGES-1 in vitro with IC50 values of 0.3 μM, 1 μM, 0.3 μM, 0.4 μM, 1 μM, 5 μM, and 0.2 μM, respectively [10]. Although many natural products are reported to be potent COX-2, mPGES-1 inhibitors or transcriptional suppressors, poor oral availability and non-negligible toxicity still limit their further clinical use. MFHs are easily absorbed and utilized by humans with relative nontoxicity. Additionally, there are some multitarget anti-inflammatory compounds in Chinese herbal medicine. Acrovestone, extracted from acronychia, inhibits mPGES-1 and 5-lipoxygenase (5-LOX) in vitro with IC50 values of 2.7 μM and 1.1 μM, respectively [11]. Cannflavin A, extracted from Cannabis sativa, inhibits mPGES-1 and 5-LOX in vitro with IC50 values of 1.8 μM and 0.9 μM, respectively. The above multi-target inhibitors can reduce the risk of drug–drug interactions and contribute to efficient anti-inflammatory therapy [12]. Therefore, MFHs are capable of providing insight into the discovery of novel, safe, and effective COX-2 and mPGES-1 inhibitors.
There has been a growing interest in the application of machine learning (ML) and deep learning (DL) in the field of biomedicine. Various emerging supervised and unsupervised learning algorithms have shown potential in providing valuable insights and predictions for drug discovery, drug repurposing, diagnostics, and pharmaceutical production [13]. Currently, virtual screening (VS) combined with experimental verification is the standard protocol for drug discovery. Some ML-facilitating VS approaches have been widely used to mine novel molecular substances from in-house or commercial compound libraries [14]. An ML-based Quantitative structure–activity relationship (QSAR) model has become a favorable tool in VS due to its ability to rapidly output predictions based on input datasets and its high hit rate [15]. Recent studies have focused on implementing VS on natural compound libraries; this is because active ingredients derived from natural products are considered an ideal starting point for designing anti-inflammatory agents [16]. Wang’s group [17] built several classification models and employed these models for identifying potential P-glycoprotein inhibitors in the Traditional Chinese Medicine Systems Pharmacology (TCMSP) database. Sattar et al. [18] utilized molecular docking to screen potential COX inhibitors on compounds isolated from Eucalyptus maculata resin, and they found that 1,6-dicinnamoyl-O-α-D-glucopyranoside exhibited an COX-2 inhibitory effect.
Currently, obscure functions and unclear therapeutic targets limit the clinical practice of medicine and food homology (MFH) substances. In this study, we constructed the first comprehensive MFH database aiming at further investigating the complex active ingredients of MFH and deciphering its biology and functional activities (detailed workflow is shown in Figure 2). Additionally, we constructed the updated comprehensive dataset of COX-2 and mPGES-1 inhibitors. Leveraging these datasets, we employed supervised machine learning algorithms and ensemble learning techniques to develop a series of QSAR models for predicting the inhibitory efficacy of COX-2 and mPGES-1. Furthermore, utilizing the well-performed machine learning models combined with unsupervised learning and scaffold analysis, we performed virtual screening on the MFH database to identify potential COX-2 and mPGES-1 inhibitors which could contribute to the efficient and safe anti-inflammatory therapy.

2. Results and Discussion

2.1. Chemical Space and Scaffold Analysis of the MFH Database

We constructed a database containing 27,319 active ingredient molecules derived from 503 kinds of MFH substances by collecting data from the Traditional Chinese Medicines Integrated Database (TCMID) [19] and Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP) [20], as well as 92 pieces of literature. Details of the constructed MFH database can be seen in Table S1. The distribution of two physicochemical properties of the MFH database is shown in Figure 3. The molecular weight (MW) of molecules derived from TCMID, TCMSP, and the literature ranges from 32 to 1900, 30 to 1466, and 59 to 1263, respectively. The octanol–water partition coefficient (LogP) of molecules derived from TCMID, TCMSP, and the literature ranges from −24.9 to 24, −11.9 to 24.5, and −8.8 to 12.6, respectively. The distribution of basic physicochemical properties across two orders of magnitude demonstrates the extensive chemical space of our MFH database.
To further evaluate the structural diversity of our MFH database, we calculated the Tanimoto coefficients (TCs) [21] on molecules (represented by 1024 bits ECFP4 fingerprints). The average TC value for the entire MFH database was 0.216, with 95.92% of molecule pairs having TC values of less than 0.6, indicating significant structural differences within our MFH database.
In addition, we extracted Murcko scaffolds of the MFH database and performed a K-means clustering [22] analysis. As a result, the MFH database was divided into 11 clusters. The Flavonoids cluster was the most abundant, comprising approximately one-fifth of the database, while Fatty Acids, Saponins, and Steroids each accounted for roughly one-tenth. Moreover, there were Lignans, Alkaloids, Triterpenoids, Sesquiterpenes, Diterpenes, and Stilbenes in the MFH database (seen in Figure 4A). These are commonly isolated components from natural products, and some Flavonoids, Steroids, Alkaloids, and Triterpenoids have been reported to play important roles in the regulation of immune responses [23]. We further summarized the top 20 Murcko scaffolds in the MFH database, as shown in Figure 4B; the scaffolds ranged from simple aromatic natural products with a single ring to complex skeletons with 7–8 ring systems. These observations demonstrate that our MFH database is comprehensive and exhibits high structural diversity.

2.2. Performances Evaluation and Comparison of Developed Models

In this study, we developed classification models using four datasets: Dataset 1 and 2, which contain COX-2 inhibitors, and Dataset 4 and 5, which are composed of mPGES-1 inhibitors. To generate training and test sets, we split the datasets for classification 10 times randomly. We then characterized these sets using three types of fingerprints: Avalon, ECFP4, and MACCS. Four machine learning methods were applied to predict high or weak inhibition on COX-2 and mPGES-1 separately. As a result, we constructed a total of 480 classification models. Subsequently, we integrated the predicted probabilities of classification models (built with the same algorithm but with the different molecular characterization) by a stacked generalization approach. This resulted in the construction of an additional 160 ensemble classification models: 80 for COX-2 inhibitors and 80 for mPGES-1 inhibitors. Additionally, we utilized two datasets, Dataset 3 and Dataset 6, represented by two types of descriptors to build 160 regression models following a similar process of developing classification models.

2.2.1. Performances of Classification and Ensemble Models on Dataset 1

The Matthews correlation coefficients (MCCs) on the test sets were employed to evaluate the predictive ability and stability of constructed classification models, and the MCC on external validation sets was applied to assess the generalization ability of models. Table 1 lists the overall model performances based on Dataset 1 (1640 COX-2 inhibitors), and the detailed performances of all models on Dataset 1 can be seen in Table S1. Figure 5A depicts the MCC values of 10 randomly split test sets for Dataset 1.
When comparing the performance of different algorithms, SVM and XGB slightly outperformed other machine learning methods. When using Avalon fingerprints to characterize the dataset, the mean MCC values on 10 randomly split test sets for SVM and XGB models were 0.783 and 0.744, respectively. When using ECFP fingerprints, the mean MCC values for the SVM and XGBoost models were 0.763 and 0.768, respectively. When using MACCS fingerprints, the mean MCC values for the SVM and XGBoost models were 0.746 and 0.753, respectively. Meanwhile, the RF algorithm had a slightly inferior performance with the average MCC values on 10 randomly split test sets of 0.721, 0.704, and 0.701 when using Avalon, ECFP, and MACCS fingerprints to characterize the dataset, respectively.
When comparing the performance of different descriptors, all three fingerprints (Avalon, ECFP, and MACCS) were equally effective in characterizing the dataset with no significant differences. The MCC values on test sets for all models exceeded 0.66. Compared to the base classification models, the ensemble models that integrated all three fingerprints showed improved performance with mean MCC values of 0.735, 0.782, 0.747, and 0.774 on test sets when using RF, SVM, DNN, and XGBoost models, respectively. The ensemble model for each algorithm had an average MCC improvement of 0.02 on the test set compared to the models constructed with a single type of fingerprints.

2.2.2. Performances of Classification and Ensemble Models on Dataset 2

The overall performances of developed classifiers on Dataset 2 were summarized in Table 2, the detailed performances of all models on Dataset 2 (containing 2925 COX-2 inhibitors) can be seen in Table S2. The detailed MCC values of 10 randomly split test sets for Dataset 2 were visualized in Figure 5B. The MCC values of the test sets of the models on dataset 2 all exceeded 0.52, which was a decrease in performance compared to the classification models constructed with Dataset 1. This was mainly because the moderately active inhibitors were removed from dataset 1, and the structural differences between the highly/weakly inhibitors were easily distinguished. Similar to the performances on Dataset 1, the SVM and XGBoost algorithms also exhibited excellent performance for Dataset 2, with MCC values of 0.632 ± 0.018, 0.63 ± 0.022, and 0.621 ± 0.025 for the SVM model on the test set of the dataset characterized using Avalon, ECFP, and MACCS, respectively. The MCC values of the XGBoost model on the test set of the dataset represented by Avalon, ECFP, and MACCS were 0.644 ± 0.014, 0.625 ± 0.028, and 0.626 ± 0.012, respectively. The MCC performances of the ensemble RF, SVM, DNN, and XGBoost models on the test sets were 0.576 ± 0.02, 0.642 ± 0.016, 0.604 ± 0.015, and 0.643 ± 0.012, respectively. Compared with models constructed with one type of fingerprints, the ensemble models of each algorithm increased the average MCC on the test set by 0.01.

2.2.3. Performances of Regression Models on Dataset 3

Table 3 displays the overall performances of the models using two criteria (R2 and RMSE) based on Dataset 3, which consisted of 1511 COX-2 inhibitors. Detailed results of all 80 QSAR models are listed in Table S3. The RMSE based on the test sets and external validation sets were utilized as the main metric to evaluate the performance of the regression models. Figure 6 showed the specific RMSE values of 10 randomly split test sets on Dataset 3.
The RMSE values for the test set of models constructed using the G_3D descriptors were all below 0.56. From the perspective of modeling algorithms, the SVM algorithm performed the best with an average RMSE of 0.465 on the test sets. Among the models constructed using RDKit descriptors, the SVM algorithm also showed the highest predictive power with a mean RMSE of 0.436 ± 0.019, followed by the comparable DNN and XGBoost algorithms with mean RMSE values of 0.471 ± 0.019, respectively. Comparing from the descriptor perspective, the RMSEs of the test sets of the models constructed using RDKit descriptors were all below 0.54. Based on the same algorithm, RDKit-based models performed slightly better than the G_3D-based models, which indicates that the RDKit 2D descriptor is more suitable for characterizing the COX-2 inhibitors collected in this study.

2.2.4. Performances of Classification and Ensemble Models on Dataset 4

Table 4 summarizes the overall model performances of Dataset 4, which included 3179 mPGES-1 inhibitors. The detailed performances of all models on Dataset 4 can be seen in Table S4 and Figure 7A. The XGBoost algorithm performed the best, with mean MCC values of 0.74, 0.83, and 0.765 for 10 randomly partitioned test sets based on the Avalon, ECFP4, and MACCS characterized datasets, respectively. Comparing from the perspective of dataset characterization, ECFP4 fingerprints slightly outperformed the other two fingerprints. The highest MCC values were obtained on the test sets of the SVM, DNN, and XGBoost models using the dataset characterized by ECFP4 fingerprints. Among the ensemble classification models constructed based on the three fingerprints, the ensemble SVM and XGBoost models achieved excellent MCC performances of 0.783 ± 0.017 and 0.788 ± 0.006 on the test sets, respectively.

2.2.5. Performances of Classification and Ensemble Models on Dataset 5

Compared to the classification models from Dataset 4 (Table 4), the classification models constructed based on Dataset 5, containing 3455 mPGES-1 inhibitors, were slightly less discriminative for highly/weakly active mPGES-1 inhibitors, shown in Table 5, Table S5, and Figure 7B. In terms of machine learning algorithms, the models established by the RF algorithm had the worst performance, with all the average MCC values below 0.6 on the test set based on the datasets characterized by three types of fingerprints. Conversely, the SVM and XGBoost algorithms exhibited better performance. XGBoost models achieved the best performance on the test set based on the datasets characterized by Avalon and MACCS. On the test set based on the ECFP4 characterization, the SVM and XGBoost models achieved excellent performance with average MCC values of 0.773 and 0.772, respectively. In addition, the ensemble SVM and XGBoost classification models obtained average MCC values of 0.738 ± 0.017 and 0.749 ± 0.011 on the test sets, respectively.

2.2.6. Performances of Regression Models on Dataset 6

Based on Dataset 6 consisting of 735 mPGES-1 inhibitors, a total of 80 regression models (Table 6) were constructed using the four algorithms to predict the bioactivities of the inhibitors; the details of all model performances are listed in Table S6. As seen in Table 6 and Figure 8, the RMSE values of the test sets of the models constructed using the G_3D descriptors were all below 0.42, with the SVM algorithm performing the best, with an average RMSE of 0.329 on the test sets. This was followed by models constructed by the XGBoost algorithm, with an average RMSE of 0.339 on the test sets. The RMSE values of the test sets of the models constructed using the 2D RDKit descriptors were all below 0.35. From the descriptors’ perspective, RDKit-based models performed slightly better than G_3D-based models for the same algorithm, with an average RMSE reduction of 0.03 in the test sets.

2.2.7. Performances of the External Validation Sets

In addition to evaluating the predictive accuracy and stability of the models through the performance of the five-fold cross-validation and test sets, four external validation sets were also employed to verify the generalization ability of the constructed classification and regression models. Only models that performed well on the external validation sets were applied to further implement virtual screening on the MFH database. The ensemble SVM models and the ensemble XGBoost models performed best in the classification models constructed based on Dataset 1, with predicted MCC values of 0.551 ± 0.012 and 0.546 ± 0.014 for external validation set A1, respectively. Similarly, the ensemble SVM and the ensemble XGBoost models outperformed others in the classification models developed based on Dataset 2, with predicted MCC values of 0.544 ± 0.011 and 0.531 ± 0.01, respectively. For the regression model built on Dataset 3, the regression models constructed using the XGBoost algorithm and the RDkit descriptors had the most accurate prediction RMSE values of 0.584 ± 0.035 for external validation set A2. In the classification model constructed based on Dataset 4, models constructed by the XGBoost with the ECFP4 fingerprints had the highest prediction for external validation set B1 with an MCC of 0.547 ± 0.015. For classification models established based on Dataset 5, the model developed using SVM and XGBoost algorithms with ECFP4 fingerprints had the best performance; the predicted MCC values for external validation set B1 reached 0.538 ± 0.016 and 0.537 ± 0.009, respectively. Models constructed through the SVM algorithm combined with RDkit descriptors outperformed the regression models built based on Dataset 6, with a predicted RMSE of 0.424 ± 0.019 for external validation set B2. The specific model results and parameters applied for virtual screening of MFH database are summarized in Table S7.

2.3. Virtual Screening on the MFH Database

In this study, we employed a ligand-based cascade virtual screening strategy to predict potential COX-2 and mPGES-1 inhibitors from the constructed MFH database containing 27,319 molecules. Our goal was to advance the development of safer anti-inflammatory therapies. The virtual screening workflow proceeded as follows, taking the screening of COX-2 inhibitors as an example. Firstly, molecules of the MFH database were predicted using 39 classification models that exhibited excellent performance on external validation sets (seen in Table S8). These 39 classification models comprised both single classifiers and ensemble classifiers constructed based on Dataset 1, as well as basic classifiers and ensemble classifiers constructed based on Dataset 2. By combining different datasets with multiple fingerprints and employing various supervised machine learning algorithms and the ensemble algorithm, the classification models used to predict the COX-2 inhibitory activity of molecules in the MFH database demonstrated strong generalization ability and robustness, thus providing reliability to the prediction results to some extent. Subsequently, the molecules predicted as highly active inhibitors by the 39 classification models were further input into 17 regression models to predict their specific inhibitory values against COX-2 (Detailed performances are listed in Table S8). The predicted values from the 17 regression models were averaged to mitigate occasional errors introduced by individual models, further enhancing the credibility of the prediction results. Based on the ranked mean predicted bioactivities from the 17 regression models, molecules with an average predicted IC50 values below 10 μM (a common threshold for inhibitory potency) were retained. Additionally, we applied the self-organized map (SOM), an unsupervised algorithm, to predict the molecules of the MFH database. The positioning of molecules on a two-dimensional neural network was used to determine their high or weak COX-2 inhibitory activity. This approach enabled the retention of MFH molecules that resided in the same location as the majority of highly active molecules (over 80%). Finally, COX-2 inhibitor candidates were selected through an extensive literature review and molecular scaffold analysis. The screening process for mPGES-1 inhibitors followed a similar procedure, as described above.

2.3.1. Potential COX-2 Inhibitors in the MFH Database

Through the cascade virtual screening described above, 10 potential COX-2 inhibitors were selected from a pool of 27,319 MFH molecules. The structures of these molecules, along with their positions on the two-dimensional grid mapped by SOM, were presented in Figure 9. From Figure 8, it can be observed that the majority of the MFH molecules are distributed in different grid cells compared to known COX-2 inhibitor molecules. This discrepancy arises due to significant structural differences between natural MFH molecules and the synthesized or modified COX-2 inhibitors. Our focus primarily lies on the MFH molecules that share the same position as reported natural COX-2 inhibitors and those located in the grids predominantly occupied by highly active inhibitors (rendered in red in Figure 9). HP represents the proportion of the grid occupied by highly active inhibitors. It is the proportion of the number of highly active inhibitors in the grid to all molecules in the grid. The larger the calculated HP of one grid (close to 1), the warmer the color tends to be (red). On the contrary, the closer the HP value of a grid is to 0, the cooler the color tends to be (blue). The origin, predicted bioactivity values, and reported effects of the candidate molecules are summarized in Table 7.
Candidate cmp_A3, also known as Humulene, is derived from Panax Ginseng and belongs to the monocyclic sesquiterpene. It has been reported to exhibit inhibitory activity against COX-2, suppressing the expression of COX-2 in mice and reducing the production of prostaglandin E2 (PGE2) [24]. Another candidate, cmp_A1 (Dehydrotanshinone II A), derived from Radix Salviae, is the benzofuran derivatives. Kwon et al. isolated Dehydrotanshinone II A from Salvia miltiorrhiza Bunge and validated its inhibitory effect on COX-2 using a platelet activation model [25]. Candidate cmp_A6 (β-Sesquiphellandrene) is commonly found in Atractylodes macrocephala and Zingiber officinale. Zingiber officinale oil has been reported to possess anti-inflammatory activity, particularly by inhibiting lipoxygenase [26]. cmp_A7, derived from Lycii Fructus, is an apocarotenoid that significantly inhibits the expression of IL-1β in vitro, thus exerting an anti-inflammatory effect [27]. Cmp_A5, named Scutianine C, is derived from in Jujubae Fructus and has been reported as a biologically active alkaloid with antimicrobial properties [28]. Candidate cmp_A10 (Hispaglabridin B) isolated from Glycyrrhiza glabra L. belongs to the isoflavone derivatives and is located in the same grid as the reported natural COX-2 inhibitors, which has demonstrated antioxidant effects [29]. Peptidomimetic candidates, cmp_A2 and cmp_A8, originate from Corneum Gigeriae Galli Endothelium and Fagopyrum esculentum, respectively. They were predicted by SOM to have a probability of 0.821 as high-activity COX-2 inhibitors, with average predicted IC50 values of 2.69 μM and 5.2 μM, respectively. Candidate cmp_A4, classified as a sesquiterpene, is derived from Angelica sinensis Radix and was predicted by SOM as a 100% high-activity COX-2 inhibitor, with an average predicted IC50 value of 3.69 μM. cmp_A9, found in Mori Follum, is a folinic acid with an average predicted IC50 value of 5.92 μM.
Table 7. The COX-2 inhibitors candidates from the MFH database found by virtual screening using the related ML models built in this work.
Table 7. The COX-2 inhibitors candidates from the MFH database found by virtual screening using the related ML models built in this work.
CandidatesOriginsIC50 μM aSOM HP bEffects
cmp_A1Radix Salviae1.74 0.91Anti-inflammation [25]
cmp_A2Gigeriae Galli Endothelium2.69 0.82
cmp_A3Panax Ginseng3.47 1COX-2 inhibition [24]
cmp_A4Angelica sinensis Radix3.69 1
cmp_A5Jujubae Fructus4.24 1Antibacterial [28]
cmp_A6Atractylodes macrocephala4.37 0.91Anti-inflammation [26]
cmp_A7Lycii Fructus4.60 0.91Anti-inflammation [27]
cmp_A8Fagopyrum esculentum5.20 0.82
cmp_A9Mori Follum5.92 0.82
cmp_A10Glycyrrhiza glabra L.6.97 1Anti-oxidation [29]
a: the average predicted IC50 values by a series of optimal regression models. b: the percentage of highly active inhibitors mapped to the target molecules in the same position during the self-organized map (SOM) training. It can also be taken as the probability of being predicted by an unsupervised algorithm as a highly active COX-2 inhibitor.

2.3.2. Potential mPGES-1 Inhibitors in the MFH Database

Through the virtual screening process described above, 15 potential mPGES-1 inhibitors were screened from the MFH database. As seen in Figure 10 and Table 8, Candidate cmp_B1, a cannabinolic acid derived from Cannabis Sativa L., was mapped by SOM into a grid inhabited by all highly active mPGES-1 inhibitors. The average predicted IC50 value from the regression models for cmp_B1 was 0.88 μM. Candidate cmp_B1 has been reported to decrease COX enzyme activity, although selectivity toward COX1/2 was still unknown [30]. Candidate cmp_B2, isolated from Ramulus Mori, is a prenylated flavanone. It was predicted by the SOM to be a highly active mPGES-1 inhibitor with a probability of 0.93 and a mean IC50 of 0.25 μM by the regression models. Candidate cmp_B3, an active ingredient of Amomum longiligularg, is an isopentenyl flavonoid with a predicted IC50 of 0.34 μM, and it has been reported to inhibit the growth of breast cancer cells in vitro and in vivo [31]. Candidate cmp_B4 is also known as Kanzonol C, which is a flavonoid-like active ingredient of Glycyrrhiza glabra L. It has been demonstrated as a PTP1B inhibitor [32] and has been found to exhibit inhibitory activity against nitric oxide (NO), making it a potential anti-inflammatory agent [33]. Candidate cmp_B5, derived from Gardeniae Fructus, is a caffeoylquinic acid with a predicted IC50 value of 0.18 μM. It possesses the ability to inhibit lipoxygenase [34]. Candidate cmp_B6 is the active ingredient of Schisandra chinensis with an average predicted IC50 value of 0.37 μM, which has been reported to inhibit UDP-glucuronosyltransferase [35] and Oxidized low-density lipoprotein (OxLDL) [36]. Candidate cmp_B7, also known as Garcinone B, derived from Rhizoma Dioscoreae (Chinese yam), has a predicted IC50 value of 0.55 μM. It has been found to reduce the production of prostaglandin E2, although the specific metabolic pathway of its action has yet to be demonstrated [37]. Candidate cmp_B8 (Kuwanon M) was isolated from Ramulus Mori. It was predicted by the SOM to be a highly active mPGES-1 inhibitor with a probability of 1 and was predicted by the regression models to have a mean IC50 of 0.18 μM. Candidate cmp_B9, derived from Mori Cortex, is an isopropenylated phenol derivative that has been reported to show an inhibitory effect on Tyrosinase [38]. Candidate cmp_B10, derived from Glehniae Radix, is a coumarin analogue with a predicted IC50 value of 0.2 μM. Candidate cmp_B11, also known as Xanthochymol, is a component of Colla and belongs to the polycyclic phloroglucinol. It has been found to regulate inflammation by downregulating the expression of several major histocompatibility complex (MHC) molecules [39]. Candidate cmp_B12, derived from Epimrdii Herba, is a flavonoid with a predicted IC50 value of 0.5 μM. It has been reported to inhibit CYP3A4, thus exhibiting anti-inflammatory effects [40]. Candidate cmp_B13, isolated from Coicis Semen, is a steroid with an average predicted IC50 value of 0.27 μM. Lee et al. [41] verified the inhibitory effect of γ-oryzanol, composed of cmp_B13 and various ferulate molecules, on inflammation-related diseases. Candidate cmp_B14, also known as Kaikasaponin III, derived from Radix Puerariae, is a triterpenoid with a predicted IC50 value of 0.2 μM. It possesses antioxidant effects [42] and has been shown to have therapeutic effects on colitis [43]. Candidate cmp_B15, located in the same grid as cmp_B14, is also a triterpenoid primarily found in Alisma Orientale. It has a predicted IC50 value of 0.29 μM and has been reported to inhibit COX-2 expression [44].

2.3.3. Molecular Docking on the Potential COX-2 and mPGES-1 Inhibitors

Through a series of ligand-based virtual screening processes, we have identified potential inhibitors of COX-2 and mPGES-1 from the established MFH database, and we further filtered these candidate molecules by the pan assay interference compounds (PAINS) rule [45]. To further validate the reliability and validity of the virtual screening, we conducted molecular docking (a widely utilized structure-based virtual screening approach) on the candidate MFH molecules. This procedure aimed to examine the binding modes between the candidate molecules and the proteins, while also assessing the ability of the candidate molecules to interact with the key amino acids of the target protein.

Molecular Docking Analysis on Potential COX-2 Inhibitors

The active site of COX-2 is demarcated from the initial substrate binding site by a constriction formed by three residues: Arg120, Tyr355, and Glu524. This structural constriction necessitates dilation to enable the entry or exit of substrates to and from the active site [46]. Ser530 and Tyr385 are also vital during the catalytic process of COX-2. Several typical binding modes of non-steroidal anti-inflammatory drugs (NSAIDs) interacting with COX-2 have been reported. These modes include the hydrophobic region of NSAIDs interacting with Tyr385, Trp387, and neighboring residues; the polar region of NSAIDs interacting with residues located above Tyr355; NSAIDs with negative charges interacting with Arg120 [47].
We calculated the structural similarity (measured by Tanimoto coefficient) between the candidate MFH molecules and the ligand molecules in the published crystal structures, we selected the protein crystal structure with the bound ligand having the highest structural similarity to the candidate MFH molecules for docking, the details of which are shown in Table S10. As shown in Table 9, candidate compounds cmp_A3, cmp_A4, cmp_A6, and cmp_A7 bind within the active cavity of COX-2 (PDB: 4PH9), with binding affinities of −7.37, −8.53, −8.68, and −7.69 kcal/mol, respectively. The co-crystallized ligand, ibuprofen, in the complex 4PH9, under the same docking conditions, exhibits a binding affinity of −9.41 kcal/mol. While the four candidate molecules displayed affinities lower than the original ligand within the protein, they all established interactions with key amino acid residues in the catalytic domain of COX-2: hydrophobic interactions with Trp388 and polar interactions involving Tyr356, Ser531, and Tyr386. Candidate cmp_A7 generated a hydrogen bond interaction with Arg121, indicating a robust binding force. Candidate cmp_A1 bound within the active pocket of COX-2 (PDB: 5KIR), exhibiting a binding affinity of −7.57 kcal/mol. Meanwhile, the co-crystallized ligand within complex 5KIR achieved a binding affinity of −9.80 kcal/mol under identical docking conditions. Candidate cmp_A3 engaged in polar interactions with Tyr355 and formed a π-H stacking interaction with Ser533. Candidate cmp_A2 and cmp_A10 occupied the active site of the protein with PDB index 6BL3, demonstrating binding affinities of −11.75 and −8.92 kcal/mol, respectively. The co-crystallized ligand within complex 6BL3 had a binding affinity of −12.15 kcal/mol under the same docking procedure. Candidate cmp_A2 generated hydrogen bond interactions with Ser530, Lys83, and Glu524, engaged in hydrophobic interactions with Trp387, and established polar interactions with Tyr355 and Tyr385. Candidate cmp_A10 formed hydrogen bond interactions with Glu524, had hydrophobic interactions with Trp387, and engaged in polar interactions with Ser530 and Tyr355. Candidates cmp_A5, cmp_A8, and cmp_A9 bound to the active cavity of COX-2 (PDB: 6BL4) with the affinities of −7.25, −11.24, and −9.61 kcal/mol, respectively. The affinity of the complex 6BL4 after docking with the original ligand was −12.27 kcal/mol. Candidate cmp_A5 engaged in hydrophobic interactions with Trp100 and formed polar interactions with Tyr115, Lys79, and Lys83. Candidate cmp_A8 established hydrogen bond interactions with Tyr355, Arg120, Glu524, Ser530, and Lys83, which significantly contributed to its enhanced protein affinity. Candidate cmp_A9 had hydrogen bond interactions with Ser119, Met522, Glu524, and Ser530.
In summary, the candidate MFH molecules could effectively bind to the active pocket of COX-2 and interact with key amino acids involved in COX-2 catalysis. These interactions were consistent with the classical binding interactions reported between NSAIDs and COX-2. Most candidate MFH molecules exhibited strong hydrogen bonding interactions with the key amino acid residues of COX-2. These observations support the reliability of the potential COX-2 inhibitor molecules discovered through the ligand-based virtual screening on the MFH database.

Molecular Docking Analysis on Potential mPGES-1 Inhibitors

mPGES-1 is a homotrimer, with each subunit consisting of four transmembrane helices. The mPGES-1 trimer contains three active site cavities, which are formed collectively by transmembrane helices 1, 2, and 4 along with neighboring monomers [48]. Among the currently resolved co-crystal complexes, key amino acids include AlaA123, ProA124, SerA127, ValA128, TyrA130, ThrA131, GlnA134, TyrB28, IleB32, AsnB36, ArgB38, LeuB39, PheB44, ArgB52, and HisB53, where A and B representing different monomers in the mPGES-1 trimer, respectively.
The selection of mPGES-1 crystals for docking was carried out using a similar methodology as the COX-2 crystal selection, with the detailed results presented in Table S10. The binding affinities and interactions between candidate MFH molecules and mPGES-1 are summarized in Table 10. Candidate cmp_B4 and cmp_B5 bound within the active pocket of the protein 4AL1, exhibiting binding affinities of −8.08 and −11.25 kcal/mol, respectively. The binding affinity of the original ligand of protein 4AL1 is −12.12 kcal/mol. Specifically, Candidate cmp_B4 formed hydrogen bond interactions with AspB49, AsnB46, and ThrA131. Candidate cmp_B5 established hydrogen bond interactions with AsnB46, ArgB52, GluA77, and TyrA117, with ionic bonds between ArgB38 and ArgA126, and π-H stacking with Ile B32. Candidate cmp_B1, cmp_B6, cmp_B7, cmp_B11, cmp_B12, cmp_B13, cmp_B14, and cmp_B15 were docked within the active site of the protein 4YL0, exhibiting binding affinities of −7.93, −7.25, −7.26, −6.89, −8.73, −7.19, −10.71, and −7.34 kcal/mol, respectively. The binding affinity of protein 4YL0 after docking with its ligand is −7.17 kcal/mol, and except for cmp_B11, the docking affinities of the remaining candidate molecules surpass that of the original ligand of protein 4YL0 under the same docking conditions. Notably, candidate cmp_B14 demonstrated significantly improved docking affinity compared to the other candidates. Candidate cmp_B1 had hydrogen bond interactions with ArgB38 and AsnB46. Hydrogen bond interactions were established between Candidate cmp_B6 and AspB49. Candidate cmp_B7 generated hydrogen bond interactions with GluA77 and SerA127, while also forming a π-π stacking interaction with TyrA130. It is noteworthy that Candidate cmp_B7 has been reported to inhibit the generation of PGE2 [37]. Coupled with its hydrogen bond interactions with key amino acids at the active site of mPGES-1, these findings underscore the potential of cmp_B7 as an inhibitor for mPGES-1. Candidate cmp_B2 and cmp_B3 bound to the active pocket of protein 4YL1 with binding affinities of −9.19 and −8.82 kcal/mol, both of which were superior to the docking affinity of the original ligand of protein 4YL1 (−8.71 kcal/mol). Candidate cmp_B2 formed hydrogen bond interactions with AspB49 and ThrA131. Candidate cmp_B3 generated a hydrogen bond interaction with AspB49. Candidate cmp_B10 docked within the active domain of protein 5K0I, displaying a binding affinity of −10.31 kcal/mol, which outperformed the binding affinity of the original ligand of protein 5K0I under identical docking conditions (−8.57 kcal/mol). Candidate cmp_B10 engaged in hydrogen bond interactions with AspB49, SerA127, GluA77, and TyrA117. Candidates cmp_B8 and cmp_B9 bound within the active cavity of mPGES-1 (PDB: 5TL9), with binding affinities of −9.55 and −9.2 kcal/mol, respectively. The binding affinity of the original ligand within complex 5TL9 was −8.63 kcal/mol, indicating that cmp_B8 and cmp_B9 exhibited superior binding performances under the same docking parameters. Candidate cmp_B8 formed a hydrogen bond interaction with SerA127. Candidate cmp_B9 had hydrogen bond interactions with SerA127, GluA77, and TyrA117.
The majority of candidate molecules exhibited the ability to form hydrogen bond interactions with key amino acids of mPGES-1. Notably, the docking affinities of several candidate molecules even surpass those of the original ligands within the co-crystal complexes. These findings collectively contribute to bolstering the plausibility of the potential mPGES-1 inhibitors identified through our virtual screening process.

3. Materials and Methods

3.1. Construction of the Catalogue for MFH Substances

The culture of traditional Chinese medicine has a long history, the “Compendium of Materia Medica”, compiled in the 16th century, records more than 300 kinds of medicine and food homology (MFH) substances. In 2002, Chinese National Health Commission issued a list of “Chinese medicine that can be used in health food” (including 114 kinds of substances), both as traditional Chinese medicines and as food [49]. In 2018, Chinese National Health Commission published a list of 110 MFH substances [50]. Additionally, the 2020 edition of the “Pharmacopoeia” records 86 substances as “healthy food” and 193 substances as “therapeutic food”, both of which are in line with the concept of medicine and food homology [51]. Therefore, by examining and integrating the above four lists, a total of 503 MFH substances were collected in the catalogue of this study, the active ingredients were further collected based on this catalog with the aim of constructing a comprehensive database of MFH substances.

3.2. Collection and Preparation of Active Ingredients from MFH Substances

Based on the integrated and constructed catalogue of MFH substances, we collected the active ingredients of 503 MFH substances from the following public Chinese Medicine database: Traditional Chinese Medicines Integrated Database (TCMID) [19], and Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP) [20]. However, some of the MFH substances were not included in the above databases, or the collected active ingredients were inadequate. Therefore, we also conducted a literature search (92 pieces of Chinese literature from the China National Knowledge Infrastructure (CNKI) [52]) and manually collected the active ingredients of MFH substances to ensure the integrity and diversity of our MFH substances database. For the collected active ingredient molecules, we further disconnected group metals in simple salts, removed minor components, screened the duplicated molecules, and retrieved the isomeric SMILES of molecules to upgrade the data quality of our MFH substances database. As a result, we obtained 15,362 active ingredients of MFH substances from TCMID, 11,154 active ingredients from TCMSP, and an additional 803 active ingredients from 92 pieces of literature.

3.3. Chemical Space Analysis on the MFH Substances Database

Chemical space analysis is a widely utilized approach for exploring, comprehending, and optimizing a multitude of potential molecules [53]. Herein, we calculated molecular weight (MW) and octanol–water partition coefficient (LogP) on the collected active ingredient molecules. These two physicochemical properties enable us to measure and demonstrate the breadth of the chemical space distribution of the MFH database. MW and LogP were calculated via Python-based RDKit packages [54], and the chemical space distribution was visualized with Matplotlib [55]. For measuring the structural similarity of our MFH database, we computed the Tanimoto coefficients (TCs) [21] based on 1024 bits ECFP4 fingerprints by using RDKit packages [54]. To thoroughly analyze the diversity of molecular structures within the MFH database, we computed Murcko scaffolds of the molecules and conducted a clustering analysis using the K-means algorithm. Prior to clustering, we characterized the calculated Murcko scaffolds by ECFP4 fingerprints and employed T-distributed Stochastic Neighbor Embedding (t-SNE) [56] to reduce the high-dimensional data into two dimensions. ECFP4 fingerprints were computed via K-means clustering and t-SNE were implemented with scikit-learn [57].

3.4. Construction of Datasets for Building Classification and Regression Models

3.4.1. Datasets for Modeling on COX-2 Inhibitors

Dataset 1 and 2 (Table 11) for developing classification models to classify highly/weakly active COX-2 inhibitors were identical to the datasets in our previous work [58]. Different thresholds were employed to label highly/weakly active inhibitors, which enabled the constructed classification models to cover different chemical spaces. This strategy would further contribute to the hit rate of VS on MFH substances. A total of 1511 molecules (Dataset 3) were collected from ChEMBL [59], Reaxys [60], and SciFinder [61], their IC50 values were tested by enzyme-linked immunoassay. The pIC50 (−log10IC50) values ranged from 5.06 to 9.52. Dataset 3 was utilized to develop regression models with the aim of accurately predicting the bioactivities of MFH substances. External validation sets A1 and A2 were collected from the newly published literature and used to evaluate the generalizability of the constructed classification and regression models, respectively. The IC50 values of molecules in the External validation set A2 were all tested by enzyme-linked immunoassay.

3.4.2. Datasets for Modeling on mPGES-1 Inhibitors

Numerous mPGES-1 inhibitors with diverse structures were collected from ChEMBL, Reaxys, and SciFinder. Dataset 4 was composed of 3179 mPGES-1 inhibitors, their IC50 values vary from 0.0001 to 20,000 μM. Molecules with IC50 > 10 μM were weakly active inhibitors; with IC50 < 0.6 μM were highly active inhibitors. Dataset 5 comprised 3455 inhibitors with IC50 values ranging from 0.0001 to 20,000 μM. Molecules with IC50 ≥ 10 μM and IC50 < 10 μM are weakly and highly active inhibitors, respectively. Datasets 4 and 5 were employed to construct classification models on mPGES-1 inhibitors (shown in Table 9). Dataset 6 containing 735 inhibitors was derived from our previous work [62], pIC50 values of inhibitors in Dataset 6 ranged from 5.54 to 9 (all tested by homogeneous time-resolved fluorescence assay). External validation sets B1 and B2 were used to evaluate the generalizability of the constructed classification and regression models, respectively.

3.4.3. Splitting Strategy for Generating the Training/Test Set

The datasets for modeling on COX-2 inhibitors were divided into training/test sets at the ratio of 4:1. The datasets for modeling on mPGES-1 inhibitors were randomly divided into training/test sets at the ratio of 3:1. The datasets of both COX-2 and mPGES-1 inhibitors were randomly split 10 times for generating the training/test sets to avoid the random error. The random splitting was conducted by using the function StratifiedSplit of the Python toolkit scikit-learn 0.22.1 [57].

3.5. Characterization of Datasets

3.5.1. Binary Fingerprints for Classification Models

Three types of well-known fingerprints were employed to comprehensively and multifacetedly characterize the structural features of molecules within our MFH database. 166 bits MACCS fingerprints (belongs to dictionary-based fingerprints), 1024 bits Avalon fingerprints (topological-based fingerprints) [63], and 1024 bits ECFP4 fingerprints (circular fingerprints) [64] were computed with RDKit [54] packages. To avoid the inclusion of redundant information, the calculated fingerprints were then filtered by the variance, and fingerprints with variance in the bottom quartile were excluded from the construction of classification models.

3.5.2. Physicochemical Molecular Descriptors for Regression Models

Two types of physicochemical molecular descriptors were utilized to represent molecules for further developing QSAR models. A total of 22 global molecular descriptors and 96 3D property-weighted autocorrelation from CORINA were calculated by the CORINA Symphony software V1.0 [65]. A total of 115 2D physicochemical molecular descriptors and 85 FragmentCount descriptors from RDKit were computed via MayaChemTools [66]. The calculated descriptors were further screened with Pearson correlation coefficient (PCC), and listed in descending order by recursive feature elimination with the random forest estimator (RF-RFE) before modeling; more details can be seen in the previous work [62]. Additionally, all the reserved calculated descriptors were auto-scaled to the same range from 0.1 to 0.9.

3.6. Supervised Machine Learning Algorithms for Modeling

Machine learning (ML) algorithms are capable of discerning relationships within large datasets and devising optimal approaches for their analysis without prior specification. Four supervised ML algorithms, including support vector machine (SVM) [67], random forest (RF) [68], deep neural networks (DNNs) [69], and eXtreme Gradient Boosting (XGBoost) [70], were utilized to provide predictions on COX-2 and mPGES-1 inhibitors.

3.6.1. Modeling with SVM, RF, and XGBoost

The SVM algorithm equipped with radial basis function (RBF) kernel was involved in developing both classification and regression models, penalty parameter (C) and γ were hyperparameters to be confirmed during the optimization of the classification models. Except for the two parameters mentioned above, the insensitive parameter (ε) needed to be confirmed in the regression model’s optimization process. For RF models, the number of trees (n_estimators) and the maximum leaf nodes (max_leaf_nodes) were determined by the grid search. In modeling with XGBoost, the number of trees (n_estimators), the maximum depth of a tree (max_depth), the subsample ratio of the training instances (subsample), and the subsample ratio of columns when constructing each tree (colsample_bytree) were optimized by grid search. Other parameters not mentioned were set as their default values.
In addition to the common hyperparameters of SVM, RF, and XGBoost algorithms, the number of descriptors also served as a hyperparameter in the grid-based optimization of regression models. During the grid search, 5-fold cross-validation was repeated 10 times, and the optimal hyperparameters were determined based on the smallest mean squared error (MSE) of the 5-fold cross-validation [71]. Detailed ranges of those mentioned hyperparameters were listed in Table S7 in Supplementary Materials.

3.6.2. Modeling with DNN

Fully connected feed-forward neural networks with four hidden layers were constructed to develop classification and regression models. Neurons within the hidden layers were activated using the Relu function [72] and further compiled using the Adam optimizer with a learning rate of 0.0001. The training epoch was determined through repetitive 5-fold cross-validation training combined with early stopping [73]. For classification models, early stopping monitored changes in the predicted accuracy of the validation set. The validation set was generated through 5-fold cross-validation and was part of the training set. For regression models, early stopping focused on changes in the MSE values of the validation set. Training was halted when the accuracy or MSE of the validation set ceased to change within 50 epochs. To mitigate the potential for contingency arising from early stopping based on a single validation set, each cross-validation was conducted through repetitive training 50 times.

3.7. Ensemble Learning Based on Developed Classification Models

Ensemble learning integrates the predictions of multiple machine learning models to improve the robustness and accuracy of predictions made by individual base models [74]. In this study, the stacked generalization [75] was employed to combine the predicted probabilities of classification models built using the same algorithm but on datasets characterized by different fingerprints. These probabilities were then input into a logistic regression algorithm to generate predicted values of ensemble models with the same algorithm.

3.8. Unsupervised Machine Learning on MFH Substances

Unsupervised learning algorithms are trained on unlabeled data to discover hidden patterns or relationships among the data, as the hidden patterns and relationships can serve as a foundation for exploratory analysis [76]. The self-organized map (SOM) is a type of unsupervised artificial neural network; during the training of SOM, when the input layer receives data with similar vectors, these vectors are mapped onto the same neuron or neurons that are close together in the two-dimensional grid [77]. In this study, SOM was applied to perform clustering analysis on molecules of the MFH database, and further predict their inhibitory effects on COX-2 and mPGES-1.

3.9. Evaluation of Model Performances

The predicted accuracy (Q) and the Matthews correlation coefficient (MCC) were utilized as indicators of the performances of classification models. The coefficient of determination (R2) and root mean squared error (RMSE) were applied to evaluate the performances of regression models. These criteria mentioned above were calculated by the following equations:
Q = TP + TN   TP + TN + FP + FN ×   100 %  
MCC = TP × TN FP × FN TP + FP TP + FN TN + FP ( TN + FN )
where true positives (TPs) and true negatives (TNs) represent the number of “1” and “0” that were correctly predicted, respectively. False positives (FPs) and false negatives (FNs) represent the number of “1” and “0” that were wrongly predicted, respectively.
R 2 y , y ^ = 1 i = 1 n y i y ^ i 2 i = 1 n y i y - 2
RMSE y , y ^ = 1 n i = 1 n y i y ^ i 2
where n represents the total number of compounds; y represents an observed value of a compound; y ^ represents predicted value of a compound; y ¯ represents the average of y .

3.10. Pan Assay Interference Compounds (PAINS) Screening

The pan assay interference compounds (PAINS) screening has evolved into a pivotal element within drug design. The PAINS rule is introduced to identify false positive compounds (frequent hitters) during biological screening initiatives. We obtained a list of known aggregators (12,645 molecules were shown in Table S9) from Aggregator Advisor [78], which represented the known aggregator molecules and our screened MFH candidates with MACCS fingerprints. We employed RDKit to match the structure of each screened MFH candidate molecule with 12,645 known aggregator molecules. As a result, none of MFH candidates (10 potential COX-2 inhibitors and 15 potential mPGES-1 inhibitors) appeared in the known aggregators list.

3.11. Molecular Docking

In this study, molecular docking was conducted using the latest release of the widely used open-source program AutoDock Vina 1.2.0 [79]. Given the availability of multiple COX-2 and mPGES-1 crystal structures, the selection of an appropriate receptor with a low resolution is a prerequisite for reliable docking computations. Therefore, we chose COX-2 and mPGES-1 co-complex crystal structures from the PDB database that exhibited similar bound ligand to the screened medicine and food homologous (MFH) candidates. Ligands within the complex crystal structures and the screened MFH candidates were characterized using MACCS fingerprints. Subsequently, RDKit functions were employed to calculate Tanimoto similarity between the molecular structures of these entities. For each screened MFH candidate molecule, docking was performed with the crystal structure that exhibited the highest structural similarity (details were listed in Table S10). Before formal docking, the ligands within the original complex crystal structures underwent re-docking. The better the alignment between the re-docked ligand and the experimentally determined ligand, the more optimal the parameter settings and system preparation of the docking calculation. Results of the re-docking of ligands within the complex crystal structures are presented in Figure S1 of the Supplementary Materials. The protein preparation process involved the removal of water and other solvents, repair of missing residue sections, addition of hydrogen atoms to heavy atoms, and subsequent pre-docking energy minimization of the entire protein. The ligand preparation process included the addition of hydrogen atoms, computation of Gasteiger charges for all atoms, definition of rotatable bonds, and energy optimization. The grid box was adjusted based on the spatial center of the ligand within the crystal structure. Vina force field was employed during docking, with the exhaustiveness parameter set to 32.

4. Conclusions

In this study, we constructed a comprehensive database of 27,319 active ingredient molecules from 503 different types of medicine and food homologous (MFH). Analysis of the distribution of molecular weight (MW) and octanol–water partition coefficient (LogP) showed a wide range of values, indicating that our MFH database covers a wide chemical space. Structural diversity was assessed using Tanimoto coefficients (TCs), showing significant structural differences between molecules, with 95.92% of molecule pairs having TC values below 0.6. In addition, we performed Murcko scaffold analysis and K-means clustering, resulting in the identification of 11 different clusters in the MFH database. Among them, flavonoid clusters were the most abundant, followed by fatty acids, saponins, and sterols. The database was further enriched by the presence of lignans, alkaloids, triterpenoids, sesquiterpenoids, diterpenoids, and stilbenoids. Furthermore, we summarized the top 20 Murcko scaffolds, revealing diverse structures ranging from simple aromatic compounds with a single ring to complex systems with 7–8 rings. These findings collectively demonstrate the comprehensiveness and high structural diversity of our MFH database. Our MFH database will serve as a foundation for future studies, as it could facilitate the assessment of the effects of MFH on health and identifies potential mechanisms to accelerate the development of MFH-inspired products with nutritional and therapeutic value.
Based on datasets with different distributions of bioactivities, we employed four supervised learning algorithms (RF, SVM, DNN, and XGBoost), incorporating various fingerprints and physicochemical descriptors for modeling. As a result, a total of 240 classification models and 80 QSAR models were constructed for COX-2 and mPGES-1 inhibitors, respectively. Additionally, we also utilized ensemble learning to develop classification models. Based on the constructed single classifiers, another 80 integrated classification models were constructed for COX-2 and mPGES-1 inhibitors, respectively. For COX-2 inhibitors, ModelA_ensemble_RF_1, built on Dataset 1, demonstrated the best classification performance with MCC values of 0.802 and 0.603 on the test set and external validation set, respectively. ModelB_MACCS_SVM_6, constructed using Dataset 2, achieved the highest performance with MCC values of 0.657 and 0.572 on the test set and external validation set, respectively. ModelC_RDKIT_SVM_2, the optimal regression model based on Dataset 3, yielded RMSE values of 0.419 and 0.513 on the test set and external validation set, respectively. For mPGES-1 inhibitors, ModelD_ECFP_SVM_4 emerged as the top-performing classification model on Dataset 4, exhibiting MCC values of 0.832 and 0.584 on the test set and external validation set, respectively. ModelE_ECFP_SVM_1, the best classification model for Dataset 5, achieved MCC values of 0.799 and 0.579 on the test set and external validation set, respectively. ModelF_3D_SVM_1, based on Dataset 6, served as the optimal regression model with RMSE values of 0.253 and 0.35 on the test set and external validation set, respectively. These well-performing machine learning models can serve as powerful tools for virtual screening of the constructed MFH database, aiming to identify potential COX-2 and mPGES-1 inhibitors from MFH substances. Moreover, these models can be employed to predict the inhibitory capabilities of unknown compounds against COX-2 and mPGES-1, thus facilitating the discovery of novel anti-inflammatory drugs.
Finally, by means of a cascade ligand-based virtual screening strategy and a PAINS screening rule, we identified 10 potential COX-2 inhibitors and 15 potential mPGES-1 inhibitors from the MFH database. We verified candidates by molecular docking, investigated the interaction of the candidate molecules upon binding to COX-2 or mPGES-1. It is worth mentioning that some of these molecules have been previously reported to exhibit COX-2 inhibitory or anti-inflammatory activities. This demonstrates the effectiveness of the cascaded ligand-based virtual screening strategy employed in this study and provides design and modification ideas for the development of new effective anti-inflammatory drugs targeting COX-2 and mPGES-1.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28196782/s1, Detailed performances of all classification models built on Dataset 1 can be seen in Table S1. Performances of all classification models of Dataset 2 are listed in Table S2. Detailed results of all 80 QSAR models built on Dataset 3 are listed in Table S3. Detailed performances of all classification models built on Dataset 4 can be seen in Table S4. Performances of all classification models of Dataset 5 are listed in Table S5. Detailed results of all regression models built on Dataset 6 are listed in Table S6. The specific model results and parameters applied for virtual screening of MFH database are summarized in Table S7. All the datasets utilized in this study are summarized in Table S8. A list of known aggregators (12,645 molecules) is shown in Table S9. The PDB complex crystals utilized in this work are listed in Table S10. The alignment results of the re-docking of ligands within the complex crystal structures are presented in Figure S1.

Author Contributions

Y.T. Tian conceived the original idea of this work, performed all the machine learning models, carried out the virtual screening, analyzed the results, and took the lead in writing the manuscript. Z.Z. contributed to the construction of the medicine and food homology database, verified the developed models. A.Y. supervised the findings of this work, provided critical feedback and helped shape the research, analysis, and manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the codes used in this work can be found at: tyj-19951029/medicine-and-food-homologous-virtual-screnning (github.com, accessed on 3 June 2023). The MFH database established in this study can be obtained by sending an email to the author, and commercial use should be avoided.

Acknowledgments

The authors thank Molecular Networks GmbH, Nuremberg, Germany for providing the programs Corina and SONNIA for this work.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Not available.

References

  1. Tabas, I.; Glass, C.K. Anti-Inflammatory Therapy in Chronic Disease: Challenges and Opportunities. Science 2013, 339, 166–172. [Google Scholar] [CrossRef] [PubMed]
  2. Dennis, E.A.; Norris, P.C. Eicosanoid Storm in Infection and Inflammation. Nat. Rev. Immunol. 2015, 15, 511–523. [Google Scholar] [CrossRef] [PubMed]
  3. Ju, Z.; Li, M.; Xu, J.; Howell, D.C.; Li, Z.; Chen, F.-E. Recent Development on COX-2 Inhibitors as Promising Anti-Inflammatory Agents: The Past 10 Years. Acta Pharm. Sin. B 2022, 12, 2790–2807. [Google Scholar] [CrossRef] [PubMed]
  4. Zhang, Y.-Y.; Yao, Y.-D.; Luo, J.-F.; Liu, Z.-Q.; Huang, Y.-M.; Wu, F.-C.; Sun, Q.-H.; Liu, J.-X.; Zhou, H. Microsomal Prostaglandin E2 Synthase-1 and Its Inhibitors: Molecular Mechanisms and Therapeutic Significance. Pharmacol. Res. 2022, 175, 105977. [Google Scholar] [CrossRef]
  5. Xue, H.; Wang, W.; Bian, J.; Gao, Y.; Hao, Z.; Tan, J. Recent Advances in Medicinal and Edible Homologous Polysaccharides: Extraction, Purification, Structure, Modification, and Biological Activities. Int. J. Biol. Macromol. 2022, 222, 1110–1126. [Google Scholar] [CrossRef]
  6. Hu, Y.; Chen, X.; Hu, M.; Zhang, D.; Yuan, S.; Li, P.; Feng, L. Medicinal and Edible Plants in the Treatment of Dyslipidemia: Advances and Prospects. Chin. Med. J. 2022, 17, 113. [Google Scholar] [CrossRef]
  7. Yang, M.; Yan, T.; Yu, M.; Kang, J.; Gao, R.; Wang, P.; Zhang, Y.; Zhang, H.; Shi, L. Advances in Understanding of Health-promoting Benefits of Medicine and Food Homology Using Analysis of Gut Microbiota and Metabolomics. Food Front. 2020, 1, 398–419. [Google Scholar] [CrossRef]
  8. Ju, Z.; Lin, X.; Lu, X.; Tu, Z.; Wang, J.; Kaliyaperumal, K.; Liu, J.; Tian, Y.; Xu, S.; Liu, Y. Botryoisocoumarin A, a New COX-2 Inhibitor from the Mangrove Kandelia Candel Endophytic Fungus Botryosphaeria Sp. KcF6. J. Antibiot. 2015, 68, 653–656. [Google Scholar] [CrossRef]
  9. Ambati, G.G.; Jachak, S.M. Natural Product Inhibitors of Cyclooxygenase (COX) Enzyme: A Review on Current Status and Future Perspectives. Curr. Med. Chem. 2021, 28, 1877–1905. [Google Scholar] [CrossRef]
  10. Khan, H.; Rengasamy, K.R.R.; Pervaiz, A.; Nabavi, S.M.; Atanasov, A.G.; Kamal, M.A. Plant-Derived MPGES-1 Inhibitors or Suppressors: A New Emerging Trend in the Search for Small Molecules to Combat Inflammation. Eur. J. Med. Chem. 2018, 153, 2–28. [Google Scholar] [CrossRef]
  11. Schaible, A.M.; Traber, H.; Temml, V.; Noha, S.M.; Filosa, R.; Peduto, A.; Weinigel, C.; Barz, D.; Schuster, D.; Werz, O. Potent Inhibition of Human 5-Lipoxygenase and Microsomal Prostaglandin E2 Synthase-1 by the Anti-Carcinogenic and Anti-Inflammatory Agent Embelin. Biochem. Pharmacol. 2013, 86, 476–486. [Google Scholar] [CrossRef] [PubMed]
  12. Hwang, S.H.; Wecksler, A.T.; Wagner, K.; Hammock, B.D. Rationally Designed Multitarget Agents Against Inflammation and Pain. Curr. Med. Chem. 2013, 20, 1783–1799. [Google Scholar] [CrossRef] [PubMed]
  13. Badwan, B.A.; Liaropoulos, G.; Kyrodimos, E.; Skaltsas, D.; Tsirigos, A.; Gorgoulis, V.G. Machine Learning Approaches to Predict Drug Efficacy and Toxicity in Oncology. Cell Rep. Methods 2023, 3, 100413. [Google Scholar] [CrossRef]
  14. Wang, M.; Wang, J.; Weng, G.; Kang, Y.; Pan, P.; Li, D.; Deng, Y.; Li, H.; Hsieh, C.Y.; Hou, T. ReMODE: A Deep Learning-Based Web Server for Target-Specific Drug Design. J. Cheminform. 2022, 14, 84. [Google Scholar] [CrossRef] [PubMed]
  15. Kwon, Y.; Park, S.; Lee, J.; Kang, J.; Lee, H.J.; Kim, W. BEAR: A Novel Virtual Screening Method Based on Large-Scale Bioactivity Data. J. Chem. Inf. Model. 2023, 63, 1429–1437. [Google Scholar] [CrossRef]
  16. Zhang, L.-X.; Dong, J.; Wei, H.; Shi, S.-H.; Lu, A.-P.; Deng, G.-M.; Cao, D.-S. TCMSID: A Simplified Integrated Database for Drug Discovery from Traditional Chinese Medicine. J. Cheminform. 2022, 14, 89. [Google Scholar] [CrossRef]
  17. Yang, M.; Chen, J.; Shi, X.; Xu, L.; Xi, Z.; You, L.; An, R.; Wang, X. Development of in Silico Models for Predicting P-Glycoprotein Inhibitors Based on a Two-Step Approach for Feature Selection and Its Application to Chinese Herbal Medicine Screening. Mol. Pharm. 2015, 12, 3691–3713. [Google Scholar] [CrossRef]
  18. Ali, D.E.; Gedaily, R.A.E.; Ezzat, S.M.; Sawy, M.A.E.; Meselhy, M.R.; Abdel-Sattar, E. In Silico and in Vitro Anti-Inflammatory Study of Phenolic Compounds Isolated from Eucalyptus Maculata Resin. Sci. Rep. 2023, 13, 2093. [Google Scholar] [CrossRef]
  19. Traditional Chinese Medicine Integrated Database. Available online: http://47.100.169.139/tcmid/ (accessed on 3 June 2023).
  20. TCMSP—Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform. Available online: https://tcmsp-e.com/tcmsp.php (accessed on 3 June 2023).
  21. Godden, J.W.; Xue, L.; Bajorath, J. Combinatorial Preferences Affect Molecular Similarity/Diversity Calculations Using Binary Fingerprints and Tanimoto Coefficients. J. Chem. Inf. Comput. Sci. 2000, 40, 163–166. [Google Scholar] [CrossRef]
  22. Sinaga, K.P.; Yang, M.S. Unsupervised K-Means Clustering Algorithm. IEEE Access 2020, 8, 80716–80727. [Google Scholar] [CrossRef]
  23. Zhang, R.; Ren, S.; Dai, Q.; Shen, T.; Li, X.; Li, J.; Xiao, W. InflamNat: Web-Based Database and Predictor of Anti-Inflammatory Natural Products. J. Cheminform. 2022, 14, 30. [Google Scholar] [CrossRef] [PubMed]
  24. Fernandes, E.S.; Passos, G.F.; Medeiros, R.; Da Cunha, F.M.; Ferreira, J.; Campos, M.M.; Pianowski, L.F.; Calixto, J.B. Anti-Inflammatory Effects of Compounds Alpha-Humulene and (−)-Trans-Caryophyllene Isolated from the Essential Oil of Cordia Verbenacea. Eur. J. Pharmacol. 2007, 569, 228–236. [Google Scholar] [CrossRef] [PubMed]
  25. Lim, Y.; Lee, S.; Kim, J.Y.; Shin, J.H.; Kwon, O. A Phellinus Baumii–Based Supplement Containing Salvia Miltiorrhiza Bunge Improves Atherothrombotic Profiles through Endothelial Nitric Oxide Synthase and Cyclooxygenase Pathways in Vitro and in Vivo. J. Funct. Foods 2016, 24, 231–243. [Google Scholar] [CrossRef]
  26. Bayala, B.; Bassole, I.H.N.; Gnoula, C.; Nebie, R.; Yonli, A.; Morel, L.; Figueredo, G.; Nikiema, J.-B.; Lobaccaro, J.-M.A.; Simpore, J. Chemical Composition, Antioxidant, Anti-Inflammatory and Anti-Proliferative Activities of Essential Oils of Plants from Burkina Faso. PLoS ONE 2014, 9, e92122. [Google Scholar] [CrossRef] [PubMed]
  27. Tan, K.C.; Pham, T.X.; Lee, Y.; Lee, J.Y.; Balunas, M.J. Identification of Apocarotenoids as Chemical Markers of in Vitro Anti-Inflammatory Activity for Spirulina Supplements. J. Agric. Food Chem. 2021, 69, 12674–12685. [Google Scholar] [CrossRef]
  28. Morel, A.F.; Araujo, C.A. Antibacterial Cyclopeptide Alkaloids from the Bark of Condalia buxifolia. Phytochemistry 2002, 61, 561–566. [Google Scholar] [CrossRef] [PubMed]
  29. Chin, Y.W.; Jung, H.A.; Liu, Y.; Su, B.N.; Castoro, J.A.; Keller, W.J.; Pereira, M.A.; Kinghorn, A.D. Anti-Oxidant Constituents of the Roots and Stolons of Licorice (Glycyrrhiza glabra). J. Agric. Food Chem. 2007, 55, 4691–4697. [Google Scholar] [CrossRef]
  30. Vitale, P.; Panella, A.; Scilimati, A.; Perrone, M.G. COX-1 Inhibitors: Beyond Structure Toward Therapy: Selective COX-1 inhibitors. Med. Res. Rev. 2016, 36, 641–671. [Google Scholar] [CrossRef]
  31. Guo, F.; Feng, L.; Huang, C.; Ding, H.; Zhang, X.; Wang, Z.; Li, Y. Prenylflavone Derivatives from Broussonetia Papyrifera, Inhibit the Growth of Breast Cancer Cells in Vitro and in Vivo. Phytochem. Lett. 2013, 6, 331–336. [Google Scholar] [CrossRef]
  32. Lin, Y.; Kuang, Y.; Li, K.; Wang, S.; Song, W.; Qiao, X.; Sabir, G.; Ye, M. Screening for Bioactive Natural Products from a 67-Compound Library of Glycyrrhiza Inflata. Bioorg. Med. Chem. 2017, 25, 3706–3713. [Google Scholar] [CrossRef]
  33. Damodar, K.; Kim, J.K.; Jun, J.G. Synthesis and Pharmacological Properties of Naturally Occurring Prenylated and Pyranochalcones as Potent Anti-Inflammatory Agents. Chin. Chem. Lett. 2016, 27, 698–702. [Google Scholar] [CrossRef]
  34. Nishizawa, M.; Izuhara, R.; Kaneko, K.; Fujimoto, Y. 3-caffeoyl-4-sinapoylquinic acid, a novel lipdxygenase inhibitor from gardeniae fructus. Chem. Pharm. Bull. 1987, 35, 2133–2135. [Google Scholar] [CrossRef] [PubMed]
  35. Song, J.H.; Cui, L.; An, L.B.; Li, W.T.; Fang, Z.Z.; Zhang, Y.Y.; Dong, P.-P.; Wu, X.; Wang, L.-X.; Gonzalez, F.J.; et al. Inhibition of UDP-Glucuronosyltransferases (UGTs) Activity by Constituents of Schisandra Chinensis: UGTs’ Inhibition by Schisandra chinensis. Phytother. Res. 2015, 29, 1658–1664. [Google Scholar] [CrossRef] [PubMed]
  36. Chiu, T.; Ku, C.; Ho, T.; Tsai, K.; Yang, Y.; Ou, H.; Chen, H. Schisanhenol Ameliorates oxLDL -caused Endothelial Dysfunction by Inhibiting LOX -1 Signaling. Environ. Toxicol. 2023, 38, tox.23788. [Google Scholar] [CrossRef]
  37. Yamakuni, T.; Aoki, K.; Nakatani, K.; Kondo, N.; Oku, H.; Ishiguro, K.; Ohizumi, Y. Garcinone B Reduces Prostaglandin E2 Release and NF-ΚB-Mediated Transcription in C6 Rat Glioma Cells. Neurosci. Lett. 2006, 394, 206–210. [Google Scholar] [CrossRef]
  38. Hu, X.; Yu, M.H.; Yan, G.R.; Wang, H.Y.; Hou, A.J.; Lei, C. Isoprenylated Phenolic Compounds with Tyrosinase Inhibition from Morus nigra. J. Asian Nat. Prod. Res. 2018, 20, 488–493. [Google Scholar] [CrossRef]
  39. Coste, C.; Gérard, N.; Dinh, C.P.; Bruguière, A.; Rouger, C.; Leong, S.T.; Awang, K.; Richomme, P.; Derbré, S.; Charreau, B. Targeting MHC Regulation Using Polycyclic Polyprenylated Acylphloroglucinols Isolated from Garcinia Bancana. Biomolecules 2020, 10, 1266. [Google Scholar] [CrossRef]
  40. Li, K.; Yu, X.H.; Maskey, A.R.; Musa, I.; Wang, Z.Z.; Garcia, V.; Guo, A.; Yang, N.; Srivastava, K.; Dunkin, D.; et al. Cytochrome P450 3A4 Suppression by Epimedium and Active Compound Kaempferol Leads to Synergistic Anti-Inflammatory Effect with Corticosteroid. Front. Pharmacol. 2023, 13, 1042756. [Google Scholar] [CrossRef]
  41. Kim, H.W.; Kim, J.B.; Shanmugavelan, P.; Kim, S.N.; Cho, Y.S.; Kim, H.R.; Lee, J.T.; Jeon, W.-T.; Lee, D.J. Evaluation of γ-Oryzanol Content and Composition from the Grains of Pigmented Rice-Germplasms by LC-DAD-ESI/MS. BMC Res. Notes 2013, 6, 149–159. [Google Scholar] [CrossRef]
  42. Lee, K.-T.; Sohn, I.-C.; Kim, D.-H.; Choi, J.-W.; Kwon, S.-H.; Park, H.-J. Hypoglycemic and Hypolipidemic Effects of Tectorigenin and Kaikasaponin III in the Streptozotocin-Lnduced Diabetic Rat and Their Antioxidant Activity in Vitro. Arch. Pharm. Res. 2000, 23, 461–466. [Google Scholar] [CrossRef]
  43. Chen, X.; Zhang, J.; Li, R.; Zhang, H.; Sun, Y.; Jiang, L.; Wang, X.; Xiong, Y. Flos Puerariae-Semen Hoveniae Medicinal Pair Extract Ameliorates DSS-Induced Inflammatory Bowel Disease through Regulating MAPK Signaling and Modulating Gut Microbiota Composition. Front. Pharmacol. 2022, 13, 1034031. [Google Scholar] [CrossRef] [PubMed]
  44. Sun, Y.; Long, J.; Chen, W.; Sun, Y.; Zhou, L.; Zhang, L.; Zeng, H.; Yuan, D. Alisol B 23-Acetate, a New Promoter for Cholesterol Efflux from Dendritic Cells, Alleviates Dyslipidemia and Inflammation in Advanced Atherosclerotic Mice. Int. Immunopharmacol. 2021, 99, 107956. [Google Scholar] [CrossRef] [PubMed]
  45. Baell, J.B.; Holloway, G.A. New Substructure Filters for Removal of Pan Assay Interference Compounds (PAINS) from Screening Libraries and for Their Exclusion in Bioassays. J. Med. Chem. 2010, 53, 2719–2740. [Google Scholar] [CrossRef] [PubMed]
  46. Marnett, L.J. Recent Developments in Cyclooxygenase Inhibition. Prostaglandins Other Lipid Mediat. 2002, 68–69, 153–164. [Google Scholar] [CrossRef]
  47. Romero-Estrada, A.; Boto, A.; González-Christen, J.; Romero-Estudillo, I.; Garduño-Ramírez, M.L.; Razo-Hernández, R.S.; Marquina, S.; Maldonado-Magaña, A.; Columba-Palomares, M.C.; Sánchez-Carranza, J.N.; et al. Synthesis, Biological Evaluation, and Molecular Docking Study of 3-Amino and 3-Hydroxy- Seco A Derivatives of α-Amyrin and 3-Epilupeol as Inhibitors of COX-2 Activity and NF-KB Activation. J. Nat. Prod. 2022, 85, 787–803. [Google Scholar] [CrossRef]
  48. Shekfeh, S.; Çalışkan, B.; Fischer, K.; Yalçın, T.; Garscha, U.; Werz, O.; Banoglu, E. A Multi-step Virtual Screening Protocol for Identification of Novel Non-acidic Microsomal Prostaglandin E2 Synthase-1 (MPGES-1) Inhibitors. ChemMedChem 2018, 14, 273–281. [Google Scholar] [CrossRef]
  49. National Health Commission of P. R. China. Management Approach of Food and Chinese Medicine Homologous Catalogue According to Tradition. Available online: http://www.nhc.gov.cn/wjw/yjzj/201411/67ac54fb05ed46929adc63f2db31d4bf.shtml (accessed on 3 June 2023).
  50. National Health Commission of P. R. China. List of 9 Pilot Food and Chinese Medicine Homologous Substances According to Tradition. Available online: http://www.nhc.gov.cn/sps/s7885/202001/1ec2cca04146450d9b14acc2499d854f.shtml (accessed on 3 June 2023).
  51. Pharmacopoeia Committee of P. R. China. Pharmacopoeia of People’s Republic of China, 11th ed.; China Medical Science and Technology Press: Beijing, China, 2020. [Google Scholar]
  52. CNKI. Available online: https://www.cnki.net/index/ (accessed on 3 June 2023).
  53. Schaduangrat, N.; Anuwongcharoen, N.; Charoenkwan, P.; Shoombuatong, W. DeepAR: A Novel Deep Learning-Based Hybrid Framework for the Interpretable Prediction of Androgen Receptor Antagonists. J. Cheminform. 2023, 15, 50. [Google Scholar] [CrossRef]
  54. RDKit: Open-Source Cheminformatics Software. Available online: https://rdkit.org/ (accessed on 3 June 2023).
  55. Matplotlib: Visualization with Python. Available online: https://matplotlib.org/ (accessed on 3 June 2023).
  56. Gisbrecht, A.; Schulz, A.; Hammer, B. Parametric Nonlinear Dimensionality Reduction Using Kernel T-SNE. Neurocomputing 2015, 147, 71–82. [Google Scholar] [CrossRef]
  57. Pedregosa, F.; Varoquaux, G.; Gramfort, A.; Michel, V.; Thirion, B.; Grisel, O.; Blondel, M.; Prettenhofer, P.; Weiss, R.; Dubourg, V.; et al. Scikit-Learn: Machine Learning in Python. J. Mach. Learn. Res. 2011, 12, 2825–2830. [Google Scholar]
  58. Qin, Z.; Xi, Y.; Zhang, S.; Tu, G.; Yan, A. Classification of Cyclooxygenase-2 Inhibitors Using Support Vector Machine and Random Forest Methods. J. Chem. Inf. Model. 2019, 59, 1988–2008. [Google Scholar] [CrossRef]
  59. European Bioinformatics Institute of European Molecular Biology Laboratory ChEMBL. Available online: https://www.ebi.ac.uk/chembl/g/ (accessed on 3 June 2023).
  60. Elsevier Reaxys. Available online: https://www.reaxys.com/#/search/quick (accessed on 3 June 2023).
  61. Chemical Abstracts Service of America SciFinder. Available online: https://scifinder-n.cas.org/ (accessed on 3 June 2023).
  62. Tian, Y.; Yang, Z.; Wang, H.; Yan, A. Prediction of Bioactivities of Microsomal Prostaglandin E2 Synthase-1 Inhibitors by Machine Learning Algorithms. Chem. Biol. Drug Des. 2023, 101, 1307–1321. [Google Scholar] [CrossRef]
  63. Kc, G.B.; Bocci, G.; Verma, S.; Hassan, M.M.; Holmes, J.; Yang, J.J.; Sirimulla, S.; Oprea, T.I. A Machine Learning Platform to Estimate Anti-SARS-CoV-2 Activities. Nat. Mach. Intell. 2021, 3, 527–535. [Google Scholar] [CrossRef]
  64. Riniker, S.; Landrum, G.A. Open-Source Platform to Benchmark Fingerprints for Ligand-Based Virtual Screening. J. Cheminform. 2013, 5, 26. [Google Scholar] [CrossRef] [PubMed]
  65. CORINA Symphony, version 1.0; Molecular Networks GmbH: Germany, Altamira; LLC: New York, NY, USA, 2018; Available online: https://mn-am.com/products/corinasymphony/ (accessed on 3 June 2023).
  66. MayaChemTools: A Growing Collection of Perl and Python Scripts, Modules, and Classes to Support a Variety of Day-to-Day Computational Discovery Needs. Available online: http://www.mayachemtools.org/ (accessed on 4 June 2023).
  67. Noble, W.S. What Is a Support Vector Machine? Nat. Biotechnol. 2006, 24, 1565–1567. [Google Scholar] [CrossRef] [PubMed]
  68. Belgiu, M.; Drăguţ, L. Random Forest in Remote Sensing: A Review of Applications and Future Directions. ISPRS J. Photogramm. Remote Sens. 2016, 114, 24–31. [Google Scholar] [CrossRef]
  69. LeCun, Y.; Bengio, Y.; Hinton, G. Deep Learning. Nature 2015, 521, 436–444. [Google Scholar] [CrossRef]
  70. Wang, L.; Ding, J.; Shi, P.; Fu, L.; Pan, L.; Tian, J.; Cao, D.; Jiang, H.; Ding, X. Ensemble Machine Learning to Evaluate the in Vivo Acute Oral Toxicity and in Vitro Human Acetylcholinesterase Inhibitory Activity of Organophosphates. Arch. Toxicol. 2021, 95, 2443–2457. [Google Scholar] [CrossRef]
  71. Krstajic, D.; Buturovic, L.J.; Leahy, D.E.; Thomas, S. Cross-Validation Pitfalls When Selecting and Assessing Regression and Classification Models. J. Cheminform. 2014, 6, 10. [Google Scholar] [CrossRef]
  72. Petersen, P.; Voigtlaender, F. Optimal Approximation of Piecewise Smooth Functions Using Deep ReLU Neural Networks. Neural Netw. 2018, 108, 296–330. [Google Scholar] [CrossRef]
  73. Vilares Ferro, M.; Doval Mosquera, Y.; Ribadas Pena, F.J.; Darriba Bilbao, V.M. Early Stopping by Correlating Online Indicators in Neural Networks. Neural Netw. 2023, 159, 109–124. [Google Scholar] [CrossRef]
  74. Chen, C.H.; Tanaka, K.; Kotera, M.; Funatsu, K. Comparison and Improvement of the Predictability and Interpretability with Ensemble Learning Models in QSPR Applications. J. Cheminform. 2020, 12, 19. [Google Scholar] [CrossRef] [PubMed]
  75. Lou, C.; Yang, H.; Deng, H.; Huang, M.; Li, W.; Liu, G.; Lee, P.W.; Tang, Y. Chemical Rules for Optimization of Chemical Mutagenicity via Matched Molecular Pairs Analysis and Machine Learning Methods. J. Cheminform. 2023, 15, 35. [Google Scholar] [CrossRef] [PubMed]
  76. Glielmo, A.; Husic, B.E.; Rodriguez, A.; Clementi, C.; Noé, F.; Laio, A. Unsupervised Learning Methods for Molecular Simulation Data. Chem. Rev. 2021, 121, 9722–9758. [Google Scholar] [CrossRef] [PubMed]
  77. Tian, Y.; Zhang, S.; Yin, H.; Yan, A. Quantitative Structure-Activity Relationship (QSAR) Models and Their Applicability Domain Analysis on HIV-1 Protease Inhibitors by Machine Learning Methods. Chemom. Intell. Lab. Syst. 2020, 196, 103888. [Google Scholar] [CrossRef]
  78. Irwin, J.J.; Duan, D.; Torosyan, H.; Doak, A.K.; Ziebart, K.T.; Sterling, T.; Tumanian, G.; Shoichet, B.K. An Aggregation Advisor for Ligand Discovery. J. Med. Chem. 2015, 58, 7076–7087. [Google Scholar] [CrossRef]
  79. Eberhardt, J.; Santos-Martins, D.; Tillack, A.F.; Forli, S. AutoDock Vina 1.2.0: New Docking Methods, Expanded Force Field, and Python Bindings. J. Chem. Inf. Model. 2021, 61, 3891–3898. [Google Scholar] [CrossRef]
Figure 1. Active ingredients of natural products with anti-inflammatory effects. (A) Molecules with COX-2 inhibition; (B) Molecules with mPGES-1 inhibition.
Figure 1. Active ingredients of natural products with anti-inflammatory effects. (A) Molecules with COX-2 inhibition; (B) Molecules with mPGES-1 inhibition.
Molecules 28 06782 g001
Figure 2. The workflow of this study. The pink circle in the interaction diagram represents polar amino acids, the green represents hydrophobic amino acids.
Figure 2. The workflow of this study. The pink circle in the interaction diagram represents polar amino acids, the green represents hydrophobic amino acids.
Molecules 28 06782 g002
Figure 3. The distribution of molecular weight (MW) and LogP on the MFH substances database. The dark red line/dot, the navy-blue line/dot, and the green line/dot represent molecules derived from TCMID database, TCMSP database, and the literature, respectively.
Figure 3. The distribution of molecular weight (MW) and LogP on the MFH substances database. The dark red line/dot, the navy-blue line/dot, and the green line/dot represent molecules derived from TCMID database, TCMSP database, and the literature, respectively.
Molecules 28 06782 g003
Figure 4. The overview of the chemical structures in the medicine and food homology (MFH) database. (A) The 11 classes of molecules of our MFH database after clustering analysis; (B) Top 20 Murcko scaffolds in the MFH database; percentages of molecules containing the Murcko scaffolds are labeled in the bottom-right corner.
Figure 4. The overview of the chemical structures in the medicine and food homology (MFH) database. (A) The 11 classes of molecules of our MFH database after clustering analysis; (B) Top 20 Murcko scaffolds in the MFH database; percentages of molecules containing the Murcko scaffolds are labeled in the bottom-right corner.
Molecules 28 06782 g004
Figure 5. (A) MCC values of single classification models and ensemble models constructed using three fingerprints and four machine learning algorithms for the test set based on Dataset 1. (B) MCC values of single classification models and ensemble models constructed using three fingerprints and four machine learning algorithms for Dataset 2 for the test set. The solid red line in the middle of the box chart represents the median, and the circle represents outliers.
Figure 5. (A) MCC values of single classification models and ensemble models constructed using three fingerprints and four machine learning algorithms for the test set based on Dataset 1. (B) MCC values of single classification models and ensemble models constructed using three fingerprints and four machine learning algorithms for Dataset 2 for the test set. The solid red line in the middle of the box chart represents the median, and the circle represents outliers.
Molecules 28 06782 g005
Figure 6. (A) RMSE values on the test sets of the regression models constructed using Corina 3D (represented as G_3D) on Dataset 3; (B) RMSE values on the test sets of the regression models constructed using RDKit descriptors for the test set based on Dataset 3.
Figure 6. (A) RMSE values on the test sets of the regression models constructed using Corina 3D (represented as G_3D) on Dataset 3; (B) RMSE values on the test sets of the regression models constructed using RDKit descriptors for the test set based on Dataset 3.
Molecules 28 06782 g006
Figure 7. (A) MCC values of single classification models and ensemble models constructed using three fingerprints and four machine learning algorithms for Dataset 4 for the test set. (B) MCC values of single classification models and ensemble models constructed using three fingerprints and four machine learning algorithms for the test set based on Dataset 5. The solid red line in the middle of the box chart represents the median, and the circle represents outliers.
Figure 7. (A) MCC values of single classification models and ensemble models constructed using three fingerprints and four machine learning algorithms for Dataset 4 for the test set. (B) MCC values of single classification models and ensemble models constructed using three fingerprints and four machine learning algorithms for the test set based on Dataset 5. The solid red line in the middle of the box chart represents the median, and the circle represents outliers.
Molecules 28 06782 g007
Figure 8. (A) RMSE values on the test sets of the regression models constructed using Corina 3D (represented as G_3D) on Dataset 6; (B) RMSE values on the test sets of the regression models constructed using RDKit descriptors for the test set based on Dataset 6.
Figure 8. (A) RMSE values on the test sets of the regression models constructed using Corina 3D (represented as G_3D) on Dataset 6; (B) RMSE values on the test sets of the regression models constructed using RDKit descriptors for the test set based on Dataset 6.
Molecules 28 06782 g008
Figure 9. The mapping positions of COX-2 inhibitor candidates in self-organizing map. HP represents the proportion of the number of highly active inhibitors in the grid to all molecules in the grid. The larger the calculated HP of one grid (close to 1), the warmer the color tends to be (red). On the contrary, the closer the HP value of a grid is to 0, the cooler the color tends to be (blue).
Figure 9. The mapping positions of COX-2 inhibitor candidates in self-organizing map. HP represents the proportion of the number of highly active inhibitors in the grid to all molecules in the grid. The larger the calculated HP of one grid (close to 1), the warmer the color tends to be (red). On the contrary, the closer the HP value of a grid is to 0, the cooler the color tends to be (blue).
Molecules 28 06782 g009
Figure 10. The mapping positions of mPGES-1 inhibitors candidates in self-organizing map. HP represents the proportion of the number of highly active inhibitors in the grid to all molecules in the grid. The larger the calculated HP of one grid (close to 1), the warmer the color tends to be (red). On the contrary, the closer the HP value of a grid is to 0, the cooler the color tends to be (blue).
Figure 10. The mapping positions of mPGES-1 inhibitors candidates in self-organizing map. HP represents the proportion of the number of highly active inhibitors in the grid to all molecules in the grid. The larger the calculated HP of one grid (close to 1), the warmer the color tends to be (red). On the contrary, the closer the HP value of a grid is to 0, the cooler the color tends to be (blue).
Molecules 28 06782 g010
Table 1. The overall performances of classification models on Dataset 1.
Table 1. The overall performances of classification models on Dataset 1.
FingerprintsAlgorithmTraining Set5-Fold Cross-ValidationTest SetExternal Validation Set A1
Q (ave_std)MCC (ave_std)Q (ave_std)MCC (ave_std)Q (ave_std)MCC (ave_std)Q (ave_std)MCC (ave_std)
AvalonRF0.911 ± 0.0050.817 ± 0.0090.853 ± 0.0070.7 ± 0.0150.863 ± 0.0140.721 ± 0.0260.777 ± 0.0270.524 ± 0.008
AvalonSVM0.97 ± 0.0070.939 ± 0.0140.892 ± 0.0060.779 ± 0.0120.893 ± 0.0110.783 ± 0.0230.786 ± 0.0370.531 ± 0.017
AvalonDNN0.996 ± 0.0020.991 ± 0.0050.972 ± 0.0060.942 ± 0.0120.883 ± 0.0110.759 ± 0.0240.778 ± 0.0330.532 ± 0.011
AvalonXGBoost0.975 ± 0.0240.949 ± 0.0540.8893 ± 0.0310.794 ± 0.0610.89 ± 0.0110.774 ± 0.0230.792 ± 0.0360.541 ± 0.02
ECFP4RF0.91 ± 0.0080.814 ± 0.0170.85 ± 0.0050.691 ± 0.010.857 ± 0.0150.704 ± 0.0310.763 ± 0.020.508 ± 0.024
ECFP4SVM0.988 ± 0.0070.975 ± 0.0140.885 ± 0.0070.764 ± 0.0150.885 ± 0.0160.763 ± 0.0330.794 ± 0.0230.537 ± 0.022
ECFP4DNN0.994 ± 0.0010.988 ± 0.0020.975 ± 0.0050.948 ± 0.0110.87 ± 0.0190.732 ± 0.0380.773 ± 0.0210.515 ± 0.017
ECFP4XGBoost0.99 ± 0.0030.981 ± 0.0050.883 ± 0.0020.761 ± 0.0040.887 ± 0.0110.768 ± 0.0230.767 ± 0.020.531 ± 0.013
MACCSRF0.897 ± 0.0040.79 ± 0.0080.844 ± 0.0050.68 ± 0.010.854 ± 0.0120.701 ± 0.0250.737 ± 0.0180.519 ± 0.011
MACCSSVM0.945 ± 0.0130.887 ± 0.0270.874 ± 0.0070.742 ± 0.0140.876 ± 0.0130.746 ± 0.0260.758 ± 0.0240.532 ± 0.011
MACCSDNN0.973 ± 0.0050.945 ± 0.010.939 ± 0.0090.874 ± 0.0190.862 ± 0.0120.716 ± 0.0230.741 ± 0.0150.526 ± 0.01
MACCSXGBoost0.968 ± 0.0060.931 ± 0.0140.875 ± 0.0050.745 ± 0.0110.881 ± 0.0120.753 ± 0.0250.775 ± 0.020.54 ± 0.023
EnsembleRF0.933 ± 0.0180.835 ± 0.0160.876 ± 0.0210.717 ± 0.020.886 ± 0.0240.735 ± 0.0290.787 ± 0.0250.535 ± 0.023
EnsembleSVM0.986 ± 0.0070.952 ± 0.0130.901 ± 0.0060.78 ± 0.010.903 ± 0.0090.782 ± 0.0160.797 ± 0.0170.551 ± 0.012
EnsembleDNN0.996 ± 0.0020.986 ± 0.0070.971 ± 0.0060.931 ± 0.0120.883 ± 0.0140.747 ± 0.0250.775 ± 0.0130.535 ± 0.009
EnsembleXGBoost0.987 ± 0.0080.963 ± 0.0180.893 ± 0.0090.777 ± 0.0180.896 ± 0.0090.774 ± 0.0190.787 ± 0.0170.546 ± 0.014
Dataset 1 contains 1630 COX-2 inhibitors (the number of inhibitors in the training/test set = 1304:326); external validation set A1 contains 368 COX-2 inhibitors. Q: Accuracy; MCC: Matthews Correlation Coefficient, ave_std: average values and standard deviation values, Ensemble: combined the predicted probabilities of classification models built with Avalon, ECFP4, and MACCS fingerprints.
Table 2. The overall performances of classification models on Dataset 2.
Table 2. The overall performances of classification models on Dataset 2.
FingerprintsAlgorithmTraining Set5-Fold Cross-ValidationTest SetExternal Validation Set A1
Q (ave_std)MCC (ave_std)Q (ave_std)MCC (ave_std)Q (ave_std)MCC (ave_std)Q (ave_std)MCC (ave_std)
AvalonRF0.838 ± 0.0040.67 ± 0.0080.786 ± 0.0060.566 ± 0.0120.795 ± 0.0140.586 ± 0.0290.747 ± 0.0170.505 ± 0.01
AvalonSVM0.932 ± 0.0170.863 ± 0.0350.814 ± 0.0060.625 ± 0.0130.817 ± 0.0080.632 ± 0.0180.785 ± 0.0140.529 ± 0.016
AvalonDNN0.992 ± 0.0020.984 ± 0.0030.943 ± 0.0070.883 ± 0.0150.802 ± 0.0120.602 ± 0.0250.747 ± 0.210.507 ± 0.01
AvalonXGBoost0.97 ± 0.0190.938 ± 0.0350.818 ± 0.0070.628 ± 0.0130.825 ± 0.0060.644 ± 0.0140.762 ± 0.0260.523 ± 0.015
ECFP4RF0.828 ± 0.0050.65 ± 0.0110.778 ± 0.0050.548 ± 0.0110.784 ± 0.0140.561 ± 0.0280.713 ± 0.0260.497 ± 0.013
ECFP4SVM0.948 ± 0.0250.894 ± 0.0520.811 ± 0.0030.616 ± 0.0060.817 ± 0.010.63 ± 0.0220.757 ± 0.0230.532 ± 0.018
ECFP4DNN0.989 ± 0.0020.977 ± 0.0030.956 ± 0.0050.91 ± 0.010.8 ± 0.0120.595 ± 0.0240.725 ± 0.0120.509 ± 0.009
ECFP4XGBoost0.947 ± 0.0230.896 ± 0.0430.811 ± 0.0060.617 ± 0.0120.813 ± 0.0130.625 ± 0.0280.755 ± 0.0460.523 ± 0.022
MACCSRF0.811 ± 0.0060.617 ± 0.0120.77 ± 0.0060.534 ± 0.0130.773 ± 0.0110.541 ± 0.0220.713 ± 0.0130.497 ± 0.009
MACCSSVM0.897 ± 0.0160.792 ± 0.0330.809 ± 0.0040.614 ± 0.010.812 ± 0.0120.621 ± 0.0250.758 ± 0.0270.529 ± 0.022
MACCSDNN0.946 ± 0.010.893 ± 0.0170.887 ± 0.010.772 ± 0.020.79 ± 0.010.575 ± 0.0220.732 ± 0.0110.506 ± 0.008
MACCSXGBoost0.909 ± 0.0120.816 ± 0.0250.807 ± 0.0040.611 ± 0.0070.817 ± 0.0090.626 ± 0.0120.753 ± 0.0130.515 ± 0.01
EnsembleRF0.84 ± 0.0060.66 ± 0.0080.792 ± 0.0070.563 ± 0.0110.798 ± 0.0110.576 ± 0.020.739 ± 0.0160.513 ± 0.01
EnsembleSVM0.94 ± 0.0090.864 ± 0.0190.827 ± 0.0080.633 ± 0.0110.83 ± 0.0070.642 ± 0.0160.78 ± 0.0140.544 ± 0.011
EnsembleDNN0.989 ± 0.0050.965 ± 0.0070.942 ± 0.0060.868 ± 0.0120.812 ± 0.0080.604 ± 0.0150.748 ± 0.0040.52 ± 0.003
EnsembleXGBoost0.953 ± 0.0140.895 ± 0.0270.823 ± 0.0050.63 ± 0.0080.83 ± 0.0060.643 ± 0.0120.767 ± 0.0180.531 ± 0.01
Dataset 2 contains 2925 COX-2 inhibitors (the number of inhibitors in the training/test set = 2340:585); external validation set A1 contains 368 COX-2 inhibitors. Q: Accuracy; MCC: Matthews Correlation Coefficient, ave_std: average values and standard deviation values, Ensemble: combined the predicted probabilities of classification models built with Avalon, ECFP4, and MACCS fingerprints.
Table 3. The overall performances of regression models on Dataset 3.
Table 3. The overall performances of regression models on Dataset 3.
DescriptorsAlgorithmTraining Set5-Fold Cross-ValidationTest SetExternal Validation Set A2
R2 (ave_std)RMSE (ave_std)R2 (ave_std)RMSE (ave_std)R2 (ave_std)RMSE (ave_std)R2 (ave_std)RMSE (ave_std)
G_3DRF0.847 ± 0.0040.346 ± 0.0040.633 ± 0.010.537 ± 0.0030.651 ± 0.0070.517 ± 0.0210.539 ± 0.0630.634 ± 0.036
G_3DSVM0.941 ± 0.0110.214 ± 0.0230.713 ± 0.0070.475 ± 0.0060.725 ± 0.0170.465 ± 0.0140.553 ± 0.050.612 ± 0.034
G_3DDNN0.951 ± 0.0210.191 ± 0.0410.833 ± 0.060.354 ± 0.0620.691 ± 0.0250.488 ± 0.0180.536 ± 0.0660.63 ± 0.036
G_3DXGBoost0.907 ± 0.0220.268 ± 0.0350.817 ± 0.0240.375 ± 0.0260.684 ± 0.0310.492 ± 0.020.539 ± 0.050.634 ± 0.032
RDKitRF0.858 ± 0.0070.334 ± 0.0090.679 ± 0.0080.503 ± 0.0070.688 ± 0.0250.485 ± 0.0280.532 ± 0.0410.607 ± 0.021
RDKitSVM0.959 ± 0.0120.179 ± 0.0270.744 ± 0.010.449 ± 0.010.746 ± 0.0220.436 ± 0.0190.581 ± 0.050.596 ± 0.038
RDKitDNN0.965 ± 0.0420.143 ± 0.0760.862 ± 0.0960.307 ± 0.1020.72 ± 0.0260.471 ± 0.0310.594 ± 0.0370.591 ± 0.03
RDKitXGBoost0.988 ± 0.0020.095 ± 0.010.923 ± 0.0120.243 ± 0.0210.703 ± 0.0250.477 ± 0.0150.594 ± 0.0380.584 ± 0.035
Dataset 3 contains 1511 COX-2 inhibitors with IC50 values tested in vitro by enzyme-linked immunoassay (the number of inhibitors in the training/test set = 1209:302) external validation set A2 contains 114 COX-2 inhibitors, R2: the coefficient of determination; RMSE: root mean squared error, ave_std: average values and standard deviation values.
Table 4. The overall performances of classification models on Dataset 4.
Table 4. The overall performances of classification models on Dataset 4.
FingerprintsAlgorithmTraining Set5-Fold Cross-ValidationTest SetExternal Validation Set B1
Q (ave_std)MCC (ave_std)Q (ave_std)MCC (ave_std)Q (ave_std)MCC (ave_std)Q (ave_std)MCC (ave_std)
AvalonRF0.914 ± 0.0030.722 ± 0.0120.882 ± 0.0040.605 ± 0.0140.888 ± 0.0060.627 ± 0.0220.744 ± 0.0130.5 ± 0.013
AvalonSVM0.978 ± 0.0030.931 ± 0.010.916 ± 0.0040.736 ± 0.0130.917 ± 0.0060.739 ± 0.0210.76 ± 0.0060.513 ± 0.009
AvalonDNN0.992 ± 0.0010.977 ± 0.0050.966 ± 0.0060.894 ± 0.0190.908 ± 0.0060.711 ± 0.0190.747 ± 0.010.501 ± 0.011
AvalonXGBoost0.99 ± 0.0040.971 ± 0.0110.913 ± 0.0040.73 ± 0.0080.918 ± 0.0050.74 ± 0.0140.747 ± 0.0170.508 ± 0.006
ECFP4RF0.912 ± 0.0030.713 ± 0.0120.878 ± 0.0040.587 ± 0.0160.882 ± 0.0060.603 ± 0.0240.726 ± 0.010.501 ± 0.013
ECFP4SVM0.993 ± 0.0040.977 ± 0.0120.94 ± 0.0030.812 ± 0.010.941 ± 0.0070.817 ± 0.0220.77 ± 0.0120.543 ± 0.02
ECFP4DNN0.999 ± 0.0010.998 ± 0.0010.985 ± 0.0030.954 ± 0.0080.925 ± 0.0080.767 ± 0.0240.768 ± 0.0120.533 ± 0.01
ECFP4XGBoost0.995 ± 0.0040.987 ± 0.010.933 ± 0.0040.796 ± 0.0080.946 ± 0.0060.83 ± 0.0090.778 ± 0.0050.547 ± 0.015
MACCSRF0.907 ± 0.0050.698 ± 0.0190.88 ± 0.0030.596 ± 0.0110.885 ± 0.0050.616 ± 0.0170.704 ± 0.0120.504 ± 0.01
MACCSSVM0.97 ± 0.0080.906 ± 0.0250.921 ± 0.0030.75 ± 0.0110.924 ± 0.0080.76 ± 0.0250.752 ± 0.0160.536 ± 0.011
MACCSDNN0.985 ± 0.0020.955 ± 0.0070.959 ± 0.0040.872 ± 0.0120.912 ± 0.0090.726 ± 0.0280.729 ± 0.0080.529 ± 0.012
MACCSXGBoost0.977 ± 0.0050.933 ± 0.0130.92 ± 0.0050.749 ± 0.0110.925 ± 0.0050.765 ± 0.0140.757 ± 0.0170.536 ± 0.014
EnsembleRF0.923 ± 0.0050.723 ± 0.0110.892 ± 0.0050.608 ± 0.010.897 ± 0.0070.627 ± 0.0190.736 ± 0.0110.514 ± 0.01
EnsembleSVM0.99 ± 0.0040.949 ± 0.0140.936 ± 0.0040.777 ± 0.0070.938 ± 0.0060.783 ± 0.0170.772 ± 0.010.542 ± 0.011
EnsembleDNN0.997 ± 0.0010.985 ± 0.0080.977 ± 0.0130.905 ± 0.0410.927 ± 0.0030.75 ± 0.0080.76 ± 0.0040.533 ± 0.005
EnsembleXGBoost0.997 ± 0.0010.974 ± 0.0050.932 ± 0.0020.769 ± 0.0050.939 ± 0.0040.788 ± 0.0060.771 ± 0.0080.54 ± 0.008
Dataset 4 contains 3179 mPGES-1 inhibitors (the number of inhibitors in the training/test set = 2384:795); external validation set B1 contains 217 mPGES-1 inhibitors. Q: the coefficient of determination; MCC: root mean squared error, ave_std: average values and standard deviation values, Ensemble: combined the predicted probabilities of classification models built with Avalon, ECFP4, and MACCS fingerprints.
Table 5. The overall performances of classification models on Dataset 5.
Table 5. The overall performances of classification models on Dataset 5.
FingerprintsAlgorithmTraining Set5-Fold Cross-ValidationTest SetExternal Validation Set B1
Q (ave_std)MCC (ave_std)Q (ave_std)MCC (ave_std)Q (ave_std)MCC (ave_std)Q (ave_std)MCC (ave_std)
AvalonRF0.913 ± 0.0030.694 ± 0.0130.885 ± 0.0030.583 ± 0.0140.882 ± 0.0080.575 ± 0.0320.699 ± 0.0120.495 ± 0.007
AvalonSVM0.967 ± 0.010.892 ± 0.0350.911 ± 0.0030.703 ± 0.0130.904 ± 0.0060.679 ± 0.0220.712 ± 0.0080.507 ± 0.007
AvalonDNN0.99 ± 0.0070.968 ± 0.0220.965 ± 0.0070.884 ± 0.0210.902 ± 0.0090.678 ± 0.030.718 ± 0.0090.509 ± 0.01
AvalonXGBoost0.987 ± 0.0050.962 ± 0.0080.911 ± 0.0050.708 ± 0.0130.908 ± 0.0060.702 ± 0.0140.713 ± 0.0120.511 ± 0.007
ECFP4RF0.904 ± 0.0040.658 ± 0.0150.878 ± 0.0030.552 ± 0.0110.879 ± 0.0080.559 ± 0.0330.718 ± 0.0050.495 ± 0.007
ECFP4SVM0.99 ± 0.0050.969 ± 0.0180.934 ± 0.0030.78 ± 0.0090.932 ± 0.0040.773 ± 0.0160.767 ± 0.0110.538 ± 0.016
ECFP4DNN0.998 ± 0.0020.994 ± 0.0070.986 ± 0.0050.954 ± 0.0160.917 ± 0.0050.722 ± 0.0190.728 ± 0.0080.511 ± 0.008
ECFP4XGBoost0.994 ± 0.0030.981 ± 0.010.933 ± 0.0030.769 ± 0.0110.934 ± 0.0040.772 ± 0.0120.766 ± 0.010.537 ± 0.009
MACCSRF0.9 ± 0.0030.646 ± 0.0130.876 ± 0.0030.543 ± 0.0130.882 ± 0.0070.571 ± 0.030.697 ± 0.0070.489 ± 0.008
MACCSSVM0.969 ± 0.0060.896 ± 0.0190.92 ± 0.0040.731 ± 0.0120.921 ± 0.0080.733 ± 0.0280.736 ± 0.0170.512 ± 0.006
MACCSDNN0.983 ± 0.0020.945 ± 0.0060.959 ± 0.0030.861 ± 0.0110.9 ± 0.010.67 ± 0.0280.71 ± 0.0090.498 ± 0.008
MACCSXGBoost0.976 ± 0.0040.923 ± 0.0130.918 ± 0.0050.727 ± 0.010.924 ± 0.0080.745 ± 0.0250.737 ± 0.0120.517 ± 0.011
EnsembleRF0.916 ± 0.0030.676 ± 0.010.89 ± 0.0020.569 ± 0.0090.891 ± 0.0050.578 ± 0.0220.715 ± 0.0050.504 ± 0.004
EnsembleSVM0.986 ± 0.0050.929 ± 0.0170.932 ± 0.0020.748 ± 0.0070.929 ± 0.0040.738 ± 0.0170.748 ± 0.0080.529 ± 0.007
EnsembleDNN0.998 ± 0.0020.979 ± 0.0080.98 ± 0.0040.91 ± 0.0140.917 ± 0.0060.7 ± 0.020.729 ± 0.0050.516 ± 0.007
EnsembleXGBoost0.994 ± 0.0040.964 ± 0.0080.929 ± 0.0040.744 ± 0.010.931 ± 0.0040.749 ± 0.0110.748 ± 0.010.531 ± 0.009
Dataset 5 contains 3455 mPGES-1 inhibitors (the number of inhibitors in the training/test set = 2591:864); external validation set B1 contains 217 mPGES-1 inhibitors. Q: the coefficient of determination; MCC: root mean squared error, ave_std: average values and standard deviation values, Ensemble: combined the predicted probabilities of classification models built with Avalon, ECFP4, and MACCS fingerprints.
Table 6. The overall performances of regression models on Dataset 6.
Table 6. The overall performances of regression models on Dataset 6.
DescriptorsAlgorithmTraining Set5-Fold Cross-ValidationTest SetExternal Validation Set B2
R2 (ave_std)RMSE (ave_std)R2 (ave_std)RMSE (ave_std)R2 (ave_std)RMSE (ave_std)R2 (ave_std)RMSE (ave_std)
G_3DRF0.929 ± 0.0230.175 ± 0.0270.718 ± 0.210.351 ± 0.0160.731 ± 0.0520.342 ± 0.0460.615 ± 0.0520.46 ± 0.047
G_3DSVM0.864 ± 0.0410.241 ± 0.0370.737 ± 0.0190.339 ± 0.0160.75 ± 0.0480.329 ± 0.0450.634 ± 0.0480.444 ± 0.052
G_3DDNN0.927 ± 0.0290.174 ± 0.030.745 ± 0.010.331 ± 0.0080.727 ± 0.0420.354 ± 0.0320.611 ± 0.0420.472 ± 0.034
G_3DXGBoost0.92 ± 0.0350.183 ± 0.0370.746 ± 0.0280.332 ± 0.0220.741 ± 0.0440.339 ± 0.0390.625 ± 0.0440.455 ± 0.045
RDKitRF0.925 ± 0.0230.181 ± 0.0260.745 ± 0.0150.335 ± 0.0090.776 ± 0.0510.306 ± 0.0290.619 ± 0.0160.449 ± 0.037
RDKitSVM0.874 ± 0.0350.232 ± 0.0330.751 ± 0.0190.33 ± 0.0140.778 ± 0.0380.304 ± 0.0290.63 ± 0.0230.424 ± 0.019
RDKitDNN0.939 ± 0.0250.161 ± 0.0330.834 ± 0.0380.265 ± 0.0350.787 ± 0.0160.301 ± 0.0160.625 ± 0.0090.451 ± 0.031
RDKitXGBoost0.98 ± 0.0180.079 ± 0.0480.75 ± 0.0170.331 ± 0.0110.773 ± 0.0280.311 ± 0.0180.625 ± 0.0120.445 ± 0.024
Dataset 6 contains 735 mPGES-1 inhibitors with IC50 values tested in vitro by homogeneous time-resolved fluorescence assay; the number of inhibitors in the training/test set of Dataset 6 = 551:184), external validation set B2 contains 60 mPGES-1 inhibitors, R2: the coefficient of determination; RMSE: root mean squared error, ave_std: average values and standard deviation values.
Table 8. The mPGES-1 inhibitors candidates from the MFH database by virtual screening.
Table 8. The mPGES-1 inhibitors candidates from the MFH database by virtual screening.
CandidatesOriginsIC50 μM aSOM HP bEffects
cmp_B1Cannabis Sativa L. 0.88 1COX inhibition [30]
cmp_B2Ramulus Mori0.25 0.93
cmp_B3Amomum longiligularg0.34 0.93Anti-breast cancer [31]
cmp_B4Glycyrrhiza glabra L.0.80 0.93PTP1B [32] and nitric oxide (NO) inhibition [33]
cmp_B5Gardeniae Fructus0.18 0.82LOX inhibition [34]
cmp_B6Schisandra chinensis0.37 0.82UDP-glucuronosyltransferase [35], and oxLDL inhibition [36]
cmp_B7Rhizoma Dioscoreae0.55 0.82prostaglandin E2 reduction [37]
cmp_B8Ramulus Mori0.18 1
cmp_B9Mori Cortex0.19 1Tyrosinase inhibition [38]
cmp_B10Glehniae Radix0.20 1
cmp_B11Colla0.37 1Anti-inflammation [39]
cmp_B12Epimrdii Herba0.50 1Anti-inflammation [40]
cmp_B13Coicis Semen0.27 1Anti-inflammation [41]
cmp_B14Radix Puerariae0.20 1Anti-oxidation [42], andenteritis treatment [43]
cmp_B15Alisma Orientale0.29 1COX-2 inhibition [44]
a: the average predicted IC50 values by a series of optimal regression models. b: the percentage of highly active inhibitors mapped to the target molecules in the same position during the self-organized map (SOM) training. It can also be taken as the probability of being predicted by an unsupervised algorithm as a highly active mPGES-1 inhibitor.
Table 9. The interactions between potential COX-2 inhibitors and COX-2 protein obtained by docking computations.
Table 9. The interactions between potential COX-2 inhibitors and COX-2 protein obtained by docking computations.
Molecules 28 06782 i001Molecules 28 06782 i002
cmp_A1 & 5KIR a (affinity = −7.57 kcal/mol)cmp_A2 & 6BL3 (affinity = −11.75 kcal/mol)
Molecules 28 06782 i003Molecules 28 06782 i004
cmp_A3 & 4PH9 (affinity = −7.37 kcal/mol)cmp_A4 & 4PH9 (affinity = −8.53 kcal/mol)
Molecules 28 06782 i005Molecules 28 06782 i006
cmp_A5 & 6BL4 (affinity= −7.25 kcal/mol)cmp_A6 & 4PH9 (affinity = −8.68 kcal/mol)
Molecules 28 06782 i007Molecules 28 06782 i008
cmp_A7 & 4PH9 (affinity = −7.69 kcal/mol)cmp_A8 & 6BL4 (affinity = −11.24 kcal/mol)
Molecules 28 06782 i009Molecules 28 06782 i010
cmp_A9 & 6BL4 (affinity = −9.61 kcal/mol)cmp_A10 & 6BL3 (affinity = −8.92 kcal/mol)
a: The PDB index of the receptor protein during docking.
Table 10. The interactions between potential mPGES-1 inhibitors and mPGES-1 protein obtained by molecular docking.
Table 10. The interactions between potential mPGES-1 inhibitors and mPGES-1 protein obtained by molecular docking.
Molecules 28 06782 i011Molecules 28 06782 i012
cmp_B1 & 4YL0 a (affinity = −7.93 kcal/mol)cmp_B2 & 4YL1 (affinity = −9.19 kcal/mol)
Molecules 28 06782 i013Molecules 28 06782 i014
cmp_B3 & 4YL1 (affinity = −8.82 kcal/mol)cmp_B4 & 4AL1 (affinity = −8.08 kcal/mol)
Molecules 28 06782 i015Molecules 28 06782 i016
cmp_B5 & 4AL1 (affinity= −11.25 kcal/mol)cmp_B6 & 4YL0 (affinity = −7.25 kcal/mol)
Molecules 28 06782 i017Molecules 28 06782 i018
cmp_B7 & 4YL0 (affinity = −7.26 kcal/mol)cmp_B8 & 5TL9 (affinity = −9.55 kcal/mol)
Molecules 28 06782 i019Molecules 28 06782 i020
cmp_B9 & 5TL9 (affinity = −9.2 kcal/mol)cmp_B10 & 5K0I (affinity = −10.31 kcal/mol)
Molecules 28 06782 i021Molecules 28 06782 i022
cmp_B11 & 4YL0 (affinity = −6.89 kcal/mol)cmp_B12 & 4YL0 (affinity = −8.73 kcal/mol)
Molecules 28 06782 i023Molecules 28 06782 i024
cmp_B13 & 4YL0 (affinity = −7.19 kcal/mol)cmp_B14 & 4YL0 (affinity = −10.71 kcal/mol)
Molecules 28 06782 i025
cmp_B15 & 4YL0 (affinity = −7.34 kcal/mol)
a: The PDB index of the receptor protein during docking.
Table 11. Datasets used in this study.
Table 11. Datasets used in this study.
DatasetsTargets NDescriptions
Dataset 1COX-21630COX-2 inhibitors for constructing classification models, model results are shown in Table 1 and Table S1, molecules with IC50 > 10 μM are weakly active inhibitors; with IC50 < 0.1 μM are highly active inhibitors
Dataset 2COX-22925COX-2 inhibitors for constructing classification models, model results are shown in Table 2 and Table S2, molecules with IC50 > 1 μM are weakly active inhibitors; with IC50 ≤ 1 μM are highly active inhibitors
Dataset 3COX-21511COX-2 inhibitors for constructing QSAR models, model results are shown in Table 3 and Table S3, molecules with IC50 values which were tested in vitro by enzyme-linked immunoassay
External validation set A1COX-2368for evaluating the constructed classification models on COX-2 inhibitors
External validation set A2COX-2114for evaluating the constructed regression models on COX-2 inhibitors
Dataset 4mPGES-13179mPGES-1 inhibitors for building classification models, model results are shown in Table 4 and Table S4, molecules with IC50 > 10 μM are weakly active inhibitors; with IC50 < 0.6 μM are highly active inhibitors
Dataset 5mPGES-13455mPGES-1 inhibitors for building classification models, model results are shown in Table 5 and Table S5, molecules with IC50 ≥ 10 μM are weakly active inhibitors; with IC50 < 10 μM are highly active inhibitors
Dataset 6mPGES-1735mPGES-1 inhibitors for constructing QSAR models, model results are shown in Table 6 and Table S6, molecules with IC50 values which were tested in vitro by homogeneous time-resolved fluorescence assay
External validation set B1mPGES-1217for evaluating the constructed classification models on mPGES-1 inhibitors
External validation set B2mPGES-160for evaluating the constructed regression models on mPGES-1 inhibitors
N: the number of molecules in the datasets, COX-2: cyclooxygenase-2; mPGES-1: microsomal prostaglandin E2 synthase.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tian, Y.; Zhang, Z.; Yan, A. Discovering the Active Ingredients of Medicine and Food Homologous Substances for Inhibiting the Cyclooxygenase-2 Metabolic Pathway by Machine Learning Algorithms. Molecules 2023, 28, 6782. https://doi.org/10.3390/molecules28196782

AMA Style

Tian Y, Zhang Z, Yan A. Discovering the Active Ingredients of Medicine and Food Homologous Substances for Inhibiting the Cyclooxygenase-2 Metabolic Pathway by Machine Learning Algorithms. Molecules. 2023; 28(19):6782. https://doi.org/10.3390/molecules28196782

Chicago/Turabian Style

Tian, Yujia, Zhixing Zhang, and Aixia Yan. 2023. "Discovering the Active Ingredients of Medicine and Food Homologous Substances for Inhibiting the Cyclooxygenase-2 Metabolic Pathway by Machine Learning Algorithms" Molecules 28, no. 19: 6782. https://doi.org/10.3390/molecules28196782

APA Style

Tian, Y., Zhang, Z., & Yan, A. (2023). Discovering the Active Ingredients of Medicine and Food Homologous Substances for Inhibiting the Cyclooxygenase-2 Metabolic Pathway by Machine Learning Algorithms. Molecules, 28(19), 6782. https://doi.org/10.3390/molecules28196782

Article Metrics

Back to TopTop