Next Article in Journal
A Review on Data-Driven Process Monitoring Methods: Characterization and Mining of Industrial Data
Next Article in Special Issue
Investigating the Molecular Mechanism of Qianghuo Shengshi Decoction in the Treatment of Ankylosing Spondylitis Based on Network Pharmacology and Molecular Docking Analysis
Previous Article in Journal
Recent Advances on Waste Plastic Thermal Pyrolysis: A Critical Overview
Previous Article in Special Issue
Implementation of System Pharmacology and Molecular Docking Approaches to Explore Active Compounds and Mechanism of Ocimum Sanctum against Tuberculosis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Promising Mechanisms of Low Molecular Weight Compounds of Panax Ginseng C.A. Meyer in Alleviating COVID-19: A Network Pharmacology Analysis

Department of Bio-Health Convergence, College of Biomedical Science, Kangwon National University, Chuncheon 24341, Korea
*
Author to whom correspondence should be addressed.
Processes 2022, 10(2), 333; https://doi.org/10.3390/pr10020333
Submission received: 11 January 2022 / Revised: 1 February 2022 / Accepted: 7 February 2022 / Published: 9 February 2022
(This article belongs to the Special Issue Network Pharmacology Modelling for Drug Discovery)

Abstract

:
Panax Ginseng C.A. Meyer (PGCAM) is a well-known phytomedicine, but most of its compounds, such as ginsenoside derivatives, have poor absorption and bioavailability profile due to high molecular weight (≥500 Daltons), which is the major hurdle for their clinical application. Hence, this research explored the efficiency of low molecular weight compounds (LMWCs) (<500 Daltons) screened from PGCAM and their anti-COVID-19 mechanisms through network pharmacology. Molecular compounds from PGCAM were identified using public databases and filtered out by the drug-likeness evaluation. Genes interacted with these filtered compounds, and COVID-19-related genes were extracted from public databases. In addition, overlapping genes between compounds and interactive genes were identified using the Venn diagram. In parallel, the networking between compounds and overlapping genes was analyzed by RStudio. The pathway enrichment analysis of overlapping genes was determined by STRING. Finally, the key bioactive compounds were documented through virtual screening. The bubble chart suggested that the mechanisms of PGCAM against COVID-19 were related to 28 signaling pathways. The key molecular anti-COVID-19 mechanisms might be the anti-inflammation, anti-permeability, and pro-apoptosis by inactivating the PI3K-Akt signaling pathway. The six key genes and the five compounds related to the PI3K-Akt signaling pathway were RELA-paeonol, NFKB1-frutinone A, IL6-nepetin, MCL1-ramalic acid, VEGFA-trifolirhizin, and IL2-trifolirhizin. The docking between these key genes and compounds demonstrated promising binding affinity with a good binding score. Overall, our proposed LMWCs from PGCAM provide a fundamental basis with noteworthy pharmacological evidence to support the therapeutic efficacy of PGCAM in relieving the main symptoms of COVID-19.

1. Introduction

Due to the intense efforts made by clinicians and researchers, only one human infection, “smallpox”, has been successfully eradicated from the world [1]. However, the recent outburst of a new type of pneumonia caused by a severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) [2] is not likely to disappear from the globe completely, even after the quick development of a sustainable vaccine. This indicates one certain thing: that the world needs novel treatment plans and possible therapeutic interventions to combat the highly infectious disease (COVID-19), no matter what else happens next.
In late December 2019, a life-threatening contagious disease with unknown etiology was first reported in Wuhan, Hubei province, China [3]. Although infectious disease experts identified the new coronavirus (SARS-CoV-2) spread from the animal host, concerns remained over how the first person might have been infected [4]. Viral genomics and epidemiological studies evidenced that the animal-to-human shifted SARS-CoV-2 can rapidly transmit from one person to another through interaction or respiratory droplets [5,6]. Hence, the World Health Organization (WHO) initially declared a global emergency on 22 January 2020, as the Wuhan coronavirus spread. On 20 August 2020, the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University reported a health crisis with intensifying global case numbers (22,260,000) and fatalities [7,8]. COVID-19 predominantly attacks the respiratory system, even though other organ systems are compromised. Lower respiratory tract infection (LRTI) associated with fever, shortness of breath, and dry cough were reported as the symptoms with the first case from Wuhan, China [9]. Additionally, muscle aches, chills, sore throat, runny nose, headache, and chest pain were observed [10]. Currently, COVID-19 symptoms are similar to pneumonia, though a significant number of COVID-19 infected cases are asymptomatic, which can affect others as silent spread [11]. Importantly, patients with underlying diseases, such as asthma, diabetes, and hypertension, are more susceptible to COVID-19 [12,13].
A report suggests that around 40% of those infected with SARS-CoV-2 might remain asymptomatic and thus, the transmission might carry a greater risk than expected [11]. At present, endeavors to develop medications for COVID-19 by repositioning existing drugs have been pursued, as well as vaccine development. However, the efforts have been hampered by the lack of scientific evidence on how COVID-19 attacks human hosts [14]. On the other hand, phytochemicals have always played a pivotal role in providing potential therapy against diverse diseases [15], including fighting a mutant coronavirus [16]. From folk remedy times, medicinal plants provided many outstanding bioactive compounds, which led to great drug discoveries, such as aspirin from Salix, taxol from Taxus brevifolia, and artemisinin from Artemisia annua [17,18,19]. During the 2003 SARS outbreak, a series of Chinese herbal treatments were effectively implemented to control the pandemic situation [20], and even now, the Chinese government impulses traditional remedies aimed at COVID-19 prevention. According to the state administration of traditional Chinese medicine (TMC), 85.2% of the country’s COVID-19 patients (total: 60,107) until 17 February 2020 had been recovered with the treatment of TMC [21]. Although researchers race to develop a vaccine, the rigorous validation of traditional medicine for identifying the wide-spectrum antiviral active components can also be a viable way to control the COVID 19 disease. Therefore, an attempt to explore novel antiviral compounds with the underlying mechanism of action should not be halted.
Korean traditional medicine literature, “Donguibogam”, established about 400 years ago, illustrates that Panax Ginseng C.A. Meyer (PGCAM) is a very potent herbal plant with multiple pharmacological activities, particularly strengthening the immune system against some pathogens [22]. According to an experiment, ginseng interrupted cellular oxidative damage induced by a respiratory syncytial virus (RSV) and downregulated pro-inflammatory gene expression level induced by RSV in the human alveolar epithelial cell [23]. Additionally, potential bioactive compounds of ginseng extracts, including ginsenoside derivatives (Rg1, Rg3, Rb1, Rb2, Rb3, Re, Rd, Rh2) and compound K, are renowned for therapeutic efficacy against several infectious diseases [24]. However, the major drawback of these compounds is their high molecular weight and low bioavailability; even the absorption and transportation rate of intact ginsenosides and their metabolites from the intestines are very poor [25]. In addition, some reports studied via Caco-2 permeability assay indicated the low permeability of ginsenoside derivatives (<2.59 ± 0.17·10−7 cm/s) and compound K (<8.65·10−7 cm/s) [26,27]. It is evident that low permeable compounds have poor bioavailability profiles, which do not reach the intended biological system; thus, therapeutic efficacy is not well consistent [28]. Therefore, these findings suggest that the pharmacological value of low molecular weight compounds (LMWCs) in ginseng should be strengthened to provide therapeutic evidence against immune dysfunction, especially in treating COVID-19 infection during an urgent pandemic outbreak.
Network pharmacology (NP), a new paradigm and combined analytical system, can efficiently explore the interaction networking of diverse factors, such as ligands, protein targets, diseases, and genes [29]. NP can unveil the target prediction and mechanism of drug action with a multidimensional perspective, which highlights shifting of the approach from “one target, one ligand” to “multiple targets, multivariate therapeutics” [30]. Importantly, it widely contributes to drug discovery by providing 40% of the recent active candidates with successful clinical evidence [31]. Currently, network pharmacology has been utilized broadly to decipher bioactive compounds and mechanisms of phytochemicals against diverse diseases [32]. In this research, we implemented network pharmacology to analyze LMWCs and mechanisms of PGCAM against COVID-19. Firstly, LMWCs in PGCAM were identified using browsing literature and public database and then filtered through Lipinski’s rule for drug-likeness compounds. Afterwards, genes that interacted with selected compounds or COVID-19-related genes were screened using public databases, and the overlapping genes between compounds and COVID-19 target genes were identified. Thirdly, KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analysis of overlapping genes was employed to explore the molecular mechanisms of PGCAM against COVID-19. Finally, the potential bioactive compounds and key genes of PGCAM against COVID-19 were selected by analyzing the binding affinity energy in virtual mode. The overall workflow is presented in Figure 1.

2. Results

2.1. Potential Low Molecular Weight Bioactive Compounds from PGCAM

The public databases identified a total of 121 compounds in PGCAM, and the drug-likeness properties of 107 compounds were selected by Lipinski’s rule, which is presented in Table 1.

2.2. Target Genes Associated with the 107 Compounds or COVID-19

As shown in Supplementary Table S1, a total number of 1197 genes (STP + SEA) related to 107 compounds was identified and, among them, the overlapped 399 genes between SEA and STP were screened (Figure 2). The 399 genes (Supplementary Table S2) were again analyzed with COVID-19-related genes (356 genes) (Supplementary Table S3). Finally, the Venn diagram (Figure 3) showed the 27 overlapping genes (Supplementary Table S4) between the 399 genes and the 356 genes. Overall, the final 27 genes were selected from 107 LMWCs of PGCAM through COVID-19-related genes.

2.3. Potential Molecular Pathways of PGCAM against COVID-19

The gene–gene networking analysis through overlapping 27 genes was visualized by STRING, which includes 24 nodes and 122 edges (Figure 4). The 3 (PTPRS, CPA3, and S1PR3) out of 27 genes were not associated with each other. The KEGG pathway enrichment analysis via STRING suggested that 27 overlapping genes were significantly enriched in 28 signaling pathways (Figure 5). Detailed descriptions of the 28 signaling pathways are listed in Table 2. Among them, both RELA and NFKB1 connected directly to 21 signaling pathways out of 28 signaling pathways and were noted as hub genes of PGCAM against the COVID-19 (Figure 6). Coincidently, both RELA and NFKB1 play vital roles in all 21 signaling pathways by the PI3K-Akt signaling pathway, indicating that the PI3K-Akt signaling pathway might be a hub signaling pathway of PGCAM against COVID-19.

2.4. Key Genes and Bioactive Compounds of PGCAM against COVID-19

Based on the PI3K-Akt signaling pathway, six key genes (RELA, NFKB1, IL6, MCL1, VEGFA, and IL2) and 24 potential key bioactive compounds were identified in PGCAM against COVID-19. The 24 compounds are associated with six key genes on PI3K-Akt signaling pathway. The 24 compounds (Figure 7) are classified as follows: eight aliphatic acyclic compounds (NN-dimethyldecanamide, linoleic acid, methyl palmitelaidate, methyl linoleate, methyl myristate, methyl stearate, methyl margarate, and methyl pentadecanoate), eight aromatic heteropolycyclic compounds (frutinone A, nepetin, L-adenosine, adenosine, N9-formylharman, inermin, trifolirhizin, and kaempferol), five aromatic homomonocyclic compounds (paeonol, p-hydroxycinnamic acid, dibuthyl phthalate, salicylic acid, and ramalic acid), one aliphatic heteromonocyclic compound (beta-D-mannose), one aliphatic heteropolycyclic compound (darutoside), and one aromatic heteromonocyclic compound (p-glucosyloxymandelonitrile). The interactions between each gene and potential active compound were visualized with RStudio (Supplementary Figure S1).

2.5. Affinity Score on Six Key Genes of Potential Bioactive Compounds in PGCAM against COVID-19

The potential bioactive compounds were docked against six genes (RELA, NFKB1, IL6, MCL1, VEGFA, IL2) to measure the binding energy. The molecular docking figure is presented in Figure 8.
Docking simulation of the affinity between R1-R2 and RELA protein in the “Homo sapiens” setting was analyzed. Based on the docking score, the order of priority of binding energy is given as: R2, R1. The two-binding energy of R1-RELA and R2-RELA demonstrated −3.9 and −5.3 kcal/mol, respectively. Interaction analysis of the best-docked compound, namely “Paeonol”, resulted in four hydrogen bonds (Arg-174, Val-163, Thr-164, Arg-95) and seven hydrophobic bonds (Pro-176, Leu-175, Pro-177, Gln-162, Arg-73). The detailed information is listed in Table 3. Docking simulation of the affinity between N1-N4 and NFKB1 protein in the “Homo sapiens” setting revealed promising binding affinity, and the order of priority of binding energy is as follows: N1, N3, N2, N4. The binding energy of N1-NFKB1, N2-NFKB1, N3-NFKB1, and N4-NFKB1 exposed −7.2, −5.2, −5.6, and −4.6 kcal/mol, respectively. Interaction analysis of the best-docked compound, namely “Frutinone A”, resulted in one hydrogen bond (Ser-113) and seven hydrophobic bonds (Val-145, Val-61, Leu-143, Arg-157, Ala-156, Thr-153, Lys-149). The detailed information is listed in Table 4. Docking simulations of the affinity between L1-L3 and IL6 protein in the “Homo sapiens” setting showed good binding energy, and the order of priority is given as: L3, L2, L1. The five-binding energy of L1-IL6, L2-IL6, and L3-IL6 were −5.0, −5.7, and −8.0 kcal/mol, respectively. Interaction analysis of the best-docked compound, namely “Nepetin”, resulted in two hydrogen bonds (Glu-165, Gln-166) and seven hydrophobic bonds (Lys-103, Phe-83, Glu-81, Ser-168, Pro-80, Ile-106, Pro-40). The detailed information is listed in Table 5. Docking simulation of the affinity between M1-M8 and MCL1 protein in the “Homo sapiens” setting was investigated, and the order of priority of binding energy is as follows: M7, M1, M2 = M8, M3, M6, M5, M4. The eight-binding energy of M1-MCL1, M2-MCL1, M3-MCL1, M4-MCL1, M5-MCL1, M6-MCL1, M7-MCL1, and M8-MCL1 unveiled −6.8, −6.5, −5.3, −4.2, −4.4, −4.7, −6.9, and −6.5 kcal/mol, respectively. Interaction analysis of the best-docked compound, namely “Ramalic acid”, resulted in three hydrogen bonds (Asp-313, Arg-201, Arg-310) and five hydrophobic bonds (Gly-311, Glu-284, Thr-191, Glu-288, Arg-187). The detailed information is listed in Table 6.
Docking simulation of the affinity between V1-V11 and VEGFA protein in the “Homo sapiens” setting was explored, and the order of priority of binding energy is as follows: V4, V3, V2 = V6 = V7, V5, V9, V1, V11, V8 = V10. The eleven-binding energy of V1-VEGFA, V2-VEGFA, V3-VEGFA, V4-VEGFA, V5-VEGFA, V6-VEGFA, V7-VEGFA, V8-VEGFA, V9-VEGFA, V10-VEGFA, and V11-VEGFA revealed −5.0, −7.1, −7.5, −8.1, −6.8, −7.1, −7.1, −3.9, −5.1, −3.9, and −4.4 kcal/mol, respectively. Interaction analysis of the best-docked compound, namely “Trifolirhizin”, resulted in two hydrogen bonds (Asp-41, Glu-42) and five hydrophobic bonds (Asn-75, Glu-38, Tyr-39, Ser-95, Leu-97). The detailed information is listed in Table 7.
Finally, docking simulation of the affinity between P1-P6 and IL2 protein in the “Homo sapiens” setting was investigated, and the order of priority of binding energy is as follows: P2, P5, P1 = P4, P3, P6. The six-binding energy of P1-IL2, P2-IL2, P3-IL2, P4-IL2, P5-IL2, and P6-IL2 indicated −7.0, −8.5, −6.3, −7.0, −7.4, and −5.2 kcal/mol, respectively. Interaction analysis of best-docked compound, namely “Trifolirhizin”, resulted in six hydrogen bonds (Asp-68, Ser-129, Ile-128, Ala-131, Arg-15, Gln-70) and six hydrophobic bonds (Ser-69, Asp-84, His-79, Leu-80, Ser-132, His-133). The detailed information is listed in Table 8.
This result illustrated that the five bioactive compounds might be anti-inflammatory agents against COVID-19 (Figure 9).

3. Discussion

Compounds–genes network indicated that the clinical efficacy of PGCAM on COVID-19-associated genes was directly involved in 27 genes. Of these, both RELA and NFKB1 genes were more significant genes than any other kind of genes. On top of that, based on the pathway enrichment, the PI3K-Akt signaling pathway is found as the topmost mechanism of PGCAM against COVID-19. Accordingly, the docking score through virtual screening on the six genes related to the PI3K-Akt signaling pathway suggested that five bioactive compounds (paeonol, frutinone A, nepetin, ramalic acid, trifolirhizin) were considered the most significant compounds of PGCAM against COVID-19. Meanwhile, the consequence of the KEGG pathway enrichment analysis of 27 genes might be an anti-COVID mechanism. The relationships of 28 signaling pathways with anti-virus properties were discussed as follows.
PPAR signaling pathway: Triggered PPAR-α regulates immunity of cytoplasmic DNA, blocks interferon production, and boosts vulnerability to viral infection [33].
MAPK (Mitogen-Activated Protein Kinase) signaling pathway: Virus-infected cells cause severe reactions that activate the MAPK signaling cascades to shield from virus attack. They can also be used by viruses to facilitate viral replication [34]. Furthermore, the p38 MAPK signaling pathway is a potent inflammatory pathway in lung and heart damage in COVID-19 patients [35].
ErbB signaling pathway: A great number of prime oncogenic virus infections are associated with ErbB signaling pathway. The Hepatitis B Virus(HBV) also induces overexpression of the Epidermal Growth Factor Receptor (EGFR) gene and protein [34,36,37,38].
RAS signaling pathway: Inhibition of RAS can curtail tissue damage severity in COVID-19 patients. Additionally, ACE (Angiotensin-converting enzyme) blockers cause a reduction in the response of the RAS system [39].
cAMP signaling pathway: cAMPs might have a viral-resistant effect by curbing viral replication and interrupting viral entry. On the other hand, they might have an inimical role by decreasing HIV antiviral immune responses, thus halting the removal of the virus and causing T cell malfunction [40].
Chemokine signaling pathway: Elevation of chemokines (CXCL10 and CXCL8) is clinically characteristic of COVID-19 infection. Moreover, its expression level indicated different clinical acuteness among the COVID-19 patients [41].
NF-κB signaling pathway: Activation of NF-κB signaling pathway shows acute inflammation induced after SARS-CoV attack, and NF-κB blockers are favorable antivirals in the disorder caused by COVID-19 [42].
HIF-1 signaling pathway: HIF-1α expression supports SARS-CoV-2 replication and sustainability, and it intensifies the cytokine production level in monocytes [43].
FOXO signaling pathway: Reduction of FOXO3 in T cells interrupted apoptosis, elevated multifunction of CD8 cells, and enhanced viral control [44].
Sphingolipid signaling pathway: Sphingolipids play a major protective role in the lungs against pulmonary and lung damage, and the control of their pathways may provide effective and improved therapeutic tactics [45]. Regulation of sphingolipids can be very advantageous and can have an effect on anti-inflammation, preservation of neuronal integrity, and anti-coagulant effects [46,47]. All of these favorable features might be utilized to counterbalance the involved disorders of COVID-19 [48].
PI3K-Akt signaling pathway: Activating of PI3K-Akt signaling is an important mechanism to prolong viral replication in both acute and constant infection, thus, virus-infected Vero E6 cells are activated by PI3K-Akt pathway [49,50].
AMPK signaling pathway: In both HBV and HBC viral infections, the AMPK activation is advantageous, but in others, such as the Ebola virus, dengue virus, and even human cytomegaloviral infections, AMPK plays a counteracting role [51].
Wnt signaling pathway: Controlling Wnt signaling could be a critical process in affecting the sustainability of viral pathogenesis [52].
Toll-like receptor signaling pathway: TLR-7 agonists may inhibit the onset of acute COVID-19 in indicative patients and synergy effects with antiviral treatment [53].
NOD-like receptor signaling pathway: NLRs may regulate inflammatory response via NF-κB inhibition and block the viral defense system via interaction with mitochondrial antiviral-signaling protein (MAVS) [54].
RIG-I-like receptor signaling pathway: A report shows that SARS-CoV-2 ORF6 blocks RIG-I-like receptor signaling pathway, thus leading to the interruption of the host’s innate immune system [55].
JAK/STAT signaling pathway: JAK/STAT pathway is linked to the induction of multiple molecular immune pathways. The interruption of this pathway may lead to the blockade of several cellular responses [56].
IL-17 signaling pathway: IL-17 blockers are immunologically probable as a strategy to inhibit severe respiratory symptoms in COVID-19 patients [57].
T cell receptor signaling pathway: It remains uncertain whether T cell responses are beneficial or detrimental in COVID-19, and whether T cell responses have minimal or maximal effects in COVID-19 infection [58].
B cell receptor signaling pathway: B cells have a crucial role in fighting viral infections; however, B cells are triggered after viral infection and before the production of immunoglobulin G (IgG) [59].
TNF signaling pathway: Anti-TNFα monoclonal antibodies are likely to weaken inflammatory disorders occurring in COVID-19 infection, lessening the release of other inflammatory-aggravating triggers [60].
Neutrophin signaling pathway: A report shows that neurofilament light chain (NfL)—a biomarker of astrocytic and neuronal damage—was elevated in the blood of patients with COVID-19 [61]. It is evident that COVID-19 infection might dampen the neurological signaling pathway.
GnRH signaling pathway: When viral infections disturb the Blood–Brain Barrier (BBB), circulating immune cells, such as B and T cells, monocytes, and granulocytes, can enter the brain parenchyma to produce cytokines, and thus, it might be able to negatively alter the functions of GnRH neurons [62].
Prolactin signaling pathway: Patients with HIV have a higher prolactin amount compared to the healthy. Additionally, prolactin is considered as a cytokine in its response in the immune system [63,64].
Adipocytokine signaling pathway: Adipose tissue can react to proinflammatory inducement triggered in the lung through systemic circulation adipocytokines and other inflammatory factors [65].
Relaxin signaling pathway: The dysfunction of the relaxin signaling pathway might bring respiratory disorders recognized in COVID-19 patients. Moreover, SARS-CoV-2 NSP7 protein renders relaxin receptors very much weaker [66].
AGE-RAGE signaling pathway in diabetic complications: AGE-RAGE can also act as an innate immune sensor, including in respiratory viruses. As a consequence, the blocking of RAGE curtails the inflammatory level and progression of cardiovascular diseases [67].
Epithelial cell signaling in Helicobacter pylori infection: Helicobacter pylori downregulates T and B cell signaling to initiate the immune system [68]. It is obvious that COVID-19 patients with Helicobacter pylori might be vulnerable to inflammatory responses.
Based on the pathway enrichment analysis, both RELA and NFKB1 were considered as hub genes in PGCAM against COVID-19. Two genes of each were directly enriched in 21 out of 28 signaling pathways by the PI3K-Akt signaling pathway, indicating that PI3K-Akt signaling pathway might be a hub signaling pathway in PGCAM against COVID-19. The other four genes (IL6, MCL1, VEGFA, and IL2) directly related to the PI3K-Akt signaling pathway might be the significant genes to exert synergistic effects against COVID-19. Liu et al. [69] showed that inhibition of both RELA and/or NFKB1 reduces cytokine secretion and thus alleviates inflammation severity. A report specified that MCL-1 inhibition in mitochondrial membrane stimulates BAX, Cytochrome C, Caspase-9, and Caspase-3, and thus it leads to cell death [70,71]. Noticeably, a report indicated that vascular endothelial growth factor A (VEGFA) was inhibited by angiotensin-converting enzyme 2 (ACE2) and upregulated by the attack of COVID-19, since COVID-19 interrupts the expression of ACE2. Consequently, VEGFA elevates vascular permeability and aggravation of endothelial damage [72]. Most recently, it has been reported that a higher level of both IL-2 and IL-6 was detected in COVID-19 patients [73]. Endothelial cell inflammation is also a major severity in COVID-19 infection, and its unmanageable cytokine production in tissues and cells aggravates severe immune reaction, which is defined as “cytokine storm”, resulting in worsening pneumonia. On top of that, the inhibition of six genes (RELA, NFKB1, IL6, MCL1, VEGFA, and IL2) associated with the PI3K-Akt signaling pathway contributes to the anti-proinflammation, anti-vascular permeability, and pro-apoptosis against COVID-19.
Therefore, the key mechanism of PGCAM against COVID-19 might be to block inflammation and vascular permeability, and trigger pro-apoptosis in tissues and/or cells by inactivating the PI3K-Akt signaling pathway.

4. Materials and Methods

4.1. Selective Compound Construction and Drug-Likeness Evaluation

The LMWCs information of PGCAM was collected by TCMSP (https://tcmspw.com/tcmsp.php), literature, SMILES, and Pubchem CID. The molecular formula of selective compounds was identified using ChemSpider (https://www.chemspider.com/StructureSearch.aspx, accessed on 13 October 2021) or Pubchem (https://pubchem.ncbi.nlm.nih.gov, accessed on 14 October 2021). The drug-likeness properties of the identified compounds were selected through Lipinski’s rule (molecular weight ≤ 500 g/mol; Moriguchi octanol–water partition coefficient ≤ 4.15; the number of nitrogen or oxygen ≤ 10; and the number of NH or OH ≤ 5) in SwissADME (http://www.swissadme.ch, accessed on 14 October 2021). The compound structures were drawn by PubChem Sketcher V2.4 (https://pubchem.ncbi.nlm.nih.gov/edit3/index.html accessed on 14 October 2021).

4.2. Target Genes Related to Selected Compounds or COVID-19

Based on SMILES, target genes associated with screened compounds were selected through both SEA (http://sea.bkslab.org, accessed on 15 October 2021) and STP (http://www.swisstargetprediction.ch, accessed on 15 October 2021) with “Homo Sapiens” mode. The COVID-19-related genes were identified through PubChem (https://pubchem.ncbi.nlm.nih.gov, accessed on 17 October 2021). The overlapping genes between compounds and COVID-19 target genes were identified and visualized by VENNY 2.1 (https://bioinfogp.cnb.csic.es/tools/venny, accessed on 18 October 2021).

4.3. Bubble Chart of Enrichment Pathway Analysis of Overlapping Genes

Genes–genes interaction figure was drawn by STRING (https://string-db.org, accessed on 22 October 2021). RStudio plotted the bubble chart of KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analysis of overlapping genes. The results of KEGG pathway enrichment provide a hint at the promising molecular mechanism of PGCAM against COVID-19.

4.4. Binding Affinity Energy Value of Key Bioactive Compounds on Genes In Silico

The binding affinity energy measurement of the uttermost potential bioactive compounds on key genes was determined by Autodock (http://autodock.scripps.edu, accessed on 25 October 2021), Vina (http://vina.scripps.edu/ accessed on 25 October 2021), Pymol (https://pymol.org/2/ accessed on 25 October 2021).

5. Conclusions

The efficiency of LMWCs extracted from PGCAM and their anti-COVID-19 mechanisms were firstly analyzed by using network pharmacology. The findings of this research indicated that potential bioactive compounds are paenol, frutinone A, ramalic acid, nepetin, and trifolirhizin, and the target genes are RELA, NFKB1, MCL1, VEGFA, IL2, and IL6 of PGCAM against COVID-19. The mechanism(s) of PGCAM against COVID-19 might be inhibited inflammation, vascular permeability, and induction of cell apoptosis against COVID-19 by inactivating the PI3K-Akt signaling pathway. This research provides a scientific indication to support the therapeutic effect of PGCAM on COVID-19, and thus, the proper application of five bioactive compounds against COVID-19 might have potential synergistic effects, such as anti-inflammation, anti-vascular permeability, and pro-apoptosis of the infected cells to prevent COVID-19 vulnerability.
Still, there are also limitations to our analysis. The incompleteness of the herbal natural products database, COVID-19-related genes and protein–protein interaction network might bring bias about target gene findings. Moreover, our approach is dependent on GSEA analysis and molecular docking simulation, either agonist or antagonist. Thus, our predicted result would need to be further improved, through in vitro and in vivo.
The current study did not verify the LMWCs of PGCAM, whose efficacy needs to be validated via clinical test. However, our approach is a holistic perspective and integrated gene–compound interaction, which might be good and promising hints against COVID-19.
Finally, our analysis did not consider the target gene expression level practically after treatment of the selected compounds, which should be implemented in the future.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/pr10020333/s1, Table S1: 1:197 genes (STP + SEA) related to 107 compounds. Table S2: The overlapping genes of SEA and STP: 399 Genes. Table S3: COVID 19-related genes: 356 genes. Table S4: Final 27 overlapping genes related to COVID-19. Figure S1: (A) Interactions between two bioactive compounds and RELA (PDB ID: 1NFI). (B) Interactions between four bioactive compounds and NFKB1 (PDB ID: 1SVC). (C) Interactions between three bioactive compounds and IL6 (PDB ID: 4CNI). (D) Interactions between nine bioactive compounds and MCL1 (PDB ID: 5JSB). (E) Interactions between eleven bioactive compounds and VEGFA (PDB ID: 4KZN). (F) Interactions between six bioactive compounds and IL2 (PDB ID: 2ERJ).

Author Contributions

Conceptualization, methodology, formal analysis, investigation, data curation, writing—original draft, K.-K.O.; software, investigation, data curation, K.-K.O. and M.A.; validation, writing—review and editing, M.A.; supervision, project administration, D.-H.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated or analyzed during this study are included in this published article (and its Supplementary Information files).

Acknowledgments

This study has been worked with the support of a research grant of Kangwon National University in 2022.

Conflicts of Interest

There are no conflicts of interest declared.

Abbreviations

ACE2Angiotensin-Converting Enzyme 2
AGEAdvanced Glycation Endproducts
AMPKAMP-activated Protein Kinase
Caco-2 cellHuman intestinal epithelial cell
cAMPCyclic Adenosine MonoPhosphate
cAMPsCyclic Adenosine MonoPhosphates
CD8Cluster of Differentiation 8
COVID-19SARS-CoV-2 virus
EGFREpidermal Growth Factor Receptor
ErbBErythroblastic Leukemia Viral Oncogene Homolog
FOXOFOrkhead boX protein O
GnRHGonadotropin-Releasing Hormone
HBCHepatitis B Core antigen
HBVHepatitis B Virus
HIF-1Hypoxia-Inducible Factor-1
HIF-1αHypoxia-Inducible Factor-1 alpha
HIVHuman Immunodeficiency Virus
IL2Interleukin 2
IL6Interleukin 6
KEGGKyoto Encyclopedia of Genes and Genomes
LMWCsLow Molecular Weight Compounds (<500 Daltons; accepted by Lipinski’s rule)
LRTILower Respiratory Tract Infection
MAPKMitogen-Activated Protein Kinase
MAVSMitochondrial AntiViral-Signaling protein
MCL1Myeloid Cell Leukemia 1
NFKB1Nuclear Factor Kappa B Subunit 1
NF-KBNuclear Factor Kappa-light-chain-enhancer of activated B cells
NfLNeurofilament Light chain
NODNucleotide-binding Oligomerization Domain
NPNetwork Pharmacology
ORF6Open Reading Frame 6
PGCAMPanax Ginseng C.A. Meyer
PI3K-AktPhosphatidylinositol 3-kinase—protein kinase B
PPARPeroxisome Proliferator-Activated Receptors
PPAR-αPeroxisome Proliferator-Activated Receptor-alpha
RAGEReceptor for Advanced Glycation Endproducts
RELATranscription factor p65
RASRenin-Angiotensin System
RIG-IRetinoic acid-Inducible Gene I
RSVRespiratory Syncytial Virus
SEASimilarity Ensemble Approach
SMILESSimplified Molecular Input Line Entry System
STPSwiss Target Prediction
TCMSPTraditional Chinese Medicine Systems
TLRToll-Like Receptor
TNFTumor Necrosis Factor
VEGFAVascular Endothelial Growth Factor A
WHOWorld Health Organization

References

  1. Belongia, E.A.; Naleway, A.L. Smallpox Vaccine: The Good, the Bad, and the Ugly. Clin. Med. Res. 2003, 1, 87. [Google Scholar] [CrossRef] [Green Version]
  2. Keni, R.; Alexander, A.; Nayak, P.G.; Mudgal, J.; Nandakumar, K. COVID-19: Emergence, Spread, Possible Treatments, and Global Burden. Front. Public Health 2020, 8, 1–13. [Google Scholar] [CrossRef] [PubMed]
  3. Wu, Y.C.; Chen, C.S.; Chan, Y.J. The outbreak of COVID-19: An overview. J. Chin. Med. Assoc. 2020, 83, 217–220. [Google Scholar] [CrossRef] [PubMed]
  4. Zheng, J. SARS-CoV-2: An emerging coronavirus that causes a global threat. Int. J. Biol. Sci. 2020, 16, 1678–1685. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Oh, K.K.; Adnan, M.; Cho, D.H. Network pharmacology approach to decipher signaling pathways associated with target proteins of NSAIDs against COVID-19. Sci. Rep. 2021, 11, 9606. [Google Scholar] [CrossRef] [PubMed]
  6. Oh, K.K.; Adnan, M.; Cho, D.H. Drug-repurposing against COVID-19 by targeting a key signaling pathway: An in silico study. Med. Hypotheses 2021, 155, 110656. [Google Scholar] [CrossRef] [PubMed]
  7. WHO Declares Coronavirus Outbreak a Global Health Emergency—STAT; WHO: Geneva, Switzerland, 2020.
  8. COVID-19 Map—Johns Hopkins Coronavirus Resource Center. Available online: https://coronavirus.jhu.edu/map.html (accessed on 20 December 2021).
  9. Huang, C.; Wang, Y.; Li, X.; Ren, L.; Zhao, J.; Hu, Y.; Zhang, L.; Fan, G.; Xu, J.; Gu, X.; et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet 2020, 395, 497–506. [Google Scholar] [CrossRef] [Green Version]
  10. Speth, M.M.; Singer-Cornelius, T.; Oberle, M.; Gengler, I.; Brockmeier, S.J.; Sedaghat, A.R. Olfactory Dysfunction and Sinonasal Symptomatology in COVID-19: Prevalence, Severity, Timing, and Associated Characteristics. Otolaryngol. Head Neck Surg. 2020, 163, 114–120. [Google Scholar] [CrossRef] [PubMed]
  11. Oran, D.P.; Topol, E.J. Prevalence of Asymptomatic SARS-CoV-2 Infection. Ann. Intern. Med. 2020, 173, 362–367. [Google Scholar] [CrossRef]
  12. Azman, S.; Sekar, M.; Wahidin, S.; Gan, S.H.; Vaijanathappa, J.; Bonam, S.R.; Alvala, M.; Lum, P.T.; Thakur, V.; Beladiya, J.V.; et al. Embelin Alleviates Severe Airway Inflammation in OVA-LPS-Induced Rat Model of Allergic Asthma. J. Asthma Allergy 2021, 14, 1511–1525. [Google Scholar] [CrossRef]
  13. Sanyaolu, A.; Okorie, C.; Marinkovic, A.; Patidar, R.; Younis, K.; Desai, P.; Hosein, Z.; Padda, I.; Mangat, J.; Altaf, M. Comorbidity and its Impact on Patients with COVID-19. SN Compr. Clin. Med. 2020, 2, 1069–1076. [Google Scholar] [CrossRef] [PubMed]
  14. Gordon, D.E.; Jang, G.M.; Bouhaddou, M.; Xu, J.; Obernier, K.; White, K.M.; O’Meara, M.J.; Rezelj, V.V.; Guo, J.Z.; Swaney, D.L.; et al. A SARS-CoV-2 protein interaction map reveals targets for drug repurposing. Nature 2020, 583, 459–468. [Google Scholar] [CrossRef] [PubMed]
  15. Adnan, M.; Chy, M.N.U.; Kamal, A.T.M.M.; Chowdhury, K.A.A.; Rahman, M.A.; Reza, A.S.M.A.; Moniruzzaman, M.; Rony, S.R.; Nasrin, M.S.; Azad, M.O.K.; et al. Intervention in Neuropsychiatric Disorders by Suppressing Inflammatory and Oxidative Stress Signal and Exploration of In Silico Studies for Potential Lead Compounds from Holigarna caustica (Dennst.) Oken leaves. Biomolecules 2020, 10, 561. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Islam, M.T.; Sarkar, C.; El-Kersh, D.M.; Jamaddar, S.; Uddin, S.J.; Shilpi, J.A.; Mubarak, M.S. Natural products and their derivatives against coronavirus: A review of the non-clinical and pre-clinical data. Phyther. Res. 2020, 34, 2471–2492. [Google Scholar] [CrossRef] [PubMed]
  17. Desborough, M.J.R.; Keeling, D.M. The aspirin story—From willow to wonder drug. Br. J. Haematol. 2017, 177, 674–683. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Liu, N.Q.; Schuehly, W.; von Freyhold, M.; Van der Kooy, F. A novel purification method of artemisinin from Artemisia annua. Ind. Crops Prod. 2011, 34, 1084–1088. [Google Scholar] [CrossRef]
  19. Rao, K.V. Taxol and Related Taxanes. I. Taxanes of Taxus brevifolia Bark. Pharm. Res. An Off. J. Am. Assoc. Pharm. Sci. 1993, 10, 521–524. [Google Scholar] [CrossRef]
  20. Chen, Z.; Nakamura, T. Statistical evidence for the usefulness of Chinese medicine in the treatment of SARS. Phyther. Res. 2004, 18, 592–594. [Google Scholar] [CrossRef] [PubMed]
  21. Zhang, D.; Zhang, X.; Peng, B.; Deng, S.; Wang, Y.; Yang, L.; Zhang, K.; Ling, C.; Wu, K. Network pharmacology suggests biochemical rationale for treating COVID-19 symptoms with a Traditional Chinese Medicine. Commun. Biol. 2020, 3, 466. [Google Scholar] [CrossRef] [PubMed]
  22. Wee, J.J.; Park, K.M.; Chung, A.S. Biological activities of ginseng and its application to human health. In Herbal Medicine: Biomolecular and Clinical Aspects, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2011; pp. 157–174. ISBN 9781439807163. [Google Scholar]
  23. Lee, J.S.; Ko, E.J.; Hwang, H.S.; Lee, Y.N.; Kwon, Y.M.; Kim, M.C.; Kang, S.M. Antiviral activity of ginseng extract against respiratory syncytial virus infection. Int. J. Mol. Med. 2014, 34, 183–190. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Nguyen, N.H.; Nguyen, C.T. Pharmacological effects of ginseng on infectious diseases. Inflammopharmacology 2019, 27, 871–883. [Google Scholar] [CrossRef]
  25. Leung, K.W.; Wong, A.S. Pharmacology of ginsenosides: A literature review. Chin. Med. 2010, 5, 20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Han, M.; Fang, X.-L. Difference in oral absorption of ginsenoside Rg 1 between in vitro and in vivo models. Acta Pharmacol. Sin. 2006, 27, 499–505. [Google Scholar] [CrossRef] [PubMed]
  27. Zhang, B.; Zhu, X.M.; Hu, J.N.; Ye, H.; Luo, T.; Liu, X.R.; Li, H.Y.; Li, W.; Zheng, Y.N.; Deng, Z.Y. Absorption mechanism of ginsenoside compound K and its butyl and octyl ester prodrugs in caco-2 cells. J. Agric. Food Chem. 2012, 60, 10278–10284. [Google Scholar] [CrossRef] [PubMed]
  28. Bennion, B.J.; Be, N.A.; Mcnerney, M.W.; Lao, V.; Carlson, E.M.; Valdez, C.A.; Malfatti, M.A.; Enright, H.A.; Nguyen, T.H.; Lightstone, F.C.; et al. Predicting a Drug’s Membrane Permeability: A Computational Model Validated with In Vitro Permeability Assay Data. J. Phys. Chem. B 2017, 121, 5228–5237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Oh, K.; Adnan, M.; Cho, D. Uncovering Mechanisms of Zanthoxylum piperitum Fruits for the Alleviation of Rheumatoid Arthritis Based on Network Pharmacology. Biology 2021, 10, 703. [Google Scholar] [CrossRef] [PubMed]
  30. Hopkins, A.L. Network pharmacology: The next paradigm in drug discovery. Nat. Chem. Biol. 2008, 4, 682–690. [Google Scholar] [CrossRef] [PubMed]
  31. Xu, L.; Zhang, Y.; Zhang, P.; Dai, X.; Gao, Y.; Lv, Y.; Qin, S.; Xu, F. Integrated Metabolomics and Network Pharmacology Strategy-Driven Active Traditional Chinese Medicine Ingredients Discovery for the Alleviation of Cisplatin Nephrotoxicity. Chem. Res. Toxicol. 2019, 32, 2411–2421. [Google Scholar] [CrossRef] [PubMed]
  32. Chowdhury, H.U.; Adnan, M.; Oh, K.K.; Cho, D.H. Decrypting molecular mechanism insight of Phyllanthus emblica L. fruit in the treatment of type 2 diabetes mellitus by network pharmacology. Phytomedicine Plus 2021, 1, 100144. [Google Scholar] [CrossRef]
  33. Tao, L.; Lowe, A.; Wang, G.; Dozmorov, I.; Chang, T.; Yan, N.; Reese, T. Metabolic Control of Viral Infection through PPAR-α Regulation of STING Signaling. bioRxiv 2019, 731208. [Google Scholar] [CrossRef] [Green Version]
  34. Shi, W.; Hou, X.; Li, X.; Peng, H.; Shi, M.; Jiang, Q.; Liu, X.; Ji, Y.; Yao, Y.; He, C.; et al. Differential gene expressions of the MAPK signaling pathway in enterovirus 71-infected rhabdomyosarcoma cells. Brazilian J. Infect. Dis. 2013, 17, 410–417. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Grimes, J.M.; Grimes, K.V. p38 MAPK inhibition: A promising therapeutic approach for COVID-19. J. Mol. Cell. Cardiol. 2020, 144, 63–65. [Google Scholar] [CrossRef] [PubMed]
  36. Ho, J.; Moyes, D.L.; Tavassoli, M.; Naglik, J.R. The Role of ErbB Receptors in Infection. Trends Microbiol. 2017, 25, 942–952. [Google Scholar] [CrossRef] [PubMed]
  37. Menzo, S.; Clementi, M.; Alfani, E.; Bagnarelli, P.; Iacovacci, S.; Manzin, A.; Dandri, M.; Natoli, G.; Levrero, M.; Carloni, G. Trans-Activation of Epidermal Growth Factor Receptor Gene by the Hepatitis B Virus X-Gene Product. Virology 1993, 196, 878–882. [Google Scholar] [CrossRef] [PubMed]
  38. Miyaki, M.; Sato, C.; Sakai, K.; Konishi, M.; Tanaka, K.; Muraoka, M.; Kikuchi-Yanoshita, R.; Nadaoka, Y.; Kanda, H.; Kitagawa, T. Malignant transformation and EGFR activation of immortalized mouse liver epithelial cells caused by HBV enhancer-X from a human hepatocellular carcinoma. Int. J. Cancer 2000, 85, 518–522. [Google Scholar] [CrossRef]
  39. Battagello, D.S.; Dragunas, G.; Klein, M.O.; Ayub, A.L.P.; Velloso, F.J.; Correa, R.G. Unpuzzling COVID-19: Tissue-related signaling pathways associated with SARS-CoV-2 infection and transmission. Clin. Sci. 2020, 134, 2137–2160. [Google Scholar] [CrossRef]
  40. Moreno-Fernandez, M.E.; Rueda, C.M.; Velilla, P.A.; Rugeles, M.T.; Chougnet, C.A. cAMP during HIV infection: Friend or foe? AIDS Res. Hum. Retrovir. 2012, 28, 49–53. [Google Scholar] [CrossRef] [Green Version]
  41. Coperchini, F.; Chiovato, L.; Croce, L.; Magri, F.; Rotondi, M. The cytokine storm in COVID-19: An overview of the involvement of the chemokine/chemokine-receptor system. Cytokine Growth Factor Rev. 2020, 53, 25–32. [Google Scholar] [CrossRef] [PubMed]
  42. DeDiego, M.L.; Nieto-Torres, J.L.; Regla-Nava, J.A.; Jimenez-Guardeno, J.M.; Fernandez-Delgado, R.; Fett, C.; Castano-Rodriguez, C.; Perlman, S.; Enjuanes, L. Inhibition of NF- B-Mediated Inflammation in Severe Acute Respiratory Syndrome Coronavirus-Infected Mice Increases Survival. J. Virol. 2014, 88, 913–924. [Google Scholar] [CrossRef] [Green Version]
  43. Codo, A.C.; Davanzo, G.G.; de Monteiro, L.B.; de Souza, G.F.; Muraro, S.P.; Virgilio-da-Silva, J.V.; Prodonoff, J.S.; Carregari, V.C.; de Biagi Junior, C.A.O.; Crunfli, F.; et al. Elevated Glucose Levels Favor SARS-CoV-2 Infection and Monocyte Response through a HIF-1α/Glycolysis-Dependent Axis. Cell Metab. 2020, 32, 437–446. [Google Scholar] [CrossRef]
  44. Sullivan, J.A.; Kim, E.H.; Plisch, E.H.; Suresh, M. FOXO3 Regulates the CD8 T Cell Response to a Chronic Viral Infection. J. Virol. 2012, 86, 9025–9034. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Abu-Farha, M.; Thanaraj, T.A.; Qaddoumi, M.G.; Hashem, A.; Abubaker, J.; Al-Mulla, F. The role of lipid metabolism in COVID-19 virus infection and as a drug target. Int. J. Mol. Sci. 2020, 21, 3544. [Google Scholar] [CrossRef] [PubMed]
  46. Hannun, Y.A.; Obeid, L.M. Sphingolipids and their metabolism in physiology and disease. Nat. Rev. Mol. Cell Biol. 2018, 19, 175–191. [Google Scholar] [CrossRef] [PubMed]
  47. Huwiler, A.; Zangemeister-Wittke, U. The sphingosine 1-phosphate receptor modulator fingolimod as a therapeutic agent: Recent findings and new perspectives. Pharmacol. Ther. 2018, 185, 34–49. [Google Scholar] [CrossRef] [PubMed]
  48. Lythgoe, M.P.; Middleton, P. Ongoing Clinical Trials for the Management of the COVID-19 Pandemic. Trends Pharmacol. Sci. 2020, 41, 363–382. [Google Scholar] [CrossRef]
  49. Ji, W.-T.; Liu, H. PI3K-Akt Signaling and Viral Infection. Recent Pat. Biotechnol. 2008, 2, 218–226. [Google Scholar] [CrossRef] [PubMed]
  50. Mizutani, T.; Fukushi, S.; Saijo, M.; Kurane, I.; Morikawa, S. Importance of Akt signaling pathway for apoptosis in SARS-CoV-infected Vero E6 cells. Virology 2004, 327, 169–174. [Google Scholar] [CrossRef]
  51. Silwal, P.; Kim, J.K.; Yuk, J.M.; Jo, E.K. AMP-Activated protein kinase and host defense against infection. Int. J. Mol. Sci. 2018, 19, 3495. [Google Scholar] [CrossRef] [Green Version]
  52. van Zuylen, W.J.; Rawlinson, W.D.; Ford, C.E. The Wnt pathway: A key network in cell signalling dysregulated by viruses. Rev. Med. Virol. 2016, 26, 340–355. [Google Scholar] [CrossRef] [PubMed]
  53. Onofrio, L.; Caraglia, M.; Facchini, G.; Margherita, V.; De Placido, S.; Buonerba, C. Toll-like receptors and COVID-19: A two-faced story with an exciting ending. Futur. Sci. OA 2020, 6, FSO605. [Google Scholar] [CrossRef]
  54. Shaw, P.J.; Lamkanfi, M.; Kanneganti, T.D. NOD-like receptor (NLR) signaling beyond the inflammasome. Eur. J. Immunol. 2010, 40, 624–627. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Lei, X.; Dong, X.; Ma, R.; Wang, W.; Xiao, X.; Tian, Z.; Wang, C.; Wang, Y.; Li, L.; Ren, L.; et al. Activation and evasion of type I interferon responses by SARS-CoV-2. Nat. Commun. 2020, 11, 1–12. [Google Scholar] [CrossRef] [PubMed]
  56. Goker Bagca, B.; Biray Avci, C. The potential of JAK/STAT pathway inhibition by ruxolitinib in the treatment of COVID-19. Cytokine Growth Factor Rev. 2020, 54, 51–61. [Google Scholar] [CrossRef] [PubMed]
  57. Pacha, O.; Sallman, M.A.; Evans, S.E. COVID-19: A case for inhibiting IL-17? Nat. Rev. Immunol. 2020, 20, 345–346. [Google Scholar] [CrossRef] [PubMed]
  58. Chen, Z.; John Wherry, E. T cell responses in patients with COVID-19. Nat. Rev. Immunol. 2020, 20, 529–536. [Google Scholar] [CrossRef] [PubMed]
  59. Purtha, W.E.; Chachu, K.A.; Virgin, H.W.; Diamond, M.S. Early B-Cell Activation after West Nile Virus Infection Requires Alpha/Beta Interferon but Not Antigen Receptor Signaling. J. Virol. 2008, 82, 10964–10974. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Catanzaro, M.; Fagiani, F.; Racchi, M.; Corsini, E.; Govoni, S.; Lanni, C. Immune response in COVID-19: Addressing a pharmacological challenge by targeting pathways triggered by SARS-CoV-2. Signal Transduct. Target. Ther. 2020, 5, 1–10. [Google Scholar] [CrossRef] [PubMed]
  61. Wood, H. New insights into the neurological effects of COVID-19. Nat. Rev. Neurol. 2020, 16, 403. [Google Scholar] [CrossRef] [PubMed]
  62. Barabás, K.; Szabó-Meleg, E.; Ábrahám, I.M. Effect of inflammation on female gonadotropin-releasing hormone (GnRH) neurons: Mechanisms and consequences. Int. J. Mol. Sci. 2020, 21, 529. [Google Scholar] [CrossRef] [Green Version]
  63. Acharya, S.; Fernando, J.J.R.; Gama, R. Gynaecomastia, hyperprolactinaemia and HIV infection. Ann. Clin. Biochem. 2005, 42, 301–303. [Google Scholar] [CrossRef] [Green Version]
  64. Borba, V.V.; Zandman-Goddard, G.; Shoenfeld, Y. Prolactin and autoimmunity: The hormone as an inflammatory cytokine. Best Pract. Res. Clin. Endocrinol. Metab. 2019, 33, 101324. [Google Scholar] [CrossRef] [PubMed]
  65. Mancuso, P. Pulmonary Physiology and Pathophysiology in Obesity Obesity and lung inflammation. J. Appl. Physiol. 2010, 108, 722–728. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Khan, M.A.-A.-K.; Islam, A.B.M.M.K. SARS-CoV-2 proteins exploit host’s genetic and epigenetic mediators for the annexation of key host signaling pathways that confers its immune evasion and disease pathophysiology. bioRxiv 2020, 050260. [Google Scholar] [CrossRef]
  67. Hudson, B.I.; Lippman, M.E. Targeting RAGE Signaling in Inflammatory Disease. Annu. Rev. Med. 2018, 69, 349–364. [Google Scholar] [CrossRef] [PubMed]
  68. Reyes, V.E.; Peniche, A.G. Helicobacter pylori deregulates T and B cell signaling to trigger immune evasion. In Current Topics in Microbiology and Immunology; Springer: Berlin/Heidelberg, Germany, 2019; Volume 421, pp. 229–265. [Google Scholar]
  69. Liu, T.; Zhang, L.; Joo, D.; Sun, S.C. NF-κB signaling in inflammation. Signal Transduct. Target. Ther. 2017, 2, 17023. [Google Scholar] [CrossRef] [Green Version]
  70. Tanaka, N. The anti-apoptotic protein MCL1, a novel target of lung cancer therapy. J. Cancer Treat. Diagn. 2017, 2, 54–58. [Google Scholar] [CrossRef]
  71. Chen, L.; Li, X.; Cheng, M.; Wang, S.; Zheng, Q.; Liu, Q. Iso-pencillixanthone A from a marine-derived fungus reverses multidrug resistance in cervical cancer cells through down-regulating P-gp and re-activating apoptosis. RSC Adv. 2018, 8, 41192–41206. [Google Scholar] [CrossRef] [Green Version]
  72. Turkia, M. COVID-19 as an Endothelial Disease and Its Relationship to Vascular Endothelial Growth Factor (VEGF) and Iodide. SSRN Electron. J. 2020, 1–14. [Google Scholar] [CrossRef]
  73. Rahmati, M. Cytokine-targeted therapy in severely ill COVID-19 patients: Options and cautions. Eurasian J. Med. Oncol. 2020, 4, 179–181. [Google Scholar] [CrossRef]
Figure 1. Workflow of network pharmacology analysis of PGCAM against COVID-19.
Figure 1. Workflow of network pharmacology analysis of PGCAM against COVID-19.
Processes 10 00333 g001
Figure 2. Overlapping genes (399 genes) between SEA (672 genes) and STP (924 genes).
Figure 2. Overlapping genes (399 genes) between SEA (672 genes) and STP (924 genes).
Processes 10 00333 g002
Figure 3. Overlapping genes between 399 overlapping genes from two databases (A) and COVID-19-related genes (356 genes) (B). * Both human and mouse ACE2 gene included in Venn diagram, thus indicating a total of 355 genes (B).
Figure 3. Overlapping genes between 399 overlapping genes from two databases (A) and COVID-19-related genes (356 genes) (B). * Both human and mouse ACE2 gene included in Venn diagram, thus indicating a total of 355 genes (B).
Processes 10 00333 g003
Figure 4. Gene–gene interaction with 27 nodes and 117 edges in PGCAM against COVID-19.
Figure 4. Gene–gene interaction with 27 nodes and 117 edges in PGCAM against COVID-19.
Processes 10 00333 g004
Figure 5. Bubble chart of 28 signaling pathways in PGCAM against COVID-19.
Figure 5. Bubble chart of 28 signaling pathways in PGCAM against COVID-19.
Processes 10 00333 g005
Figure 6. Degree values of 17 genes in 28 signaling pathways.
Figure 6. Degree values of 17 genes in 28 signaling pathways.
Processes 10 00333 g006
Figure 7. Chemical structures of potential bioactive compounds in PGCAM against COVID-19.
Figure 7. Chemical structures of potential bioactive compounds in PGCAM against COVID-19.
Processes 10 00333 g007
Figure 8. Molecular docking interaction between best docked compounds from PGCAM and target proteins. (A) Paeonol on 1NFI (B) Frutinone A on 1SVC (C) Nepetin on 4CNI (D) Ramalic acid on 5JSB (E) Trifolirhizin on 4KZN (F) Trifolirhizin on 2ERJ.
Figure 8. Molecular docking interaction between best docked compounds from PGCAM and target proteins. (A) Paeonol on 1NFI (B) Frutinone A on 1SVC (C) Nepetin on 4CNI (D) Ramalic acid on 5JSB (E) Trifolirhizin on 4KZN (F) Trifolirhizin on 2ERJ.
Processes 10 00333 g008aProcesses 10 00333 g008b
Figure 9. Illustration of key findings in this study.
Figure 9. Illustration of key findings in this study.
Processes 10 00333 g009
Table 1. A final list of the selected 107 compounds in PGCAM for network analysis based on Lipinski’s rule.
Table 1. A final list of the selected 107 compounds in PGCAM for network analysis based on Lipinski’s rule.
No.CompoundMolecular Weight (g/mol)Pubchem ID No.CompoundMolecular Weight (g/mol)Pubchem ID
1pisol186.33819355protopanaxadiol460.79920281
23691-11-0204.3594275563,5-dimethyl-p-anisic acid180.288944
3ε-cadinene204.3512302130573-methylheptane114.2311519
4δ-elemene204.3512309449584-methyldodecane184.36521958
5(+)-maalioxide222.376452862595-(E)-heptadec-12-enyl benzene-1,3-diol346.55318012
6(20S)-protopanaxatriol476.711468733607-tetradecyne194.36141979
73,4-dimethylheptane128.251353461nepetin316.265317284
83-ethyl-3-methylheptane142.2814021362ZINC8680004456.77048528
9acetal118.177765633′-biotin244.315315463
10MAV194.1444614564celabenzine379.5442847
11γ- selinene204.3552133465(R)-(+)-citronellal154.2575427
12ginsenoyne A258.35531763266dianthramine289.24441562
13ginsenoyne B294.8531763367ditertbutyl phthalate278.34121712
14ginsenoyne C276.4531763468arachidonate304.5444899
15ginsenoyne D262.4531763569frutinone A264.23441965
16panaxytriol278.49348470folienetriol460.711048822
17methylselenocysteine182.0914700471(20S)-protopanaxadiol460.711213350
18vulgarin264.329425372ginsenoyne E258.355320336
19palmitic acid256.4298573girinimbin263.396943
20linoleic acid280.4528045074L-erythro-isocitric acid192.12439238
21tetradecane198.391238975D-erythro-Isocitric acid192.12447805
22α-gurjunene204.351556027676malvic acid280.410416
23α-cadinol222.371039865677beta-D-mannose180.16439680
24loxanol V214.39820978methyl tricosanoate368.675519
259-hexadecenoic acid254.41466879NN-dimethyldecanamide199.3326690
26methyl myristate242.43128480stearyl acetate312.569968
27pentadecyclic acid242.41384981panaxadiol460.773498
28methyl linoleate294.5528442182panaxatriol476.773599
292,6-dimethyl-3,7-octadiene-2,6-diol170.25535245183pancratistatin325.27441597
30hydron-propanedioate104.0623511544842-formylpyrrole95.113854
31methyl stearate298.5820185ramalic acid346.35320886
32methyl palmitelaidate268.463830386suchilactone368.4132350840
33methyl pentadecanoate256.422351887trifolirhizin446.4442827
34panaxynol244.37528114988pangamic acid476.657346693
35bicyclogermacrene204.351389453789D-pantothenic Acid219.236613
36neocnidilide194.27308385790β-elemene204.356918391
37OCT114.2335691β-santalol220.356857681
38L-adenosine267.2444837492dammarane414.79548714
39elemicin208.251024893methyl margarate284.515609
40diop390.63393494oleanane412.79548717
41β-humulene204.35531810295p-glucosyloxymandelonitrile311.29441465
42stigmasterol412.7528079496darutoside484.644715524
43β-sitosterol414.722228497dibutyl phthalate278.343026
44β-caryophyllene204.35528151598fumarine353.44970
45inermin284.269151099paeonol166.1711092
46methyl (Z)-icos-11-enoate324.55463047100methyl palmitate270.58181
47kaempferol286.245280863101N9-formylharman212.25129650345
48α-humulene epoxide220.3514038843102adenosine267.2460961
49TDA214.3412530103maltol (3-hydroxy-2-methyl-4-pyrone)126.118369
50n-heptadecanol256.515076104salicylic acid138.12338
511-hexadecyne222.4112396105p-hydroxycinnamic acid164.16637542
5213-tetradecenyl acetate254.41521718106gomisin N400.5158103
532,6,10,15-tetramethylheptadecane296.641209107gomisin A416.53001662
542-methyltridecane198.3915269
Table 2. Target genes in 28 signaling pathways enrichment related to COVID-19.
Table 2. Target genes in 28 signaling pathways enrichment related to COVID-19.
KEGG IDDescriptionTarget Genes
hsa04933AGE-RAGE signaling pathway in diabetic complicationsRELA, MAPK8, NFKB1, IL1B, JUN, VEGFA, CXCL8
hsa04657IL-17 signaling pathwayRELA, MAPK8, NFKB1, IL1B, JUN, IL6, CXCL8, MAP3K7
hsa04620Toll-like receptor signaling pathwayRELA, MAPK8, NFKB1, IL1B, JUN, IL6, CXCL8, MAP3K7
hsa04621NOD-like receptor signaling pathwayRELA, MAPK8, NFKB1, IL1B, JUN, IL6, CXCL8, MAP3K7
hsa04668TNF signaling pathwayRELA, MAPK8, NFKB1, IL1B, JUN, IL6, MAP3K7
hsa05120Epithelial cell signaling in Helicobacter pylori infectionRELA, JUN, MAPK8, NFKB1, CXCL8
hsa04622RIG-I-like receptor signaling pathwayRELA, NFKB1, MAPK8, MAP3K7, CXCL8
hsa04010MAPK signaling pathwayRELA, MAPK8, NFKB1, IL1B, JUN, VEGFA, MAP3K7
hsa04064NF-kappa B signaling pathwayRELA, NFKB1, IL1B, CXCL8, MAP3K7
hsa04660T cell receptor signaling pathwayRELA, NFKB1, JUN, IL2,MAP3K7
hsa04071Sphigolipid signaling pathwayRELA, NFKB1,MAPK8,S1PR1,S1PR3
hsa04926Relaxin signaling pathwayRELA, NFKB1,JUN, MAPK8, VEGFA
hsa04920Adipocytokine signaling pathwayRELA, NFKB1, MAPK8, PPARA
hsa04151PI3K-Akt signaling pathwayRELA, NFKB1, IL6, MCL1, VEGFA, IL2
hsa04024cAMP signaling pathwayRELA, NFKB1, JUN, MAPK8, PPARA
hsa04066HIF-1 signaling pathwayRELA, NFKB1, IL6, VEGFA
hsa04722Neurotrophin signaling pathwayRELA, NFKB1,JUN, MAPK8
hsa04917Prolactin signaling pathwayRELA, NFKB1, MAPK8
hsa04662B cell receptor signaling pathwayRELA, NFKB1,JUN
hsa03320PPAR signaling pathwayPPARA, PPARG, FABP2
hsa04014RAS signaling pathwayRELA, NFKB1, MAPK8, VEGFA
hsa04068FoxO signaling pathwayIL6, MAPK8, S1PR1
hsa04310Wnt signaling pathwayJUN, MAPK8, MAP3K7
hsa04630Jak-STAT signaling pathwayIL2, IL6, MCL1
hsa04062Chemokine signaling pathwayRELA, NFKB1, CXCL8
hsa04012ErbB signaling pathwayJUN, MAPK8
hsa04912GnRH signaling pathwayJUN, MAPK8
hsa04152AMPK signaling pathwayPPARG, MAP3K7
Table 3. Binding energy of potential active compounds on RELA (PDB ID:1NFI).
Table 3. Binding energy of potential active compounds on RELA (PDB ID:1NFI).
Hydrogen Bond InteractionsHydrophobic Interactions
ProteinLigandPubChem IDSymbolBinding Energy (kcal/mol)Amino Acid ResidueAmino Acid Residue
1NFINN-dimethyldecanamide26690R1−3.9Thr-67Gly-68, Pro-69
Lys-314, Arg-166
Pro-317, Tyr-66
Asp-167
Paeonol11092R2−5.3Arg-174, Val-163Pro-176, Leu-175
Thr-164, Arg-95Pro-177, Gln-162
Arg-73
Table 4. Binding energy of potential active compounds on NFKB1 (PDB ID: 1SVC).
Table 4. Binding energy of potential active compounds on NFKB1 (PDB ID: 1SVC).
Hydrogen Bond InteractionsHydrophobic Interactions
ProteinLigandPubChem IDSymbolBinding Energy (kcal/mol)Amino Acid ResidueAmino Acid Residue
1SVCFrutinone A441965N1−7.2Ser-113Val-145, Val-61
Leu-143, Arg-157
Ala-156, Thr-153
Lys-149
Paeonol11092N2−5.2Ser-113, Arg-157His-144, Lys-149
Thr-153Ala-111, His-112
Leu-143, Val-145
Val-61
P-hydroxycinnamic acid637542N3−5.6Ala-111, Ser-113His-144, Lys-149
Arg-157Leu-143, Thr-146
Thr-153, His-112
Val-145, Val-61
Dibuthyl phthalate3026N4−4.6Lys-337Lys-278, Phe-301
Ile-281, Pro-303
Gly-299, Phe-298
Gln-282
Table 5. Binding energy of potential active compounds on IL6 (PDB ID: 4CNI).
Table 5. Binding energy of potential active compounds on IL6 (PDB ID: 4CNI).
Hydrogen Bond InteractionsHydrophobic Interactions
ProteinLigandPubChem IDSymbolBinding Energy (kcal/mol)Amino Acid ResidueAmino Acid Residue
4CNILinoleic acid5280450L1−5.0Lys-39Glu-81, Pro-80
Ser-168, Phe-83
Gln-105, Leu-104
Gln-166, Lys-103
Glu-165, Pro-40
β-d-mannose439680L2−5.7Val-170, Asp-167His-171, Lys-169
Asn-138, Asp-170
Ser-168, Gly-169
Nepetin5317284L3−8.0Glu-165, Gln-166Lys-103, Phe-83
Glu-81, Ser-168
Pro-80, Ile-106
Pro-40
Table 6. Binding energy of potential active compounds on MCL1 (PDB ID: 5JSB).
Table 6. Binding energy of potential active compounds on MCL1 (PDB ID: 5JSB).
Hydrogen Bond InteractionsHydrophobic Interactions
ProteinLigandPubChem IDSymbolBinding Energy (kcal/mol)Amino Acid ResidueAmino Acid Residue
5JSBL-adenosine448374M1−6.8Gly-192, Ser-202Arg-201, Glu-115
Tyr-116, Gln-189Arg-222, Gln-221
Glu-188Lys-276, Asp-218
Adenosine60961M2−6.5Ser-202, Gly-192Lys-276, Arg-222
Glu-188, Tyr-116
Asp-218, Gln-189
Salicylic acid338M3−5.3Ala-193, Thr-191Glu-173, Asp-195
Arg-184, Arg-187
Glu-188
Methyl palmitelaidate638303M4−4.2Lys-102Tyr-106, Arg-101
Ala-99, Glu-98
Lys-11, Ala-14
Val-18, Leu-95
Glu-15
Methyl pentadecanoate23518M5−4.4Tyr-106Val-18, Ala-99
Glu-98, Arg-101
Lys-11, Glu-15
Ala-14, Leu-95
Lys-102
Methyl linoleate 5284421M6−4.7Lys-102, Tyr-106Glu-15, Ala-99
Leu-95, Glu-98
Lys-11, Ala-14
Arg-101, Val-18
Ramalic acid5320886M7−6.9Asp-313, Arg-201Gly-311, Glu-284
Arg-310Thr-191, Glu-288
Arg-187
N9-formylharman129650345M8−6.5n/aArg-187, Asp-195
Glu-188, Arg-184
Ala-193, Glu-173
Arg-201, Thr-191
Table 7. Binding energy of potential active compounds on VEGFA (PDB ID: 4KZN).
Table 7. Binding energy of potential active compounds on VEGFA (PDB ID: 4KZN).
Hydrogen Bond InteractionsHydrophobic Interactions
ProteinLigandPubChem IDSymbolBinding Energy (kcal/mol)Amino Acid ResidueAmino Acid Residue
4KZNβ-d-mannose439680V1−5.0Asn-75, Ser-95Phe-96, Leu-97
Tyr-39, Glu-38
Darutoside44715524V2−7.1Asp-41, Ser-95Glu-42, Ile-43
Glu-38, Tyr-39
Phe-96
Inermin91510V3−7.5Glu-38, Leu-97Phe-96, Ser-95
Asp-41Pro-40, Tyr-39
Trifolirhizin442827V4−8.1Asp-41, Glu-42Asn-75, Glu-38
Tyr-39, Ser-95
Leu-97
P-glucosyloxymandelonitrile441465V5−6.8Leu-97, Asn-75Pro-40, Glu-42
Asp-41, Glu-38Ser-95, Phe-96
Tyr-39
Nepetin 5317284V6−7.1n/aAsp-41, Tyr-39
Phe-96, Leu-97
Kaempferol5280863V7−7.1Asn-75Phe-96, Ser-95
Tyr-39, Asp-41
Methyl myristate31284V8−3.9n/aAsn-75, Glu-38
Tyr-39, Ser-95
Leu-97
Methyl stearate8201V9−5.1n/aGlu-42, Asp-41
Glu-38, Asn-75
Leu-97, Ser-95
Tyr-39
Methyl pentadecanoate23518V10−3.9n/aTyr-39, Glu-38
Leu-97, Ser-95
Asn-75
Methyl margarate15609V11−4.4n/aGlu-73, Tyr-39
Asn-75, Leu-97
Table 8. Binding energy of potential active compounds on IL2 (PDB ID: 2ERJ).
Table 8. Binding energy of potential active compounds on IL2 (PDB ID: 2ERJ).
Hydrogen Bond InteractionsHydrophobic Interactions
ProteinLigandPubChem IDSymbolBinding Energy (kcal/mol)Amino Acid ResidueAmino Acid Residue
2ERJDarutoside44715524P1−7.0Ile-63, Gln-46Gln-34, Cys-48
Val-88, Arg-85
Trp-90,Thr-47
Asn-61
Trifolirhizin442827P2−8.5Asp-68, Ser-129Ser-69, Asp-84
Ile-128, Ala-131His-79, Leu-80
Arg-15, Gln-70Ser-132, His-133
P-glucosyloxymandelonitrile441465P3−6.3Asn-61, Arg-85Cys-48, Ile-63
Trp-90, Glu-49Gln-34
Nepetin 5317284P4−7.0Pro-109, Gln-199Ile-110, Leu-112
Met-107, Pro-196
Trp-197, Ala-108
Ser-198
Kaempferol5280863P5−7.4Asp-20, Arg-15Leu-85, Asp-84
His-133Ser-69, Gln-70
Arg-81, His-79
Ala-66, Asp-68
Met-23
Methyl linoleate 5284421P6−3.7n/aCys-60, Cys-48
Leu-51, Glu-49
Ile-63, Asn-61
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Oh, K.-K.; Adnan, M.; Cho, D.-H. The Promising Mechanisms of Low Molecular Weight Compounds of Panax Ginseng C.A. Meyer in Alleviating COVID-19: A Network Pharmacology Analysis. Processes 2022, 10, 333. https://doi.org/10.3390/pr10020333

AMA Style

Oh K-K, Adnan M, Cho D-H. The Promising Mechanisms of Low Molecular Weight Compounds of Panax Ginseng C.A. Meyer in Alleviating COVID-19: A Network Pharmacology Analysis. Processes. 2022; 10(2):333. https://doi.org/10.3390/pr10020333

Chicago/Turabian Style

Oh, Ki-Kwang, Md. Adnan, and Dong-Ha Cho. 2022. "The Promising Mechanisms of Low Molecular Weight Compounds of Panax Ginseng C.A. Meyer in Alleviating COVID-19: A Network Pharmacology Analysis" Processes 10, no. 2: 333. https://doi.org/10.3390/pr10020333

APA Style

Oh, K. -K., Adnan, M., & Cho, D. -H. (2022). The Promising Mechanisms of Low Molecular Weight Compounds of Panax Ginseng C.A. Meyer in Alleviating COVID-19: A Network Pharmacology Analysis. Processes, 10(2), 333. https://doi.org/10.3390/pr10020333

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop