Exploration on Shufeng Jiedu Capsule for Treatment of COVID-19 Based on Network Pharmacology and Molecular Docking ()
1. Introduction
Coronavirus disease 2019 (COVID-19) is caused by a novel coronavirus and is highly infectious. It was named by the World Health Organization on January 12, 2020, and its pathogen was named SARS-CoV-2 [1]. From the traditional Chinese medicine angle, COVID-19 is a blight [2] [3]. Chinese medicine therapy is recommended from the third edition of “Novel Coronavirus Pneumonia Diagnosis and Treatment Scheme” issued by National Health Commission of the People’s Republic of Chian and National Administration of Traditional Chinese Medicine, and achieved good results in treatment [4] [5] [6] [7]. Shufeng Jiedu Capsule was recommended from the sixth edition of “Novel Coronavirus Pneumonia Diagnosis and Treatment Scheme”; the formula was from “Qudu San”, which is made from Polygoni Cuspidati Rhizoma Et Radix (Huzhang), Forsythiae Fructus (Lianqiao), Isatidis Radix (Banlangen), Radix Bupleuri (Chaihu), Herba Patriniae (Baijiangcao), Verbenae Herb (Mabiancao), Phragmitis Rhizoma (Lugen), and licorice (Gancao). The capsule has the functions of expelling wind, clearing heat and detoxification. Clinical studies showed that the capsule had antibacterial, antiviral and immunomodulatory effects, and had a good effect in the treatment of acute upper respiratory tract infection and community-acquired pneumonia [8] [9] [10] [11]. In this study, we selected relevant targets and pathological mechanisms of Shufeng Jiedu Capsule for the treatment of COVID-19 by network pharmacology and molecular docking [12], and then provided the basis for further study.
2. Materials and Methods
2.1. The Main Active Components and Corresponding Targets of Shufeng Jiedu Capsule
Searched names of 8 the traditional Chinese medicine (TCM) in Shufeng Jiedu Capsule as Polygoni Cuspidati Rhizoma Et Radix (Huzhang), Forsythiae Fructus (Lianqiao), Isatidis Radix (Banlangen), Radix Bupleuri (Chaihu), Herba Patriniae (Baijiangcao), Verbenae Herb (Mabiancao), Phragmitis Rhizoma (Lugen), licorice (Gancao), selected the active components and corresponding targets with filter conditions as oral bioavailability (OB) ≥ 30, drug-like (DL) ≥ 0.18.
2.2. Constructed Drug-Compound-Target Network
The gene symbols of the active components of Shufeng Jiedu Capsules were extracted from UniProtKB (http://www.uniprot.org), and the target was limited to Homo sapiens. Then construct PPI network with filter conditions as median confidence > 0.4. Then construct drug-compound-target network by software Cytoscape 3.6.1, then select potential active components of Shufeng Jiedu Capsule against COVID-19 according to the node degree method [13].
2.3. Gene Ontology Analysis and Pathway Analysis
The targets of active components of Shufeng Jiedu Capsule were analyzed by Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. GO analysis included Biological Process (BP), Molecular Function (MF), Cellular Component (CC). Two analyses were worked on software R 5.3.2 with package of “org.Hs.eg.db”, “clusterProfiler”, “enrichplot” and “ggplot2” [14], adjust P < 0.05.
2.4. Molecular Docking
Download the active ingredients of Shufeng Jiedu Capsule from TCMSP and saved them in mol2 format, then saved them in PDBQT format through software Auto Dock. The 3D structure (PDB ID: 6LU7) of SARS-CoV-2 3CL hydrolase and angiotensin converting enzyme II (ACE2) download from RCSB database (http://www.rcsb.org/), then saved in PDB format, removed water and other small molecules, added hydrogenation through software Auto Dock, then saved in PDBQT format. Made molecular docking by software Vina (http://vina.scripps.edu/index.html), then selected the minimum binding energy of drug and protein as the optimal site, viewed and saved images by software Discovery Studio 2016.
3. Results
3.1. The Main Active Components and Corresponding Targets of Shufeng Jiedu Capsule
Selected 207 components of Shufeng Jiedu Capsule through TCMSP platform, included 10 of Polygoni Cuspidati Rhizoma Et Radix (Huzhang), 23 of Forsythiae Fructus (Lianqiao), 39 of Isatidis Radix (Banlangen), 17 of Radix Bupleuri (Chaihu), 13 of Herba Patriniae (Baijiangcao), 12 of Verbenae Herb (Mabiancao), 1 of Phragmitis Rhizoma (Lugen), 92 of Phragmitis Rhizoma (Gancao).
3.2. Constructed Drug-Compound-Target Network
Drug-compound-target network of Shufeng Jiedu Capsule includes 400 nodes, 7756 edges, 8 TCM, 155 active compounds, 236 targets. As shown in Figure 1, green—TCM, yellow—compounds, red—targets. Our study also selected key nods according to degree of topological analysis, including 3 TCM, 45 targets, 26 compounds (shown in Table 1).
Figure 1. Drug-compound-target network of Shufeng Jiedu Capsule.
Table 1. Main active components of Shufeng Jiedu Capsule.
3.3. Gene Ontology Analysis and Pathway Analysis
We got 2363 GO items through Gene ontology analysis, including 2334 BP items, 103 CC items, 198 MF items. We also got 177 pathways through KEGG analysis, including Human cytomegalovirus infection, PI3K-Akt signaling pathway and Hepatitis B, etc. The results are shown in Figure 2.
3.4. Molecular Docking
Molecular docking with SARS-CoV-2 3CL hydrolase was carried out between the 26 main active components of Shufeng Jiedu Capsule selected above. The lower the binding energy, the greater the possibility of interaction between the active compound and the hydrolase. The binding energy of 26 compounds is from −32.21 to −20.94 kJ·mol−1 (shown in Table 2). Bicuculline (MOL000791) had the greatest effect with hydrolase (−32.21 kJ·mol−1), which is the active compound of Forsythiae Fructus (Lianqiao), the second was MOL004978 (−30.54 kJ·mol−1), the third was MOL004966 (−30.12 kJ·mol−1), the results are shown in Figure 3.
4. Discussion
Currently, there is no specific drug treatment for COVID-19, and existing chemical drugs can only inhibit some clinical symptoms within a limited range, but within higher cure rate and lower mortality rate of confirmed patients with TCM treatment [15]. The sixth edition of “Novel Coronavirus Pneumonia Diagnosis and Treatment Scheme” advised different TCM treatments according to TCM dialectics in treatment in patients in different pathological stages, and Shufeng Jiedu Capsule was recommended for prevention and treatment during medical observation stage.
Table 2. The binding energy of 26 compounds and SARS-CoV-2 3CL hydrolase.
Note: 1) HZ—Polygoni Cuspidati Rhizoma Et Radix, LQ—Forsythiae Fructus, Isatidis BLG—Radix, CH—Radix Bupleuri, BJCb—Herba Patriniae, MBC—Verbenae Herb, LG—Phragmitis Rhizoma, licorice—GC. 2) √—contain compound, — —free from compound.
GO KEGG
Figure 2. Gene ontology analysis and KEGG pathway analysis of Shufeng Jiedu Capsule. (The picture above is gene ontology analysis, the picture below is KEGG pathway analysis. In the figure, the gradient from blue to red represents that the smaller the P value is, and the gradient from small to large represents the larger the number of enrichment targets.)
Figure 3. Molecular docking mode of 3 active compounds in Shufeng Jiedu Capsule and SARS-CoV-2 3CL hydrolase.
In this study, through network pharmacology and molecular docking, we got 26 active compounds of Shufeng Jiedu Capsule, and binding energy of compounds and SARS-CoV-2 3CL hydrolase were from −32.21 to −20.94 kJ·mol−1. Quercetin (−25.52 kJ·mol−1) in our study is similar to the published research [13], whose binding energy with SARS-CoV-2 3CL hydrolase was −26.36 kJ·mol−1. The results of molecular docking are shown that combination of Shufeng Jiedu Capsule and hydrolase includes hydrogen bond, pi-pi conjugate and van der Waals force. Moreover, the capsule contains many active compounds, then it can be predicted that Shufeng Jiedu Capsule has a high possibility of anti-COVID-19.
Our study selected the main compounds and targets of Shufeng Jiedu Capsule by network pharmacology, then analyzed the potential molecular mechanism on its binding energy with SARS-CoV-2 3CL hydrolase by molecular docking. The results showed that Shufeng Jiedu Capsule acts on hydrolase by multiple compounds, targets and pathways. The method in this study provides the basis in mechanism of action of Chinese herbal compound preparation. However, this study still has some limitations, such as insufficient understanding of the pathogenesis of COVID-19 and insufficient research on the key targets of the viral infection process. In this study, the pharmacodynamics, evaluation and clinical efficacy of COVID-19 were studied and discussed at a theoretical level. In the later stage, the pharmacodynamics, evaluation and clinical efficacy should be further studied to provide a basis for further development and utilization of COVID-19 therapeutic drugs.
Funding
This study was funded by the Study and Development Fund for Sciences and Technology in Chengde City (No. 201701A087).