Skip to main content

ORIGINAL RESEARCH article

Front. Physiol., 22 January 2021
Sec. Systems Biology Archive

CSF3 Is a Potential Drug Target for the Treatment of COVID-19

\r\nChao Fang,&#x;Chao Fang1,2†Jie Mei&#x;Jie Mei1†Huixiang TianHuixiang Tian3Yu-Ligh LiouYu-Ligh Liou1Dingchao RongDingchao Rong4Wei Zhang*Wei Zhang1*Qianjin Liao*Qianjin Liao2*Nayiyuan Wu,*Nayiyuan Wu1,2*
  • 1Hunan Key Laboratory of Pharmacogenetics, Department of Clinical Pharmacology, Xiangya Hospital, Institute of Clinical Pharmacology, Central South University, Changsha, China
  • 2Hunan Cancer Hospital, The Affiliated Cancer Hospital of Xiangya School of Medicine, Central South University, Changsha, China
  • 3Department of Pharmacy, Xiangya Hospital, Central South University, Changsha, China
  • 4Department of Spine Surgery, Xiangya Hospital, Central South University, Changsha, China

Coronavirus Disease 2019 (COVID-19) is an acute respiratory infectious disease that appeared at the end of 2019. As of July 2020, the cumulative number of infections and deaths have exceeded 15 million and 630,000, respectively. And new cases are increasing. There are still many difficulties surrounding research on the mechanism and development of therapeutic vaccines. It is urgent to explore the pathogenic mechanism of viruses to help prevent and treat COVID-19. In our study, we downloaded two datasets related to COVID-19 (GSE150819 and GSE147507). By analyzing the high-throughput expression matrix of uninfected human bronchial organoids and infected human bronchial organoids in the GSE150819, 456 differentially expressed genes (DEGs) were identified, which were mainly enriched in the cytokine–cytokine receptor interaction pathway and so on. We also constructed the protein–protein interaction (PPI) network of DEGs to identify the hub genes. Then we analyzed GSE147507, which contained lung adenocarcinoma cell lines (A549 and Calu3) and the primary bronchial epithelial cell line (NHBE), obtaining 799, 460, and 46 DEGs, respectively. The results showed that in human bronchial organoids, A549, Calu3, and NHBE samples infected with SARS-CoV-2, only one upregulated gene CSF3 was identified. Interestingly, CSF3 is one of the hub genes we previously screened in GSE150819, suggesting that CSF3 may be a potential drug target. Further, we screened potential drugs targeting CSF3 by MOE; the top 50 drugs were screened by flexible docking and rigid docking, with 37 intersections. Two antiviral drugs (Elbasvir and Ritonavir) were included; Elbasvir and Ritonavir formed van der Waals (VDW) interactions with surrounding residues to bind with CSF3, and Elbasvir and Ritonavir significantly inhibited CSF3 protein expression.

Introduction

Coronavirus Disease 2019, first detected in Wuhan, China, in December 2019, has been raging around the world from its detection through to current times (Kannan et al., 2020; Thompson, 2020). On January 9, 2020, the World Health Organization (WHO) announced the discovery of the new coronavirus. This coronavirus was originally called 2019-nCoV and then officially named SARS-CoV-2, which had never been found among humans. On February 11, 2020, the respiratory disease originating from SARS-CoV-2 infection was called COVID-19 (Lu et al., 2020; Mahase, 2020). Although the virus was first discovered in Wuhan, there is still controversy about its specific origin (Andersen et al., 2020; Forster et al., 2020). WHO announced that COVID-19 had become a global health problem because it has spread in more than 140 countries around the world (Zhai et al., 2020).

Fever, sore throat, fatigue, cough, or dyspnea are the typical symptoms of COVID-19 (Kim et al., 2020; Ozma et al., 2020). As of 10 August 2020, confirmed COVID-19 cases have reached 19,718,030, including 728,013 deaths in the world1. The government shared the genome sequence of SARS-CoV-2 with the public in a timely manner. This was effective in helping scientists deal with the emergency (Harcourt et al., 2020; Wu et al., 2020). Although most countries in the world have effectively controlled COVID-19, it is still far from being completely over. This is because research on the treatment and prevention of COVID-19 is still insufficient, and the vaccine under research has not yet been fully applied in clinic (Dhama et al., 2020).

At present, the treatment of COVID-19 is still limited, and it has been reported that some drugs may have a therapeutic effect on COVID-19. As an antimalarial agent, chloroquine and its derivative, hydroxychloroquine, have shown to be active against COVID-19 in vitro (Wang M. et al., 2020). However, there is not enough clinical data to support the therapeutic effect of chloroquine or hydroxychloroquine on COVID-19 (Chen et al., 2020; Huang et al., 2020).

In addition, remdesivi (Jean et al., 2020), lopinavir, and emetine (Choy et al., 2020) have been reported to inhibit SARS-CoV-2 in vitro. But these drugs are not currently recommended for clinical use, and are only recommended for use in clinical trials (Kotwani and Gandra, 2020). Moreover, research on the screening of anti-SARS-CoV-2 drugs and their mechanism is mainly focused on angiotensin-converting enzyme 2 (ACE2) (Williamson and Kerimi, 2020). Whether there are other drug targets for SARS-CoV-2 remains an unanswered question for researchers.

In our study, we analyzed the high-throughput sequencing data of uninfected human bronchial organoids and SARS-CoV2 infected human bronchial organoids in the GEO database to explore the effects of SARS-CoV2 (Blanco-Melo et al., 2020). Then we verified the high-throughput data of A549, Calu3, and NHBE cell lines, finding that only CSF3 (Tamada et al., 2006) was up-regulated in the four sets of samples. It has been reported that the granulocyte colony-stimulating factor (CSF3) was the most upregulated gene after SARS-CoV2-infected genes (Nunnari et al., 2020), suggesting that CSF3 is a significant target after infection and it may also be a potential target for drug therapy. Therefore, we have screened drugs that may potentially target CSF3 for the treatment of COVID-19 from the FDA-approved drug library.

Materials and Methods

Data Download

The Gene Expression Omnibus (GEO)2 database was used to downloaded datasets about COVID-19. With “COVID-19 OR SARS-CoV-2, Homo sapiens” as the screening criteria, we selected GSE150819 and GSE147507 for analysis. In GSE150819, human bronchial organoids (hBO) from commercially available cryopreserved primary human bronchial epithelial cells (hBEpC) were constructed for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) research. We chose three uninfected samples and three SARS-CoV-2 infected samples to analyze. In GSE147507, we selected primary human lung epithelium cells (NHBE) and lung adenocarcinoma cells (A549, Calu-3), which were mock treated or infected with SARS-CoV-2.

Identification of DEGs

In our study, | log fold change (FC)| > 1.5 and p < 0.05 were used as the cut-off criteria to find DEGs in uninfected and SARS-CoV-2-infected organoids or cell lines. Then, these DEGs were visualized by volcano plots and heatmaps, which were performed by R package in RStudio. VENNY 2.1.03 online software was used to identify overlapping DEGs in the four sets of samples.

GO and KEGG Enrichment Analysis

Gene Ontology (GO) analysis classified the genes in the difference tables according to their functions, which included biological process (BP), cellular component (CC), and molecular function (MF). The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed the key pathways of clustered genes. GO and KEGG enrichment analysis of DEGs in GSE150819 were performed by the Database for Annotation, Visualization, and Integrated Discovery (Huang et al., 2009) (DAVID)4 v6.8. And a p < 0.05 was considered to be statistically significant.

Protein–Protein Interaction (PPI) Network Construction, Module Analysis, and Identification of Hub Genes

In order to facilitate understanding of the interaction between the proteins corresponding to these DEGs, a PPI network was constructed by Search Tool for the Retrieval of Interacting Gene (STRING)5 database (Szklarczyk et al., 2015). In the network, core modules and top 10 hub genes were identified by ClusterONE and CytoHubba in Cytoscape software, suggesting the key genes and modules which might play a significant role in COVID-19.

Gene Set Enrichment Analysis

Gene set enrichment analysis (GSEA) is a software which can analyze the positive and negative regulatory pathway of a hub gene or all DEGs. To better understand the role of CSF3 in COVID-19, we chose KEGG function of CSF3 in GSE10819 to analyze, obtaining the key positive and negative regulatory pathways by NOM p-value ranking (p < 0.05).

Molecular Docking

The X-ray structure of CSF3 was downloaded from RCSB Protein Data Bank (PDB ID: 2D9Q)6. The FDA-approved drug library was obtained from DrugBank7. The 2D structures of the FDA-approved drug molecules were drawn in ChemBioDraw 2014 and converted to 3D in MOE through energy minimization. MOE dock was used for homology modeling and molecular docking simulations of the binding affinity between FDA-approved drugs and the target protein. Prior to docking, the force field of AMBER10:EHT and the implicit solvation model of Reaction Field (R-field) were selected. The protonation state of the protein and the orientation of the hydrogens were optimized by LigX, at a PH of 7 and temperature of 300 K. The docking workflow followed the “induced fit” protocol, in which the side chains of the receptor pocket were allowed to move according to ligand conformations, with a constraint on their positions. The weight used for tethering side chain atoms to their original positions was 10. All docked poses were ranked by London dG scoring first, then a force field refinement was carried out on the top 10 poses followed by a rescoring of GBVI/WSA dG. The conformation with the lowest free energy of binding was selected as the best (probable) binding mode. Molecular graphics were generated by PyMOL.

Western Blotting

Crude cellular proteins were separated by SDS-PAGE and transferred to nitrocellulose membranes. Membranes were blocked in 5% non-fat milk and probed with primary antibodies (1:1,000 dilution) against CSF3 (Abcam, ab181053) overnight at 4°C. Membranes were incubated with horseradish peroxidase (HRP)-conjugated secondary antibodies for 1.5 h and immunolabeling detected using a Bio-Rad imaging system.

Results

Identification of DEGs in Datasets

In GSE150819, we analyzed three samples of uninfected human bronchial organoids and three samples of SARS-CoV-2-infected human bronchial organoids, obtaining 456 DEGs, of which 205 were upregulated DEGs and 251 downregulated DEGs (Table 1). To visualize the expression of the DEGs in GSE150819, a heatmap and volcano plot were constructed, as shown in Figure 1.

TABLE 1
www.frontiersin.org

Table 1. The consensual DEGs were identified in GSE150819.

FIGURE 1
www.frontiersin.org

Figure 1. Visual expression of DEGs in GSE150819. (A) The heatmap of DEGSs. (B) The volcano plot of DEGs. Red and green represent relative upregulation and downregulation of gene expression, respectively.

Enrichment Analysis of DEGs

To identify the significant pathways associated with the DEGs between uninfected human bronchial organoids and infected human bronchial organoids in GSE150819, we performed the GO and KEGG pathway enrichment analysis by DAVID.

The result showed that in GO enrichment analysis (Table 2), genes were enriched in biological processes (BP) of inflammatory response, immune response, and cell chemotaxis. As for cellular component (CC), these genes showed enrichment in extracellular space, extracellular region, and integral component of plasma membrane. The molecular function (MF) results showed cytokine activity, chemokine activity, and structural molecule activity. KEGG pathway enrichment analysis was shown in Table 2. DEGs were mainly enriched in cytokine receptor interaction, TNF signaling pathway, and African trypanosomiasis.

TABLE 2
www.frontiersin.org

Table 2. GO and KEGG enrichment analysis of DEGs.

Construction of PPI Network and Identification of Hub Genes

STRING online software constructed the PPI network of 456 DEGs, composed of 282 nodes and 930 edges after excluding the isolated nodes, as shown in Figure 2A. ClusterONE identified the top three modules in the PPI network, suggesting that these submodules might be highly significant in the process of COVID-19. Module 1 (Figure 2B) and module 2 (Figure 2C), consisting of 22 and 12 DEGs, respectively, were both enriched in KEGG pathway, including the Cytokine–cytokine receptor interaction and TNF signaling pathway. Module 3 (Figure 2D) contained 14 nodes and 45 edges.

FIGURE 2
www.frontiersin.org

Figure 2. PPI network and core modules analysis of DEGs in GSE150819. (A) The PPI network of DEGs was constructed by STRING. (B–D) ClusterONE identified the core modules in the PPI network.

The top 10 hub genes were also calculated by the Degree algorithm in CytoHubba. They were MMP9, CCL20, CSF3, CCL5, CXCL10, IL1B, SAA1, ICAM1, IL10, and CXCL8 (Figure 2D). Interestingly, all of them were upregulated after being infected with SARS-CoV-2.

Identification of Overlapping Genes in Different Datasets and GSEA Analysis

Similarly, we also analyzed the three cell lines samples, which were NHBE, A549, and Calu3, and their SARS-CoV-2-infected cells in GSE147507. After screening, the results suggested that there were 45 DEGs in the NHBE group, 158 DEGs in the A549 group, and 460 DEGs in the Calu3 group. VENNY 2.1 online software was used to find the overlapping DEGs in the four sets of samples, and finally only one overlapping DEG was found, as shown in Figures 3A,B. Only CSF3 showed a common upregulated gene in infected human bronchial organoids and cell lines.

FIGURE 3
www.frontiersin.org

Figure 3. Identification and analysis of overlapping DEGs. VENNY software identified the upregulated DEGs (A) and downregulated DEGs (B) of the four sets of samples in GSE150819 and GSE147507. (C) The overlapping gene CSF3 was analyzed by GSEA. There were four positive regulatory pathways of CSF3, which were proteasome, Parkinson’s disease, oxidative phosphorylation, and Graft vs. host disease. The following were four negative regulatory pathways of CSF3: mismatch repair, DNA replication, homologous recombination, and inositol phosphate metabolism.

For GSEA, we analyzed CSF3 enrichment in GSE150819 to explore its potential molecular functions, obtaining the positive and negative regulatory pathways of CSF3. The results were shown in Figure 3C. The results showed that the top four positive regulatory pathways of CSF3 were proteasome, Parkinson’s disease, oxidative phosphorylation, and Graft vs. host disease. While the top four negative regulatory pathways of CSF3 were mismatch repair, DNA replication, homologous recombination, and inositol phosphate metabolism.

Using CSF3 as the Target to Find Potential Targeted Drugs

Giuseppe Nunnari discovered the potential value of CSF3 in COVID-19 research (Nunnari et al., 2020). We obtained the X-ray structure of CSF3 (Tamada et al., 2006) from the RCSB PDB and aimed to screen out the 2,470 FDA-approved drugs to potentially target CSF3. We investigated the binding mode of FDA-approved drugs with CSF3 by rigid and flexible docking simulation. The scores of the top 50 in flexible docking are shown in Table 3. Among them, ubiquinol (DB11340), an active antioxidant (Ernster and Forsmark-Andrée, 1993; Sies, 1997), got the highest score of -8.0513706 kcal/mol. In addition, two antiviral drugs, Elbasvir (DB11574) (Myers et al., 2015) with scores of -6.5620656 kcal/mol and Ritonavir (DB00503) (Cameron et al., 1998) with scores of -6.259393 kcal/mol, have attracted our attention.

TABLE 3
www.frontiersin.org

Table 3. KEGG pathway analysis of DEGs.

The binding mode of Elbasvir with CSF3 were illustrated in Figure 4A. Elbasvir has formed a suitable steric complementarity with the binding site of CSF3. The oxygen atom of carbonyl in Elbasvir, regarded as hydrogen bond acceptors, form the hydrogen bonds with Lys 16 and Lys 23 in CSF3. The binding mode of Ritonavir with CSF3 were illustrated in Figure 4B. The nitrogen atom in Ritonavir, regarded as a hydrogen bond donor, forms hydrogen bonds with Gln 119 in CSF3. Moreover, Elbasvir and Ritonavir formed van der Waals (VDW) interactions with surrounding residues, which mainly contribute to the binding energy between Elbasvir and Ritonavir with CSF3.

FIGURE 4
www.frontiersin.org

Figure 4. Docked pose of CSF3 with Elbasvir and Ritonavir. (A) Overall structural view of CSF3 with Elbasvir, wherein Elbasvir are chimeric into protein; Detailed interaction view between Elbasvir and CSF3 wherein hydrogen bonding is shown as a red dashed line; 2D interaction view between Elbasvir and CSF3, with green arrows representing hydrogen bonding interactions. (B) Overall structural view of CSF3 with Ritonavir, wherein Ritonavir are chimeric into protein; Detailed interaction view between Ritonavir and CSF3 wherein hydrogen bonding is shown as red dashed line; 2D interaction view between Ritonavir and CSF3, with green arrows representing hydrogen bonding interactions. (C,D) Both Elbasvir and Ritonavir inhibited CSF3 expression in A549 cells when treated with 0.1 and 1 μM of Elbasvir/Ritonavir for 24 h.

We further verified the effect of Elbasvir and Ritonavir on CSF3; WB results show that both Elbasvir and Ritonavir significantly inhibited the CSF3 protein expression in A549 cells (Figures 4C,D).

Discussion

In the current global pandemic, COVID-19 has infected more than 19 million people and caused more than 720,000 deaths. It has a very high infection rate (Li R. et al., 2020), and the case fatality rate has been reported to be as high as 5% (Li L.Q. et al., 2020). Clinical data in the United States showed that among patients hospitalized in the intensive care unit, case fatality was up to 40% (Wiersinga et al., 2020). The virus has also bene shown to spread among people of all ages8. Studies have shown that the risk of death is highly correlated with the age gradient, and younger age may be associated with lower mortality (Verity et al., 2020).

This is an epidemic that has affected all mankind and requires our joint efforts. In the winter when the virus was raging (Smit et al., 2020), China quickly quarantined Wuhan, promptly and effectively suppressing the spread of COVID-19. The Chinese people came together with the government to combat the disease. As of March 20209, the daily number of local cases in mainland China has dropped to two digits, and the infected are mainly from people returning from abroad.

But the virus still exists, and we still have not completely defeated it. The existing obstacles are the limitations of treatment methods, the lack of effective drugs (Triggle et al., 2020), the difficulty of vaccine development (Uddin et al., 2020), and the inability to completely block the spread of the virus (Rothan and Byrareddy, 2020). All of these issues require continuous efforts by researchers, and research into treatment and prevention should be given equal attention. On the one hand, reliable treatments and therapeutic drugs need to be found. On the other hand, the spread of the virus needs to be blocked and vaccines need to be researched as quickly as possible.

In our study, we hope to find potential drug therapeutic targets from the high-throughput data of virus-infected organoids and cells and explore whether there are clinical drugs in use that can target the infectious protein. We found 456 DEGs in the analysis of uninfected and SARS-CoV-2-infected human bronchial organoids. KEGG enrichment analysis showed that these genes were mainly enriched in Cytokine–cytokine receptor interaction and TNF signaling pathway. Through the verification of high-throughput data of A549, Calu3, and NHBE cell lines, we have found an important target protein CSF3, which, along with its receptor CSF3R, participate in the regulation of granulopoiesis, neutrophil function, and hematopoietic stem cell mobilization (Zhang et al., 2018). It has been reported that CSF3 mutations have carcinogenic effects (Hollmén et al., 2016; Jing et al., 2016). Virus-related research has found that human respiratory syncytial virus (hRSV) infection disrupted the polarity of the pediatric respiratory epithelial secretome and was associated with immune modulating proteins (CXCL6, CXCL16, and CSF3) never before linked with this virus before (Touzelet et al., 2020). The latest reports showed that CSF3 is the most modulated gene in NHBE cells infected with SARS-CoV-2 (Nunnari et al., 2020), indicating the potential value of CSF3.

Regarding CSF3 as a drug target, we use molecular docking simulation technology by MOE to screen out potential drugs from the FDA-approved drug library, which were downloaded from DrugBank. These drugs were scored separately by rigid docking and flexible docking. Among the top 50 drugs, a total of 37 drugs were screened by both methods. The screened drugs included anti-HCV drugs, such as Ombitasvir (DB09296) (Ascione et al., 2018), Daclatasvir (DB09102) (Lee, 2013), and Elbasvir (DB11574) (Myers et al., 2015), and the anti-HIV drugs, Cobicistat (DB09065) (Angione et al., 2018) and Ritonavir (DB00503) (Cameron et al., 1998). The drug with the highest score is Ubiquinol (DB11340), which is an active antioxidant drug (Ernster and Forsmark-Andrée, 1993). No antiviral drugs have been proven to be effective against COVID-19. Therefore, we give priority to antiviral drugs such as Elbasvir to explore whether it has a potential effect to COVID-19 in our study. Similarly, Junmei Wang screened out Elbasvir, carfilzomib, eravacycline, valrubicin, and lopinavir as potential inhibitors of SARS-CoV-2’s main protease (Wang, 2020). In Meenakshisundaram Balasubramaniam’s study, Elbasvir was predicted to bind multiple SARS-CoV-2 proteins, blocking virus replication, and had a good antiviral effect (Balasubramaniam and Reis, 2020). As an FDA-approved drug for anti-HIV and HCV, Ritonavir (Feld et al., 2014) was tried early in the treatment of COVID-19 (Lim et al., 2020; Wang Y. et al., 2020), although the results were not satisfactory (Griffin, 2020). At present, these predicted antiviral drugs are still not recommended for clinical use, and are only supported for clinical trials. In a word, further mechanism exploration and clinical trials are needed for the application of these antiviral drugs.

Conclusion

In conclusion, our study identified the DEGs and underlying pathways in uninfected and SARS-CoV-2-infected human bronchial organoids and cell lines, found the only upregulated gene was CSF3, and explored its potential value. We regarded CSF3 as a potential drug target, and used molecular docking and virtual drug screening technology to find drugs that may interact with it to find potential treatments for COVID-19. Elbasvir, Ritonavir, and other antiviral drugs were screened through both rigid docking and flexible docking, and they were potential anti-COVID-19 drugs that have been studied. Targeting CSF3 may be a potential therapeutic mechanism of these drugs, but a more in-depth exploration is needed.

Data Availability Statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.

Author Contributions

CF, NW, WZ, and Y-LL: conception and design. NW and QL: administrative support. JM, CF, and HT: collection and assembly of data. CF, DR, and JM: data analysis and interpretation. NW, JM, and CF: manuscript writing. All authors finally approved the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (82003050 and 81972836), the Natural Science Foundation of Hunan Province (2020JJ5338 and 2020JJ8020), the Foundation from Social Development Science and Technology Division (2018SK2121, kq2004138, and kq1801104), the Scientific Research Project of Hunan Health Commission (20201487), and the Hunan Provincial Health Commission (20201487).

Conflict of Interest

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.

Acknowledgments

We would like to thank all the members of the WZ Laboratory.

Footnotes

  1. ^ https://www.who.int/emergencies/diseases/novel-coronavirus-2019
  2. ^ https://www.ncbi.nlm.nih.gov/geo/
  3. ^ https://bioinfogp.cnb.csic.es/tools/venny/
  4. ^ https://david.ncifcrf.gov/
  5. ^ https://string-db.org/cgi/input.pl
  6. ^ http://www.rcsb.org/
  7. ^ https://www.drugbank.ca/
  8. ^ https://www.who.int/publications-detail/report-of-the-who-china-joint-mission-on-coronavirus-disease-2019-(covid-19)
  9. ^ https://www.who.int/docs/default-source/coronaviruse/situation-reports/20200331-sitrep-71-covid-19.pdf?sfvrsn=4360e92b_8

References

Andersen, K. G., Rambaut, A., Lipkin, W. I., Holmes, E. C., and Garry, R. F. (2020). The proximal origin of SARS-CoV-2. Nat. Med. 26, 450–452. doi: 10.1038/s41591-020-0820-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Angione, S. A., Cherian, S. M., and özdener, A. E. (2018). A review of the efficacy and safety of genvoya® (Elvitegravir, Cobicistat, Emtricitabine, and Tenofovir Alafenamide) in the management of HIV-1 infection. J. Pharm. Pract. 31, 216–221. doi: 10.1177/0897190017710519

PubMed Abstract | CrossRef Full Text | Google Scholar

Ascione, A., De Luca, M., Melazzini, M., Montilla, S., Trotta, M. P., Petta, S., et al. (2018). Safety and efficacy of ombitasvir/paritaprevir/ritonavir/dasabuvir plus ribavirin in patients over 65 years with HCV genotype 1 cirrhosis. Infection 46, 607–615. doi: 10.1007/s15010-018-1157-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Balasubramaniam, M., and Reis, R. (2020). Computational target-based drug repurposing of elbasvir, an antiviral drug predicted to bind multiple SARS-CoV-2 proteins. ChemRxiv [preprint]. doi: 10.26434/chemrxiv.12084822

PubMed Abstract | CrossRef Full Text | Google Scholar

Blanco-Melo, D., Nilsson-Payant, B. E., Liu, W. C., Uhl, S., Hoagland, D., Møller, R., et al. (2020). Imbalanced host response to SARS-CoV-2 drives development of COVID-19. Cell 181, 1036–1045. doi: 10.1016/j.cell.2020.04.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Cameron, D. W., Heath-Chiozzi, M., Danner, S., Cohen, C., Kravcik, S., Maurath, C., et al. (1998). Randomised placebo-controlled trial of ritonavir in advanced HIV-1 disease. the advanced HIV disease ritonavir study group. Lancet 351, 543–549. doi: 10.1016/s0140-6736(97)04161-5

CrossRef Full Text | Google Scholar

Chen, J., Liu, D., Liu, L., Liu, P., Xu, Q., Xia, L., et al. (2020). [A pilot study of hydroxychloroquine in treatment of patients with moderate COVID-19]. Zhejiang Da Xue Xue Bao Yi Xue Ban 49, 215–219. doi: 10.3785/j.issn.1008-9292.2020.03.03

PubMed Abstract | CrossRef Full Text | Google Scholar

Choy, K. T., Wong, A. Y., Kaewpreedee, P., Sia, S. F., Chen, D., Hui, K., et al. (2020). Remdesivir, lopinavir, emetine, and homoharringtonine inhibit SARS-CoV-2 replication in vitro. Antiviral. Res. 178:104786. doi: 10.1016/j.antiviral.2020.104786

PubMed Abstract | CrossRef Full Text | Google Scholar

Dhama, K., Sharun, K., Tiwari, R., Dadar, M., Malik, Y. S., Singh, K. P., et al. (2020). COVID-19, an emerging coronavirus infection: advances and prospects in designing and developing vaccines, immunotherapeutics, and therapeutics. Hum. Vaccin. Immunother 16, 1232–1238. doi: 10.1080/21645515.2020.1735227

PubMed Abstract | CrossRef Full Text | Google Scholar

Ernster, L., and Forsmark-Andrée, P. (1993). Ubiquinol: an endogenous antioxidant in aerobic organisms. Clin. Investig. 71, S60–S65. doi: 10.1007/BF00226842

PubMed Abstract | CrossRef Full Text | Google Scholar

Feld, J. J., Kowdley, K. V., Coakley, E., Sigal, S., Nelson, D. R., Crawford, D., et al. (2014). Treatment of HCV with ABT-450/r-ombitasvir and dasabuvir with ribavirin. N. Engl. J. Med. 370, 1594–1603. doi: 10.1056/NEJMoa1315722

PubMed Abstract | CrossRef Full Text | Google Scholar

Forster, P., Forster, L., Renfrew, C., and Forster, M. (2020). Phylogenetic network analysis of SARS-CoV-2 genomes. Proc. Natl. Acad. Sci. U S A. 117, 9241–9243. doi: 10.1073/pnas.2004999117

PubMed Abstract | CrossRef Full Text | Google Scholar

Griffin, S. (2020). Covid-19: Lopinavir-ritonavir does not benefit hospitalised patients. UK trial finds. BMJ 370:m2650. doi: 10.1136/bmj.m2650

PubMed Abstract | CrossRef Full Text | Google Scholar

Harcourt, J., Tamin, A., Lu, X., Kamili, S., Sakthivel, S. K., Murray, J., et al. (2020). Severe acute respiratory syndrome coronavirus 2 from patient with coronavirus disease, United States. Emerg. Infect. Dis. 26, 1266–1273. doi: 10.3201/eid2606.200516

PubMed Abstract | CrossRef Full Text | Google Scholar

Hollmén, M., Karaman, S., Schwager, S., Lisibach, A., Christiansen, A. J., Maksimow, M., et al. (2016). G-CSF regulates macrophage phenotype and associates with poor overall survival in human triple-negative breast cancer. Oncoimmunology 5:e1115177. doi: 10.1080/2162402x.2015.1115177

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, D. W., Sherman, B. T., and Lempicki, R. A. (2009). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4, 44–57. doi: 10.1038/nprot.2008.211

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, M., Tang, T., Pang, P., Li, M., Ma, R., Lu, J., et al. (2020). Treating COVID-19 with Chloroquine. J. Mol. Cell Biol. 12, 322–325. doi: 10.1093/jmcb/mjaa014

PubMed Abstract | CrossRef Full Text | Google Scholar

Jean, S. S., Lee, P. I., and Hsueh, P. R. (2020). Treatment options for COVID-19: the reality and challenges. J. Microbiol. Immunol. Infect. 53, 436–443. doi: 10.1016/j.jmii.2020.03.034

PubMed Abstract | CrossRef Full Text | Google Scholar

Jing, X., Sonoki, T., Miyajima, M., Sawada, T., Terada, N., Takemura, S., et al. (2016). EphA4-deleted microenvironment regulates cancer development and leukemoid reaction of the isografted 4T1 murine breast cancer via reduction of an IGF1 signal. Cancer Med. 5, 1214–1227. doi: 10.1002/cam4.670

PubMed Abstract | CrossRef Full Text | Google Scholar

Kannan, S., Shaik, S. A. P., Sheeza, A., and Hemalatha, K. (2020). COVID-19 (Novel Coronavirus 2019) - recent trends. Eur. Rev. Med. Pharmacol. Sci. 24, 2006–2011. doi: 10.26355/eurrev_202002_20378

CrossRef Full Text | Google Scholar

Kim, E. S., Chin, B. S., Kang, C. K., Kim, N. J., Kang, Y. M., Choi, J. P., et al. (2020). clinical course and outcomes of patients with severe acute respiratory syndrome coronavirus 2 infection: a preliminary report of the first 28 patients from the korean cohort study on COVID-19. J. Korean Med. Sci. 35:e142. doi: 10.3346/jkms.2020.35.e142

PubMed Abstract | CrossRef Full Text | Google Scholar

Kotwani, A., and Gandra, S. (2020). Potential pharmacological agents for COVID-19. Ind. J. Public Health 64, S112–S116. doi: 10.4103/ijph.IJPH_456_20

CrossRef Full Text | Google Scholar

Lee, C. (2013). Daclatasvir: potential role in hepatitis C. Drug Des. Dev. Ther. 7, 1223–1233. doi: 10.2147/DDDT.S40310

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, L. Q., Huang, T., Wang, Y. Q., Wang, Z. P., Liang, Y., Huang, T. B., et al. (2020a). COVID-19 patients’ clinical characteristics, discharge rate, and fatality rate of meta-analysis. J. Med. Virol. 92, 577–583. doi: 10.1002/jmv.25757

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, R., Pei, S., Chen, B., Song, Y., Zhang, T., Yang, W., et al. (2020b). Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV-2). Science 368, 489–493. doi: 10.1126/science.abb3221

PubMed Abstract | CrossRef Full Text | Google Scholar

Lim, J., Jeon, S., Shin, H. Y., Kim, M. J., Seong, Y. M., Lee, W. J., et al. (2020). Case of the index patient who caused tertiary transmission of COVID-19 infection in korea: the application of lopinavir/ritonavir for the treatment of COVID-19 infected pneumonia monitored by quantitative RT-PCR. J. Korean Med. Sci. 35:e79. doi: 10.3346/jkms.2020.35.e79

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, R., Zhao, X., Li, J., Niu, P., Yang, B., Wu, H., et al. (2020). Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding. Lancet 395, 565–574. doi: 10.1016/S0140-6736(20)30251-8

CrossRef Full Text | Google Scholar

Mahase, E. (2020). China coronavirus: WHO declares international emergency as death toll exceeds 200. BMJ 368:m408. doi: 10.1136/bmj.m408

PubMed Abstract | CrossRef Full Text | Google Scholar

Myers, R. P., Shah, H., Burak, K. W., Cooper, C., and Feld, J. J. (2015). An update on the management of chronic hepatitis C: 2015 consensus guidelines from the Canadian association for the study of the liver. Can. J. Gastroenterol. Hepatol. 29, 19–34. doi: 10.1155/2015/692408

PubMed Abstract | CrossRef Full Text | Google Scholar

Nunnari, G., Sanfilippo, C., Castrogiovanni, P., Imbesi, R., Volti, G. L., Barbagallo, I., et al. (2020). Network perturbation analysis in human bronchial epithelial cells following SARS-CoV2 infection. Exp. Cell Res. 395:112204. doi: 10.1016/j.yexcr.2020.112204

PubMed Abstract | CrossRef Full Text | Google Scholar

Ozma, M. A., Maroufi, P., Khodadadi, E., Köse, ş, Esposito, I., Ganbarov, K., et al. (2020). Clinical manifestation, diagnosis, prevention and control of SARS-CoV-2 (COVID-19) during the outbreak period. Infez Med. 28, 153–165.

Google Scholar

Rothan, H. A., and Byrareddy, S. N. (2020). The epidemiology and pathogenesis of coronavirus disease (COVID-19) outbreak. J. Autoimmun. 109:102433. doi: 10.1016/j.jaut.2020.102433

PubMed Abstract | CrossRef Full Text | Google Scholar

Sies, H. (1997). Oxidative stress: oxidants and antioxidants. Exp. Physiol. 82, 291–295. doi: 10.1113/expphysiol.1997.sp004024

PubMed Abstract | CrossRef Full Text | Google Scholar

Smit, A. J., Fitchett, J. M., Engelbrecht, F. A., Scholes, R. J., Dzhivhuho, G., and Sweijd, N. A. (2020). Winter is coming: a southern hemisphere perspective of the environmental drivers of SARS-CoV-2 and the potential seasonality of COVID-19. Int. J. Environ. Res. Public Health 17:5634. doi: 10.3390/ijerph17165634

PubMed Abstract | CrossRef Full Text | Google Scholar

Szklarczyk, D., Franceschini, A., Wyder, S., Forslund, K., Heller, D., Huerta-Cepas, J., et al. (2015). STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 43, D447–D452. doi: 10.1093/nar/gku1003

PubMed Abstract | CrossRef Full Text | Google Scholar

Tamada, T., Honjo, E., Maeda, Y., Okamoto, T., Ishibashi, M., Tokunaga, M., et al. (2006). Homodimeric cross-over structure of the human granulocyte colony-stimulating factor (GCSF) receptor signaling complex. Proc. Natl. Acad. Sci. U S A. 103, 3135–3140. doi: 10.1073/pnas.0511264103

PubMed Abstract | CrossRef Full Text | Google Scholar

Thompson, R. (2020). Pandemic potential of 2019-nCoV. Lancet Infect. Dis. 20:280. doi: 10.1016/S1473-3099(20)30068-2

CrossRef Full Text | Google Scholar

Touzelet, O., Broadbent, L., Armstrong, S. D., Aljabr, W., Cloutman-Green, E., Power, U. F., et al. (2020). the secretome profiling of a pediatric airway epithelium infected with hRSV identified aberrant apical/basolateral trafficking and novel immune modulating (CXCL6, CXCL16, CSF3) and antiviral (CEACAM1) proteins. Mol. Cell. Proteom. 19, 793–807. doi: 10.1074/mcp.ra119.001546

CrossRef Full Text | Google Scholar

Triggle, C. R., Bansal, D., Farag, E., Ding, H., and Sultan, A. A. (2020). COVID-19: learning from lessons to guide treatment and prevention interventions. mSphere 5, e317–e320.

Google Scholar

Uddin, M., Mustafa, F., Rizvi, T. A., Loney, T., Suwaidi, H. A., Al-Marzouqi, A., et al. (2020). SARS-CoV-2/COVID-19: viral genomics, epidemiology, vaccines, and therapeutic interventions. Viruses 12:526. doi: 10.3390/v12050526

PubMed Abstract | CrossRef Full Text | Google Scholar

Verity, R., Okell, L. C., Dorigatti, I., Winskill, P., Whittaker, C., Imai, N., et al. (2020). Estimates of the severity of coronavirus disease 2019: a model-based analysis. Lancet Infect. Dis. 20, 669–677. doi: 10.1016/S1473-3099(20)30243-7

CrossRef Full Text | Google Scholar

Wang, J. (2020). Fast identification of possible drug treatment of coronavirus disease-19 (COVID-19) through computational drug repurposing study. J. Chem. Inf. Model. 60, 3277–3286. doi: 10.1021/acs.jcim.0c00179

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, M., Cao, R., Zhang, L., Yang, X., Liu, J., Xu, M., et al. (2020a). Remdesivir and chloroquine effectively inhibit the recently emerged novel coronavirus (2019-nCoV) in vitro. Cell Res. 30, 269–271. doi: 10.1038/s41422-020-0282-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Zhang, D., Du, G., Du, R., Zhao, J., Jin, Y., et al. (2020b). Remdesivir in adults with severe COVID-19: a randomised, double-blind, placebo-controlled, multicentre trial. Lancet 395, 1569–1578. doi: 10.1016/S0140-6736(20)31022-9

CrossRef Full Text | Google Scholar

Wiersinga, W. J., Rhodes, A., Cheng, A. C., Peacock, S. J., and Prescott, H. C. (2020). Pathophysiology, transmission, diagnosis, and treatment of coronavirus disease 2019 (COVID-19): a review. JAMA 324, 782–793. doi: 10.1001/jama.2020.12839

PubMed Abstract | CrossRef Full Text | Google Scholar

Williamson, G., and Kerimi, A. (2020). Testing of natural products in clinical trials targeting the SARS-CoV-2 (Covid-19) viral spike protein-angiotensin converting enzyme-2 (ACE2) interaction. Biochem. Pharmacol. 178:114123. doi: 10.1016/j.bcp.2020.114123

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, F., Zhao, S., Yu, B., Chen, Y. M., Wang, W., Song, Z. G., et al. (2020). A new coronavirus associated with human respiratory disease in China. Nature 579, 265–269. doi: 10.1038/s41586-020-2008-3

CrossRef Full Text | Google Scholar

Zhai, P., Ding, Y., Wu, X., Long, J., Zhong, Y., and Li, Y. (2020). The epidemiology, diagnosis and treatment of COVID-19. Int. J. Antimicrob Agents 55:105955. doi: 10.1016/j.ijantimicag.2020.105955

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H., Coblentz, C., Watanabe-Smith, K., Means, S., Means, J., Maxson, J. E., et al. (2018). Gain-of-function mutations in granulocyte colony-stimulating factor receptor (CSF3R) reveal distinct mechanisms of CSF3R activation. J. Biol. Chem. 293, 7387–7396. doi: 10.1074/jbc.ra118.002417

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: COVID-19, MOE, ritonavir, elbasvir, CSF3

Citation: Fang C, Mei J, Tian H, Liou Y-L, Rong D, Zhang W, Liao Q and Wu N (2021) CSF3 Is a Potential Drug Target for the Treatment of COVID-19. Front. Physiol. 11:605792. doi: 10.3389/fphys.2020.605792

Received: 22 September 2020; Accepted: 22 December 2020;
Published: 22 January 2021.

Edited by:

James J. Cai, Texas A&M University, United States

Reviewed by:

Joyce Chen, University of Chicago, United States
Ali Salehzadeh-Yazdi, University of Rostock, Germany

Copyright © 2021 Fang, Mei, Tian, Liou, Rong, Zhang, Liao and Wu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Nayiyuan Wu, wunayiyuan@163.com; Qianjin Liao, march-on@126.com; Wei Zhang, csuzhangwei@csu.edu.cn

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.