key: cord-0841179-gy35gsl2 authors: Zhang, Deng-hai; Wu, Kun-lun; Zhang, Xue; Deng, Sheng-qiong; Peng, Bin title: In silico screening of Chinese herbal medicines with the potential to directly inhibit 2019 novel coronavirus date: 2020-02-20 journal: J Integr Med DOI: 10.1016/j.joim.2020.02.005 sha: d3fb8500862db1ce9cc299443a6d31209fab5602 doc_id: 841179 cord_uid: gy35gsl2 OBJECTIVE: In this study we execute a rational screen to identify Chinese medical herbs that are commonly used in treating viral respiratory infections and also contain compounds that might directly inhibit 2019 novel coronavirus (2019-nCoV), an ongoing novel coronavirus that causes pneumonia. METHODS: There were two main steps in the screening process. In the first step we conducted a literature search for natural compounds that had been biologically confirmed as against sever acute respiratory syndrome coronavirus or Middle East respiratory syndrome coronavirus. Resulting compounds were cross-checked for listing in the Traditional Chinese Medicine Systems Pharmacology Database. Compounds meeting both requirements were subjected to absorption, distribution, metabolism and excretion (ADME) evaluation to verify that oral administration would be effective. Next, a docking analysis was used to test whether the compound had the potential for direct 2019-nCoV protein interaction. In the second step we searched Chinese herbal databases to identify plants containing the selected compounds. Plants containing 2 or more of the compounds identified in our screen were then checked against the catalogue for classic herbal usage. Finally, network pharmacology analysis was used to predict the general in vivo effects of each selected herb. RESULTS: Of the natural compounds screened, 13 that exist in traditional Chinese medicines were also found to have potential anti-2019-nCoV activity. Further, 125 Chinese herbs were found to contain 2 or more of these 13 compounds. Of these 125 herbs, 26 are classically catalogued as treating viral respiratory infections. Network pharmacology analysis predicted that the general in vivo roles of these 26 herbal plants were related to regulating viral infection, immune/inflammation reactions and hypoxia response. CONCLUSION: Chinese herbal treatments classically used for treating viral respiratory infection might contain direct anti-2019-nCoV compounds. Toward the end of December 2019, a novel coronavirus (2019-nCoV) with human-to-human transmission and severe human infection, originating in Wuhan, China, was identified [1] . This virus has affected many persons in China and spread to other countries in a very short time. On January 30, 2020, the Director-General of the World Health Organization declared that the outbreak of 2019-nCoV constitutes a public health emergency of international concern and issued temporary recommendations under the International Health Regulations [2] . According to the Daily Report of China National Health Commission, as of this writing on February 2, 2020, 14,488 cases, including 304 deaths, have been confirmed in China; 146 cases, including 1 death, have also been reported among 23 other countries. This pandemic is still ongoing, so it is urgent to find new preventive and therapeutic agents as soon as possible. In addition, commensurate with the risk, strong measures for early detection, isolation and treatment of cases, as well as minimization of transmission through social interaction must be implemented. While specific vaccines and antiviral agents are the most effective methods to prevent and treat viral infection, there are not yet effective treatments that target the 2019-nCoV. Development of these treatments may require months or years, meaning that a more immediate treatment or control mechanism should be found if possible. Herbs used in traditional Chinese medicine present a potentially valuable resource to this end. The effectiveness of herbal treatment to control contagious disease was demonstrated during the 2003 severe acute respiratory syndrome (SARS) outbreak [3] . As such, the Chinese government is encouraging the use of herbal plants in fighting this new viral pneumonia. However, the application of herbal treatment is mainly guided by the type of herb (based on the catalogue of classic literature on herbs) and the patient's symptoms or signs. There is often not enough information to predetermine whether the herbs in question can directly target the viral cause, in other words, herbal usage is generally not guided by viral pathology. We think more detailed knowledge about the direct antiviral effects of different plants would be greatly helpful to the doctors selecting them. In fact, after the outbreak of SARS, many groups dedicated themselves to finding anti-coronavirus agents, including some natural compounds that exist in traditional Chinese herbal medicines [4] [5] [6] [7] [8] [9] [10] [11] [12] . The coronavirus encodes more than one dozen proteins, some of which are essential to viral entry and replication. Among these proteins, the most well-studied are papain-like protease (PLpro), 3C-like protease (3CLpro) and spike protein. Coronavirus PLpro not only processes the viral polypeptide onto functional proteins but is also a deubiquitinating enzyme that can dampen host anti-viral response by hijacking the ubiquitin (Ub) system. For example, SARS PLpro cleaves ISG15, a two-domain Ub-like protein, and Lys48-linked polyUb chains, releasing diUbLys48 products [13, 14] . SARS-3CLpro is a cysteine protease indispensable to the viral life cycle [15] . Coronavirus spike protein uses angiotensinconverting enzyme 2 as a receptor to help the virus enter cells [16] . These three proteins make attractive targets for drug development. Through in silico and biological processing, a series of small molecules, including those from natural compounds, have been screened and confirmed to directly inhibit these important proteins in SARS or Middle East respiratory syndrome (MERS) coronavirus [17] [18] [19] [20] [21] [22] [23] . The gene sequence of 2019-nCoV has been released, which suggests high similarities between the main proteins in this virus and those previously identified in SARS-Cov or MERS-Cov [24, 25] . In this sense, previously reported anti-SARS-Cov or anti-MERS-Cov natural compounds may become a valuable guide to finding anti-coronavirus (2019-nCoV) herbal plants among the traditional Chinese herbs used to treat viral pneumonia. It is a challenge to screen out the herbs containing anticoronavirus (2019-nCoV) compounds from the large number of those possibly being used for patients infected with this pathogen, especially in very short time. Here, we propose two principles to guide such work: oral effectiveness and traditional usage compatibility. The first principle refers to the fact that most Chinese herbal plants are orally ingested after boiling with water, meaning that the anti-coronavirus (2019-nCoV) ingredients in selected plants should be absorbable via oral preparation. The second principle recognizes that candidate plants should be consistent with the type classifications for traditional herbal usage, since type-guided applications are integral to herbal use, as mentioned above. Following these two principles, we used a 6-step selection process (3 for each principle), including drug-likeness, evaluation of oral bioavailability, molecular docking, network pharmacology analysis and other methods to identify herbs that have both a high possibil-ity of containing effective anti-coronavirus (2019-nCoV) compounds and are classified as treating virus-caused respiratory infection. PubMed literature concerning natural compounds against SARS or MERS coronavirus activity was selected using the query ''coronavirus AND inhibitor AND (SARS OR MERS OR SARS-CoV OR MERS-CoV)." After careful reading of the studies returned by this search, the natural compounds that had biologically confirmed antiviral activities were compared with the Traditional Chinese Medicine Systems Pharmacology database (TCMSP, http://www.tcmspw.com/browse.php?qc=herbs), the Encyclopedia of Traditional Chinese Medicine (ETCM, http://www.nrc.ac.cn:9090/ETCM/) and SymMap (https://www.symmap.org/). Natural compounds both associated with antiviral activity and contained in herbs were examined in the next step of our study. Since Chinese herbal treatments are always taken orally after boiling with water, an in silico integrative model of absorption, distribution, metabolism and excretion (ADME) was used to screen for natural compounds that may be bioactive via oral administration. The indices used for the screening include evaluation of oral bioavailability, Caco-2 permeability, drug-like value, and drug half-life. The threshold values indicating effectiveness for these four indices were > 30%, > À0.4, > 0.18 and > 3 h, respectively, as recommended by Hu et al [26] . The values of these four indices can be obtained from the TCMSP database. We used molecular docking software AutoDock 4 to perform protein compound docking analysis, according to the following procedure: (1) We built three-dimensional (3D) structure files of the proteins of interest. We used the online server SWISS-MODEL (https://swissmodel.expasy.org/) to build the 3D structures of the proteins of interest by template-based modeling, these template structures being the reported 3D structures of the corresponding proteins from SARS-CoV. The models built were of Protein Data Bank (PDB) format. (2) To retrieve the required 3D structure files of compounds, the structure data file (SDF) format of compounds were retrieved from the PubChem website and then converted to PDB format by Discovery Studio. (3) AutoDock 4.2 was used to prepare PDBQT format files for target and ligand screening (Target. pdbqt and Ligand.pdbqt) and grid and docking parameter files (a. gpf and a.dpf). (4) Molecular docking was performed using Auto-Dock in Cygwin and finally the results were analyzed. The process and parameters used were detailed by Rizvi et al [27] . Herbs were selected through three steps. (1) Primary selection: molecules chosen from the above steps were used as input for the TCMSP, ETCM and SymMap to search for plants containing that input and the plants were filtered by the numbers of antiviral compounds they contain. Those containing 2 or more antiviral compounds were selected for the next step. (2) Classic usage catalogue cross-reference: only herbs traditionally used to treat viral respiratory infection were retained for further study. (3) Predication of general effects in vivo with network pharmacology analysis, which is detailed as follows. The TCMSP provided the main components of each herb and the protein targets for each component. We identified the reported chemical constituents for each plant in the final analysis and used the ADME indices listed above to find the orally absorbable and drug-like compounds for the plant. The protein targets of these compounds were downloaded from the TCMSP database. All protein targets for each individual plant were used as input for the String online server (https://string-db.org/) to perform proteinprotein interaction analysis and pathway enrichment. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enriched (with P < 0.01) by the input were downloaded. All data were processed using the statistical language R (3.6.2), unless otherwise specified. We received 261 hits from conducting our search in the PubMed database. After careful evaluation of the abstracts from these citations we downloaded and carefully analyzed the full text of 23 highly relevant papers. The natural compounds reported to have biologically confirmed anti-coronavirus activity were identified and then compared to the ingredients listed in three Chinese herbal databases. The result was 115 overlapping ingredients, which we used for further testing (Fig. 1 ). The antiviral activities of the 115 natural compounds were reportedly confirmed with enzyme-based (cell-free) or cell-based experimental systems. To be utilized as a Chinese herbal medicine, they must be absorbable via oral prescription. Therefore, we performed ADME screening for the 115 natural compounds, reducing the number of candidates to 13. To perform the docking analysis, the 3D structure files of 2019-nCoV PLpro, 3CLpro and spike proteins were built based on the corresponding SARS-CoV templates, i.e., PDB 5e6j, 1uj1 and 6cad, respectively. Then, molecule-protein docking was carried out between the molecules and their reported targets. If the molecules were reported to inhibit viral entry, they were docked with spike proteins (Table 1) . Each separate analysis returned positive results ( Table 1 , Fig. 2 and online supplementary Fig. S1 ), indicating the natural compounds we selected might directly inhibit 2019-nCoV. The molecules selected to target PLpro (M2, M3, M7, M9, M10, M11 and M13) mainly bound in the region between the thumb and palm domains, which might interfere with substrate entering this enzyme's active sites, located at the bottom of the two domains [29] . The molecules reported to inhibit 3CLpro (M1, M2, M3, M4, M5, M7, M8, M10, M11, M12 and M13) mainly entered the region between domains 2 and 3, and this region is important for 3CLpro to form a dimmer [30] . M6 was reported to inhibit viral entry, accordingly it bound the fusion cone of spike protein; this cone structure is important for viral membrane fusion [31] (Fig. 2 and online supplementary Fig. S1 ). The 13 molecules passing the three-round selection process were then compared to the three Chinese herbal databases, and we found 230 herbs containing these molecules (online supplementary file Table S1 ). We then evaluated these herbs for those containing 2 or more of the 13 natural compounds, leaving 125 results. We cross-referenced the 125 results with the classic categorizations for herbal usage in the TCMSP database, finally choosing 11 types that are traditionally used to treat viral respiratory infections. There are 26 herbal plants within the 11 types. The timeframe during the course of a viral infection that each of these 26 herbal plants ( Step 1 Step 2 Step 3 Step 2 Step 3 Step 1 Part 2 Traditional usage-compatible herbal selection qi-reinforcing drugs, antitussive antiasthmatics, pungent cool diaphoretics and phlegm-resolving medicines may all be used throughout the course of infection, whereas drugs belonging to the interior warming group may be best utilized in prevention. Each of the potentially effective herbal remedies contains many ingredients in addition to the antiviral ones found here. Thus, the general effects of each plant should be examined by combining the effects of all of the orally absorbable and biologically active ingredients in it. To evaluate the possible general in vivo effects of each of our identified herbs, we used the ADME indices listed above to examine each of the orally absorbable and drug-like ingredients recorded in the TCMSP database for each plant. We then extracted the target proteins for each ingredient which had passed the screening process. All proteins belonging to a single plant were combined as input on the online protein-protein interaction analysis server, String, to find the pathway enrichment. For the 26 herbs, about 1/3 of the top 30 KEGG-enriched pathways (mean = 11) were related to regulating viral infection, immune/inflammatory reactions and hypoxia response, indicating that they are potentially effective treatments for viral respiratory infection ( Fig. 3 and online supplementary Fig. S2 ). Note that some of the herbal plants selected here had been reported to be effective for SARS-CoV infection in 2003 (online supplementary Table S2 ). In this work, we undertook a multiple step selection process and screened out 26 herbal plants with a high probability of directly inhibiting the novel coronavirus (2019-nCoV), possibly providing instant help in the prevention and treatment of the pneumonia that it can cause. While mainly in China at this point, viral spread is ongoing and has affect persons worldwide. Two principles guided our screening work. The first is that the anti-coronavirus (2019-nCoV) components contained in the source plants should be absorbable via oral prescription. This principle requires that the herbs selected should contain biologically proven anti-coronavirus (2019-nCoV) ingredients, and that these natural compounds should pass the drug-likeness and oral bioavailability evaluations. Therefore, we conducted a three-step screening process. First, we extracted natural compounds verified in PubMed as being effective in treating SARS or MERS coronavirus and then cross-checked these compounds in the Chinese herbal databases. There were 115 overlapping compounds. This method was an expeditious way to identify natural components both preexisting in Chinese herbal treatment and having a high possibility of anti-coronavirus (2019-nCoV) activity. This is important, as the anti-coronavirus effects of the selected compounds have been biologically confirmed, and the genetic similarities between coronavirus (2019-nCoV) and SARS or MERS coronavirus are high [24, 25] . The anti-coronavirus effects of the natural compounds screened by the above method have been mainly confirmed in vitro by direct loading onto cultured cells, thus it does not guarantee their effectiveness in vivo, especially with oral preparation-the principal way in which Chinese herbals are administered. Therefore, to meet the first principle, we ran ADME filters on the natural compounds selected by 4 indices, as used by Hu et al. [26] . Among the 115 compounds highlighted by our first step, only 13 passed this screening, showing the necessity of such a test. The novel coronavirus has some mutations when compared to SARS or MERS coronavirus, so the natural compounds effectiveness against the two previous coronaviruses might not be present in the new virus. To reduce this risk, and as the third step of our first principle, we reconstructed the 3D structure of the new coronavirus using the reported structures of SARS and MERS coronavirus proteins as a guide, and then used molecular docking technology to simulate whether the 13 natural compounds selected could combine with the structures we constructed for the new coronavirus proteins. All 13 compounds could bind to the proteins as predicted for the new coronavirus. We believe that the high success rate of our docking screening was due to the high genetic similarity between the new coronavirus and the SARS or MERS virus [24, 25] . Our second principle for screening should also be emphasized and elaborated upon. It states that the selected herbal plants must conform to traditional usages. There are many kinds of Chinese herbs that have been used for thousands of years. Based on this rich history and experience, Chinese herbal medicines are divided into different types, each type dedicated to certain kinds of diseases. Ignoring these grouping guidelines can lead to serious side effects. Therefore, as a further condition for the medicine screened here, we verified that they have been routinely used to treat viral pneumonia. To meet this principle, we conducted another threestep screening process for the herbal plants. First, we searched the Chinese medicine database for herbs containing the 13 natural compounds identified. Herbs containing at least 2 of these potentially useful compounds were selected, and a total of 125 herbal plants were identified. The second step in targeted plant selection was based on type classification. Of the 125 results, only 26 herbs were found to be routinely used in treating viral respiratory infection. Finally, network pharmacological analysis was performed to predict the possible therapeutic effects of these 26 plants. Because Chinese herbal medicines contain many ingredients, and multiple absorbable ingredients might exert their effects on the body, the general effects of herbs may be dictated by all of the absorbable ingredients they contain. With this consideration in mind, we extracted the recorded ingredients of each of the plants selected from the Chinese medicine database and screened these ingredients for drug-likeness and oral availability (via ADME filter) [26, 32] . The target proteins of all ingredients passing ADME selection were used for network enrichment to predict the general effects of the herbal plants. For all the plants analyzed, nearly half of the top 30 pathways enriched in KEGG are related to antiviral, immune/inflammatory responses and hypoxia response indicating that these herbs are suitable for anti-viral usage. In fact, some of the herbal plants selected here had been reported effective in against SARS-CoV infection in 2003 (online supplementary Table S2 ). We thought that the general antiviral and immune/inflammation effects predicted for the 26 plants are correlated with the fact that these plants were selected according to Chinese herbal type classifications. Of course, it should be pointed out that Chinese herbs that have not been identified through this screening process may still have beneficial effects. Further, considering that the biologically validated natural compounds reported in the literature cannot cover all antiviral natural compounds, and the natural compounds included in the Chinese medicine database are not complete, the process that we have followed may have excluded herbs that would be well suited to this treatment. Nevertheless, the purpose of this screening was to provide a rational approach for selecting Chinese herbal medicines with a high potential efficacy in treating 2019-nCoV and related viruses. The specific dosage and usage of each herb should be determined based on patients' manifestations. Finally, the key step in this screening was molecular docking. The 3D structures of the proteins used here are based on reported gene sequences. If the virus mutates during transmission, a new screening is recommended. In conclusion, this work has identified several Chinese medicinal plants classified as antiviral/pneumonia-effective that might directly inhibit the novel coronavirus, 2019-nCoV. Additionally, we propose screening principles and methods which may provide guidance in screening antiviral drugs from other natural drug databases. DZ conceived the study, participated in its design, coordination, and all the work processes. KW participated in herbal selection. XZ participated in data collection and network pharmacology analysis. SD helped to collect data. BP helped to draft the manuscript. Full-genome evolutionary analysis of the novel corona virus (2019-nCoV) rejects the hypothesis of emergence as a result of a recent recombination event World Health Organization. Statement on the second meeting of the International Health Regulations (2005) Emergency Committee regarding the outbreak of novel coronavirus (2019-nCoV) Statistical evidence for the usefulness of Chinese medicine in the treatment of SARS Quaternary structure, substrate selectivity and inhibitor design for SARS 3C-like proteinase Virtual screening for finding natural inhibitor against cathepsin-L for SARS therapy Synthesis of novel test compounds for antiviral chemotherapy of severe acute respiratory syndrome (SARS) Small molecules targeting severe acute respiratory syndrome human coronavirus SARS-CoV protease inhibitors design using virtual screening method from natural products libraries Antiviral activity of glycyrrhizic acid derivatives against SARS-coronavirus Identification of natural compounds with antiviral activities against SARS-associated coronavirus Binding interaction of quercetin-3-b-galactoside and its synthetic derivatives with SARS-CoV 3CL(pro): structure-activity relationship studies reveal salient pharmacophore features Evaluation of polyphenols from Broussonetia papyrifera as coronavirus protease inhibitors Ubiquitylation in innate and adaptive immunity Ubiquitination, ubiquitin-like modifiers, and deubiquitination in viral infection Inhibitors of SARS-3CLpro: virtual screening, biological evaluation, and molecular dynamics simulation studies Angiotensinconverting enzyme 2 is a functional receptor for the SARS coronavirus Specific plant terpenoids and lignoids possess potent antiviral activities against severe acute respiratory syndrome coronavirus SARS-CoV 3CLpro inhibitory effects of quinone-methide triterpenes from Tripterygium regelii Tanshinones as selective and slow-binding inhibitors for SARS-CoV cysteine proteases Dieckol, a SARS-CoV 3CL(pro) inhibitor, isolated from the edible brown algae Ecklonia cava Papain-like protease (PLpro) inhibitory effects of cinnamic amides from Tribulus terrestris fruits Chalcones isolated from Angelica keiskei inhibit cysteine proteases of SARS-CoV High-throughput screening and identification of potent broad-spectrum inhibitors of coronaviruses A pneumonia outbreak associated with a new coronavirus of probable bat origin 2019-nCoV: new challenges from coronavirus A network pharmacology study on the active ingredients and potential targets of Tripterygium wilfordii Hook for treatment of rheumatoid arthritis A simple click by click protocol to perform docking: AutoDock 4.2 made easy for non-bioinformaticians Safe, high-throughput screening of natural compounds of MERS-CoV entry inhibitors using a pseudovirus expressing MERS-CoV spike protein Severe acute respiratory syndrome coronavirus papain-like protease: structure of a viral deubiquitinating enzyme Coronavirus main proteinase (3CLpro) structure: basis for design of anti-SARS drugs Crystal structure of severe acute respiratory syndrome coronavirus spike protein fusion core Systems pharmacology for investigation of the mechanisms of action of traditional Chinese medicine in drug discovery Thanks to Mr. Steven Schwab for tidying up the English expressions and grammar. The authors declare no competing interests. Supplementary data to this article can be found online at https://doi.org/10.1016/j.joim.2020.02.005.