3.1 Retrieval of chemical compositions, effective components and targets
There were a large number of pharmaceutical chemical components in SGMD that were retrieved by TCMSP. A total of 194 effective ingredients were obtained after the second retrieval based on the screening conditions of bioavailability (OB) > 30% and drug-like property (DL) > 0.18.Among them, 13 in Pinelliae Rhizoma, 5 in Citri Reticulatae Pericarpium, 21 in Codonopsis Radix, 20 in Astragali Radix, 7 in Atractylodis Macrocephalae Rhizoma, 15 in Poria, 11 in Pogostemonis Herba, 10 in Amomi Fructus and 92 in Glycyrrhizae Radix et Rhizoma were found. Due to the multi-component and multi-target characteristics of traditional Chinese medicine, there exists the phenomenon that multiple traditional Chinese medicines contain the same component. 173 effective ingredients and 253 corresponding action targets were obtained after removing duplicates. The basic information of some of the effective ingredients of SGMD is shown in table 1.
Table 1 Basic information of active ingredients of SGMD
Herbs
|
MolID
|
MolName
|
OB
|
DL
|
Atractylodis Macrocephalae Rhizoma
|
MOL000020
|
12-senecioyl-2E,8E,10E-atractylentriol
|
62.4
|
0.22
|
|
MOL000021
|
14-acetyl-12-senecioyl-2E,8E,10E-atractylentriol
|
60.31
|
0.31
|
|
MOL000022
|
14-acetyl-12-senecioyl-2E,8Z,10E-atractylentriol
|
63.37
|
0.3
|
Pinelliae Rhizoma
|
MOL001755
|
24-Ethylcholest-4-en-3-one
|
36.08
|
0.76
|
|
MOL002670
|
Cavidine
|
35.64
|
0.81
|
|
MOL002714
|
baicalein
|
33.52
|
0.21
|
|
MOL002776
|
Baicalin
|
40.12
|
0.75
|
|
MOL000358
|
beta-sitosterol
|
36.91
|
0.75
|
|
MOL000449
|
Stigmasterol
|
43.83
|
0.76
|
Citri Reticulatae Pericarpium
|
MOL005100
|
5,7-dihydroxy-2-(3-hydroxy-4-methoxyphenyl) chroman-4-one
|
47.74
|
0.27
|
|
MOL005815
|
Citromitin
|
86.9
|
0.51
|
|
MOL005828
|
nobiletin
|
61.67
|
0.52
|
Codonopsis Radix
|
MOL002140
|
Perlolyrine
|
65.95
|
0.27
|
|
MOL005321
|
Frutinone A
|
65.9
|
0.34
|
|
MOL000006
|
luteolin
|
36.16
|
0.25
|
|
MOL007059
|
3-beta-Hydroxymethyllenetanshiquinone
|
32.16
|
0.41
|
|
MOL008400
|
glycitein
|
50.48
|
0.24
|
Poria
|
MOL000296
|
hederagenin
|
36.91
|
0.75
|
Pogostemonis Herba
|
MOL005573
|
Genkwanin
|
37.13
|
0.24
|
|
MOL005911
|
5-Hydroxy-7,4'-dimethoxyflavanon
|
51.54
|
0.27
|
|
MOL005916
|
irisolidone
|
37.78
|
0.3
|
|
MOL005918
|
phenanthrone
|
38.7
|
0.33
|
|
MOL005921
|
quercetin 7-O-β-D-glucoside
|
49.57
|
0.27
|
Astragali Radix
|
MOL000239
|
Jaranol
|
50.83
|
0.29
|
|
MOL000354
|
isorhamnetin
|
49.6
|
0.31
|
|
MOL000371
|
3,9-di-O-methylnissolin
|
53.74
|
0.48
|
|
MOL000378
|
7-O-methylisomucronulatol
|
74.69
|
0.3
|
|
MOL000379
|
9,10-dimethoxypterocarpan-3-O-β-D-glucoside
|
36.74
|
0.92
|
|
MOL000380
|
(6aR,11aR)-9,10-dimethoxy-6a,11a-dihydro-6H-benzofurano [3,2-c] chromen-3-ol
|
64.26
|
0.42
|
|
MOL000387
|
Bifendate
|
31.1
|
0.67
|
|
MOL000392
|
formononetin
|
69.67
|
0.21
|
|
MOL000417
|
Calycosin
|
47.75
|
0.24
|
|
MOL000422
|
kaempferol
|
41.88
|
0.24
|
|
MOL000442
|
1,7-Dihydroxy-3,9-dimethoxy pterocarpene
|
39.05
|
0.48
|
Amomi Fructus
|
MOL001755
|
24-Ethylcholest-4-en-3-one
|
36.08
|
0.76
|
|
MOL001771
|
poriferast-5-en-3beta-ol
|
36.91
|
0.75
|
Glycyrrhizae Radix et Rhizoma
|
MOL004941
|
(2R)-7-hydroxy-2-(4-hydroxyphenyl) chroman-4-one
|
71.12
|
0.18
|
|
MOL001792
|
DFV
|
32.76
|
0.18
|
|
MOL004835
|
Glypallichalcone
|
61.6
|
0.19
|
|
MOL004841
|
Licochalcone B
|
76.76
|
0.19
|
|
MOL004985
|
icos-5-enoic acid
|
30.7
|
0.2
|
|
MOL004996
|
gadelaidic acid
|
30.7
|
0.2
|
|
MOL003896
|
7-Methoxy-2-methyl isoflavone
|
42.56
|
0.2
|
|
MOL000500
|
Vestitol
|
74.66
|
0.21
|
|
MOL004957
|
HMO
|
38.37
|
0.21
|
|
MOL004328
|
naringenin
|
59.29
|
0.21
|
|
MOL000392
|
formononetin
|
69.67
|
0.21
|
|
MOL000422
|
kaempferol
|
41.88
|
0.24
|
|
MOL000417
|
Calycosin
|
47.75
|
0.24
|
|
MOL000098
|
quercetin
|
46.43
|
0.28
|
3.2 Retrieval of COVID-19 related targets
The GeneCards and NCBI databases were used to retrieve the keyword "novel coronavirus". In total, 346 related targets were found in the GeneCards database, and 48 related targets were found in the NCBI database, among which 46 were common targets. After removing common targets, 348 targets related to COVID-19 were obtained.
3.3 Potential targets of COVID-19 in drug therapy
The Venny mapping software was used to perform intersection alignment of 253 effective ingredients regulated targets retrieved from the TCMSP database, and a total of 348 targets related to COVID-19 identified from the GeneCards and NCBI databases. 50 cross-targets were obtained, which were potential targets for the treatment of COVID-19 by SGMD, as shown in fig.1.
3.4 Analysis of regulatory network of TCM compounds
After the cross-targets were imported into Cytoscape 3.7.1 software (http://www. cytoscape.org/), the TCM compound regulatory network was obtained, with 168 nodes and 588 edges, as shown in fig. 2. The 118 effective ingredients in SGMD acted on 50 cross-targets of drugs and diseases to treat COVID-19. In fig. 2, the ginger yellow nodes represent the effective ingredients of Pinelliae Rhizoma, the grey nodes represent the effective ingredients of Citri Reticulatae Pericarpium, and the purple nodes represent the effective ingredients of Codonopsis Radix. The light-yellow nodes represent the effective ingredients of Astragali Radix, and the yellow node represents the effective ingredients of Atractylodis Macrocephalae Rhizoma, while the brown node represents the effective ingredients of Poria, and the pink node represents the effective ingredients of Pogostemonis Herba. The blue node represents the effective ingredients of Amomi Fructus, and the green node represents the effective ingredients of Glycyrrhizae Radix et Rhizoma, while the red nodes represent the chemical components that exist in multi-herb Chinese medicine and act on COVID-19. The sky-blue nodes represent the potential targets of drug treatment, and each edge represents the interaction between the effective ingredients and the potential therapeutic targets. There were 8 red nodes in total. MOL000239-Jaranol, MOL000354-isorhamnetin, MOL000392-formononetin, MOL000417-claycosin, and MOL000422-kaempferol were the common components of Astragali Radix and Glycyrrhizae Radix et Rhizoma; MOL004328- naringin was the common component of Citri Reticulatae Pericarpium and Glycyrrhizae Radix et Rhizoma; MOL003896-7-methoxy-2-methyl isoflavone was the common component of Codonopsis Radix and Glycyrrhizae Radix et Rhizoma; and MOL000098- quercetin was the common component of Pogostemonis Herba, Astragali Radix and Glycyrrhizae Radix et Rhizoma. These common components are all flavonoids. The first 15 active ingredients, which were sorted in descending order of degree, were quercetin, luteolin, kaempferol, naringin, licochalcone A, nobiletin, irisolidone, baicalein, formononetin, isorhamnetin, 7-O-methylisomucronulatol, Odoratin, Vestitol, 7-Acetoxy-2-methylisoflavone and Quercetin der./3,3'-di-O-methylquercetin. These compounds are also flavonoids, which are key components in the treatment of COVID-19.
3.5 PPI network
Fifty drug and disease cross-targets were imported into the Bisogenet plug-in, and the protein–protein interaction (PPI) network was obtained after the network analysis, as shown in fig. 3. A total of 2,133 target proteins were loaded.
3.6 Network topology analysis
The CytoNCA plug-in was used to carry out the topological analysis of target proteins in the PPI network. In order of degree, 101 target proteins and sub-network 1 were obtained by filtering with DC>61 as the screening condition (fig. 4b). The resulting target proteins were filtered twice in betweenness order, and 28 target proteins and sub-network 2 were obtained by filtering with BC>100 as the screening condition (fig. 4c). The 28 targets were TP53, EGFR, SRC, AR, ABL1, GRB2, ABL1, GRB2, AKT1, MAPK1, HSP90AB1, BRCA1, CTNNB1, RAF1, MAP3K3, ESR1, SMAD2, HSP90AA1, YWHAZ, UBC, LYN, STAT3, HSPA8, PIK3R1, RB1, SMAD3, FYN, CASP3, IKBKG and SP1. As the core targets in the network, they are the most closely related to the treatment of COVID-19 by SGMD.
3.7 GO function analysis and KEGG pathway analysis
SGMD treatment targets related to COVID-19 were entered into Bioconductor for the GO functional analysis and KEGG enrichment analysis. P value <0.05 was set as the default option, and the first 20 lines of the screening results were plotted into a histogram. The GO functional analysis includes the biological process (BP), molecular function (MF) and cell composition (CC). According to panels A, B and C in fig. 5, the relevant targets for the treatment of COVID-19 were mainly concentrated in biological processes such as responses to metal ions, molecules of bacterial origin, lipopolysaccharides, toxic substances and oxidative stress. The main molecular functions were concentrated in cytokine receptor binding, phosphatase binding, BH domain binding, cytokine activity and protein phosphatase binding. The main cellular components were concentrated in membrane raft, membrane microdomain and the membrane region. According to fig. 5d, the relevant targets for the treatment of COVID-19 were produced by treating the effect mainly through the AGE−RAGE signaling pathway in diabetic complications as well as the TNF and IL-17 signaling pathways.
3.8 Target -KEGG pathway network
Cytoscape 3.7.1 software was used to map the 20 pathways and path-related targets obtained from the KEGG enrichment analysis in order to understand the pathway mechanism. The hsa04933-AGE−RAGE signaling pathway in diabetic complications as well as the hsa04668-TNF and hsa04657-IL−17 signaling pathways are of great importance in the enriched pathways. MAPK1, MAPK3, RELA, IL6, MAPK14 and other targets also play important roles in this network (fig. 6).
3.9 Molecular docking
15 key components, including quercetin, luteolin, kaempferol, naringin, licochalcone A, nobiletin and irisolidone were connected with the SARS-CoV-2 3CL hydrolase and angiotensin converting enzyme 2 (ACE2) using the AutoDock 4.2.6 software. It is generally believed that the lower the binding energy of the ligand receptor, the easier the binding. According to the binding energy, the first three core components with the lowest molecular weight in connection with 6LU7 and 1R42 were identified as formononetin, Quercetin der., 7-Acetoxy-2-methylisoflavone and 7-Acetoxy-2-methylisoflavone, irisolidone and formononetin, respectively. Formononetin and 7-Acetoxy-2-methylisoflavone are bound by hydrogen bonds with SER-10, ALA-7, and LYS-5. Quercetin der. is bound by hydrogen bonding between SER-10, ALA-7 and 6LU7. 7-Acetoxy-2-methylisoflavone is bound by hydrogen bonds formed between IL-256, MET-249 and 1R42, while irisolidone and formononetin is formed between SER-254 and ASN-250. The docking mode is shown in fig. 7.