# Using a System Pharmacology Method to Search for the Potential Targets and Pathways of Yinqiaosan against COVID-19

Authors: Li-hua Cao; Xing-yuan Jia; Hong-juan He; Zhen-zhen Wang; Yi-ying Zhao; Xue Yang; Xiu-min Li; Ming-san Miao
Journal: Journal of Healthcare Engineering
 (2022)

---

## Abstract

The first reported case of coronavirus disease 2019 (COVID-19) occurred in Wuhan, Hubei, China. Thereafter, it spread through China and worldwide in only a few months, reaching a pandemic level. It can cause severe respiratory illnesses such as pneumonia and lung failure. Since the onset of the disease, the rapid response and intervention of traditional Chinese medicine (TCM) have played a significant role in the effective control of the epidemic. Yinqiaosan (YQS) was used to treat COVID-19 pneumonia, with good curative effects. However, a systematic overview of its active compounds and the therapeutic mechanisms underlying its action has yet to be performed. The purpose of the current study is to explore the compounds and mechanism of YQS in treating COVID-19 pneumonia using system pharmacology. A system pharmacology method involving drug-likeness assessment, oral bioavailability forecasting, virtual docking, and network analysis was applied to estimate the active compounds, hub targets, and key pathways of YQS in the treatment of COVID-19 pneumonia. With this method, 117 active compounds were successfully identified in YQS, and 77 potential targets were obtained from the targets of 95 compounds and COVID-19 pneumonia. The results show that YQS may act in treating COVID-19 pneumonia and its complications (atherosclerosis and nephropathy) through Kaposi sarcoma-related herpesvirus infection and the AGE-RAGE signaling pathway in diabetic complications and pathways in cancer. We distinguished the hub molecular targets within pathways such as TNF, GAPDH, MAPK3, MAPK1, EGFR, CASP3, MAPK8, mTOR, IL-2, and MAPK14. Five of the more highly active compounds (acacetin, kaempferol, luteolin, naringenin, and quercetin) have anti-inflammatory and antioxidative properties. In summary, by introducing a systematic network pharmacology method, our research perfectly forecasts the active compounds, potential targets, and key pathways of YQS applied to COVID-19 and helps to comprehensively clarify its mechanism of action.


---

## Body

## 1. Introduction

Coronavirus infectious disease 2019 (COVID-19) is caused by 2019-nCoV virus and it seriously endangers human health. On January 12, 2020, the World Health Organization (WHO) confirmed and named its pathogen 2019-nCoV [1], which belongs to the novel coronavirus pneumonia of genus β. On January 31, 2020, it was designated a Public Health Emergency of International Concern (PHEIC). According to the latest statistics of the World Health Organization, as of October 28, 2021, there were approximately 244385444 confirmed cases and more than 4961489 deaths worldwide, including over 7053806 confirmed cases and over 108334 deaths in the United States.The main clinical features of these cases are a high body temperature, cough without mucous, breathing difficulties (dyspnea), headache, and pneumonia. Disease attack may lead to progressive respiratory failure due to alveolar damage (which has been observed by transverse chest computerized tomography images) and cytokine storm syndrome and it could also develop a severe form with sometimes fatal outcomes [2, 3]. The cause of the disease was identified by clinicians as virus-induced pneumonia in line with the clinical signs and other standards, including increased body temperature, decreased lymphocytes and white blood cells (although the levels of the latter seemed normal sometimes), new pulmonary infiltrates on chest radiography, and no evident improvement after treatment with antibiotics for three days [4].In traditional Chinese medicine (TCM), COVID-19 falls under the category of “pestilences,” which occur because of direct or indirect contact with epidemic pathogens. Since the onset of the disease, the rapid response and intervention of TCM have played a vital role in effectively controlling the epidemic, and the classic prescription has become one of the highlights in preventing and controlling the epidemic [5, 6]. Yinqiaosan (YQS) is a traditional prescription for treating the influenza virus from the “Wen Bing Tiao Bian.” YQS is a herbal formula comprised of 10 medicinal herbs (Jinyinhua (JYH, Lonicera japonica Thunb.), Lianqiao (LQ, Forsythia suspensa (Thunb.)), Jiegeng (JG, Platycodon grandiflorum (Jacq.) A. DC), Bohe (BH, Mentha haplocalyx Briq.), Dantouchi (DDC, Glycine max (L.) Merr.), Danzhuye (DZY, Lophatherum gracile Brongn.), Niubangzi (NBZ, Arctium lappa L.), Jingjie (JJ, Schizonepeta tenuifolia Eriq.), Lugen (LG, Phragmites communis Trin.), and Gancao (GC, Glycyrrhiza uralensis Fisch., Glycyrrhiza inflata Bat., or Glycyrrhiza glabra L.). The National Health Commission of the People’s Republic of China and the National Administration of TCM issued the “Diagnosis and Treatment Program for COVID-19 (Trial Version 3).” They began to recommend TCM therapy, among which the recommended prescriptions are Maxing Shigan decoction and YQS.Tang et al. [7], based on the TCMATCOV platform, analyzed the potential effects of commonly used traditional formulations for treating COVID-19 and discovered that YQS is a frequently used therapy for pneumonia. For the past few years, systems pharmacology [8, 9] means have been widely introduced in the prescription and mechanism of TCM. System pharmacology is an emerging systematic methodology and is guided by the “multicomponent multitarget network” theory, which conforms to the integrity and systematicity of TCM prescriptions. In addition, it combines oral bioavailability screening and multiple drug target prediction, providing an overall method for exploring the targets and potential mechanisms of TCM, and it can be used as an alternative strategy to analyze the therapeutic effects of the active ingredients [10, 11]. Wang et al. [12] discussed the targets and signaling pathways of Maxing Shigan decoction for treating pneumonia through systems pharmacology and found that the possible targets were IL-6, TNF, MAPK8, etc. Xia et al. found that Akt1 was a potential target of Lianhua Qingwen to treat and prevent COVID-19 using network pharmacology and molecular docking analyses [13].Aiming to explore the possibility and mechanism of YQS in treating COVID-19, this paper uses systems pharmacology techniques and approaches to forecast the possible targets and signaling pathways of YQS and provides data for subsequent experimental research.


## 2. Materials and Methods

### 2.1. Finding Potential Targets of YQS

Taking the herbs (Jinyinhua, Lianqiao, Jiegeng, Bohe, Dantouchi, Danzhuye, Niubangzi, Jingjie, Lugen, and Gancao) in YQS as keywords, the TCMSP database (http://tcmspw.com/tcmsp.php) was used to retrieve their relevant chemical compositions, and oral bioavailability (OB) ≥ 30%, drug-likeness (DL) ≥ 0.18, and drug half-life (HL) ≥ 4 h were adopted as screening conditions to obtain active ingredients that meet the conditions [14–16]. Among them, the active ingredients of Danzhuye were obtained through a combination of a literature search [17] and the TCMSP database.Through the TCMSP database and the PubChem database (https://pubchem.ncbi.nlm.nih.gov/) combined with the Swiss Target Prediction database, the targets of the active ingredients were obtained. UniProt KB (http://www.UniProt.org/), which prevents, for example, overannotation of semblable proteins such as paralogs and putative products of pseudogenes, was utilized to normalize gene names and organisms. Through the UniProt database (https://www.UniProt.org/), the species was set as “Homo sapiens (Human)” to standardize the obtained drug targets.


### 2.2. Screening of COVID-19-Associated Genes

Through the GeneCards database (http://www.genecards.org/) and OMIM (https://www.omim.org/) database, using “novel coronavirus pneumonia” as the keyword, COVID-19-related disease targets were retrieved and obtained. Only “Homo sapiens” proteins linked to COVID-19 were selected. The UniProt database was applied to standardize the obtained disease targets in Homo sapiens.


### 2.3. Screening of Candidate Targets

The acting targets of the active ingredients of YQS and the disease targets of COVID-19 were processed through Venny2.1 software (http://bioinfogp.cnb.csic.es/tools/venny/index.html) to acquire their common targets, which were taken as alternative targets of YQS for treating COVID-19, and its Venny map was also acquired.


### 2.4. Protein-Protein Interaction (PPI)

The retrieved YQS active ingredient targets were linked to COVID-19 targets by the STRING (https://string-db.org/, ver.11.0) database [18]. The requirement was restrained to “Homo sapiens.” In this paper, a high confidence score with a correlation degree ≥0.400 [19] as the cutoff value was set to obtain the network.


### 2.5. GO and KEGG Pathway Enrichment Analysis

Gene Ontology (GO) enrichment analysis was conducted using the OmicShare tools (http://www.omicshare.com/tools). The KOBAS 3.0 database was applied to analyze the KEGG pathway enrichment of the candidate targets, and the RStudio 3.6.3 and ggplot2 packages were applied to illustrate the results.


### 2.6. Networks Construction and Analysis

Cytoscape 3.5.1, an open-source software platform to visualize complex networks, was employed to illustrate the networks. With this software, the network building was handled as follows: (1) compound-COVID19 PPI network, (2) compound-targets network, and (3) YQS (Yinqiaosan)-compound-targets-pathway network. In addition, four topological features of the hub network, such as “Degree” and “Betweenness Centrality,” were determined to choose YQS alternative targets with topological significance [20, 21]. To determine the relationship between network clusters and to identify high-connectivity hub genes, “Cytohubba” (a plug-in of Cytoscape) was introduced to check the node element. The top 10 nodes were ranked by EPC [22].


### 2.7. Molecular Docking

For the sake of validation of the drug-targeted correlations, the molecular docking simulation further proceeded on every drug docking with their targets. Each drug molecule, such as luteolin, naringenin, acacetin, kaempferol, and quercetin, was downloaded from the TCMSP database and converted to a PDBQT file by Chimera (version 1.10.2).The structures of the proteins in each alternative target were downloaded from the RCBS Protein Data Bank (http://www.rcsb.org/pdb), and every protein file was opened with ADT (version 1.5.6), a free GUI for AutoDock. In every file, the water molecules were canceled, the polar hydrogen atoms were increased, and then they were written to a PDBQT file.The intersection of possible alternative targets and drug targets of PCOS was chosen for in-depth investigations. To validate the interaction and observe the docking sites between drugs and targets, molecular docking was performed by an open-source program called Autodock vina 1.1.2. The docking results were observed with ADT.The foundation for developing medications for the clinical therapy of COVID-19 is shown in Figure1.Figure 1
Process overview. To understand the active ingredients of each single drug in Yinqiaosan (YQS) through the TCMSP (Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform) and TCMID (Traditional Chinese Medicine Integrated Database), the targets of the active ingredients were identified by Swiss TargetPrediction and SEA (similarity ensemble approach). Similarly, COVID-19 targets were obtained through the GeneCards and OMIM (Online Mendelian Inheritance in Man) databases. Intersecting targets were then assessed by molecular docking and GO and KEGG enrichment analysis. GO (Gene Ontology) enrichment analysis was carried out with the OmicShare tools. The KOBAS 3.0 database was used for KEGG enrichment analysis, and RStudio was used to visualize the results. Molecular docking was conducted by an open-source program named Autodock vina 1.1.2. Networks were constructed to provide a visual view by Cytoscape 3.5.1 software.


## 2.1. Finding Potential Targets of YQS

Taking the herbs (Jinyinhua, Lianqiao, Jiegeng, Bohe, Dantouchi, Danzhuye, Niubangzi, Jingjie, Lugen, and Gancao) in YQS as keywords, the TCMSP database (http://tcmspw.com/tcmsp.php) was used to retrieve their relevant chemical compositions, and oral bioavailability (OB) ≥ 30%, drug-likeness (DL) ≥ 0.18, and drug half-life (HL) ≥ 4 h were adopted as screening conditions to obtain active ingredients that meet the conditions [14–16]. Among them, the active ingredients of Danzhuye were obtained through a combination of a literature search [17] and the TCMSP database.Through the TCMSP database and the PubChem database (https://pubchem.ncbi.nlm.nih.gov/) combined with the Swiss Target Prediction database, the targets of the active ingredients were obtained. UniProt KB (http://www.UniProt.org/), which prevents, for example, overannotation of semblable proteins such as paralogs and putative products of pseudogenes, was utilized to normalize gene names and organisms. Through the UniProt database (https://www.UniProt.org/), the species was set as “Homo sapiens (Human)” to standardize the obtained drug targets.


## 2.2. Screening of COVID-19-Associated Genes

Through the GeneCards database (http://www.genecards.org/) and OMIM (https://www.omim.org/) database, using “novel coronavirus pneumonia” as the keyword, COVID-19-related disease targets were retrieved and obtained. Only “Homo sapiens” proteins linked to COVID-19 were selected. The UniProt database was applied to standardize the obtained disease targets in Homo sapiens.


## 2.3. Screening of Candidate Targets

The acting targets of the active ingredients of YQS and the disease targets of COVID-19 were processed through Venny2.1 software (http://bioinfogp.cnb.csic.es/tools/venny/index.html) to acquire their common targets, which were taken as alternative targets of YQS for treating COVID-19, and its Venny map was also acquired.


## 2.4. Protein-Protein Interaction (PPI)

The retrieved YQS active ingredient targets were linked to COVID-19 targets by the STRING (https://string-db.org/, ver.11.0) database [18]. The requirement was restrained to “Homo sapiens.” In this paper, a high confidence score with a correlation degree ≥0.400 [19] as the cutoff value was set to obtain the network.


## 2.5. GO and KEGG Pathway Enrichment Analysis

Gene Ontology (GO) enrichment analysis was conducted using the OmicShare tools (http://www.omicshare.com/tools). The KOBAS 3.0 database was applied to analyze the KEGG pathway enrichment of the candidate targets, and the RStudio 3.6.3 and ggplot2 packages were applied to illustrate the results.


## 2.6. Networks Construction and Analysis

Cytoscape 3.5.1, an open-source software platform to visualize complex networks, was employed to illustrate the networks. With this software, the network building was handled as follows: (1) compound-COVID19 PPI network, (2) compound-targets network, and (3) YQS (Yinqiaosan)-compound-targets-pathway network. In addition, four topological features of the hub network, such as “Degree” and “Betweenness Centrality,” were determined to choose YQS alternative targets with topological significance [20, 21]. To determine the relationship between network clusters and to identify high-connectivity hub genes, “Cytohubba” (a plug-in of Cytoscape) was introduced to check the node element. The top 10 nodes were ranked by EPC [22].


## 2.7. Molecular Docking

For the sake of validation of the drug-targeted correlations, the molecular docking simulation further proceeded on every drug docking with their targets. Each drug molecule, such as luteolin, naringenin, acacetin, kaempferol, and quercetin, was downloaded from the TCMSP database and converted to a PDBQT file by Chimera (version 1.10.2).The structures of the proteins in each alternative target were downloaded from the RCBS Protein Data Bank (http://www.rcsb.org/pdb), and every protein file was opened with ADT (version 1.5.6), a free GUI for AutoDock. In every file, the water molecules were canceled, the polar hydrogen atoms were increased, and then they were written to a PDBQT file.The intersection of possible alternative targets and drug targets of PCOS was chosen for in-depth investigations. To validate the interaction and observe the docking sites between drugs and targets, molecular docking was performed by an open-source program called Autodock vina 1.1.2. The docking results were observed with ADT.The foundation for developing medications for the clinical therapy of COVID-19 is shown in Figure1.Figure 1
Process overview. To understand the active ingredients of each single drug in Yinqiaosan (YQS) through the TCMSP (Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform) and TCMID (Traditional Chinese Medicine Integrated Database), the targets of the active ingredients were identified by Swiss TargetPrediction and SEA (similarity ensemble approach). Similarly, COVID-19 targets were obtained through the GeneCards and OMIM (Online Mendelian Inheritance in Man) databases. Intersecting targets were then assessed by molecular docking and GO and KEGG enrichment analysis. GO (Gene Ontology) enrichment analysis was carried out with the OmicShare tools. The KOBAS 3.0 database was used for KEGG enrichment analysis, and RStudio was used to visualize the results. Molecular docking was conducted by an open-source program named Autodock vina 1.1.2. Networks were constructed to provide a visual view by Cytoscape 3.5.1 software.


## 3. Results

### 3.1. Active Components and Target Proteins of Yinqiaosan

By a systematic search of the public databases, 141 components were retrieved in JYH (17), LQ (11), JG (7), BH (10), DDC (2), DZY (4), NBZ (4), JG (9), LG (1), and GC (76), and the specific information of these components can be found in Table1. After deleting the duplicate values, 117 active ingredients were acquired. Refer to the results in Figure S1. Through observation, it was discovered that some special compounds in this network interacted with multiple targets and took part in the regulation of multiple targets, such as luteolin, naringenin, acacetin, kaempferol, and quercetin.Table 1
Active ingredients of Yinqiaosan (YQS).
Chinese herbsMol idMolecule nameOB (%)DLHL(h)JYHMOL000006Luteolin36.160.2515.94JYHMOL000098Quercetin46.430.2814.4JYHMOL000358Beta-sitosterol36.910.755.36JYHMOL000422Kaempferol41.880.2414.74JYHMOL000449Stigmasterol43.830.765.57JYHMOL001494Mandenol420.195.39JYHMOL001495Ethyl linolenate46.10.26.2JYHMOL002773Beta-carotene37.180.584.36JYHMOL002914Eriodyctiol (flavanone)41.350.2415.88JYHMOL003006(-)-(3R,8S,9R,9aS,10aS)-9-Ethenyl-8-(beta-D-glucopyranosyloxy)-2,3,9,9a,10,10a-hexahydro-5-oxo-5H,8H-pyrano[4,3-d]oxazolo[3,2-a]pyridine-3-carboxylic acid_qt87.470.235.5JYHMOL003036ZINC0397878143.830.765.79JYHMOL003044Chrysoeriol35.850.2716.31JYHMOL003059kryptoxanthin47.250.574.37JYHMOL0030624,5′-Retro-.beta.,.beta.-Carotene-3,3′-dione, 4′,5′-didehydro-31.220.555.39JYHMOL0030955-hydroxy-7-methoxy-2-(3,4,5-trimethoxyphenyl)chromone51.960.4115.98JYHMOL003111Centauroside_qt55.790.55.18JYHMOL003124XYLOSTOSIDINE43.170.649.15LQMOL000006luteolin36.160.2515.94LQMOL000098quercetin46.430.2814.4LQMOL000173wogonin30.680.2317.75LQMOL000211Mairin55.380.788.87LQMOL000358Beta-sitosterol36.910.755.36LQMOL000422Kaempferol41.880.2414.74LQMOL000791Bicuculline69.670.8815.83LQMOL00328120(S)-Dammar-24-ene-3β,20-diol-3-acetate40.230.829.14LQMOL0033153beta-Acetyl-20,25-epoxydammarane-24alpha-ol33.070.797.82LQMOL003365Lactucasterol40.990.855.53LQMOL003370Onjixanthone I79.160.314.86JGMOL000006Luteolin36.160.2515.94JGMOL001689Acacetin34.970.2417.25JGMOL004355Spinasterol42.980.765.32JGMOL004580cis-Dihydroquercetin66.440.2714.51JGMOL0059962-O-Methyl-3―O-β-D-glucopyranosyl platycogenate A45.150.256.03JGMOL006026Dimethyl 2-O-methyl-3-O-a-D-glucopyranosyl platycogenate A39.210.255.04JGMOL006070Robinin39.840.7116.67BHMOL000006Luteolin36.160.2515.94BHMOL000359Sitosterol36.910.755.37BHMOL000471Aloe-emodin83.380.2431.49BHMOL001689Acacetin34.970.2417.25BHMOL001790Linarin39.840.7116.07BHMOL002881Diosmetin31.140.2716.34BHMOL004328Naringenin59.290.2116.98BHMOL005190Eriodictyol71.790.2415.81BHMOL005573Genkwanin37.130.2416.1BHMOL011616Fortunellin35.650.7414.19DDCMOL008400Glycitein50.480.2416.32DDCMOL0116916′-O-Malonylglycitin30.40.8117.25DZYMOL000006Luteolin36.160.2515.94DZYMOL000359Sitosterol36.910.755.37DZYMOL002322Isovitexin31.290.7216.45DZYMOL003137Swertiajaponin32.120.7816.28NBZMOL000358Beta-sitosterol36.910.755.36NBZMOL000422Kaempferol41.880.2414.74NBZMOL002773Beta-carotene37.180.584.36NBZMOL010868Neoarctin A39.990.275.82JJMOL000006Luteolin36.160.2515.94JJMOL000098Quercetin46.430.2814.4JJMOL000358Beta-sitosterol36.910.755.36JJMOL000359Sitosterol36.910.755.37JJMOL000449Stigmasterol43.830.765.57JJMOL002881Diosmetin31.140.2716.34JJMOL005043Campest-5-en-3beta-ol37.580.714.43JJMOL0051005,7-Dihydroxy-2-(3-hydroxy-4-methoxyphenyl)chroman-4-one47.740.2716.51JJMOL011856Schkuhrin I54.450.525.89LGMOL000449Stigmasterol43.830.765.57GCMOL000098Quercetin46.430.2814.4GCMOL000211Mairin55.380.788.87GCMOL000239Jaranol50.830.2915.5GCMOL000354Isorhamnetin49.60.3114.34GCMOL000359Sitosterol36.910.755.37GCMOL000392Formononetin69.670.2117.04GCMOL000417Calycosin47.750.2417.1GCMOL000422Kaempferol41.880.2414.74GCMOL000497Licochalcone a40.790.2916.2GCMOL001484Inermine75.180.5411.72GCMOL001792DFV32.760.1817.89GCMOL002311Glycyrol90.780.679.85GCMOL002565Medicarpin49.220.348.46GCMOL003656Lupiwighteone51.640.3715.63GCMOL0038967-Methoxy-2-methyl isoflavone42.560.216.89GCMOL004328Naringenin59.290.2116.98GCMOL004805(2S)-2-[4-Hydroxy-3-(3-methylbut-2-enyl)phenyl]-8,8-dimethyl-2,3-dihydropyrano[2,3-f]chromen-4-one31.790.7214.82GCMOL004806Euchrenone30.290.5715.89GCMOL004808Glyasperin B65.220.4416.1GCMOL004810Glyasperin F75.840.5415.64GCMOL004814Isotrifoliol31.940.427.91GCMOL004815(E)-1-(2,4-Dihydroxyphenyl)-3-(2,2-dimethylchromen-6-yl)prop-2-en-1-one39.620.3516.16GCMOL004824(2S)-6-(2,4-dihydroxyphenyl)-2-(2-hydroxypropan-2-yl)-4-methoxy-2,3-dihydrofuro[3,2-g]chromen-7-one60.250.634.31GCMOL004827Semilicoisoflavone B48.780.5517.02GCMOL004828Glepidotin A44.720.3516.09GCMOL004829Glepidotin B64.460.3415.98GCMOL004835Glypallichalcone61.60.1917.01GCMOL0048388-(6-Hydroxy-2-benzofuranyl)-2,2-dimethyl-5-chromenol58.440.388.71GCMOL004841Licochalcone B76.760.1917.02GCMOL004848Licochalcone G49.250.3215.75GCMOL004855Licoricone63.580.4716.37GCMOL004856Gancaonin A51.080.416.82GCMOL004857Gancaonin B48.790.4516.49GCMOL004860Licorice glycoside E32.890.2725.39GCMOL0048633-(3,4-Dihydroxyphenyl)-5,7-dihydroxy-8-(3-methylbut-2-enyl)chromone66.370.4115.81GCMOL0048645,7-Dihydroxy-3-(4-methoxyphenyl)-8-(3-methylbut-2-enyl)chromone30.490.4114.99GCMOL0048662-(3,4-Dihydroxyphenyl)-5,7-dihydroxy-6-(3-methylbut-2-enyl)chromone44.150.4116.77GCMOL004882Licocoumarone33.210.369.66GCMOL004883Licoisoflavone41.610.4216.09GCMOL004884Licoisoflavone B38.930.5515.73GCMOL004885Licoisoflavanone52.470.5415.67GCMOL004891Shinpterocarpin80.30.736.5GCMOL004898(E)-3-[3,4-Dihydroxy-5-(3-methylbut-2-enyl)phenyl]-1-(2,4-dihydroxyphenyl)prop-2-en-1-one46.270.3115.24GCMOL004903Liquiritin65.690.7417.96GCMOL004907Glyzaglabrin61.070.3521.2GCMOL004910Glabranin52.90.3116.24GCMOL004912Glabrone52.510.516.09GCMOL0049131,3-Dihydroxy-9-methoxy-6-benzofurano[3,2-c]chromenone48.140.438.87GCMOL0049141,3-Dihydroxy-8,9-dimethoxy-6-benzofurano[3,2-c]chromenone62.90.539.32GCMOL004915Eurycarpin A43.280.3717.1GCMOL004917Glycyroside37.250.7914.62GCMOL004924(-)-Medicocarpin40.990.9513.2GCMOL004935Sigmoidin-B34.880.4114.49GCMOL004941(2R)-7-Hydroxy-2-(4-hydroxyphenyl)chroman-4-one71.120.1818.09GCMOL004945(2S)-7-Hydroxy-2-(4-hydroxyphenyl)-8-(3-methylbut-2-enyl)chroman-4-one36.570.3217.95GCMOL004948Isoglycyrol44.70.846.69GCMOL004949Isolicoflavonol45.170.4215.55GCMOL004957HMO38.370.2116.56GCMOL0049591-Methoxyphaseollidin69.980.649.53GCMOL004961Quercetin der.46.450.3316.61GCMOL004985Icos-5-enoic acid30.70.25.28GCMOL004988Kanzonol F32.470.899.98GCMOL0049896-Prenylated eriodictyol39.220.4116.52GCMOL0049917-Acetoxy-2-methylisoflavone38.920.2617.49GCMOL0049938-Prenylated eriodictyol53.790.415.7GCMOL004996Gadelaidic acid30.70.25.25GCMOL005000Gancaonin G60.440.3916.13GCMOL005001Gancaonin H50.10.7816.64GCMOL005003Licoagrocarpin58.810.589.45GCMOL005007Glyasperins M72.670.5915.57GCMOL005008Glycyrrhiza flavonol A41.280.613.71GCMOL005012Licoagroisoflavone57.280.4919.64GCMOL00501318α-Hydroxyglycyrrhetic acid41.160.714.96GCMOL005016Odoratin49.950.316.35GCMOL005017Phaseol78.770.589.64GCMOL005018Xambioona54.850.8714.5A total of 810 targets were obtained from public databases, but not all information on all 810 targets is shown. The YQS compounds may jointly affect all of these targets, producing pharmacological impacts on COVID-19.


### 3.2. Screening Results of COVID-19-Related Disease Targets

The GeneCards database, as an online catalog of human genes and genetic illness, can enable effective navigation of gene-disease linkages [23]. The OMIM database deals with over 15,500 gene entries and concentrates on elaborating gene-phenotype correlations [24]. Through the two databases, we obtained 435 genes in total (not all shown).


### 3.3. Screening Results of Candidate Targets

Altogether, 77 overlapping genes were acquired by searching the overlaps of the aforementioned compound targets with the 435 COVID-19 gene targets. These genes included PIK3CD, PLA2G5, and ECE1 (TableS1, Figure 2(a)).Figure 2
Overlapping targets (a) and the protein-protein network (b). Node colors change gradually from red to yellow, and node sizes are in proportion to their degree.
(a)(b)By importing compound-disease cotarget information into STRING, a compound-COVID-19 target PPI network was obtained. Taking a score ≥0.4 as the limiting condition [19], the interaction relationship of 77 targets obtained from the String network database was screened. Then, a network diagram of 77 targets protein interactions was obtained. Combined with Cytoscape 3.5.1 software, we obtained an interaction network graph of alternative targets of YQS for treating COVID-19, including 77 nodes and 719 edges. The node colors are illustrated from red to yellow in descending order of degree values (Figure 2(b)). Among the targets of YQS for treating COVID-19, the 5 targets with the highest degree were TNF, GAPDH, MAPK1, MAPK3, and EGFR.


### 3.4. Results of the Construction of the Compound-Target Network

The data of the active ingredients of YQS as well as its targets for treating COVID-19 were input into Cytoscape 3.5.1 software to obtain the network chart of active ingredient-acting targets. The active ingredients of YQS were represented by their Mol IDs. There were 173 nodes and 831 edges forming the compound-target network. The diamonds represent 95 active ingredients, and the ellipse indicates 77 targets of active ingredients (Figure3).Figure 3
The compound-target networks of Yinqiaosan (YQS). The diamond nodes denote ingredients, and the circular nodes denote targets. The node colors from red to blue and node size are proportional to its degree.


### 3.5. Enrichment Analysis Results of Biological Processes and Pathways

To carry out an in-depth recognition of the biological functions of the above 77 potential targets in biological networks, a GO enrichment analysis was implemented with the OmicShare tools, a free online platform for data analysis (http://www.omicshare.com/tools) [25]. As shown in Figure 4, the results fell into three strata: biological processes (BP), cellular components (CC), and molecular functions (MF). The limiting condition was P < 0.01, where the smaller the P value is, the closer the relationship with YQS will be in treating pneumonia. The results show that the BP analysis mainly includes the response to organic substances, the cellular response to organic substances, and the cellular response to chemical stimulus (Figure 4(a)). The CC analysis mainly included membrane microdomains, membrane rafts, and membrane regions (Figure 4(b)). The MF analysis was mainly related to enzyme binding, phosphotransferase activities, alcohol groups as acceptors, and kinase activities (Figure 4(c)).Figure 4
GO enrichment analysis of Yinqiaosan (YQS) targets. The three strata were biological processes (BP) (a), cellular components (CC) (b), and molecular functions (MF) (c).Y-axis: top 20 biological processes associated with the enriched targets. X-axis: enrichment factors (rich factors). Rich factors are proportional to the degree of enrichment. The node size shows the number of genes, and the larger the node is, the more the genes enriched in the term are. The color depth of the node demonstrates the significance level, and the redder the node is, the higher the significance will be to the term.
(a)(b)(c)Utilizing the online platform KOBAS 3.0, 190 important KEGG pathways (P < 0.01) were found on the basis of differentially expressed coding transcripts. The most important pathways included Kaposi sarcoma-related herpesvirus infection, AGE-RAGE signaling pathways in diabetic complications, and pathways in cancer. The result was constructed by KOBAS 3.0, and the top 20 enrichment pathways were visualized by RStudio (Figure 5). This indicated that YQS might be a candidate for treating COVID-19 through these pathways.Figure 5
Analysis of the 20 most enriched KEGG pathways of candidate targets of Yinqiaosan (YQS) for COVID-19. The node size denotes the number of target genes that belong to the pathway. TheX-axis denotes the enrichment factor (rich factor), that is, the number of genes in the pathway in the target gene/the number of genes contained in the pathway in the background gene set. The Y-axis denotes the pathway. The color depth of the node demonstrates the significance level, and the redder the node is, the higher the significance will be to the pathway.


### 3.6. Construction of “YQS-Compounds-Targets-Pathway” Network

The top ten channels in the enrichment analysis of the KEGG pathway and the compound and compound targets were imported into Cytoscape 3.5.1 software to build a “YQS-compounds-targets-pathway” network (Figure6). The figure includes 194 nodes and 1185 edges. The red diamonds denote YQS, yellow hexagons represent the traditional Chinese herbs in YQS, green rectangles represent the potential active ingredients, blue ellipses represent YQS candidate treatment targets for COVID-19, and purple triangles represent the possible signaling pathways for the active ingredients to function. It can be seen from the figure that YQS includes 10 herbs and 95 active ingredients. It mainly acts on TNF, GAPDH, MAPK3, MAPK1, EGFR, CASP3, MAPK8, MTOR, IL-2, and MAPK14 and uses Kaposi sarcoma-related herpesvirus infection, AGE-RAGE signaling pathways in diabetic complications, pathways in cancer, apoptosis, hepatitis B, human cytomegalovirus signaling pathways such as infection, measles, and hepatitis C to treat pneumonia.Figure 6
Yinqiaosan (YQS)-compound-target-pathway network. The red bubbles denote YQS (Yinqiaosan), the yellow bubbles denote Chinese herbs, the pink bubbles denote compounds, and the violet bubbles denote pathways.


### 3.7. Molecular Docking

Molecular docking was further utilized in our research to estimate the interaction between components and targets to reduce the complexity and increase the precision of the constituent target network [26]. The cytoHubba control panel is a useful tool to retrieve subnetworks from the whole large PPI set [27]. Hub gene recognition was processed by the cytoHubba plug-in in Cytoscape software. All active ingredients in the top 5 online pharmacology degrees of freedom were taken as the research object. To determine the valid connecting effects between YQS-main active ingredients and their predicted targets, molecular docking was assessed using the YQS-hub genes connecting the energy. Molecular docking analyzed 10 hub genes (CASP3, EGFR, GAPDH, IL-2, MAPK1, MAPK14, MAPK3, MAPK8, mTOR, and TNF) and 5 compounds (acacetin, kaempferol, luteolin, naringenin, and quercetin).Zhou et al. [3] ascertained that 2019-nCoV makes use of the same cell entry receptor, angiotensin transforming enzyme II (ACE2), as SARS-CoV. 3CL hydrolase is the core protease of single-stranded RNA virus precursor polyproteolysis and it plays a vital role when replicating the single-stranded RNA virus [28]. Therefore, in this study, the 5 compounds were also molecularly docked with SARS-CoV-2 3CL hydrolase (Mpro) and ACE2 to provide references for treating COVID-19 by TCM.The lower the binding energy, the more stable the conformation of the ligand-receptor binding [29]. When the binding energy is less than −5.0 kcal/mol, this suggests that the compound can bind to the target. When the binding energy is less than -7.0 kcal/mol, this suggests that the compound has a good ability to bind to the target. We selected the conformation with the lowest binding energy to analyze the docking binding mode and used Discovery Studio for mapping.Molecular docking results showed that all five active ingredients could be combined with 10 hub genes (TableS2), Mpro or ACE2 (Table 2, Figure 7). Therefore, the 5 compounds could efficiently act on the 10 targets, especially TNF, MAPK14, MAPK3, and MAPK8. Kaempferol (Figures 7(a) and 7(b)) and Mpro have a better binding capacity, which is consistent with the docking results of the molecules discovered by Ling [30]. Acacetin (Figures 7(c) and 7(d)), luteolin (Figures 7(e) and 7(f)), and ACE2 have a better combining ability.Table 2
Docking results of main ingredients with target proteins.
No.Ligand nameMol idBinding free energy/(kcal⋅mol−1)MproACE21AcacetinMOL001689−7.5−8.82KaempferolMOL000422−7.9−8.63LuteolinMOL000006−7.4−8.84NaringeninMOL004328−6.8−7.95QuercetinMOL000098−7.5−8.5Figure 7
Molecular docking of compounds from YQS for COVID-19 targets. Predicted binding mode of Mpro (PDB id: 6lu7) with kaempferol in the active pocket (a). 2D binding view of kaempferol with Mpro (b). Predicted binding mode of ACE2 (PDB id: 1r4l) with acacetin in the active pocket (c). 2D binding view of acacetin with ACE2 (d). Predicted binding mode of ACE2 (PDB id: 1r4l) with luteolin in the active pocket (e). 2D binding view of luteolin with ACE2 (f).


## 3.1. Active Components and Target Proteins of Yinqiaosan

By a systematic search of the public databases, 141 components were retrieved in JYH (17), LQ (11), JG (7), BH (10), DDC (2), DZY (4), NBZ (4), JG (9), LG (1), and GC (76), and the specific information of these components can be found in Table1. After deleting the duplicate values, 117 active ingredients were acquired. Refer to the results in Figure S1. Through observation, it was discovered that some special compounds in this network interacted with multiple targets and took part in the regulation of multiple targets, such as luteolin, naringenin, acacetin, kaempferol, and quercetin.Table 1
Active ingredients of Yinqiaosan (YQS).
Chinese herbsMol idMolecule nameOB (%)DLHL(h)JYHMOL000006Luteolin36.160.2515.94JYHMOL000098Quercetin46.430.2814.4JYHMOL000358Beta-sitosterol36.910.755.36JYHMOL000422Kaempferol41.880.2414.74JYHMOL000449Stigmasterol43.830.765.57JYHMOL001494Mandenol420.195.39JYHMOL001495Ethyl linolenate46.10.26.2JYHMOL002773Beta-carotene37.180.584.36JYHMOL002914Eriodyctiol (flavanone)41.350.2415.88JYHMOL003006(-)-(3R,8S,9R,9aS,10aS)-9-Ethenyl-8-(beta-D-glucopyranosyloxy)-2,3,9,9a,10,10a-hexahydro-5-oxo-5H,8H-pyrano[4,3-d]oxazolo[3,2-a]pyridine-3-carboxylic acid_qt87.470.235.5JYHMOL003036ZINC0397878143.830.765.79JYHMOL003044Chrysoeriol35.850.2716.31JYHMOL003059kryptoxanthin47.250.574.37JYHMOL0030624,5′-Retro-.beta.,.beta.-Carotene-3,3′-dione, 4′,5′-didehydro-31.220.555.39JYHMOL0030955-hydroxy-7-methoxy-2-(3,4,5-trimethoxyphenyl)chromone51.960.4115.98JYHMOL003111Centauroside_qt55.790.55.18JYHMOL003124XYLOSTOSIDINE43.170.649.15LQMOL000006luteolin36.160.2515.94LQMOL000098quercetin46.430.2814.4LQMOL000173wogonin30.680.2317.75LQMOL000211Mairin55.380.788.87LQMOL000358Beta-sitosterol36.910.755.36LQMOL000422Kaempferol41.880.2414.74LQMOL000791Bicuculline69.670.8815.83LQMOL00328120(S)-Dammar-24-ene-3β,20-diol-3-acetate40.230.829.14LQMOL0033153beta-Acetyl-20,25-epoxydammarane-24alpha-ol33.070.797.82LQMOL003365Lactucasterol40.990.855.53LQMOL003370Onjixanthone I79.160.314.86JGMOL000006Luteolin36.160.2515.94JGMOL001689Acacetin34.970.2417.25JGMOL004355Spinasterol42.980.765.32JGMOL004580cis-Dihydroquercetin66.440.2714.51JGMOL0059962-O-Methyl-3―O-β-D-glucopyranosyl platycogenate A45.150.256.03JGMOL006026Dimethyl 2-O-methyl-3-O-a-D-glucopyranosyl platycogenate A39.210.255.04JGMOL006070Robinin39.840.7116.67BHMOL000006Luteolin36.160.2515.94BHMOL000359Sitosterol36.910.755.37BHMOL000471Aloe-emodin83.380.2431.49BHMOL001689Acacetin34.970.2417.25BHMOL001790Linarin39.840.7116.07BHMOL002881Diosmetin31.140.2716.34BHMOL004328Naringenin59.290.2116.98BHMOL005190Eriodictyol71.790.2415.81BHMOL005573Genkwanin37.130.2416.1BHMOL011616Fortunellin35.650.7414.19DDCMOL008400Glycitein50.480.2416.32DDCMOL0116916′-O-Malonylglycitin30.40.8117.25DZYMOL000006Luteolin36.160.2515.94DZYMOL000359Sitosterol36.910.755.37DZYMOL002322Isovitexin31.290.7216.45DZYMOL003137Swertiajaponin32.120.7816.28NBZMOL000358Beta-sitosterol36.910.755.36NBZMOL000422Kaempferol41.880.2414.74NBZMOL002773Beta-carotene37.180.584.36NBZMOL010868Neoarctin A39.990.275.82JJMOL000006Luteolin36.160.2515.94JJMOL000098Quercetin46.430.2814.4JJMOL000358Beta-sitosterol36.910.755.36JJMOL000359Sitosterol36.910.755.37JJMOL000449Stigmasterol43.830.765.57JJMOL002881Diosmetin31.140.2716.34JJMOL005043Campest-5-en-3beta-ol37.580.714.43JJMOL0051005,7-Dihydroxy-2-(3-hydroxy-4-methoxyphenyl)chroman-4-one47.740.2716.51JJMOL011856Schkuhrin I54.450.525.89LGMOL000449Stigmasterol43.830.765.57GCMOL000098Quercetin46.430.2814.4GCMOL000211Mairin55.380.788.87GCMOL000239Jaranol50.830.2915.5GCMOL000354Isorhamnetin49.60.3114.34GCMOL000359Sitosterol36.910.755.37GCMOL000392Formononetin69.670.2117.04GCMOL000417Calycosin47.750.2417.1GCMOL000422Kaempferol41.880.2414.74GCMOL000497Licochalcone a40.790.2916.2GCMOL001484Inermine75.180.5411.72GCMOL001792DFV32.760.1817.89GCMOL002311Glycyrol90.780.679.85GCMOL002565Medicarpin49.220.348.46GCMOL003656Lupiwighteone51.640.3715.63GCMOL0038967-Methoxy-2-methyl isoflavone42.560.216.89GCMOL004328Naringenin59.290.2116.98GCMOL004805(2S)-2-[4-Hydroxy-3-(3-methylbut-2-enyl)phenyl]-8,8-dimethyl-2,3-dihydropyrano[2,3-f]chromen-4-one31.790.7214.82GCMOL004806Euchrenone30.290.5715.89GCMOL004808Glyasperin B65.220.4416.1GCMOL004810Glyasperin F75.840.5415.64GCMOL004814Isotrifoliol31.940.427.91GCMOL004815(E)-1-(2,4-Dihydroxyphenyl)-3-(2,2-dimethylchromen-6-yl)prop-2-en-1-one39.620.3516.16GCMOL004824(2S)-6-(2,4-dihydroxyphenyl)-2-(2-hydroxypropan-2-yl)-4-methoxy-2,3-dihydrofuro[3,2-g]chromen-7-one60.250.634.31GCMOL004827Semilicoisoflavone B48.780.5517.02GCMOL004828Glepidotin A44.720.3516.09GCMOL004829Glepidotin B64.460.3415.98GCMOL004835Glypallichalcone61.60.1917.01GCMOL0048388-(6-Hydroxy-2-benzofuranyl)-2,2-dimethyl-5-chromenol58.440.388.71GCMOL004841Licochalcone B76.760.1917.02GCMOL004848Licochalcone G49.250.3215.75GCMOL004855Licoricone63.580.4716.37GCMOL004856Gancaonin A51.080.416.82GCMOL004857Gancaonin B48.790.4516.49GCMOL004860Licorice glycoside E32.890.2725.39GCMOL0048633-(3,4-Dihydroxyphenyl)-5,7-dihydroxy-8-(3-methylbut-2-enyl)chromone66.370.4115.81GCMOL0048645,7-Dihydroxy-3-(4-methoxyphenyl)-8-(3-methylbut-2-enyl)chromone30.490.4114.99GCMOL0048662-(3,4-Dihydroxyphenyl)-5,7-dihydroxy-6-(3-methylbut-2-enyl)chromone44.150.4116.77GCMOL004882Licocoumarone33.210.369.66GCMOL004883Licoisoflavone41.610.4216.09GCMOL004884Licoisoflavone B38.930.5515.73GCMOL004885Licoisoflavanone52.470.5415.67GCMOL004891Shinpterocarpin80.30.736.5GCMOL004898(E)-3-[3,4-Dihydroxy-5-(3-methylbut-2-enyl)phenyl]-1-(2,4-dihydroxyphenyl)prop-2-en-1-one46.270.3115.24GCMOL004903Liquiritin65.690.7417.96GCMOL004907Glyzaglabrin61.070.3521.2GCMOL004910Glabranin52.90.3116.24GCMOL004912Glabrone52.510.516.09GCMOL0049131,3-Dihydroxy-9-methoxy-6-benzofurano[3,2-c]chromenone48.140.438.87GCMOL0049141,3-Dihydroxy-8,9-dimethoxy-6-benzofurano[3,2-c]chromenone62.90.539.32GCMOL004915Eurycarpin A43.280.3717.1GCMOL004917Glycyroside37.250.7914.62GCMOL004924(-)-Medicocarpin40.990.9513.2GCMOL004935Sigmoidin-B34.880.4114.49GCMOL004941(2R)-7-Hydroxy-2-(4-hydroxyphenyl)chroman-4-one71.120.1818.09GCMOL004945(2S)-7-Hydroxy-2-(4-hydroxyphenyl)-8-(3-methylbut-2-enyl)chroman-4-one36.570.3217.95GCMOL004948Isoglycyrol44.70.846.69GCMOL004949Isolicoflavonol45.170.4215.55GCMOL004957HMO38.370.2116.56GCMOL0049591-Methoxyphaseollidin69.980.649.53GCMOL004961Quercetin der.46.450.3316.61GCMOL004985Icos-5-enoic acid30.70.25.28GCMOL004988Kanzonol F32.470.899.98GCMOL0049896-Prenylated eriodictyol39.220.4116.52GCMOL0049917-Acetoxy-2-methylisoflavone38.920.2617.49GCMOL0049938-Prenylated eriodictyol53.790.415.7GCMOL004996Gadelaidic acid30.70.25.25GCMOL005000Gancaonin G60.440.3916.13GCMOL005001Gancaonin H50.10.7816.64GCMOL005003Licoagrocarpin58.810.589.45GCMOL005007Glyasperins M72.670.5915.57GCMOL005008Glycyrrhiza flavonol A41.280.613.71GCMOL005012Licoagroisoflavone57.280.4919.64GCMOL00501318α-Hydroxyglycyrrhetic acid41.160.714.96GCMOL005016Odoratin49.950.316.35GCMOL005017Phaseol78.770.589.64GCMOL005018Xambioona54.850.8714.5A total of 810 targets were obtained from public databases, but not all information on all 810 targets is shown. The YQS compounds may jointly affect all of these targets, producing pharmacological impacts on COVID-19.


## 3.2. Screening Results of COVID-19-Related Disease Targets

The GeneCards database, as an online catalog of human genes and genetic illness, can enable effective navigation of gene-disease linkages [23]. The OMIM database deals with over 15,500 gene entries and concentrates on elaborating gene-phenotype correlations [24]. Through the two databases, we obtained 435 genes in total (not all shown).


## 3.3. Screening Results of Candidate Targets

Altogether, 77 overlapping genes were acquired by searching the overlaps of the aforementioned compound targets with the 435 COVID-19 gene targets. These genes included PIK3CD, PLA2G5, and ECE1 (TableS1, Figure 2(a)).Figure 2
Overlapping targets (a) and the protein-protein network (b). Node colors change gradually from red to yellow, and node sizes are in proportion to their degree.
(a)(b)By importing compound-disease cotarget information into STRING, a compound-COVID-19 target PPI network was obtained. Taking a score ≥0.4 as the limiting condition [19], the interaction relationship of 77 targets obtained from the String network database was screened. Then, a network diagram of 77 targets protein interactions was obtained. Combined with Cytoscape 3.5.1 software, we obtained an interaction network graph of alternative targets of YQS for treating COVID-19, including 77 nodes and 719 edges. The node colors are illustrated from red to yellow in descending order of degree values (Figure 2(b)). Among the targets of YQS for treating COVID-19, the 5 targets with the highest degree were TNF, GAPDH, MAPK1, MAPK3, and EGFR.


## 3.4. Results of the Construction of the Compound-Target Network

The data of the active ingredients of YQS as well as its targets for treating COVID-19 were input into Cytoscape 3.5.1 software to obtain the network chart of active ingredient-acting targets. The active ingredients of YQS were represented by their Mol IDs. There were 173 nodes and 831 edges forming the compound-target network. The diamonds represent 95 active ingredients, and the ellipse indicates 77 targets of active ingredients (Figure3).Figure 3
The compound-target networks of Yinqiaosan (YQS). The diamond nodes denote ingredients, and the circular nodes denote targets. The node colors from red to blue and node size are proportional to its degree.


## 3.5. Enrichment Analysis Results of Biological Processes and Pathways

To carry out an in-depth recognition of the biological functions of the above 77 potential targets in biological networks, a GO enrichment analysis was implemented with the OmicShare tools, a free online platform for data analysis (http://www.omicshare.com/tools) [25]. As shown in Figure 4, the results fell into three strata: biological processes (BP), cellular components (CC), and molecular functions (MF). The limiting condition was P < 0.01, where the smaller the P value is, the closer the relationship with YQS will be in treating pneumonia. The results show that the BP analysis mainly includes the response to organic substances, the cellular response to organic substances, and the cellular response to chemical stimulus (Figure 4(a)). The CC analysis mainly included membrane microdomains, membrane rafts, and membrane regions (Figure 4(b)). The MF analysis was mainly related to enzyme binding, phosphotransferase activities, alcohol groups as acceptors, and kinase activities (Figure 4(c)).Figure 4
GO enrichment analysis of Yinqiaosan (YQS) targets. The three strata were biological processes (BP) (a), cellular components (CC) (b), and molecular functions (MF) (c).Y-axis: top 20 biological processes associated with the enriched targets. X-axis: enrichment factors (rich factors). Rich factors are proportional to the degree of enrichment. The node size shows the number of genes, and the larger the node is, the more the genes enriched in the term are. The color depth of the node demonstrates the significance level, and the redder the node is, the higher the significance will be to the term.
(a)(b)(c)Utilizing the online platform KOBAS 3.0, 190 important KEGG pathways (P < 0.01) were found on the basis of differentially expressed coding transcripts. The most important pathways included Kaposi sarcoma-related herpesvirus infection, AGE-RAGE signaling pathways in diabetic complications, and pathways in cancer. The result was constructed by KOBAS 3.0, and the top 20 enrichment pathways were visualized by RStudio (Figure 5). This indicated that YQS might be a candidate for treating COVID-19 through these pathways.Figure 5
Analysis of the 20 most enriched KEGG pathways of candidate targets of Yinqiaosan (YQS) for COVID-19. The node size denotes the number of target genes that belong to the pathway. TheX-axis denotes the enrichment factor (rich factor), that is, the number of genes in the pathway in the target gene/the number of genes contained in the pathway in the background gene set. The Y-axis denotes the pathway. The color depth of the node demonstrates the significance level, and the redder the node is, the higher the significance will be to the pathway.


## 3.6. Construction of “YQS-Compounds-Targets-Pathway” Network

The top ten channels in the enrichment analysis of the KEGG pathway and the compound and compound targets were imported into Cytoscape 3.5.1 software to build a “YQS-compounds-targets-pathway” network (Figure6). The figure includes 194 nodes and 1185 edges. The red diamonds denote YQS, yellow hexagons represent the traditional Chinese herbs in YQS, green rectangles represent the potential active ingredients, blue ellipses represent YQS candidate treatment targets for COVID-19, and purple triangles represent the possible signaling pathways for the active ingredients to function. It can be seen from the figure that YQS includes 10 herbs and 95 active ingredients. It mainly acts on TNF, GAPDH, MAPK3, MAPK1, EGFR, CASP3, MAPK8, MTOR, IL-2, and MAPK14 and uses Kaposi sarcoma-related herpesvirus infection, AGE-RAGE signaling pathways in diabetic complications, pathways in cancer, apoptosis, hepatitis B, human cytomegalovirus signaling pathways such as infection, measles, and hepatitis C to treat pneumonia.Figure 6
Yinqiaosan (YQS)-compound-target-pathway network. The red bubbles denote YQS (Yinqiaosan), the yellow bubbles denote Chinese herbs, the pink bubbles denote compounds, and the violet bubbles denote pathways.


## 3.7. Molecular Docking

Molecular docking was further utilized in our research to estimate the interaction between components and targets to reduce the complexity and increase the precision of the constituent target network [26]. The cytoHubba control panel is a useful tool to retrieve subnetworks from the whole large PPI set [27]. Hub gene recognition was processed by the cytoHubba plug-in in Cytoscape software. All active ingredients in the top 5 online pharmacology degrees of freedom were taken as the research object. To determine the valid connecting effects between YQS-main active ingredients and their predicted targets, molecular docking was assessed using the YQS-hub genes connecting the energy. Molecular docking analyzed 10 hub genes (CASP3, EGFR, GAPDH, IL-2, MAPK1, MAPK14, MAPK3, MAPK8, mTOR, and TNF) and 5 compounds (acacetin, kaempferol, luteolin, naringenin, and quercetin).Zhou et al. [3] ascertained that 2019-nCoV makes use of the same cell entry receptor, angiotensin transforming enzyme II (ACE2), as SARS-CoV. 3CL hydrolase is the core protease of single-stranded RNA virus precursor polyproteolysis and it plays a vital role when replicating the single-stranded RNA virus [28]. Therefore, in this study, the 5 compounds were also molecularly docked with SARS-CoV-2 3CL hydrolase (Mpro) and ACE2 to provide references for treating COVID-19 by TCM.The lower the binding energy, the more stable the conformation of the ligand-receptor binding [29]. When the binding energy is less than −5.0 kcal/mol, this suggests that the compound can bind to the target. When the binding energy is less than -7.0 kcal/mol, this suggests that the compound has a good ability to bind to the target. We selected the conformation with the lowest binding energy to analyze the docking binding mode and used Discovery Studio for mapping.Molecular docking results showed that all five active ingredients could be combined with 10 hub genes (TableS2), Mpro or ACE2 (Table 2, Figure 7). Therefore, the 5 compounds could efficiently act on the 10 targets, especially TNF, MAPK14, MAPK3, and MAPK8. Kaempferol (Figures 7(a) and 7(b)) and Mpro have a better binding capacity, which is consistent with the docking results of the molecules discovered by Ling [30]. Acacetin (Figures 7(c) and 7(d)), luteolin (Figures 7(e) and 7(f)), and ACE2 have a better combining ability.Table 2
Docking results of main ingredients with target proteins.
No.Ligand nameMol idBinding free energy/(kcal⋅mol−1)MproACE21AcacetinMOL001689−7.5−8.82KaempferolMOL000422−7.9−8.63LuteolinMOL000006−7.4−8.84NaringeninMOL004328−6.8−7.95QuercetinMOL000098−7.5−8.5Figure 7
Molecular docking of compounds from YQS for COVID-19 targets. Predicted binding mode of Mpro (PDB id: 6lu7) with kaempferol in the active pocket (a). 2D binding view of kaempferol with Mpro (b). Predicted binding mode of ACE2 (PDB id: 1r4l) with acacetin in the active pocket (c). 2D binding view of acacetin with ACE2 (d). Predicted binding mode of ACE2 (PDB id: 1r4l) with luteolin in the active pocket (e). 2D binding view of luteolin with ACE2 (f).


## 4. Discussion

At present, COVID-19 has expanded to over 200 countries and areas, and the lives and health of people around the world are under serious threat. Special drugs against COVID-19 are still under development. TCM adopts the principle of dialectical treatment with the characteristics of multiple targets and multiple pathways, and it plays a vital role in antiepidemic treatment [31–33]. The “COVID-19 Diagnosis and Treatment Program (Trial Ver. 3)” issued by the National Health Commission of the People’s Republic of China pointed out that COVID-19 is categorized as a “pestile” in TCM, which happens under direct or indirect contact with epidemic pathogens. The main features are “damp, heat, poisonous, stasis.” COVID-19 is divided into four types of TCM syndromes: cold-dampness stagnating in the lung, damp-heat accumulating in the lung, damp toxin stagnating in the lung, and internal blocking causing external collapse. The clinical symptoms of “damp-heat accumulating in the lung” are low body temperature or no fever, mild aversion to cold, fatigue, a heavy sensation in the head and body, dry coughs with scanty phlegm, sore throat, and a dry mouth with no wish to drink water. Alternatively, chest tightness, epigastric fullness, absence of sweating or inhibited sweating, vomiting, a poor appetite, and loose stools or sticky stools with a feeling of incomplete bowel movement may also be present. The tongue is pale red with a white, thick, and greasy or thin, yellow coating. The pulse is soft or slippery and rapid. The recommended prescriptions are Maxing shigantang and YQS.YQS was described by Wu in “Wen Bing Tiao Bian.” He is a famous infectious febrile disease specialist. Part of the composition of the Lianhua Qingwen capsule comes from YQS, which is particularly effective in treating pneumonia. YQS is comprised of ten Chinese herbs, JYH, LQ, JG, BH, DDC, DZY, NBZ, JJ, LG, and GC, and mainly treats the manifestations caused by diseases such as damp-heat accumulation in the lung. Bohe and Niubangzi have flavor and cool properties and can help sore throats. Jingjie and Dandouchi can help the main medicine drive away damp heat. Dan Zhuye and Lugen can clear the heat, and Jiegeng and Gancao can relieve cough. YQS is commonly used in the clinic against influenza A virus infection [34–36].This study uses network pharmacology to forecast the main active components and potential molecular mechanism of YQS against COVID-19. Through the TCMSP database, 141 active ingredients of each single medicine in Yinqiao powder were obtained, of which the five most powerful components were acacetin, kaempferol, luteolin, naringenin, and quercetin. Studies have shown that luteolin may suppress the infection of respiratory epithelial cells byEscherichia coli through immunomodulation [37]. Wu et al. established a mouse model of Staphylococcus aureus pneumonia and found that naringenin could significantly alleviate the disease [38]. Acacia has antibacterial, antioxidant, anti-inflammatory, antiviral, antipyretic, and analgesic effects [39]. It can protect against atherosclerosis by antagonizing the LysoPC-induced apoptosis of vascular smooth muscle cells and stimulating the antioxidant regulatory factor Nrf2 [40]. Kaempferol has an extensive scope of pharmacological activities, including antioxidant, anti-inflammatory, antimicrobial, etc. [41], and its derivatives are used as antiviral drugs against coronavirus 3a channel protein [42]. Quercetin and its derivatives also have strong anti-inflammatory and analgesic pharmacological activities [43]. Wang et al. [12] found that active compounds such as quercetin, kaempferol, naringenin, and luteolin also played a pivotal role when treating COVID-19 using Maxing Shigan decoction. Kong et al. [44] discovered through molecular docking that kaempferol, quercetin, luteolin, and new coronavirus (SARS-CoV-2) 3CL hydrolase have good affinity.There are 77 candidate targets for YQS in COVID-19 treatment. TNF, GAPDH, MAPK3, MAPK1, EGFR, CASP3, MAPK8, mTOR, IL-2, and MAPK14 are important targets for YQS in COVID-19 treatment. COVID-19 is similar to other viral pneumonia conditions. Some scholars have found that cytokine storms (CSs) may affect the progression of COVID-19 illness [45]. The cytokines involved in CS mainly include six categories: tumor necrosis factor (TNF), interleukin (IL), interferon (IFN), colony stimulating factor (CSF), and growth factor (GF). The cytokine storm induced by SARS-CoV-2 includes the cytokine IL-2, and the TNF-α- and H7N9-induced cytokine storms include cytokines, including IL-2 [46]. Studies have shown that, in throat swabs and sputum specimens of 4 cases of COVID-19, the housekeeping gene GAPDH of human cells showed a typical amplification signal curve [47]. Staphylococcus aureus pneumonia and virulence factor A can stimulate the AMPK signaling pathway and inhibit the mTOR pathway to induce autophagy in alveolar epithelial cells, an important mechanism during early infection of S. pneumoniae [48]. Liu et al. [49] predicted MAPK1, MAPK3, and MAPK14 through network pharmacology as important targets for Dayuanyin to treat COVID-19. Huang et al. [4] predicted that MAPK1 and MAPK14 were important targets for Huanglian Jiedu decoction to treat COVID-19. Other studies have shown that EGFR may be a key factor in the inflammation caused by Klebsiella pneumoniae infection [50].Furthermore, molecular docking results showed that the 5 compounds could efficiently act on the 10 targets, especially TNF, MAPK14, MAPK3, and MAPK8. Kaempferol and Mpro have a better binding capacity. Acacetin, luteolin, and ACE2 have a better combining ability. Whether they are docked with molecules of hub genes or docked with the proven targets of pneumonia, this shows the feasibility of using these five compounds to treat pneumonia. Reportedly, AGEs modulate COVID-19 pathogenesis and related comorbidities through their interactions with RAGE, among other molecules [51]. Interestingly, Yalcin Kehribar et al. [52] found that the RAGE pathway plays an important role in the aggravation of COVID-19. These reports provide strong evidence for our data.To clarify the role of YQS in COVID-19 treatment based on gene functions and signaling pathways, this study performed GO functional enrichment and pathway enrichment analysis of the candidate targets. The GO enrichment analysis showed that the biological procedure of YQS against COVID-19 was mainly reflected in the positive regulation of the adhesion of white blood cells and vascular endothelial cells and the negative regulation of the apoptosis process of cholangiocarcinoma cells. The cellular components were mainly reflected in cytoplasmic vesicles, membrane rafts, and organelle membranes. The molecular functions were mainly reflected in phosphotransferase activity, alcohol receptor, and cytokine receptor binding. The results of KEGG pathway enrichment analysis show that YQS may comprehensively play a role in preventing and treating pneumonia by way of PTGS2, MAPK14, CCR1, CCR3, MAPKAPK2, MAPK3, MAPK1, ICAM1, JAK1, NFKB1, MAPK8, PIK3R1, MTOR, RELA, CASP3, FGF2, TBK1, EIF2AK2, CASP8, and CDK4 targets that are related to the Kaposi sarcoma-associated herpesvirus infection signaling pathway; TGFB1, MAPK14, SERPINE1, STAT1, PIK3CA, PIK3CB, PIK3CD, MAPK3, MAPK1, PIK3R1, NFKB1, RELA, ICAM1, PRKCA, PRKCB, PRKCE, PRKCZ, MAPK8, CASP3, TNF, CDK4, and BCL2 targets that are related to the AGE-RAGE signaling pathway; and TGFB1, MAPK14, SERPINE1, STAT1, PIK3CA, PIK3CB, PIK3CD, MAPK3, MAPK1, PIK3R1, NFKB1, RELA, ICAM1, PRKCA, PRKCB, PRKCE, PRKCZ, MAPK8, CASP3, TNF, CDK4, and BCL2 targets that are related to the pathways in cancer signaling pathways.


## 5. Summary

In summary, our study found that active ingredients such as luteolin, naringenin, farnesin, kaempferol, and quercetin can act on targets such as TNF, GAPDH, MAPK3, MAPK1, EGFR, CASP3, MAPK8, MTOR, IL-2, and MAPK14 to regulate signaling pathways such as Kaposi sarcoma-associated herpesvirus infection and AGE-RAGE, finally achieving the effect of curing COVID-19.In the next step, further in vitro and in vivo experimental validation will be carried out to advance our research. The follow-up plan of our research group is to carry out more in-depth experimental research on YQS and to fully explore the mechanism of its efficacy.



---

*Source: 9248674-2022-03-15.xml*