Defining the Potential Targets for Biological Activity of Isoegomaketone Based on Network Pharmacology and Molecular Docking Methods

Isoegomaketone is a water-soluble natural ketone compound that is commonly present in Rabdosia angustifolia and Perilla frutescens. At present, it is known that isoegomaketone has a wide range of pharmacological activity, but there has been no thorough investigation of its potential targets. As a result, we examined the potential targets of isoegomaketone using the network pharmacology approach. In our study, the TCM Database@Taiwan was utilized to search for the chemical formula. The pharmacological characteristics of isoegomaketone were then evaluated in silico using the Swiss Absorption, Distribution, Metabolism, and Excretion (Swiss ADME) and Deep Learning–Acute Oral Toxicity (DL-AOT) methods, and the potential isoegomaketone target genes were identified using a literature study. Additionally, using the clusterProfiler R package 3.8.1, the Gene Ontology (GO) enrichment analysis and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of target genes were performed. In order to obtain the protein interaction network, we simultaneously submitted the targets to the STRING database. After this, we performed molecular docking with respect to targets and isoegomaketone. Finally, we created visual networks of protein–protein interactions (PPI) and examined these networks. Our results showed that isoegomaketone had good drug-likeness, bioavailability, medicinal chemistry friendliness, and acceptable toxicity. Subsequently, through the literature analysis, 48 target genes were selected. The bioinformatics analysis and network analysis found that these target genes were closely related to the biological processes of isoegomaketone, such as atherosclerotic formation, inflammation, tumor formation, cytotoxicity, bacterial infection, virus infection, and parasite infection. These findings show that isoegomaketone may interact with a wide range of proteins and biochemical processes to form a systematic pharmacological network, which has good value for the creation and use of drugs.


Introduction
Traditional Chinese medicine contains a large number of active drug molecules [1]. Through the mining of traditional Chinese medicine, the development of drugs can be greatly enriched [2]. Isoegomaketone is a water-soluble natural ketone compound [3]. It is commonly present in traditional Chinese medicines Rabdosia angustifolia and Perilla frutescens [4]. Recent studies have shown that isoegomaketone has various types of biological activity, such as anti-inflammatory, anti-tumor, anti-rheumatoid arthritis, etc. [5][6][7][8][9]. Thus, isoegomaketone is considered to be an important active component of Perilla frutescens and has great drug development potential. However, the potential targets of isoegomaketone and its potential for further development are not yet fully understood. At present, the extraction of potential targets from existing research and the use of computer methods to predict the potential development directions has become the main method [10][11][12]. Compared with traditional methods, this method provides more convenient technology and a clearer direction for drug design and development.
As a result, we applied the network pharmacology approach to comprehensively examine the pharmacological effects of isoegomaketone. Firstly, the Swiss Absorption, Distribution, Metabolism, and Excretion (Swiss ADME) and Deep Learning-Acute Oral Toxicity (DL-AOT) methods were used to evaluate in silico the drug properties of isoegomaketone. Furthermore, we provided candidate target genes through literature analysis. Additionally, Gene Ontology (GO) enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, and molecular docking were performed on these potential target genes. Finally, we thoroughly demonstrated the probable targets of the compound by creating the isoegomaketone pharmacological interaction network.

Molecular Formula and In Silico Drug Properties of Isoegomaketone
An important database resource that is free to use and readily available is the TCM Database@Taiwan (http://tcm.cmu.edu.tw/ (accessed on 25 March 2022)), which houses sources of chemical and pharmacological information on isoegomaketone [13]. The Swiss ADME (http://www.swissadme.ch/ (accessed on 25 March 2022)), which supports drug development, enables researchers to calculate the physicochemical descriptors and predict the ADME parameters, pharmacokinetic features, drug-like nature, and medicinal chemistry friendliness of one or more small compounds [14]. Thus, in our study, the in silico drug properties of isoegomaketone were analyzed using the Swiss ADME.

Target Gene Screening for Isoegomaketone
All candidate potential targets are derived from literature analysis. Literature analysis was carried out through various academic databases, with the keyword "isoegomaketone" [5,6,[16][17][18][19][20][21][22][23][24]. Finally, eleven published, peer-reviewed pharmacological studies on isoegomaketone were selected. All eleven selected pharmacological studies on isoegomaketone demonstrated high research quality and were all carried out via in vivo or in vitro methods.

Analysis of PPI Network
One of the online databases that compiles all publicly available sources of knowledge of PPI is STRING 11.0 (https://string-db.org/ (accessed on 29 March 2022)), which enables users to supplement the knowledge already known about PPI with computational predictions [25]. We uploaded 48 putative targets of isoegomaketone to the STRING database in order to build a PPI network, with the species defined as "Homo sapiens" and the minimum interaction score set to 0.4. After this, the outcomes were imported into Cytoscape 3.7.2 for visual evaluation. Cytoscape, an open-source software platform for visualizing complex networks, can calculate the parameters of each node in the network diagram, such as the degree, betweenness centrality (BC), and closeness centrality (CC) [26]. We utilized the cytoHubba plug-in topological method to identify the significant protein nodes and subnetworks in the network, after filtering the target nodes by the matching median values of degree values, BC and CC in the PPI network [27].

Gene Function and Pathway Enrichment Analysis
GO functional annotation analysis is a standard method used to carry out large-scale functional enrichment analyses of genes, and it comprises biological processes (BP), molecular functions (MF), and cellular components (CC) [28]. It is possible to assign functional meanings to genes and genomes at the molecular and higher levels by adopting the widely used Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway [29]. Using the clus-Life 2022, 12, 2115 3 of 15 terProfiler R package 3.8.1 with FDR < 0.2 (FDR, false discovery rate) and p value < 0.05, KEGG pathway analysis and GO analysis were carried out, and the most important targets engaged in the pertinent biological processes were examined.

Compound-Target Molecular Docking
The 2D structure of isoegomaketone, downloaded from PubChem (https://pubchem. ncbi.nlm.nih.gov/ (accessed on 26 March 2022)) in mol2 format, was imported into the AutoDockTools 1.5.6 software. After checking its spatial structure, adding atomic charges, assigning atomic types, and changing all flexible bonds to rotatable by default, it was saved in pdbqt format as a docking ligand. The 3D structures of target-gene-associated proteins were downloaded from the Uniprot database (https://www.uniprot.org/ (accessed on 26 March 2022)) and the pdb database (https://www.rcsb.org/ (accessed on 26 March 2022)), and they were then saved in pdb format. AutoDockTools 1.5.6 software was used to carry out the removal of origin ligands, removal of water molecules, hydrogenation, charge addition, and various optimizations. The grid box was set as the default value and the final optimized proteins were saved in pdbqt format as docking acceptors. Molecular docking was performed using AutoDock Vina 1.1.2. The protein structure was set during molecular docking to be a rigid macromolecule, and the Genetic Algorithm Parameters algorithm was used. PyMOL was used to illustrate the outcomes for the group in which each protein had the lowest binding energy.

Molecular Formula and In Silico Drug Properties of Isoegomaketone
We obtained the chemical formula for isoegomaketone using the TCM Database@Taiwan database ( Figure 1). Through Swiss ADME, we obtained essential ADME-related data on isoegomaketone, and the drug-likeness as well as the medicinal chemistry of isoegomaketone were subsequently evaluated via the classic formula. The results are provided in Table 1. Regarding the drug-likeness, lead-drug likeness, and medicinal chemistry friendliness, the pharmacokinetic characteristics of isoegomaketone were consistent with Lipinski's rule of five, the Ghose filter, Veber's Rule, and the RO (3) rule, and the synthetic accessibility was 2.9. Definitions of these rules and parameters can be found in References [30][31][32][33][34]. molecular functions (MF), and cellular components (CC) [28]. It is possible to assign functional meanings to genes and genomes at the molecular and higher levels by adopting the widely used Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway [29]. Using the clusterProfiler R package 3.8.1 with FDR < 0.2 (FDR, false discovery rate) and p value < 0.05, KEGG pathway analysis and GO analysis were carried out, and the most important targets engaged in the pertinent biological processes were examined.

Compound-Target Molecular Docking
The 2D structure of isoegomaketone, downloaded from PubChem (https://pubchem.ncbi.nlm.nih.gov/ (accessed on 26 March 2022)) in mol2 format, was imported into the AutoDockTools 1.5.6 software. After checking its spatial structure, adding atomic charges, assigning atomic types, and changing all flexible bonds to rotatable by default, it was saved in pdbqt format as a docking ligand. The 3D structures of target-gene-associated proteins were downloaded from the Uniprot database (https://www.uniprot.org/ (accessed on 26 March 2022)) and the pdb database (https://www.rcsb.org/ (accessed on 26 March 2022)), and they were then saved in pdb format. AutoDockTools 1.5.6 software was used to carry out the removal of origin ligands, removal of water molecules, hydrogenation, charge addition, and various optimizations. The grid box was set as the default value and the final optimized proteins were saved in pdbqt format as docking acceptors. Molecular docking was performed using AutoDock Vina 1.1.2. The protein structure was set during molecular docking to be a rigid macromolecule, and the Genetic Algorithm Parameters algorithm was used. PyMOL was used to illustrate the outcomes for the group in which each protein had the lowest binding energy.

Molecular Formula and in Silico Drug Properties of Isoegomaketone
We obtained the chemical formula for isoegomaketone using the TCM Database@Taiwan database ( Figure 1). Through Swiss ADME, we obtained essential ADME-related data on isoegomaketone, and the drug-likeness as well as the medicinal chemistry of isoegomaketone were subsequently evaluated via the classic formula. The results are provided in Table 1. Regarding the drug-likeness, lead-drug likeness, and medicinal chemistry friendliness, the pharmacokinetic characteristics of isoegomaketone were consistent with Lipinski's rule of five, the Ghose filter, Veber's Rule, and the RO (3) rule, and the synthetic accessibility was 2.9. Definitions of these rules and parameters can be found in References [30][31][32][33][34].   The LD50 value of isoegomaketone was predicted by DL-AOT, as shown in Table 2. The result showed that the LD50 value of isoegomaketone was 31.2 mg/kg, while the evaluation indicated "caution".

No.
Gene ID Gene Name Gene Description

Construction of the PPI Network and Screening of Key Targets
In order to obtain the PPI network association, we imported the 48 target genes into the STRING database and chose Homo sapiens as the organism. The free nodes were then eliminated, and the data were loaded into Cytoscape 3.7.2 for visualization (Figure 2a). The circle's size and color varied with the degree value, while the width and color of the edges varied with the total score. The image shows 48 targets and 485 interactions, indicating that the more projected disease-related targets there are, the more likely it is that there are successful interactions between those targets. Calculated by the Network Analyzer plug-in in Cytoscape, the average closeness centrality (CC) was 0.63, the average betweenness centrality (BC) was 0.01, and the average degree value (DV) of the node was 19.1. Twelve of them had DV, BC, and CC values that were higher than the norm, while AKT1, TP53, JUN, MAPK8, CASP3, IL6, MTOR, and MAPK14 all displayed high connection degrees that were greater than 30. These could be the key target proteins that are the focus of isoegomaketone. The subnetwork of these key targets is shown in Figure 2b. Table S1 displays the specific node properties in detail.

KEGG Enrichment Analysis
We performed KEGG enrichment analysis of these target genes concurrently. KEGG pathway enrichment analysis revealed that 48 possible target genes were enriched and that there were significant relationships between the target genes and 194 signal pathways (FDR ≤ 0.05). The top 20 pathways with the highest enrichment ratios are demonstrated in Figure 4. The pathway with the largest count and lowest p value was the lipid and atherosclerosis pathway.

Results of Molecular Docking
In general, the likelihood of interaction decreases with the energy of conformational stability of the ligand and receptor. There is yet no single standard for the target screening of active compounds, and, at present, binding energy ≤ −5.0 kcal/mol is typically used as the basis for screening. The findings of molecular docking revealed that 6 of the 12 chosen target proteins exhibited an affinity for isoegomaketone that was less than −5.0 kcal/mol. All docking results are shown in Table 4 and Figure 5.

KEGG Enrichment Analysis
We performed KEGG enrichment analysis of these target genes concurrently. KEGG pathway enrichment analysis revealed that 48 possible target genes were enriched and that there were significant relationships between the target genes and 194 signal pathways (FDR ≤ 0.05). The top 20 pathways with the highest enrichment ratios are demonstrated in Figure 4. The pathway with the largest count and lowest p value was the lipid and atherosclerosis pathway.

Results of Molecular Docking
In general, the likelihood of interaction decreases with the energy of conformational stability of the ligand and receptor. There is yet no single standard for the target screening of active compounds, and, at present, binding energy ≤ −5.0 kcal/mol is typically used as

Discussion
Currently, network pharmacology is being used with increasing frequency in the development and utilization of novel drugs. It converts the "single target, single drug" concept of drug research into a "network target, multicomponent treatment" model [35]. As a result, in this study, we built a multidimensional network through target prediction and protein interaction networks based on the network pharmacology and molecular docking research methodologies to clarify the probable targets of isoegomaketone.
Drug-likeness is recognized as a guideline for chemical library screening optimization in the area of drug development [36]. Lipinski's rule of five, the Ghose filter, and Veber's Rule are some of the most widely used rules of thumb. The main purpose of these tests is to determine whether a substance may be used as a drug or whether a substance, having pharmacological or biological activity, can be converted into an oral drug for human use. All of these formulas evaluate specific values of a given compound, including logP, TPSA, MW, MR, nHBA, nRotB, and nHBD. According to Table 1, isoegomaketone, with MW ≤ 500, logP < 5, nRotB ≤ 10, nHBD ≤ 5, nHBA ≤ 10, MR ≤ 140, and TPSA < 140 Å, perfectly passes the test of Lipinski's rule of five, the Ghose filter, and Veber's rule. Typically, a molecule that satisfies two of the three rules can be considered to fit the bioavailability of drugs. Thus, this study strongly indicates that isoegomaketone has excellent bioavailability in drugs, making it plausible that it is an active oral substance in humans. In addition, according to Table 1, isoegomaketone was predicted to have high gastrointestinal absorption, and it was also predicted to be able to permeate the blood-brain barrier. The LogS value for isoegomaketone is −2.39, which indicates its solubility in water. Moreover, isoegomaketone complies with the RO (3) rule and has good synthetic accessibility, which suggests its potential as a lead drug. The LD50 value of isoegomaketone was predicted to be 31.2 mg/kg, as shown in Table 2. According to DL-AOT, compounds are divided into four categories, namely "danger/poison", "warning", "caution", and "none required", and isoegomaketone is placed in the "caution" group, meaning that it is predicted to have no apparent toxicity to normal cells. All of these results indicate that isoegomaketone is likely to be a good candidate for drug development.
As indicated in Table 3, 48 possible isoegomaketone targets were found in the literature analysis. The STRING database was then used to examine the 48 possible targets, and there were 12 nodes whose degree value, betweenness centrality, and closeness centrality were all greater than the average, most of which (such as AKT1, TP53, JUN, MAPK8, CASP3, IL6, MTOR, MAPK14) are involved in the process of apoptosis and the stress response. These selected isoegomaketone targets suggest that isoegomaketone may have potential anti-inflammatory, anti-tumor, and cytotoxic activity. We also performed molecular docking to further investigate any potential interactions between isoegomaketone and these targets. With a decrease in the conformational energy of the ligand binding to the receptor, there is a higher chance of interaction. If −5.0 kcal/mol is utilized as the screening criterion, isoegomaketone may interact directly with these targets, since 6 out of the 12 targets had binding energies that were less than −5.0 kcal/mol [37][38][39][40].
Most of them are associated with anti-inflammatory, anti-tumor, and cytotoxic responses. They might, thus, be the main targets of isoegomaketone's pharmacological activity.
To elucidate the potential biological activity of isoegomaketone based on gene functions and signaling pathways, we further conducted GO analysis and KEGG pathway analysis of these key targets with the clusterProfiler R package 3.8.1. According to the analysis of BP items, isoegomaketone is closely associated with the following: neuron death, regulation of neuron death, cellular response to chemical stress, aging, response to radiation, cellular response to abiotic stimulus, etc. The results of KEGG pathway enrichment analysis show that the lipid and atherosclerosis pathway is highly enriched. It should be mentioned that no previous studies have clarified the relationship between isoegomaketone and atherosclerosis. The result indicates that isoegomaketone may play a comprehensive role in preventing and treating atherosclerosis, as pathways such as the lipid and atherosclerosis pathway and fluid shear stress and atherosclerosis pathway are involved. In addition, the anti-inflammatory activity of isoegomaketone may also be involved. Further in vitro and in vivo experimental validation could be carried out. Moreover, we found that 9 of the top 20 enriched KEGG pathways were related to bacterial, viral, and parasite infections, such as hepatitis B, Kaposi sarcoma-associated herpesvirus, Epstein-Barr virus, human cytomegalovirus, human immunodeficiency virus 1, tuberculosis, Chagas disease (American trypanosomiasis), measles, and toxoplasmosis. According to these findings, isoegomaketone could possess antiviral, antibacterial, and parasitic action. Some pathways are related to inflammatory diseases-for example, the AGE-RAGE signaling pathway in non-alcoholic fatty liver disease and diabetic complications-which is in accordance with a previous in vivo study [9]. In addition, the KEGG enrichment results also include the TNF signaling pathway, Toll-like receptor signaling pathway, C-type lectin receptor signaling pathway, and other immuno-inflammatory-response-related pathways. As a result, our research revealed that isoegomaketone may be involved by preventing an overactive immune response, an influx of inflammatory substances, and the growth of tumor cells. These findings are similar to some of the in vitro studies of isoegomaketone that have reported its anti-inflammatory, anti-tumor, and cytotoxic activity [41]. Therefore, on one hand, it suggests the accuracy and rationality of our GO and KEGG research. On the other hand, the result further reveals the anti-inflammatory, anti-tumor, and cytotoxic activity target of isoegomaketone, suggesting the potential of isoegomaketone as an anti-inflammatory, anti-tumor, and cytotoxic drug. Taking these results together, we created a drug-target pathway network diagram that more obviously showed that isoegomaketone could exert a variety of pharmacological effects.
However, there are limitations to our research. The most significant is that, as mentioned above, our target screening was obtained through a large number of previously published studies. These provide high-quality, real-world evidence, based, in part, on in vivo and in vitro studies [5,6,[16][17][18][19][20][21][22][23][24]. However, although our study provides some prospective evidence through the in silico method, some results have not yet been suitably empirically validated. Therefore, real-world validation using in vivo or in vitro methods based on the current prospective results is highly desirable.

Conclusions
Our research revealed that isoegomaketone has a wide range of pharmacological effects. At the same time, we investigated isoegomaketone's potential targets, which may be used to further create secure and efficient anti-atherosclerosis, anti-inflammatory, antitumor, cytotoxic, antiviral, antibacterial, and anti-parasite medications. Our work offers a new perspective on the study, creation, and clinical use of isoegomaketone.