key: cord-0704320-ee6xeha7
authors: Wu, Huaying; Gong, Ke; Qin, You; Yuan, Zhiying; Xia, Shuaishuai; Zhang, Shiying; Yang, Jingjing; Yang, Ping; Li, Liang; Xie, Mengzhou
title: In silico analysis of the potential mechanism of a preventive Chinese medicine formula on coronavirus disease 2019
date: 2021-07-15
journal: J Ethnopharmacol
DOI: 10.1016/j.jep.2021.114098
sha: 871072c8f969c9d99e7cf4755a286f9ac2c11ba9
doc_id: 704320
cord_uid: ee6xeha7

ETHNOPHARMACOLOGICAL RELEVANCE: With the spread of Coronavirus Disease (2019) (COVID-19), combination with traditional Chinese medicine (TCM) has been widely used as a prevention and therapy strategy in China. Xin guan No.1 (XG-1) prescription is a preventive formula recommended by the Hunan Provincial Administration of TCM to prevent the pandemic of COVID-19. AIM OF THE STUDY: To explore the potential preventive mechanisms of XG-1 against COVID-19 in the combination of network pharmacology approach, single-cell RNA expression profiling analysis, molecular docking and retrospective study. MATERIALS AND METHODS: Encyclopedia of Traditional Chinese Medicine (ETCM) database was used to determine the meridian tropism, active components and target genes of XG-1. Gene ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis were conducted by R Cluster Profiler package (3.14.3). Single cell RNA sequencing (scRNA-seq) data of human lung (GSE122960) was downloaded from Gene Expression Omnibus (GEO) database and analyzed by R Seurat package (3.1.2). Cytoscape (3.7.2) was used to construct the interaction network. The main ingredients in XG-1 were identified by HPLC- Q-TOF- MS and used for molecular docking with COVID-19 3CL hydrolytic enzyme and angiotensin converting enzyme II (ACE2). A retrospective study of 47 close contact participants from Dongtang Community of Hunan Province was conducted to evaluated the preventive effect of XG-1. RESULTS: According to the network pharmacology analysis, XG-1 formula was closely related to lung-, spleen- and stomach-meridians and include a total of 206 active components and 853 target genes. GO and KEGG pathway enrichment revealed that XG-1 mainly regulated cellular amino acid metabolism process and neuroactive ligand-receptors interaction. The scRNA-seq profiling showed that angiotensin converting enzyme 2 (ACE2) was principally expressed in alveolar type 2 epithelial cells (AT2). 153 genes were up-regulated in AT2 cells expressing ACE2 and 12 genes were obtained by intersecting with XG-1 target genes, of which 3 were related to immunity. Five main chemical ingredients were detected in XG-1 sample by HPLC-Q-TOF-MS. The molecular docking showed that Rutin, Liquiritin and Astragaloside Ⅳ had a good affinity with COVID-19 3CL hydrolytic enzyme and ACE2. Compared with participants who didn't take XG-1, preventive treatment with XG-1gradules resulted in a significant lower rate of testing positive for SARS-CoV-2 nucleic acid (P < 0.0001). CONCLUSION: The present study showed that XG-1 exerts a preventive effect in close contacts against COVID-19. The underlying mechanism may be related to modulate immunity response through multiple components, pathways, and several target genes co-expressed with ACE2. These findings provide preliminary evidences and methodological reference for the potential preventive mechanism of XG-1 against COVID-19.

In December 2019, an acute pneumonia caused by a new coronavirus (SARS-CoV-2) broke out, which named as coronavirus disease 2019 . The World Health Organization (WHO) declared a public health emergency of international concern on January 30, 2020 (World Health Organization, 2020a) . The COVID-19 epidemic has rapidly spread to the world and the current situation is still severe (World Health Organization, 2020b) . Clinical symptoms of COVID-19 often manifest as fever, dry cough, muscle pain or fatigue. A part of patients is accompanied by nasal congestion, runny nose, diarrhea and other upper respiratory or digestive tract symptoms. Patients usually died of acute respiratory distress syndrome, septic shock, coagulation dysfunction or multiple organ failure .

SARS-CoV-2 has been sequenced and it's reported that SARS-CoV-2 has a similar receptor binding domain (RBD) structure with SARS-CoV spike (S) protein, and with 76.5% homology in amino acid sequence Xu, X. et al., 2020b) . Human angiotensin-converting enzyme II (ACE2) has a strong binding affinity to SARS-CoV S protein and is a functional receptor for SARS-CoV where the virus enters host cells in vitro or in vivo and replicates the virus Li et al., 2003 Li et al., , 2005 . Research showed that SARS-CoV-2 can identify human ACE2 more effectively than SARS-CoV, which could regulate the capacity and efficiency of human-to-human transmission (Wan et al., 2020) . Overexpression of ACE2 in HeLa cells of different species (human, pig and civet) can cause infection and replication of SARS-CoV-2, which directly proves that SARS-CoV-2 and ACE2 binds to promote virus entry and affected replication on the surface of host cells .

Traditional Chinese medicine (TCM) has a long history of fighting against epidemic disease and established a complete theoretical system of infectious diseases intervention. Previous studies have shown that TCM therapy exerted an essential role in combating coronavirus pneumonia such as severe acute respiratory syndrome (SARS) and H1N1 influenza (Leung, 2007; Zhang and Chen, 2008) . Evidences also have demonstrated that TCM synergy therapy could alleviate clinical symptoms such as fever, control the progression of disease, reduce the dosage of hormone and the occurrence of complications in patients with COVID-19 Wu et al., 2020) . Xin guan No.1 (XG-1) is a preventive prescription proposed by Hunan Provincial Administration of TCM for the prevention of COVID-19 (Hunan Provincial Administration of Traditional Chinese Medicine, 2020). It's based on YuPing-Feng San, a classic and famous prescription with immune enhancement and inflammation inhibitory effects (Li, Y. et al., 2017b; Liu et al., 2017; Wang et al., 2016; Zhou et al., 2019) . It had been used to protect the at-risk hospital workers and showed a significant preventive effect during the SARS epidemic (Lau et al., 2005) . A recent study has demonstrated that XG-1 has a certain regulatory effect on the immune function of young and elderly healthy people (Liu, L. et al., 2020a) . XG-1 consists of 8 traditional Chinese herbs including Astragalus mongholicus Bunge (Huang Qi) (Shi Chang Pu) . At the peak of the COVID-19 epidemic in Hunan Province, close contacts with the confirmed cases were under high risks of infection. Although the XG-1 formula has been widely used in Hunan Province to prevent COVID-19 invasion, its bioactive ingredients and their underlying preventive mechanism remains to be elucidated.

Network pharmacology is a method based on systems biology and system pharmacology that analyzes the network of biological systems and selects specific signal nodes for drug molecule design. Chinese herbal compound is the main form of TCM for disease prevention and treatment, with the characteristics of multiple components, multiple targets and multiple pathways. It was demonstrated that network pharmacology has been widely applied to recognize the active ingredients and potential targets of herbals, to predict pharmacological mechanisms of herbal formulas (Wu et al., 2019; Yu et al., 2020) . Single-cell sequencing has been a hot topic in recent years and single-cell RNA sequencing (scRNA-Seq) technology is the most mature and widely used (Shalek et al., 2013) . Its purpose is to comprehensively measure the expression level of genes in cells, and to improve the accuracy of our understanding of cellular gene function (Stubbington et al., 2017) . Molecular docking is a computational tool used to predict the binding capacity and binding mode of receptor-drug molecular complexes. It's confirmed that the combination of the S protein of SARS-CoV-2 and ACE2 in the human body could help the virus invade the body and cause disease. The crystal structure of SARS-CoV-2 3CL hydrolase is reported and considered to be an effective target of SARS-CoV-2, which provides a basis for screening active ingredients against COVID-19 (Lu et al., 2020; Elmezayen et al., 2020) .

Therefore, in our study, network pharmacology was applied to screen out the active ingredients and action targets of XG-1, and a singlecell RNA expression profiling of human lung was performed to cluster analysis. The co-action targets between ACE2 expressing cells and XG-1 bioactive components were identified via computationally analysis. The main active components in XG-1 were identified by HPLC-Q-TOF-MS, and the binding energy of the key components with COVID-19 3CL hydrolase and ACE2 was predicted through molecular docking. A retrospective study was conducted to verify the preventive efficacy of XG-1 treatment in the close contacts. We hope our study can provide a preliminary evidence of the potential precautionary mechanism of XG-1 against COVID-19.

XG-1 consists of 8 herbs (Table 1) . "HUANG QI", "BAI ZHU", "LIAN QIAO", "GAN CAO", "FANG FENG", "JIN YIN HUA", "GUANG HUO XIANG" and "SHI CHANG PU" were searched in the Chinese Medicine Encyclopedia (ETCM, http://www.nrc.ac.cn:9090/ETCM/) database to obtain the contained compounds and the target genes (Xu et al., 2019) .

In order to further observe the biological function of the target, Gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of differentially expressed genes were performed using the Cluster Profiler package (3.14.3) in R platform (Yu et al., 2012) . The results with the threshold P < 0.05 and in the top 20 entries were visualized into a bubble chart.

The human lung tissue single cell transcriptome sequencing data (GSE122960) was downloaded from the GEO database, and the Seurat software package in R (3.6.2) ((http://satijalab.org/seurat/, version: 3.1.2) was used to perform cell clustering. Low-quality cells that express more than 10% of mitochondrial genes in cells with less than 200 RNAs were filtered out, then the 2000 most highly mutated genes were calculated and analyzed by PCA dimensionality reduction. For the top 10 PCs with explanatory power, Find Clusters was performed to identify each cell cluster, and t-distribution random neighbor embedding (tSNE) was used for clustering and visualization, then the initial clusters were checked and merged. Find All Markers function and Wilcoxon rank sum tests were used to analyze and screen differentially expressed marker genes. Differentially expressed genes were expressed in at least 25% of the cells in the cluster, and the fold change is greater than 0.25 (Log Scale), which was considered to be a marker gene. The tSNE curve and violin curve were generated using Seurat.

Cytoscape (http://www.cytoscape.org/, version: 3.7.2) was used to construct the network and analyze the network (Shannon et al., 2003) . In a network, nodes represent active ingredients, targets, or pathways.

HPLC-Q-TOF -MS was used to control quality and identify the main chemical constituents of XG-1 granules. Rutin (Lot: 1000801-201811, purity: 91.7%), Liquiritin (Lot:111610-201607, purity: 93.1%), Chlorogenic acid (Lot:110753-201817, purity: 96.8%) and Astragaloside IV (Lot:110781-201717, purity: 96.9%) standards were purchased from National Institutes for Food and Drug Control, China. Phillyrin (Lot:110821-200609, purity: 93.5%) standard was purchased from National Institute for the Pharmaceutical and Biological Products Control, China. All of these standards were used as quality control of XG-1.

Sample preparation: 0.3 g of XG-1 prescription was accurately weighed, then placed in a stoppered Erlenmeyer flask. 4 mL of water was precisely added and weighed then ultrasonicated for 40 min. Weighed again when it was cooled and the lost weight was supplied with water.

Shaked well and then filter through a 0.22 μm microporous membrane to obtain the sample solution.

Standard preparation: 1 mg of each standard substance was accurately weighed and placed in an EP tube then 4 mL of methanol was added to dissolve.

Chromatographic conditions: An Agilent 1290 HPLC system (Agilent, USA) with an Inertsil ODS-2 Rapid Resolution HT column (4.6 × 250 mm, 5 μm, Shimadzu, Japan) was used for HPLC analysis at the flow of 0.5 mL/min and the temperature of the column was 30 • C. The injection volume was 10 μL, with wavelength 210 nm, and full wavelength scan range was 190-400 nm. The gradient elution procedure was shown in Table 3 . Mass spectrometry conditions: All the samples were analyzed with ESI source in negative mode. The operating parameters of MS were: drying gas, 8 L/min; gas temperature, 320 • C; sheath gas temperature, 350 • C; nebulizer, 35 psig; sheath gas flow, 10 L/min; nozzle voltage, 1.0 kV; V Cap, 3.5 kV; fragmentor, 110 V; skimmer voltage, 65 V; OCT 1 RF Vpp, 750 V; MS scan range, m/z 100-1000 Da. The detection results were shown in Fig. 1 and the quantitative analysis information of active ingredients in XG-1 was shown in Table 4 .

AutoDock Vina 4.2.6 was used for molecular docking (Trott and Olson, 2010) . The key active ingredients of XG-1 as ligands and the SARS-CoV-2 3CL hydrolase and ACE2 as receptors. The PDB format file of 3D structure of SARS-CoV-2 3CL hydrolase (PDB ID: 6LU7) and ACE2 (PDB ID: 1R42) were downloaded from the RSCB PDB database (http://www.rcsb.org/). AutoDock Tools 1.5.6 was used to remove water molecules in the SARS-CoV-2 3CL hydrolase protein and ACE2, separate the protein, add non-polar hydrogen, calculate the Gasteiger charge, and save it as a pdbqt file. Three key active ingredients of XG-1, Rutin (Compound CID: 5280805), Liquiritin (Compound CID: 503737) and Astragaloside Ⅳ (Compound CID: 13943297), were downloaded from PubChem database (https://pubchem.ncbi.nlm.nih.gov) to obtain the SDF format file of 2D structure. After checking the spatial structure in Chem3D, the energy is minimized and converted to pdb format. Chem3D was used to convert them into a mol2 file and AutoDock Tools 1.5.6 program was used to load the structure, add atomic charges and assign atomic types. All flexible bonds are rotatable by default and saved in pdbqt format as a docking ligand.

The retrospective study in this research was reviewed and approved by the Ethics Committee of the Dongtang Community Health Service Center in Changsha (No. 202001), and informed consent was obtained. All information was kept anonymously and complied with national laws and institutional requirements.

Participants were selected that close contacts who resident in Dongtang community of Changsha City in Hunan Province between February 13 and March 12, 2020. The SARS-CoV-2 nucleic acid test of pharyngeal swab sample was conducted by real-time polymerase chain Table 2 Function and distribution information of 12 intersection genes between ACE2 expressed AT2 cells target genes and XG-1 target genes. reaction (RT-PCR) on the day when they're traced to be the close contact of the confirmed COVID-19 cases. Individuals only who possess negative nucleic acid results enrolled in our study. Close contacts should meet one of the following criteria: (1) contacts who live, study or work together with the COVID-19 confirmed patients;

(2) contacts who were unprotected exposed to the COVID-19 confirmed patients within 1 m within 2 days before the onset of patient's symptom. Exclusion criteria included:

(1) contacts with suspected COVID-19 cases;

(2) younger than 18 years old;

(3) pregnant woman. The diagnosis of confirmed or suspected COVID-19 cases was in accordance with the "Clinical Guidance of Diagnosis and Treatment for COVID-19 (Trial version 8)" (General Office of the National Health Commission, 2020).

All participants were subject to isolation and observation for 14 days according to relevant regulations. During the quarantine period, SARS-CoV-2 nucleic acid tests were performed by RT-PCR every 7 days. Each participant received the preventive treatment of XG-1 granules (dissolved in 150 mL boiled water with twice a day) for 14 days. The granules were purchased from the First Affiliated Hospital of Hunan University of Traditional Chinese Medicine and the quality was in accordance with the "Pharmacopoeia of the People's Republic of China (Edition, 2015)". Following indicators were observed and recorded: 1. The number of participants taking traditional Chinese medicine; 2. The outcomes of SARS-CoV-2 nucleic acid testing; 3. Adverse reactions during the treatment. Additionally, the baseline characteristics including age, gender, contact time and relationship of participants were collected.

Statistical analysis of the data was performed using the IBM Social Science Statistical Software Package (SPSS, version 24). Continuous variables were calculated with Mean ± standard deviation, and the percentage of categorical variables. Comparisons between groups were performed using Student's t-test or Mann Whitney U test (the normal distribution of the data has been verified by the Kolmogorov-Smirnov test) and the chi-square test for continuous variables (Fisher's exact test was used when appropriate).

An herb-meridian network was constructed based on the herb information of XG-1 compound, as shown in Fig. 2A . Among all 8 herbs in XG-1, the most affinity meridians were Lung (5), Spleen (5) and Stomach (5), followed by heart (4), Small intestine (1), Bladder (1) and Liver (1). Therefore, XG-1 compound mainly acted on the lungs, spleen and stomach system.

Totally 385 active ingredients and 2297 targets were found in XG-1 by searching ETCM database, including 27 active ingredients and 461 targets in Astragalus mongholicus Bunge (Huang Qi) (Guang Huo Xiang) and 18 active ingredients and 41 targets in Acorus calamus var. angustatus Besser (Shi Chang Pu) . After removing the duplicate results, a total of 206 active ingredients and 853 targets were obtained (Table 1) . By calculating the number of target genes in each Chinese medicine herb, there were 380 genes targeting 1 herb, 128 genes targeting 2 herbs, 72 genes targeting 3 herbs, 100 genes targeting 4 herbs, 60 genes targeting 5 herbs, 23 genes targeting 6 herbs, 53 genes targeting 7 herbs, and 27 genes targeting 8 herbs (Fig. 2B) . These 27 genes were as follows: CYP1A1, CYP1A2, CYP2C8, CYP2C9, CYP3A4, GABRA1, GABRA2, GABRA3, GABRA4, GABRA5, GABRA6, CYP1B1 CYP2C19, CYP2D6, GABRB1, GABRB2, GABRB3, CYP2A6, GABRD, GABRE, GABRP, GABRQ, GABRG1, GABRG2, GABRG3, ALB, and CYP2B6.

GO analysis and KEGG pathway enrichment analysis were performed based on the 853 target genes screened from the database. In the Top 20 results of GO functional enrichment analysis, the biological process (BP) mainly involving organic anion transport, cellular amino acid Table 3 The gradient elution procedure in high-performance liquid chromatography. (Fig. 3) . KEGG pathway enrichment screened 100 signal pathways (P < 0.05), and the Top 20 results showed that the targets of XG-1 were mainly related to neuroactive ligand-receptors, glutamate synapses, steroid hormone synthesis, drug metabolism, calcium signaling pathway, and PPAR signaling pathways (Fig. 4) .

The tSNE classification results showed 8 types of cells, of which 2038 were alveolar type II epithelial cells (AT2), accounting for 54%; 1055 alveolar macrophages, accounting for 28%; 191 monocytes, accounting for 5%; 126 endothelial cells, accounting for 3%; 119 club cells, accounting for 3%; 112 alveolar type I cells (AT1), accounting for 2%; 55 ciliated cells, accounting for 1%; 51 dendritic cells, accounting for 1% ( Fig. 5 A and B) . As a functional receptor for SARS-CoV S protein, ACE2 is highly expressed in human lung AT2 cells (Fig. 5C ). The heat map shows the expression of marker genes in 8 cell clusters, which can observe whether the marker gene expression is specific and determine the cell type based on the marker gene. 585, 858, 564, 826, 168, 1064, 1498 and 154 marker genes were highly expressed in cluster 0 to cluster 7 respectively (Supplementary materials). Representative high expressed marker genes in each cell cluster are displayed (Fig. 6) .

By co-expression analysis of the ACE2 overexpressed AT2 cells and 853 Chinese medicine target genes, 12 intersecting target genes were obtained: PLA2G1B, SFTPD, MAOA, SLCO4C1, SLC22A3, ATP1A1, SLC7A2, RNASE1, ASPH, SQLE, CA2, DBI, with their expression levels in different cells visualized by Violin diagrams (Fig. 7) . There're some gene-gene networks were established connecting with these interaction genes except PLA2G1B. In the networks, ABL1, ITGA5, CAD, PPT1, GABRE, LGALS3 and SLC25A32 were found to be the target gene in XG-1 (Fig. 8) . The expression distribution information of the 12 intersection target genes in AT2 cells was shown in Table 2 . Among them, PLA2G1B, SFTPD and SLCO4C were related to immunity function.

In this study, three key active compounds of Rutin, Liquiritin and Astragaloside Ⅳ were considered as ligands and the SARS-CoV-2 3CL hydrolase and ACE2 were receptors. Their interaction with surrounding key residues and their binding at the active site are shown in Figs. 9 and 10, and the energy values of the compounds shown in the docking results are shown in Table 5 . The smaller the docking energy means the more stable binding. From the results, it showed that the binding energy of Rutin, Liquiritin and Astragaloside Ⅳ to the SARS-CoV-2 3CL hydrolase (6LU7) were − 7.9, − 7.7 and − 8.6 kJ/mol respectively, and their binding energy to the ACE2 (1R42) were − 8.5, − 8.5 and − 8.8 kJ/mol respectively, indicating these compounds and the receptors may bind and form a stable conformation.

A total of 47 close contacts were included in this study and the baseline characteristics were summarized in Table 6 . Of total contacts, 18 (38.3%) contacts were males. 36 (76.6%) contacts were adults aged 20-59 years; 12 (25.5%) contacts were spouse and 21 (44.7%) are nonspouse family members of the index cases. Regarding to the infection, 7 (14.9%) contacts were identified to be infected with SARS-COV-2, which included 6 (20.7%) females. The group aged 10-19 years, 50-59 years and over 70 years shared 33.3% attack rate. In addition, the attack rate differs from relationship with index cases: 16.7% for spouse; 19.1% for non-spouse family members; 12.5% for relatives. Participants with more contact frequency has a higher attack rate:18.8% for often; 7.7% for moderate, etc. XG-1 granules were distributed to everyone, the medication taking information was recorded and the data was analyzed. There're 36 participants (76.6%) took XG-1, while 11 (23.4%) refused to take Chinese medicine treatment during the quarantine time. There's no significance in the gender between the participants with taking XG-1 group and without (P = 0.12), while we observed that younger people are more likely to take XG-1 than older people (P = 0.01). Interestingly, regarding to the SARS-CoV-2 nucleic acid testing, of 11 participants who refused to take XG-1, 7 (63.64%) participants were found to be positive, while participants who took XG-1 were all confirmed negative results (P < 0.001). During the quarantine time, few participants expressed some slight symptoms and signs such as headache and insomnia. The characteristics were shown in Table 7 . We followed up all the close contacts who took XG-1 and found nobody detected nucleic acid positive and no 

With the spread of COVID-19 epidemic, China suffered a major public health emergency. In the early stage of the outbreak, there're no specific antiviral drugs for COVID-19. TCM treatment combined with symptomatic treatment was the main strategy for patients in China. Studies referring to TCM treatment in combating COVID-19 showed significant therapeutic effects especially in early use (Ai et al., 2020; Shi et al., 2020; Wu et al., 2020) . XG-1is a recommended preventive TCM formula and has been widely used to prevent the epidemic. To explore the material basis for the preventive effect of XG-1, network pharmacology approach was integrated to investigate the bioactive ingredients and potential mechanism. According to the network pharmacology study, the meridian information of each herb in XG-1 compound prescription was obtained and a meridian-herbs network of XG-1 was constructed. The meridian theory refers to the position of the herbs (or its chemical ingredients) action in the body, which is summarized through long-term observation of the curative effect in TCM (Zhao et al., 2008) . Patients with COVID-19 manifest clinical symptoms such as fever, dry cough, fatigue, loss of appetite and nausea, as well as the abnormal radiological changes in lung tissue (Guan et al., 2020; Yang et al., 2020) , while in the meridian network, 5 of 8 herbs in XG-1 have a specific relationship with the meridian of lung, spleen and stomach, indicating that XG-1 could acts on the lungs, spleen and stomach system to alleviate symptoms. Moreover, a total of 206 active ingredients and 853 target genes were obtained by screening the ETCM database. The target genes were collected to perform enrichment analysis. GO enrichment analysis revealed multiple pathways related to biological processes, molecular functions, and cell components; the enrichment results of KEGG pathway showed that XG-1 may exert effects via multiple signaling pathways such as neuroactive ligand-receptor, glutamate synapse, drug metabolism, calcium signaling pathway, and PPAR signaling pathway. Lin et al. stated upregulated gene expressions related to the neuroactive ligand-receptor interaction pathway and elevated internal Ca 2+ were detected after H5N1 infection (Lin et al., 2015) . Studies found that calcium signaling pathway could control the production of TNF-α and IL-10, which played a key role in immunity function regulation (Planes et al., 2018) . Furthermore, calcium signaling pathway, as a mechanosensitive ion channel, could impact the cell stretch to stimulate surfactant secretion in AT II cells (Willner et al., 2003) . Another study proved that PPAR expressed in lung macrophages and could affect the development of acute host disease and the restoration of tissue homeostasis after respiratory virus infection (such as respiratory syncytial virus) by regulating lung inflammation .

ACE2 has been identified as a functional receptor that could effectively bind to the S1 domain of the SARS-CoV spike protein, helping SARS virus to enter host cells and perform replication (Li et al., 2003) . Recent studies have confirmed that SARS-CoV-2/2019-nCoV used the same cell receptor-ACE2-as SARS-CoV Zhou et al., 2020) . ACE2, as a homologue of ACE, could negatively regulate the renin-angiotensin system by balancing ACE activity, causing vasoconstriction and exerting multiple biological functions, which helps to prevent damages from acute lung injury Kuba et al., 2006) . A previous study has showed that ACE2 was abundantly distribute in the epithelia of the lung and small intestine in humans (Hamming et al., 2004) . A recent scRNA-seq analysis from normal lung tissues of 8 adult donors proved that 83% of ACE2 expressing cells were alveolar type II epithelial cells (AT2), which believed that these cells could serve as a possible route for virus entry. In addition, GO enrichment analysis found AT2 expressing ACE2 had high expression of multiple genes related to viral processes, including regulatory genes for viral processes, viral life cycle, viral assembly and viral genome replication . Therefore, AT2, as ACE2 expressing target cells, may be the most common and easily damaged cell in SARS-CoV-2 infection. In addition, host immune responses would be triggered during infection.

The immune responses include innate immune response and adaptive immune response, both of them playing an indispensable role in defending against pathogens. Owing to only a small number of monocytes-macrophages in the lung express ACE2, how the SARS-CoV-2 infects innate immune cells still unclear. A latest retrospective study reported that T lymphocyte subsets (CD3 + , CD4 + , and CD8 + T-cell) were remarkably decreased and inflammatory cytokines (SAA, CRP, IL-6, and IL-10) were significantly increased in SARS-CoV-2 infected patients, which was positively related to the deaths of in-hospital and severity of disease (Xu, B. et al., 2020a) . Moreover, a survey from lung immune microenvironment by scRNA-seq showed the presence of more CD8 + T-cells effectors and highly expanded clones CD8 + T-cells formed in the mildly COVID-19 patients, suggesting that a rapid and powerful adaptive immune response is potentially crucial for controlling COVID-19 pandemic .

To determine whether XG-1 acts on AT2 cells expressing ACE2, the scRNA-seq data of human lung was studied. It's found that AT2 cells were the highest expression type and were the ACE2-expressing cell type in the profiling. By labeling the AT2 cell genes expressed by ACE2, it found that 153 genes were up-regulated, and then co-analyzed with the target genes of XG-1 to obtain 12 intersection target genes, indicating that it may be the potential target for XG-1 to prevent COVID-19; Among them, 3 intersection target genes (PLA2G1B, SFTPD, SLCO4C1) were related to immunity, indicating that XG-1 may exert preventive effect by immune regulation.

In this study, HPLC-Q-TOF-MS was used to control quality and identify the main chemical constituents of XG-1 granules. Five chemical components in XG-1 were detected. Rutin and Chlorogenic acid (CGA) from Lonicera japonica Thunb. (Jin Yin Hua), Phillyrin (KD-1) from Forsythia suspensa (Thunb.) Vahl (Lian Qiao), Liquiritin from Glycyrrhiza uralensis Fisch. ex DC. (Gan Cao) and Astragaloside IV (AS-IV) from Astragalus mongholicus Bunge (Huang Qi) . Jin Yin Hua is commonly used in Chinese medicine herbs. Numerous compounds such as flavonoids, organic acids, iridoids, triterpenoids and volatile components have been reported from this plant (Song et al., 2014) . Rutin is a flavonoids compound and has been proved antioxidant and anti-inflammatory effects through promoting production of antioxidant enzyme and reducing the expression of proinflammatory cytokine such as NF-κB . CGA is an organic acids compound and as a quality control index component of Jin Yin Hua. Studies have shown that CGA has strong biological activities including antibacterial, antiviral, antitumor, antioxidative, anti-inflammatory and metabolic regulatory activities . KD-1 is a high content of lignan compounds and is used as an important indicator for the quality control of Lian Qiao. It's investigated that KD-1 could attenuate the inflammatory reaction and pathological changes caused by influenza A virus and LPS-infected mice (Qu et al., 2016; Zhong et al., 2013) . Furthermore, KD-1 could protect against SARS-CoV-2 and human coronavirus 229E attack via suppressing the activity of NF-кB signaling pathway . A recent study revealed that CGA and KD-1 has an effective ability to prevent the combination of SARS-CoV-2-spike protein and ACE2 through the molecular docking method . AS-IV, a major active ingredient of Huang Qi, possesses anti-inflammatory, immune regulation, antioxidant, anti-apoptotic, anti-tumor and many other bioactive effects (Li et al., 2017a, b) . Experiment showed that AS-IV has a therapeutic action on inflammation through activating the immune function of Tregs inhibited by high-mobility group box 1 (HMGB1) protein in mice . Liquiritin is a representative ingredient of Gan Cao and it has a neuroprotective, antioxidant and anti-apoptosis effects via regulating multiple pathways .

It's reported that SARS-CoV-2 3CL hydrolase was considered to be an effective target of SARS-CoV-2, which provides a basis for screening active ingredients against COVID-19 (Elmezayen et al., 2020) . The binding of ACE2 and the S protein of SARS-CoV-2 in the human body was also proved to lead the virus invading the body and cause disease . In this study, key active compounds of Chlorogenic acid, Phillyrin, Rutin, Liquiritin and Astragaloside Ⅳ were screened by network pharmacology and HPLC-Q-TOF-MS were verified. Chlorogenic acid and Phillyrin has been explored to block the combination of SARS-CoV-2-spike protein and ACE2 through the molecular docking method . Therefore, Rutin, Liquiritin and Astragaloside Ⅳ were elected to perform molecular docking with SARS-CoV-2 3CL hydrolase and ACE2. From the results, it showed all the binding energy of the three active compounds to the SARS-CoV-2 3CL hydrolase (PDB: 6LU7) and ACE2 (PDB: 1R42) was all less than − 7 kJ/mol, indicating these compounds may bind and form a stable conformation with the receptors to resist SARS-CoV-2. The epidemiological investigations found that individuals who has a close contact exposure history was more likely confirmed as COVID-19 patient, with 1 out of 4.4-16.5 close contacts being identified to be infected with SARS-COV-2, so that tracing and quarantine the close contacts timely is essential to the epidemic spread and control Liu, T. et al., 2020b; Jian et al., 2020) . XG-1is a recommended preventive TCM formula and has been widely used in Hunan Province. It's reported that no infection case was found in the population taking the prescription and immune regulatory effect has been proven on adults (Liu, L. et al., 2020a) . In order to clarify the preventive effect of XG-1 against COVID-19, a retrospective study was conducted, in which 47 close contacts with COVID-19 confirmed cases were screened. Among them, a total of 36 participants took XG-1 decoction and 11 participants refused to take. 7 participants were confirmed positive SARS-CoV-2 nucleic acid by RT-PCR and none of them took the XG-1 preventive prescription, which indicates the potential efficacy of this preparation as a COVID-19 preventive agent.

This study has some limitations. First, due to the incomplete understanding of the virus and disease pathogenic mechanism as well as the limitations of molecular docking itself, the results obtained may be biased. Second, we performed a community-based single-center retrospective study to assess the preventive effects of XG-1, however, relevant in vitro or in vivo experiments could not be carried out due to the biosecurity problem of the disease. The specific mechanism of action still needs further exploration. Third, the sample size in this study was relatively small that might influence the accuracy in statistical analysis. Table 6 General characteristics of close contacts and attack rates of COVID-19 in contacts with different characteristics. Thus, larger, multiple-centers and prospective studies are demanded to ascertain the benefit of XG-1 in COVID-19 prevention and control.

The present study showed that XG-1 exerts a potential preventive effect in close contacts against COVID-19. The underlying mechanism may relate to modulate immunity response and eliminate inflammatory process through several target genes co-expressed with ACE2. These findings provide preliminary evidences and methodological reference for the potential preventive mechanism of XG-1 against COVID-19.

LL, and MZX proposed and designed the study. ZYY, SSX and SYZ screened the data and conducted the data analysis. KG, PY and JJY collected the data of retrospective study. YQ and MZX performed the quality control. HYW wrote this manuscript, and LL revised it. All of the authors reviewed and approved the submitted version of the paper.

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Comparison of the characteristics between the groups with or without XG-1.

With XG-1 (n = 36) Fever 0 (0%) 0 (0%) Cough 0 (0%) 0 (0%) Dyspnea 0 (0%) 0 (0%) Headache 0 (0%) 1 (9.09%) Nausea 0 (0%) 0 (0%) Muscle soreness 0 (0%) 0 (0%)

Mild insomnia 1 (2.78%) 2 (18.18%)

Fei yan No. 1" as a combined treatment for COVID-19: an efficacy and potential mechanistic study

Epidemiology and transmission of COVID-19 in 391 cases and 1286 of their close contacts in Shenzhen, China: a retrospective cohort study

Effects of Rutin on wound healing in hyperglycemic rats

Drug repurposing for coronavirus (COVID-19): in silico screening of known drugs against coronavirus 3CL hydrolase and protease enzymes

General Office of the National Health Commission, 2020. Diagnosis and Treatment Plan for COVID-19 (Trial Version 8)

Clinical characteristics of coronavirus isease 2019 in China

Tissue distribution of ACE2 protein, the functional receptor for SARS coronavirus. A first step in understanding SARS pathogenesis

Clinical features of patients infected with 2019 novel coronavirus in

PPAR-gamma in macrophages limits pulmonary inflammation and promotes host recovery following respiratory viral infection

Hunan Traditional Chinese Medicine Letter (2020) No. 23 Notice of the Hunan Provincial Administration of Traditional Chinese Medicine on the promotion of the use of traditional Chinese medicine in government agencies, enterprises, institutions, and schools to prevent new coronavirus infections

Angiotensin-converting enzyme 2 protects from severe acute lung failure

Contact tracing with digital assistance in Taiwan's COVID-19 outbreak response

Angiotensin-converting enzyme 2 in lung diseases

A crucial role of angiotensin converting enzyme 2 (ACE2) in SARS coronavirus-induced lung injury

Using herbal medicine as a means of prevention experience during the SARS crisis

The efficacy of Chinese medicine for SARS: a review of Chinese publications after the crisis

Structure of SARS coronavirus spike receptor-binding domain complexed with receptor

Astragaloside IV attenuates inflammatory reaction via activating immune function of regulatory T-cells inhibited by HMGB1 in mice

Research review on the pharmacological effects of astragaloside IV

Angiotensin-converting enzyme 2 is a functional receptor for the SARS coronavirus

Liquiritin protects PC12 cells from corticosterone-induced neurotoxicity via regulation of metabolic disorders, attenuation ERK1/2-NF-kappaB pathway, activation Nrf2-Keap1 pathway, and inhibition mitochondrial apoptosis pathway

Yupingfeng Powder relieves the immune suppression induced by dexamethasone in mice

Single-cell landscape of bronchoalveolar immune cells in patients with COVID-19

Insights into human astrocyte response to H5N1 infection by microarray analysis

The effect of hunan preventing corona virus disease 2019 (COVID-19) No.1 prescription on human cellular immunity and humoral immunity

Yupingfeng san inhibits NLRP3 inflammasome to attenuate the inflammatory response in asthma mice

Risk factors associated with COVID-19 infection: a retrospective cohort study based on contacts tracing

Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding

Phillyrin (KD-1) exerts anti-viral and anti-inflammatory activities against novel coronavirus (SARS-CoV-2) and human coronavirus 229E (HCoV-229E) by suppressing the nuclear factor kappa B (NF-kappaB) signaling pathway

Combating COVID-19 with integrated traditional Chinese and Western medicine in China

HIV-1 envelope glycoproteins induce the production of TNF-alpha and IL-10 in human monocytes by activating calcium pathway

Protective effects of phillyrin against influenza A virus in vivo

Single-cell transcriptomics reveals bimodality in expression and splicing in immune cells

Cytoscape: a software environment for integrated models of biomolecular interaction networks

Traditional Chinese medicine formulation therapy in the treatment of coronavirus disease 2019 (COVID-19)

Research progress on chemical constituents from Lonicerae Flos

Single-cell transcriptomics to explore the immune system in health and disease

AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading

Receptor recognition by the novel coronavirus from wuhan: an analysis based on decade-long structural studies of SARS coronavirus

Advances in research on the pharmacological effects and mechanism of action of chlorogenic acid

Yupingfeng pulvis regulates the balance of T cell subsets in asthma mice

Secretory activity and cell cycle alteration of alveolar type II cells in the early and late phase after irradiation

Emergency Committee Regarding the Outbreak of Novel Coronavirus (2019-nCoV

World Health Organization, 2020b. Coronavirus Disease (COVID-19) Outbreak Situation

Network pharmacology based analysis on the molecular biological mechanisms of Xin Hui Tong Formula in coronary heart disease treatment

Traditional Chinese Medicine as a complementary therapy in combat with COVID-19-A review of evidence-based research and clinical practice

Suppressed T cell-mediated immunity in patients with COVID-19: a clinical retrospective study in Wuhan

ETCM: an encyclopaedia of traditional Chinese medicine

Evolution of the novel coronavirus from the ongoing Wuhan outbreak and modeling of its spike protein for risk of human transmission

Comparison of clinical characteristics of patients with asymptomatic vs symptomatic coronavirus disease 2019 in wuhan, China

clusterProfiler: an R package for comparing biological themes among gene clusters

Exploring the active compounds of traditional Mongolian medicine in intervention of novel coronavirus (COVID-19) based on molecular docking method

Anticomplementary principles of a Chinese multiherb remedy for the treatment and prevention of SARS

Single-cell RNA expression profiling of ACE2, the receptor of SARS-CoV-2

The experimental study of Cortex Eucommiae on meridian tropsim: the distribution study of aucubin in rat tissues

Phillyrin attenuates LPS-induced pulmonary inflammation via suppression of MAPK and NF-kappaB activation in acute lung injury mice

Restoration of immune suppressor function of regulatory B cells collected from patients with allergic rhinitis with Chinese medical formula Yupingfeng San

A pneumonia outbreak associated with a new coronavirus of probable bat origin