key: cord-0070698-roujrnph authors: Zhu, Kunyao; Zhang, Man; Long, Jia; Zhang, Shuqi; Luo, Huali title: Elucidating the Mechanism of Action of Salvia miltiorrhiza for the Treatment of Acute Pancreatitis Based on Network Pharmacology and Molecular Docking Technology date: 2021-11-24 journal: Comput Math Methods Med DOI: 10.1155/2021/8323661 sha: cdfdcbdfcd661af5c25d6e49598a54f08b06c8f2 doc_id: 70698 cord_uid: roujrnph Using network pharmacology and molecular docking, this study investigated the molecular mechanisms by which the active components in Salvia miltiorrhiza can alleviate acute pancreatitis. Initially, the active components of Salvia miltiorrhiza and the targets collected from the GeneCards database were screened based on the platform of systematic pharmacology analysis of traditional Chinese medicine. Subsequently, the active components were intersected with the disease targets. Also, interactions among the targets were computed using the STRING database. Biological function and pathway enrichment were analyzed using the Cluster Profiler package in the R software. Protein-protein interaction and component target pathway network were constructed using the Cytoscape software. Ultimately, the key targets and their corresponding components in the network were verified using the AutoDock Vina software. The results showed Salvia miltiorrhiza had 111 targets for acute pancreatitis. The biological process (BP) analysis showed that the active components of Salvia miltiorrhiza induced a drug response, positive regulation of transcription by RNA polymerase II promoter, signal transduction, positive regulation of cell proliferation, and negative regulation of apoptosis. Furthermore, the KEGG enrichment analysis screened 118 (P < 0.05) signaling pathways, such as the pathways related to cancer, neuroactive ligand-receptor interaction, PI3K-Akt signaling pathway, and cAMP signaling pathway, to name a few. Finally, molecular docking showed that the active components of Salvia miltiorrhiza had a good binding affinity with their corresponding target proteins. Through network pharmacology, this study predicted the potential pharmacodynamic material basis and the mechanisms by which Salvia miltiorrhiza can treat acute pancreatitis. Moreover, this study provided a scientific basis for mining the pharmacodynamic components of Salvia miltiorrhiza and expanding the scope of its clinical use. Acute pancreatitis (AP) is a common clinical inflammatory disease caused by the abnormal activation of trypsinogen. Approximately 15-25% of the patients with this abnormality might develop severe acute pancreatitis. Also, clinical treatment of AP is difficult for several reasons, including the rapid development of the disease, many complications, and high mortality. Therefore, it is necessary to find a cure for AP. Salvia miltiorrhiza, an inexpen-sive and easily available traditional medicine, may be used as a routine drug in the comprehensive treatment of AP. It promotes blood circulation, removes blood stasis, dredges channels, relieves pain, dilates microvessels, and improves microcirculation. Clinical studies reported the significant role of [1] Salvia miltiorrhiza in the treatment of AP. However, the specific mechanism by which Salvia miltiorrhiza can alleviate AP is unclear. Further studies are required to identify the specific disease targets of the components of Salvia miltiorrhiza. Network pharmacology has recently become popular globally. It provides a holistic perspective of the relationship between drugs and diseases by constructing a "drug-gene target disease" network. Due to its holistic and systematic characteristics, it has been widely used to predict the potentially active components in traditional Chinese medicine and study their mechanism of action on disease targets. Molecular docking, an essential method of virtual drug screening, is important for modern pharmacological research. It predicts the binding affinity between drugs and receptors by computing the interactions and binding energy between small-molecule ligands (drugs) and protein receptors. This study used a combination of network pharmacology and molecular docking methods to investigate the mechanism of action of the active components of Salvia miltiorrhiza for the treatment of acute pancreatitis. It provides a theoretical basis for the screening of the active components of Salvia miltiorrhiza and promotes further research on Salvia miltiorrhiza. The database, analysis platform, and software used in this study are shown in Table 1 . Salvia miltiorrhiza Bge. We used the Chinese medicine systems pharmacology analysis platform TCMSP and inserted "Danshen" as the keyword to retrieve all its chemical components. Pharmacokinetics (oral bioavailability) and druglike (DL) are often used by researchers to screen for active ingredients. OB values indicate the percentage of drugs that enter the human circulation. High OB values are usually the key indicators for evaluating the effectiveness of drugs. DL values indicate the degree of similarity between the components of traditional Chinese medicine and known chemical drugs. It is important for determining whether the components of Chinese Materia medica can exhibit pharmacological activities in organisms. Drug-likeness (DL) is based on the analysis of the physicochemical properties and/or structural characteristics of existing small organic drugs, especially those with undesirable properties, for example, compounds with poor ADMET (absorption, distribution, metabolism, excretion, and toxicity) [2] . Based on the molecular descriptors and the Tanimoto coefficients, the research team of RU [3] and others constructed a model to evaluate DL in the TCMSP database represented as TðA, BÞ = A ⋅ B/kAk 2 + kBk 2 − A ⋅ B, where A is the molecular descriptor for a given herbal component, and B represents the average value of this property for all molecules in the drug bank database [4] . DL is related to simple molecular properties such as molecular weight, physicochemical properties, and the number of convertible bonds or aromatic rings. In general, DL is related to the potential success of drug discovery in terms of pharmacokinetics and safety [5] . Because the discovery of drugs from candidate compounds is not always successful, the DL method has been widely used to screen compounds with undesirable properties, especially those with poor ADMET spectra. For TCPSP database compounds, DL ≥ 0:18 has been used as a standard for screening bioactive compounds through the traditional Chinese medicine systems pharmacology analysis [3, 4] . In conclusion, the active components of Salvia miltiorrhiza were screened under the conditions of OB ≥ 30% and DL ≥ 0:18, and the information of the qualified compounds was recorded and maintained in a tabular form using Microsoft Excel. Additionally, the quality control components of Salvia miltiorrhiza and the reported compounds with definite biological activities in the Chinese Pharmacopoeia were also analyzed as candidates. Then, the chemical structure of the above candidate components was constructed by the ChemBioDraw Ultra software and stored in MOL2 format. The final list of components was prepared after removing the duplicate active components. This list was used to search the TCMSP analysis platform for targets corresponding to the active components. Furthermore, the gene names of all the targets were obtained using the UniProt database. With "acute pancreatitis" as the keyword and the screening criterion set as median screening, disease-related targets were retrieved using the GeneCards database. "Median screening" ensured the reliability of the candidate targets. Furthermore, the names of AP-related targets were standardized using the UniProt database. Construction of the PPI Network. Venny 2.1 was used to build a Wayne map by providing Salvia miltiorrhiza and AP as the component target and disease target, respectively. The map provided information about common targets. These common targets were imported into the STRING database after specifying the "species" and "minimum required interaction score" as "Homo sapiens" and Table 1 : Information of the databases, analysis platform, and software. 3.5. Gene Enrichment Analysis. The function of these common targets in signal pathways was further investigated using the GO function and KEGG pathway enrichment (P < 0:05), which were in turn computed using the Cluster Profiler package in the R software. 3.6. Network Construction. Data comprising the active components and their corresponding targets were stored in files and imported to the Cytoscape software, where the component-target network diagram was constructed. This network was used to investigate the pharmacological mechanism of action of the active components of Salvia miltiorrhiza. Furthermore, the components and their AP targets underwent reverse screening based on the first 20 pathways enriched by KEGG. Subsequently, these screened data were used to build a network relation table of "componentstargets-pathways" using the Cytoscape software. Verification. From the constructed "Component-Target-Pathway" network, the components and their corresponding targets with moderate to high values were selected for molecular docking studies using AutoDock Vina. A series of steps were performed during molecular docking. (i) The 3D structure of the target protein was retrieved from the PDB database, and the ligand and nonprotein molecules (such as water molecules) were removed using the PyMoL software. The target protein files were saved in the PDB format. (ii) The target protein and the components of Salvia miltiorrhiza (in mol2 format) were converted into the pdbqt file format. (iii) The active site in the target protein was prepared for docking by setting the proper Grid Box size and coordinates. Then, the active components of Salvia miltiorrhiza were docked to the active site of the target protein. (iv) Finally, the components that demonstrated the lowest binding energy and highest binding affinity toward the target proteins were selected. The docked complexes of the selected components were visualized using the PyMoL software. A total of 202 active chemical components of Salvia miltiorrhiza were retrieved in this study. Using specific screening criteria, such as OB ≥ 30%, DL ≥ 0:18, searching for quality control components in the Chinese Pharmacopoeia, and the removal of nontarget components, 65 active components of Salvia miltiorrhiza were selected. These components are presented in Table 2 . Bunge. From the TCPSP database, 65 candidate targets of Salvia miltiorrhiza were selected. After removing the inactive components, 58 active components of Salvia miltiorrhiza Bunge were selected, and their corresponding 919 targets were obtained. Finally, the gene names corresponding to the target protein names were obtained from the Uni-Prot database. Analysis of Salvia miltiorrhiza Bunge. As shown in Figure 1 , the DANSHEN component-target network contained 173 nodes (58 candidate components and 115 associated targets) and 744 edges. Network topology analysis was performed using the network analyzer function in Cytoscape. The degree value (degree) is an important index for describing network nodes. The average number of target sites per compound was 12.83, and the average number of target sites per compound was 6.47. It was found that 14 targets interacted with ≥20 active compounds. The top 10 compounds were luteolin , and tanshinone I (Mol007119). Additionally, more than 25 targets interacted with ≥10 active compounds. The top 10 targets were PTGS2, NCOA1, SCN5A, OPRM1, CHRM1, ADRB2, F2, CA2, CHRNA7, and ACHE. The results suggested that a single component of Salvia miltiorrhiza affected multiple targets. Screening and PPI Network Analysis of Salvia miltiorrhiza in Treatment of Acute Pancreatitis. Using "acute pancreatitis" as the keyword, 8,319 targets were obtained from the Genecards database. After the median screening of these targets, 4,167 targets were selected. Later, the inter-action of 145 active components of Salvia miltiorrhiza with the 4,167 targets related to acute pancreatitis was investigated, and 111 common targets were found. These results indicated the active role of the components of Salvia miltiorrhiza in acute pancreatitis through their interaction with multiple targets (Figure 2 ). To further investigate the mechanism of interaction between the components of Salvia miltiorrhiza and the targets of AP, the 111 common targets were introduced into the STRING database for PPI network analysis (Figure 3 ). The results of this analysis were imported into Cytoscape 3.8.0 to construct the PPI network. There were 54 nodes and 82 edges in the PPI network, and the average degree of freedom was 9.13. In the PPI network, 65 protein degrees were higher than the average, and eight key nodes had moderate values according to their target nodes. These eight nodes, viz., PTGS2, NCOA1, SCN5A, OPRM1, CHRM1, ADRB2, F2, and PTGS1, were selected for further analysis. The proteins indicated by these key nodes interacted with a higher number of target proteins compared to the proteins in the other nodes in the PPI network. Hence, these nodes played an important role in the PPI network. In the GO enrichment analysis, 503 biological processes (defined as a series of events generated by an orderly combination of one or more molecular functions, BP), 113 molecular functions (described as activity at the molecular level, MF), and 59 cell components (each part of the cell and extracellular environment, CC) were obtained. The resultant BP, MF, and CC were sorted according to corrected P values. Then, to better describe the statistical significance of the conclusion data, the top 20 BP, MF, and CC were used to draw bubble charts and histograms, as shown in In the bubble chart, the size of the bubble represented the number of genes enriched in BP, CC, and MF entries, while different bubble colors represented different enrichment degrees of target genes in each BP, CC, and MF entry. From the results of the GO analysis, it was observed that many BP, with a high enrichment degree and a large number of enriched genes, were mainly involved in drug response, positive regulation of RNA polymerase, promoter transcription, signal transduction, positive regulation of cell proliferation, and negative regulation of apoptosis. Additionally, 118 signaling pathways (P < 0:05), screened using KEGG enrichment analysis, mainly involved pathways related to cancer, neuroactive ligand-receptor interaction signaling pathway, PI3K-Akt signaling pathway, and cAMP signaling pathway (Table 3) . Analysis of "Composition-Target-Pathway" Network of Salvia miltiorrhiza. From the KEGG enrichment analysis, the top 20 pathways with 153 targets for AP and 58 active ingredients of Salvia miltiorrhiza were selected. A "component-target-pathway" was constructed with the selected targets and compounds using the Cytoscape software ( Figure 12 ). Among the compounds, luteolin, tanshinone IIA, dihydrodanshen lactone, 2-isopropyl-8-methylphenyl-3, 4- .0 kcal/mol indicated reasonable binding affinity, good binding affinity, and strong binding affinity, respectively. To confirm the results of the network analysis, four targets with high median scores from the "component-target-pathway," which included PTGS2, RELA, TNF, and AKT1, were selected as receptors for molecular docking studies. Six active components of Salvia miltiorrhiza, designated as the six core components, were used as ligands for docking analysis using AutoDock Vina. The docking results are shown in Table 4 . All the six core components of Salvia miltiorrhiza-Paleonia vulgaris showed binding energies of ≤-5 kcal/mol with the receptors (target proteins). It was observed that 83.3% of the components 7 Computational and Mathematical Methods in Medicine and quercetin to receptor-OPRM1 were higher than the binding affinity between these receptors and their original ligands. Additionally, the docking results confirmed the reliability of the results of network pharmacological analysis (Table 4 ). Finally, the PyMoL software was used to visualize the docked complexes. The docked complexes with low binding energies are illustrated in Figure 13 . TCM deems that inflammation of internal organs may be due to accumulation of pathogenic qi, dampness, and heat. Qi stagnation causes blood stasis, and dampness leads to distention. When the pathogenic heat flames extremely, carbuncle will form. Acute pancreatitis is classified into Computational and Mathematical Methods in Medicine "abdominal pain" or "internal carbuncle" in TCM. It can be relieved by the treating principles as follows: promoting qi and blood circulation and removing pathogenic heat by purgation [6] . In this study, the best tanshinones and flavonoid were screened from the high-frequency danshen which was used to treat AP by TCPSP and PubChem database, it is suggested that tanshinones and flavones in Salvia miltiorrhiza may be the main material basis for the treatment of AP. The former includes tanshinone iia, 4-methylenemiltirone, salviolone, neocryptotanshinone ii, and isocryptotanshi-none, and the latter includes luteolin. To further identify the key compounds of Salvia miltiorrhiza in the treatment of AP, we constructed a composite target network, and the compounds in the top 10 degrees were luteolin, tanshinone iia, dihydrotanshinolide, dan-shexinkum d,4-methylenemiltirone salviolone, isocryptotanshi-none, 2-isopropyl-8-methylphenanthrene-3,4-dione, neocryptotanshinone ii, and miltionone I. It was found that [7] luteolin protected SAP mice by inducing ho-1-mediated anti-inflammatory and antioxidant activities and inhibiting the activation of the NF-ΚB pathway, which demonstrated that luteolin could significantly improve the infiltration and necrosis of inflammatory cells in the pancreas. A study by Chen et al. [8] confirmed that tanshinone IIA can inhibit oxidative stress through the NRF2ROS pathway and protect mice from acute pancreatitis. Isocryptotanshinone is a tanshinone monomer, which is the active component of Salvia miltiorrhiza. Antibacterial, anti-inflammatory, blood stasis promoting, wound healing, and other effects can play a positive therapeutic role in all aspects of the treatment of AP [9] . In conclusion, several key compounds in Salvia miltiorrhiza can not only inhibit the release of inflammatory factors and reduce the inflammatory reaction but also reduce the absorption of intestinal endotoxin and the translocation of pathogenic bacteria. They can also inhibit the growth of cancer cells and play Computational and Mathematical Methods in Medicine an important role in improving microcirculation and alleviate edema and necrosis of pancreatic tissue. Many studies [10, 11] have shown that neocryptotanshinone, which is a diterpenoid isolated from Salvia miltiorrhiza, can inhibit LPS-induced inflammation by inhibiting NF-κB and iNOS signals. In a study by Fan et al. [12] , it was found that neocryptotanshinone II, a monomer found in Salvia miltiorrhiza, plays an important role in the regulation of inflammatory cytokines and their signal pathways. TNF, IL-6, and NFKB1 are involved in the inflammation-related RB1, GSK3B, CDKN1A, GSTP1, PTGS2, ADCY8, RELA, EGFR, PIK3CG, CASP9, EDNRA, RXRA, CCND1, MYC, CASP3, ERBB2, AKT1, MAPK1, PRKACA, JUN, HSP90AA1, NOS2, MMP1, MMP2, STAT3, FOS, MMP9, VEGFA, NFKBIA, AR, IL6, CDK4, CDK2, BCL2, MDM2, BIRC5, PPARG The composition-target network of Radix Salvia miltiorrhiza reflects its anti-inflammatory and anti-injury characteristics. Based on network analysis, PTGS2, NCOA1, SCN5A, OPRM1, CHRM1, ADRB2, CHRNA7, and ACHE were selected as the key targets of AP. Studies [13] have shown that PTGS2 has an antiinflammatory effect on AP. NCOA1 can regulate IL-17, participate in the acute inflammatory process of AP, and play an important role in autoimmunity [14] . It has been reported [15] that tanshinone IIA, an active component of Salvia miltiorrhiza, can relieve the mechanical and thermal pain induced by Freund's complete adjuvant by inhibiting the expression of p-ERK and NF-κB. The levels of PTGS2 can reduce inflammation by reducing the levels and mRNA expression of IL-1β, IL-6, and tumor necrosis factor in sercytokines IL-1 β, IL-6, and TNF-α, and hence, can reduce pain. OPRM1 is an opioid receptor. The mu-opioid receptor (MOR) encoded by the OPRM1 gene is the main target of most opioid drugs. It is also an important factor that exerts analgesic effects, imparts tolerance, and provides dependence [16] . Tanshinone IIA might exert an analgesic effect through the action of OPRM1 on the inflammatory sites of the central nervous system and pancreas. CHRM1 might reduce inflammation and disease activity in various animal models of intestinal inflammation through cholinergic activation mechanisms involving the activation of the +7 NACHR subtype [17] . There is growing evidence [18, 19] from AP patients for the existence of cholinergic anti-inflammatory pathways that are completed by the efferent Vagus nerve fibers and are thus able to respond quickly and effectively to the inflammatory state. The results showed [20] that Salvia miltiorrhiza can dilate blood vessels, promote blood circulation, and remove blood stasis, and thus, inhibit platelet aggregation, anticoagulation, antithrombus, and reduce plasma viscosity. At the onset of AP, the body experiences increased activation of the sympathetic nervous system, which is marked by an increase in the heart rate and blood pressure. The sympathetic nervous system acts on the β1 and β2 receptors in the heart and blood vessels, respectively, regulating blood pressure. The ADRB2 gene regulates blood pressure [21] , and the β2 adrenergic receptor (β2-AR), expressed by ADRB2, is extensively present in immune cells. In mice and humans, stimulation of β2 AR activates NOD2 signaling, promotes T cell differentiation, and regulates the secretion of related cytokines and inflammation, according to a study [22] on dendritic cells. CHRNA7 regulates the function of immune cells in the spleen and is involved in regulating the production of proinflammatory cytokines, which are an integral part of the cholinergic anti-inflammatory pathway that regulates inflammation [23] . The mechanism by which Salvia miltiorrhiza plays an anti-inflammatory role through the cholinergic antiinflammatory pathway is not clear. By analyzing the KEGG pathway in this study, the authors speculate that CHRNA7 might be expressed in acute pancreatitis and play a key role in regulating inflammation. Salvia miltiorrhiza might play an anti-inflammatory role in the cholinergic anti-inflammatory pathway, which might elucidate the anti-inflammatory mechanisms and suggest clinical applications of Salvia miltiorrhiza. 13 Computational and Mathematical Methods in Medicine ACHE, a substance released by cholinergic neurons, is a major transmitter of the parasympathetic nervous system that rapidly hydrolyzes the neurotransmitter acetylcholine ACHE and acts as a marker of systemic inflammation. Clinical studies [24] have shown that a decrease in the serum cholinesterase level is related to the severity of the condition in patients, and a greater decrease indicates higher mortality. Serum cholinesterase levels also drop when the pancreas is severely damaged. Sheng et al. [25] found that serum cholinesterase concentration was an independent risk factor for the prognosis of acute pancreatitis. Acute pancreatitis is a complicated, acute, and severe disease of the abdomen, with many complications and high mortality, involving different metabolic pathways. The compounds in Salvia miltiorrhiza were extracted from the database. The potential targets of Salvia miltiorrhiza were predicted among these compounds, and the KEGG pathway enrichment analysis was performed. The potential pathway of Salvia miltiorrhiza was obtained based on these targets. According to the KEGG enrichment analysis of the key targets and associated studies, the acute pancreatitis pathway of Salvia miltiorrhiza is mainly involved in the PI3K-Akt signaling pathway, HIF-1 signaling pathway, cAMP signaling pathway, cholinergic synapse, FOXO signaling pathway, insulin resistance, AMPK signaling pathway, pathways associated with cancer, pancreatic cancer, neuroactive ligandreceptor interaction, etc. As an important signal transduction pathway, PI3K/AKT controls a variety of signal transduction pathways, including proliferation, apoptosis, and stress. In the AP model induced by ligation of the pancreatic duct, the expression of p-Akt and the activity of NF-κB in pancreatic tissue were significantly increased [26] . The expression of hypoxia-inducible factor-lα (hypoxia-inducible factors-1 α, HIF-1 α), Beclin1 gene, and microtubule-associated protein-1 light chain 3-II (microtubule-associated protein 1 light chain 3-II, LC3-II) are closely related to the pathological development of inflammation [27] [28] [29] . Ischemia, hypoxia, and autophagy are the main characteristics of acute pancreatitis cell injury. HIF-1α is very sensitive to the concentration of oxygen in cells and reflects the degree of hypoxia in cells. Inflammatory injury increases the expression of HIF-1α and the release of Beclin1 and LC3-II, which are markers of autophagic activity; it also increases the level of autophagy. In pancreatitis, autophagy is activated but lysosomal degradation is inhibited, and the imbalance between increased autophagy and decreased lysosomal degradation leads to autophagy flux arrest, which can cause severe acinar cell degeneration and trypsinogen activation, leading to the development of AP [30] . At the onset of AP, the acute inflammatory response of AP activates the insulin resistance pathway, triggering inflammation and immune responses. Studies [31] have shown that tanshinone monomers (such as tanshinone IIA and cryptotanshinone) can effectively increase glucose consumption, improve insulin resistance, and decrease the activity of the insulin pathway in insulin-resistant cells, which attenuates the inflammatory response and reduces the oxidative stress response in the cell. The cAMP/PKA signaling pathway plays a role in regulating oxidative stress and inflammatory response in AP. Yang et al. [32] showed that inhibition of NF-κb, which is dependent on estrogen receptor α, inhibits the production of IL-6, and thus, reduces inflammation. Saito et al. [33] showed that cholinergic synapses can induce an immune-inflammatory response that can cause pancreatic edema and bleeding. The FOXO1 signaling pathway plays an important role in cell metabolism, oxidative stress, cell proliferation, and apoptosis [34] . The oxidative stress response induced by AP can target the insulin signaling pathway via the FOXO signaling pathway to maintain homeostasis and mitigate tissue damage. Imbalance in the AMPK signal pathway is closely related to oxidative stress [35] [36] [37] [38] . The mechanism of action of Salvia miltiorrhiza Bunge on the AMPK signal pathway is not clear; however, the authors speculate that Salvia miltiorrhiza Bunge might increase the activity of AMPK, decrease the level of ROS, increase the production of SOD, and enhance the antioxidative ability of an organism. The reduction of oxidative stress in acute pancreatitis might provide important directions for the study of the mechanism of oxidative stress and clinical application of Salvia miltiorrhiza. In summary, using network pharmacology methods and molecular docking techniques, this study identified quercetin and other active components in Salvia miltiorrhiza to be effective against acute pancreatitis by targeting PTGS2, NCOA1, and other key genes. The results of the KEGG enrichment analysis indicated that Salvia miltiorrhiza might alleviate AP by regulating the PI3K-AKT signal pathway and other signal pathways. This study showed that Salvia miltiorrhiza probably affects AP through multiple targets and pathways. In this study, the possible pharmacodynamic basis and the mechanism of action of Salvia miltiorrhiza in treating acute pancreatitis were predicted by network pharmacology. The results of molecular docking also showed that the active molecule of Salvia miltiorrhiza has good binding activity with key target proteins. This study provides a scientific basis for further experimental research and clinical application. The data used to support the findings of this study are all included within the article. GO:0005886~plasma membrane GO:0045211~postsynaptic membrane GO:0005887~integral component of plasma membrane GO:0045202~synapse GO:1902711~GABA-A receptor complex GO:0005901~caveola GO:0030054~cell junction GO:0005654~nucleoplasm GO:0005829~cytosol GO:0045121~membrane raft GO:0048471~perinuclear region of cytoplasm GO:0005634~nucleus GO:0043679~axon terminus GO:0009897~external side of GO:0019899:enzyme binding GO:0008144:drug binding GO:0005230:extracellular ligand-gated ion channel activity GO:0042803:protein homodimerization activity GO:0042802:identical protein binding GO:0008134:transcription factor binding GO:0046982:protein heterodimerization activity GO:0003707:steroid hormone receptor activity GO:0005496:steroid binding GO:0016907:G-protein coupled acetylcholine receptor activity GO:0004890:GABA-A receptor activity GO:0004879:nuclear receptor activity GO:0019903:protein phosphatase binding GO:0032403:protein complex binding GO:0051378~serotonin binding GO:0003682~chromatin binding GO:00019901~protein kinase binding GO:0030235:nitric-oxide synthase regulator activity GO:0004672:protein kinase activity GO:0002020:protease binding Description Figure Management of severe acute pancreatitis: itʼs all about timing The application of _in silico_ drug-likeness predictions in pharmaceutical research TCMSP: a database of systems pharmacology for drug discovery from herbal medicines Systematic understanding the mechanisms of vitiligo pathogenesis and its treatment by Qubaibabuqi formula Considering the impact drug-like properties have on the chance of success Combined lipidomics and network pharmacology study of protective effects of Salvia miltiorrhiza against blood stasis syndrome Luteolin protects mice from severe acute pancreatitis by exerting HO-1-mediated antiinflammatory and antioxidant effects Tanshinone IIA protects against acute pancreatitis in mice by inhibiting oxidative stress via the Nrf2/ROS pathway Salvia miltiorrhiza traditional medicinal uses, chemistry, and pharmacology Neocryptotanshinone inhibits lipopolysaccharide-induced inflammation in RAW264.7 macrophages by suppression of NF-κ B and iNOS signaling pathways Chem inform abstract: two new diterpenoids from cell cultures of Salvia miltiorrhiza-Chem Inform Pharmacological network reveals the active mechanism of qi-replenishing, spleen-strengthening, phlegm-dispelling, and blood-nourishing fufang on coronary heart disease MicroRNA-101-3p inhibits fibroblast-like synoviocyte proliferation and inflammation in rheumatoid arthritis by targeting PTGS2 Biochanin A enhances RORγ activity through STAT3-mediated recruitment of NCOA1 Anti-nociceptive effects of tanshinone IIA (TIIA) in a rat model of complete Freund's adjuvant (CFA)-induced inflammatory pain Opioids and their receptors: are we there yet? The vagal innervation of the gut and immune homeostasis Determinantbased classification of acute pancreatitis severity Validation of the revised Atlanta classifica tion and the determinant-based classification for severity of acute pancreatitis The efficacy of Danshen injection as adjunctive therapy in treating angina pectoris: a systematic review and meta-analysis ADRB2 polymorphisms and dyslipidemia risk in Chinese hypertensive patients Co-delivery of hypoxia inducible factor-1α small interfering RNA and 5-fluorouracil to overcome drug resistance in gastric cancer SGC-7901 cells Deletion of muscarinic type 1 acetylcholine receptors alters splenic lymphocyte functions and splenic noradrenaline concentration Application of serum cholinesterase in the evaluation of acute pancreattis The evaluation value of the concentration of serum cholinesterase combined with acute physiology andchronic health evaluationllin predicting the condition and prognosis of patients with severe acuteeatts Bile-pancreatic juice exclusion promotes Akt/NF-κB activation and chemokine production in ligation-induced acute pancreatitis Targeting protein-protein interactions in the HIF system An essential role for chaperone-mediated autophagy in cell cycle progression Hypoxiainduced autophagy is mediated through hypoxia-inducible factor induction of BNIP3 and BNIP3L via their BH3 domains Nuclear TP53: an unraveled function as transcriptional repressor of PINK1 Tanshinone IIA improves endoplasmic reticulum stress-induced insulin resistance through AMP-activated protein kinase Peroxiredoxin 2 inhibits granulosa cell apoptosis during follicle atresia through the NFKB pathway in mice1 Epigenetics of neuroinflammation: Immune response, inflammatory response and cholinergic synaptic involvement evidenced by genome-wide DNA methylation analysis of delirious inpatients Eriodictyol inhibits survival and inflammatory responses and promotes apoptosis in rheumatoid arthritis fibroblast-like synoviocytes through AKT/ FOXO1 signaling COVID-19: recent developments in therapeutic approaches Historical and introductory aspects of carotenoids Pinolenic acid ameliorates oleic acid-induced lipogenesis and oxidative stress via AMPK/SIRT1 signaling pathway in HepG2 cells Cigarette smoke exposure induces ROS-mediated autophagy by regulating sestrin, AMPK, and mTOR level in mice We thank Huali Luo for her leadership and support of this research; Ling Yao and Yanhua Chen for assistance with data analysis. Constructive, thoughtful reviews led the way to additional insights and greatly improved the manuscript. The work was supported by grants from Scientific and Technological Research Program of Science and Technology Bureau of Jiulongpo District, Chongqing (no. 2019-02-006-Z). The authors have no conflict of interest. Conceptualization was done by HL. Data curation was done by KZ, MZ, JL, and SZ. Formal analysis was done by KZ and MZ. Funding acquisition was done by KZ. Methodology was done by KZ and HL. Project administration was done by HL. Supervision was done by KZ. Validation was done by KZ. Visualization was done by KZ and MZ. Writing-original draft was done by KZ, MZ, JL, and SZ. Writing-review and editing was done by KZ and HL.