Xiu-Yan Li1,2, Xiao-Tong Zhang3, Yi-Cheng Jiao1, Hang Chi1, Ting-Ting Xiong1, Wen-Jing Zhang1, Mi-Nan Li1, Yan-Hong Wang4. 1. Key Laboratory of Basic and Application Research of Beiyao, Heilongjiang University of Chinese Medicine, Harbin 150040, Heilongjiang Province, China. 2. College of Pharmacy, Harbin Medical University, Harbin 150040, Heilongjiang Province, China. 3. Key Laboratory of Basic and Application Research of Beiyao, Heilongjiang University of Chinese Medicine, Ministry of Education, College of Pharmacy, Harbin 150040, Heilongjiang Province, China. 4. Key Laboratory of Basic and Application Research of Beiyao, Heilongjiang University of Chinese Medicine, Harbin 150040, Heilongjiang Province, China. 799378826@qq.com.
Core Tip: Although some studies have suggested that Ruyi Jinhuang powder (RHP) has a therapeutic effect on diabetic foot ulcers (DFUs), few have used component analysis, investigated the mechanism of action, and utilized wound-healing experiments. The components of RHP were used to predict the mechanism of action, and wound healing was observed by establishing a DFU rat model to further prove the therapeutic effect of RHP on DFU to finally determine a possible mechanism of action.
INTRODUCTION
Ruyi Jinhuang powder (RHP) is included in “Authenticity of Surgery” written by Chen Shigong in the Ming Dynasty, which includes “THF [Trichosanthin (Tian Hua Fen)], DH [Rhubarb (Da Huang)], HB [Phellodendron (Huang Bai)], JH [Turmeric (Jiang Huang)], BZ [Angelica dahurica (Bai Zhi)], TNX [Arisaema (Tian Nan Xing)], CZ [Atractylodes lancea (Cang Zhu)], HP [Magnolia officinalis (Hou Po)], CP [Pericarpium Citri Reticulatae (Chen Pi)] and GC [Glycyrrhiza uralensis (Gan Cao)]”[1] and is now included in the Chinese Pharmacopoeia. This prescription is used for swelling relief and pain relief. In this prescription, Trichosanthes is the monarch medicine, and the minister medicine rhubarb with the same cold and bitter taste is used to purge fire detumescence; angelica dahurica and turmeric are the minister pungent medicines with compatibility to discharge pus pain; pericarpium citri reticulatae, Magnolia officinalis, atractylodes lancea and glycyrrhiza uralensis are combined to remove dampness and regulate Qi; and arisaema alleviates swelling pain. The above five medicines act together as adjuvants with glycyrrhiza uralensis to reconcile and detoxify the medicine. Modern clinical applications mainly include cutaneous vasculitis, gouty arthritis, herpes zoster and diabetic foot ulcer (DFU). After many studies on its pharmacological effects, it was found that RHP can inhibit bacterial infection, increase lysosomal content, enhance immune defenses and inhibit inflammation. The traditional preparation is through the addition of honey to the powder, which is then directly applied to the affected area to treat diseases; however, the formulation has been innovated using the original preparation through the continuous implementation of modern technology and made into creams, cataplasms, films, and sponges[2,3].DFU is a diabetic complication. Its pathogenic causes are often vascular and nerve lesions, resulting in lower limb infection and the formation of foot ulcers[4]. The clinical manifestations are ischemic necrosis or damage to skin tissue, incomplete skin, wound exudation, abscess generation, etc. From the perspective of traditional Chinese medicine, DFUs are caused by deficiency of Qi and Yin, weakness of pulse, and blockage of dampness-heat and blood stasis toxin. Traditional chinese medicine also has shown promising effects with regard to safety and renoprotection in some prospective, multicenter, randomized, controlled clinical trial conducted[5]. And studies have found that RHP has a good effect in the treatment of DFUs in traditional Chinese medicine. Shao et al[6] found that its combination with antibiotics can quickly reduce the swelling of patients' ulcers to shorten the course of treatment. Liu et al[7] treated 40 patients with diabetic skin abscesses with the external application of RHP and found that the rate of effective treatment in the treatment group was higher than the control group. This result shows that RHP has a therapeutic effect on DFU. Zhang[8] treated patients with RHP and Simiao Tongluo decoction and found that this prescription can promote wound healing to promote improvement and play a therapeutic role.Gas chromatography-mass spectrometry (GC–MS) and ultra-performance liquid chromatography-mass spectrometry (UPLC–MS) are widely used in the separation and identification of complex components, of which GC–MS is mainly used to separate the sample into volatile products in the instrument by pyrolysis of the sample at high temperature, with the advantages of high sensitivity, large information content, high efficiency, and small sample requirements. UPLC–MS technology allows the sample to be separated in the mobile phase after the ionization process through fragment ion mass number analysis and identification. This technology can compensate for the disadvantages of GC–MS, which cannot analyze components with features such as strong polarity, thermal instability, and difficult volatilization; UPLC-MS has the advantages of low detection limit, high automation, wide analysis range, and short analysis time. Network pharmacology research is mainly based on databases and software to obtain important target information for drug treatment of diseases and establish network connections to predict its mechanism of action. Because network pharmacology analysis shows synergy and compatibility with traditional Chinese medicine in the treatment of diseases and the therapeutic principle of syndrome differentiation and treatment, it is widely used in traditional Chinese medicine to explore the mechanism of action in the treatment of diseases from multiple targets. Wound healing is caused by a variety of molecular proteins that affect cells and remodel the tissue. Hematoxylin-eosin (HE) staining electron microscope observation is a more intuitive way to observe the healing of the wound surface. Fibroblasts can be clearly observed to evaluate drug treatment.Although it has been proven that RHP has a therapeutic effect on DFU, the specific mechanism of action is not clear. In this study, the chemical components of RHP were extracted by supercritical extraction and alcohol extraction and separated and identified by GC–MS and UPLC–MS techniques, respectively, to obtain the active ingredients of the RHP formula, and the target pathway for the treatment of DFU was studied using network pharmacological analysis. In addition, a rat model of DFU was established to verify the therapeutic effect of RHP on wound healing, providing a direction for further clinical research.
MATERIALS AND METHODS
Materials
RHP (201202056, Jilin Shuangshi Pharmaceutical Co., Ltd) was purchased from Harbin Shiyitang Chinese Herbal Medicine Co., Ltd. (Harbin, China). Methanol (purity) and acetonitrile (purity) were all obtained from Merck Co., Inc. (Germany). Formic acid (purity) was obtained from West Asia Chemical Technology Co., Ltd. All other chemicals and solvents were of analytical grade.
Animals and drug administration
Thirty healthy male Sprague Dawley rats (weights, 200 ± 20 g) were purchased from Jinan Pengyue Experimental Animal Breeding Co., Ltd. (Jinan, China) and housed for adaptive feeding for 2 d. The animal facilities and protocols were approved by the Animal Ethics Committee of Heilongjiang University of Chinese Medicine (Heilongjiang, China).Thirty rats were divided into 3 groups with 10 rats in each group, including the blank, model and RHP groups. The model and RHP groups were injected with streptozotocin/0.1 mol·L-1 citrate buffer solution (1/100, g/v). After 72 h, fasting blood glucose was measured with an electronic blood glucose meter (Sannuo Biosensor Co., Ltd, China). The modeling standard was that the random blood glucose level was greater than 12.0 mmol/L, accompanied by the typical symptoms of diabetes mellitus. In the establishment of the ulcer model, the rat hindfoot skin in each group was cut with scissors (the depth of the wound that would reach the fascia). The wound was cleaned before each treatment every day. The wound area was measured on the 5th, 7th, and 14th days, and the wound-healing rate was calculated. The formula for the healing rate was as follows: Healing rate (%) = (original wound area-unhealed area)/original wound area × 100%.After 14 d, the rats were sacrificed, and the skin around the wound was cut into 3 cm × 3 cm areas. The cut tissue was fixed with a 2.5% pentanediol solution and stored at low temperature (-80 °C).The sample tissues were placed in liquid paraffin, freeze-fixed and sliced. Dewaxing, hydration and dyeing with hematoxylin staining solution and eosin dye solution were performed. The samples were air dried, sealed and observed with transmission electron microscopy (HT7700, Hitachi, Japan).
Composition analysis of RHP by GC–MS
Supercritical fluid extraction (SFE) equipment (HA220-40-11, Nantong Huaan Supercritical Extraction Co., Ltd, China) was used to extract volatile oil. First, an appropriate amount of RHP sample was crushed and sifted through a 40 μm mesh; the medicinal material was placed into supercritical extraction equipment, and the pressure was boosted to the set parameters for extraction. The pressure of the extraction kettle was 25 MPa, and its temperature was 45 °C. The pressure of the separating kettle was 8 MPa, and its temperature was 60 °C. The pressure of the other separating kettle was 4.5 MPa, and its temperature was 37 °C. The pump frequency of the SFE was 18 Hz, and the flow rate was 60 L/h. After extraction, the materials were removed to obtain the SFE.GC–MS (5975B, Agilent Technologies, Inc., United States) was used to analyze the chemical composition of the SFE extraction. The prepared SFE extract was vortexed for 2 min, extracted with a solid-phase microextraction needle for 20 min, centrifuged for 20 min, and filtered with a membrane. Each of the samples was injected into GC–MS equipment equipped with an HP-INNOWAX (25 m × 0.20 mm, 0.40 μm) column (Agilent, United States) at 100 °C for 5 min. Then, the temperature was raised to 150 °C at 5 °C/min and then to 280 °C at 30 °C/min. The inlet temperature was 240 °C, and the flow rate of the carrier gas was 1.0 mL/min. The ion source temperature of MS with an EI source was 200 °C, and its transmission line temperature was 250 °C. The bombardment voltage was 70 eV.
Composition analysis of RHP by UPLC–MS
UPLC–MS (Ultimate 3000LC, Q Exactive HF, Thermo Fisher, United States) was used to analyze the chemical composition of RHP. First, the sample was pulverized, passed through a 40-μm mesh sieve, accurately weighed and placed in a stoppered conical flask. Hydrochloric acid/ethanol solution (1/100, v/v) was accurately added and ultrasonically treated for 40 min. The filtrate was shaken well, filtered and diluted into a 50-mL volumetric flask to obtain a sample solution. The reference substance was precisely weighed and placed in a 10-mL volumetric flask; methanol solution was added, and the solution was diluted to volume after ultrasonic treatment and shaken well to obtain the reference substance solution. Each of the samples was injected into UPLC–MS equipment equipped with a C18 chromatographic (2.1 mm × 100 mm, 1.8 μm) column (Zorbax Eclipse, Agilent, United States) at 30 °C. The flow rate was 0.3 mL/min. The mobile phase was water/formic acid (0.1/100, v/v) (A)/acetonitrile (B). The injection volume was 2 μL.Mass spectrometry conditions utilized positive and negative modes for UPLC–MS. Electrospray ionization was used in ionization mode with a sheath gas flow rate of 45 arb. The auxiliary gas flow rate was 15 arb. The purge gas flow rate was 1 arb. The electrospray voltage was 3.5 KV. The capillary temperature was 330 °C. The S-Lens RF level was 55%. The scan mode was full scan/dd-MS2 (TopN = 10) with a scanning range of 100-1500 m/z and a resolution of 120000/60000. The collision mode was high energy collision dissociation.
Building a drug-component-target-disease network relationship graph
Cytoscape (v3.7.2.) is an analysis software that shows the complex corresponding relationship of "drug-component-target-disease" in the form of a network graph. It can conveniently visualize a network relationship, perform network topology analysis, analyze the degree of connection between each node according to the relevant parameters, and thus enable researchers to draw the corresponding conclusions.
Determining the "Chemical Composition-Target" correspondence
The Traditional Chinese Medicine Systems Pharmacology (TCMSP) database (http://Lsp.nwu.edu.cn/) is a platform for the analysis and study of traditional Chinese medicines in many aspects, and it synthesizes the chemical composition and drug target data of traditional Chinese medicines[9], closely links diseases with targets and components, and explains the mechanism of action of traditional Chinese medicines[10]. The chemical components extracted from ten traditional Chinese medicines of the RHP formula were searched by entering the CAS number of chemical components in the TCMSP database, and the information related to the components could be obtained, of which the "Relatedtarget" column contained the targets corresponding to the component, and the obtained targets were integrated to determine the corresponding relationship between the components and the target.
Screening targets of RHP
Comparative Toxicogenomics Database (CTD) (http://ctdbase.org/) is a database that brings together detailed information on the intersection between genes, proteins and diseases[11], and the combination of these data with that of their pathways and functions can further elucidate the mechanism of diseases[12]. A component usually corresponds to one or more targets, and information unrelated to the treatment of DFU in the above integrated "component-target" dataset needs to be screened out. In the TCMSP database, the information related to each target in the "Relatedtarget" column above was viewed, the relevant target with the keywords "diabetes-related diseases", "pain" and "bacterial infection" was selected according to the disease type in "Relateddiseases", and its "GenecardID" was recorded. To prevent the limitations of the TCMSP database on the target and disease correspondence and make the study more accurate, the targets without keywords were searched in the CTD database to further determine whether they were related to DFU in the "Diseases" category. Two databases, TCMSP and CTD, were used to screen for and obtain the chemical compositions of targets related to DFU.
Molecular docking of ingredients and targets
DiscoveryStudio (DS v19.1.0) software is a life science molecular simulation software that can be used to establish molecular docking models. Cytoscape was used to visualize the targets and components, analyze the key targets and important components, perform molecular docking, and observe the binding effect of chemical components to target proteins. First, the structural model of the component was downloaded from the TCMSP database; the human protein number corresponding to the target from the UniProt database was queried, and the number from the PDB database was the input to obtain the three-dimensional structural model of the target protein. The component structure and protein structure were opened in DiscoveryStudio software; H2O and ligand were removed, the protein was hydrogenated, and the relationship between the molecule and the protein was established. In the end, the binding sites can be identified.
Pathway enrichment analysis
The STRING database (https://cn.string-db.org/) and the DAVID database (https://david.ncifcrf.gov) are mainly used to provide information on the individual target genes and proteins or the interaction between multiple histones and to perform GO or KEGG pathway enrichment analysis on them. The STRING database was used to observe the correlation between target proteins of RHP formula herbs, and the DAVID database was used to enrich KEGG pathways for target components and targets to obtain possible pathways for the treatment of DFU.
Statistical analyses
SPSS software (Version 24.0.0, Chicago, IL, United States) was used to analyze the wound-healing rate results of rats at different times, and all data are expressed as the mean ± SD. The differences were considered significant at P < 0.05.
RESULTS
In vivo wound healing effect
During continuous culture for 14 d, the wound areas of the rats on the first, third, seventh, and fourteenth days were recorded, and the cure rate was calculated. The results are shown in Table 1; the difference was significant. The results of the three groups of experiments showed that the cure rate of the RHP group was increased and reached 53.11%, indicating that RHP has a good therapeutic effect on the healing of DFU wounds in rats. The HE staining electron microscope results showed that after 14 d of uninterrupted administration, the ulcer tissues of 10 rats in the model group showed varying degrees of demyelination changes, unlike rats in the blank group, which showed symptoms of neuritis, as shown in Figure 1. In addition, neutrophil infiltration and granulation tissue loosening were observed. The mice treated with powder had good fibroblast function, mature granulation tissue, and restored nerve cells. The thickness of the stratum corneum and the integrity of the epidermis were good, showing a good state of recovery, and intact lymphocytes were observed, but there were a few mast cells. RHP promoted the healing of DFUs in rats by affecting fibroblasts and nerve cells.
Table 1
Statistical analysis of wound healing rate at different times
Time
0 d (%)
3 d (%)
7 d (%)
14 d (%)
Control group
0
10.41 ± 0.24
17.43 ± 0.13
33.73 ± 0.28
Experimental group
0
32.78 ± 0.46
38.14 ± 0.28
53.11 ± 0.07
Homogeneity of variance
—
√
√
√
P value
—
< 0.05
< 0.05
< 0.05
Figure 1
Hematoxylin-eosin staining electron micrograph of rat diabetic foot ulcer. A: The fibroblasts of blank group; B: The neuritis of model group; C: Mast cells of powder group; D: Nerve cells of powder group. The mice treated with powder had good fibroblasts function, mature granulation tissue, and restored nerve cells. The thickness of the stratum corneum and the integrity of the epidermis were good, showing a good state of recovery, and intact lymphocytes would be observed, but there were a few mast cells. Ruyi Jinhuang powder promoted the healing of diabetic foot ulcers in rats by affecting fibroblasts and nerve cells.
Hematoxylin-eosin staining electron micrograph of rat diabetic foot ulcer. A: The fibroblasts of blank group; B: The neuritis of model group; C: Mast cells of powder group; D: Nerve cells of powder group. The mice treated with powder had good fibroblasts function, mature granulation tissue, and restored nerve cells. The thickness of the stratum corneum and the integrity of the epidermis were good, showing a good state of recovery, and intact lymphocytes would be observed, but there were a few mast cells. Ruyi Jinhuang powder promoted the healing of diabetic foot ulcers in rats by affecting fibroblasts and nerve cells.Statistical analysis of wound healing rate at different timesThe GC–MS full spectrum identification results of RHP are shown in Table 2, and a total of 43 components were detected to obtain the total ion flow diagram (Figure 2), which mainly included alcohols, enes, esters, phenols, and other components, accounting for 27.91%, 20.93%, 20.93%, 16.28%, 4.65%, and 9.30% of the total, respectively. Among them, the higher components were ar-turmerone 24.33% (No. 23), tigerone 20.74% (No. 17), agarospirol 13.80% (No. 18), β-eudesmol 9.06% (No. 22), palmitic acid 2.31% (No. 43), and (E)-allylantone 2.23% (No. 29).
Table 2
Gas chromatography–mass spectrometry full spectrum identification results
The total ion flow diagram of gas chromatography–mass spectrometry. Among them, the higher components were ar-turmerone 24.33%, tumerone 20.74%, agarospirol 13.80%, β-eudesmol 9.06%, palmitic acid 2.31%, (E) atlantone 2.23%, β-quadruphellandrene 1.95%, methyloleate 1.87%, α-zingiberene 1.86%, and methylmitate 1.75%.
The total ion flow diagram of gas chromatography–mass spectrometry. Among them, the higher components were ar-turmerone 24.33%, tumerone 20.74%, agarospirol 13.80%, β-eudesmol 9.06%, palmitic acid 2.31%, (E) atlantone 2.23%, β-quadruphellandrene 1.95%, methyloleate 1.87%, α-zingiberene 1.86%, and methylmitate 1.75%.Gas chromatography–mass spectrometry full spectrum identification resultsThe UPLC–MS full spectrum identification results of RHP are shown in Table 3. A total of 46 components were detected. The main components include flavonoids, organic nitrogens, isopentenol esters, carboxylic acids and their derivatives, benzene and its substituted derivatives, diarylheptans, fatty acyls, anthracyclines and other components, accounting for 23.92%, 15.23%, 10.88%, 8.71%, 6.52%, 6.52%, 6.52%, 4.26% and 17.44% of the total, respectively. The positive and negative ion flow diagrams are shown in Figure 3. Among them, the components with the highest levels were α,α-trehalose (409.04%; No. 6), curcumin (307.20%; No. 37), berberine (232.25%; No. 24), (3R,5R)-1,3,5-trihydroxy-4-{[(2E)-3-(4-hydroxy-3-methoxyphenyl)-2-propenoyl]oxy}cyclohexanecarboxylic acid (201.34%; No. 16), bisdemethoxycurcumin (153.92%; No. 33), genistein (153.83%; No. 39), demethoxycurcumin (151.99%; No. 35), and citricol acid (145.29%; No. 9).
Table 3
Ultra-performance liquid chromatography full spectrum identification results
Ultra-performance liquid chromatography ion flow diagram of Ruyi Jinhuang powder. A: Positive ion flow diagram; B: Negative ion flow diagram. Among them, the components with high levels were α,α-trehalose 409.04% (No. 6), curcumin 307.20% (No. 37), berberine 232.25% (No. 24), (3R,5R)-1,3,5-trihydroxy-4-{[(2E)-3-(4-hydroxy-3-methoxyphenyl)-2-propenoyl]oxy}cyclohexanecarboxylic acid 201.34% (No. 16), bisdemethoxycurcumin 153.92% (No. 33), genistein 153.83% (No. 39), demethoxycurcumin 151.99% (No. 35), and citricol acid 145.29% (No. 9).
Ultra-performance liquid chromatography ion flow diagram of Ruyi Jinhuang powder. A: Positive ion flow diagram; B: Negative ion flow diagram. Among them, the components with high levels were α,α-trehalose 409.04% (No. 6), curcumin 307.20% (No. 37), berberine 232.25% (No. 24), (3R,5R)-1,3,5-trihydroxy-4-{[(2E)-3-(4-hydroxy-3-methoxyphenyl)-2-propenoyl]oxy}cyclohexanecarboxylic acid 201.34% (No. 16), bisdemethoxycurcumin 153.92% (No. 33), genistein 153.83% (No. 39), demethoxycurcumin 151.99% (No. 35), and citricol acid 145.29% (No. 9).Ultra-performance liquid chromatography full spectrum identification resultsRT: Radiotherapy.
Establishment and analysis of the network relationship of "Traditional Chinese Medicine-Ingredients-Target"
Eighty-nine components were obtained by GC–MS and UPLC–MS. According to the "target-disease" relationship in the database, 43 relevant components for the treatment of DFU were finally screened from the ten medicines. There were 36 related targets. The results are shown in Table 4. Cytoscape software was used to visualize the network relationship of "TCM-Ingredient-Target", which was constructed by the parameters of "Degree", "Betweenness" and "Closeness" of each node, as shown in Figure 4. The important medicinal materials of RHP for the treatment of DFU are turmeric, Magnolia officinalis, and rhubarb, and the important components are luteolin, naringin, demethoxycurcumin, gallic acid, methyl linoleate, caryophyllene, arylcurcumin, methyl palmitate, berberine and rheic acid. The important targets are posttranscriptional silencing (PTGS2), muscarinic acetylcholine receptor M2 (CHRM2), Dipeptidyl peptidase 4, topoisomerase II alpha (TOP2A), retinoic X receptor alpha (RXRA), Protein tyrosine phosphatase nonreceptor type 1, tumor necrosis factor (TNF), interleukin 6 (IL-6), etc.
The “traditional Chinese medicine-ingredient-target” network of Ruyi Jinhuang powder topology diagram. There are a total of 89 nodes and 186 edges; diamonds represent ten herbs of Ruyi Jinhuang powder, circles represent 43 related chemical components (represented in the form of CAS number), squares represent 36 target targets (represented in the form of GenecardID), nodes indicate that the greater Degree value, the darker color indicates the greater Closness value, and the thicker lines indicate the greater Betweenness value.
The “traditional Chinese medicine-ingredient-target” network of Ruyi Jinhuang powder topology diagram. There are a total of 89 nodes and 186 edges; diamonds represent ten herbs of Ruyi Jinhuang powder, circles represent 43 related chemical components (represented in the form of CAS number), squares represent 36 target targets (represented in the form of GenecardID), nodes indicate that the greater Degree value, the darker color indicates the greater Closness value, and the thicker lines indicate the greater Betweenness value.Ruyi Jinhuang powder “Chinese medicine-ingredients-target” correspondenceTNF: Tumor necrosis factor; IL: Interleukin; PTGS2: Posttranscriptional silencing; TOP2A: Topoisomerase II alpha; CHRM2: Muscarinic acetylcholine receptor M2; RXRA: Retinoic X receptor alpha; DPP4: Dipeptidyl peptidase 4; SCN5A: Sodium voltage-gated channel alpha subunit 5; HTR2A: 5-hydroxytryptamine (serotonin) receptor 2A; OPRM1: Opioid Receptor Mu 1; EGFR: Epidermal Growth Factor Recepto; PTPN1: Protein tyrosine phosphatase nonreceptor type 1; THF: Trichosanthin (Tian Hua); DH: Rhubarb (Da Huang); HB: Phellodendron (Huang Bai); JH: Turmeric (Jiang Huang); BZ: Angelica dahurica (Bai Zhi); TNX: Arisaema (Tian Nan Xing); CZ: Atractylodes lancea (Cang Zhu); HP: Magnolia officinalis (Hou Po); CP: Pericarpium Citri Reticulatae (Chen Pi); GC: Glycyrrhiza uralensis (Gan Cao).
Analysis of the docking results of important components and key target molecules
Molecular docking was performed between important components of luteolin, methyl palmitate, ar-turmeric, methyl linoleate, palmitic acid, demethoxycurcumin, and naringin and the key targets PTGS2, TOP2A, CHRM2, and RXRA. The results are shown in Table 5, in which LibDockScore indicates the docking effect, and the value indicates the docking effect. Seven components (luteolin, methyl palmitate, ar-turmeric, methyl linoleate, palmitic acid, demethoxycurcumin and naringin) had a large LibDockScore, indicating that the binding effect between these components and key targets was good. The component docking conformations shown in Figure 5 indicate the following: The complex naringin-PTGS2 was stabilized by five conventional hydrogen bonds, eight carbon-hydrogen bonds, and two bonds to the alkyl with residues including GLY, ANS, ARG, PHE, ASN, HIS, GLY, GLN, LEU; the complex methyl linoleate-RXRA was stabilized by three carbon-hydrogen bonds, two Pi bonds to the alkyl with residues including DA and DG; the complex luteolin-TOP2A was stabilized by four conventional hydrogen bonds, two carbon-hydrogen bonds, two Pi bonds to the alkyl, six Pi bonds to the benzene ring with residues including LYS, GLU, SER, and ASP; the complex naringin-TOP2A was stabilized by three conventional hydrogen bonds, four carbon-hydrogen bonds, one Pi bond to the alkyl, four Pi bonds to the benzene ring with residues including ASP, ARG, ASN, and LYS; the complex methyl palmitate-PTGS2 was stabilized by two conventional hydrogen bonds, one Pi bond to the alkyl with residues including VAL, ASN, and TRP; the complex ar-turmeric-CHRM2 was stabilized by nine Pi bonds to the alkyl, two Pi bonds to the benzene ring and Van der Waals interactions with residues including PHE, LYS, ILE, VAL, LEU, ASP, ARG, ASN, and THR; the complex methyl linoleate-PTGS2 two carbon-hydrogen bonds, one Pi bond and one bond to the alkyl with residues including VAL, ASN, ILE, and TRP; the complex ar-turmeric-PTGS2 was stabilized by one conventional hydrogen bond, one Pi bond and three bonds to the alkyl, two Pi bonds to the benzene ring with residues including GLN, PHE, VAL, PRO, and ASN; the complex palmitic acid-PTGS2 was stabilized by two conventional hydrogen bonds with residues including VAL and GLY; the complex demethoxycurcumin-PTGS2 was stabilized by four conventional hydrogen bonds, one carbon-hydrogen bond, three Pi bonds to the benzene ring with residues including GLY, VAL, ASN, and PHE.
Table 5
Docking parameters of chemical components and target protein molecules
Name
Target
LibDockScore
Luteolin
TOP2A
104.333
Methyl palmitate
PTGS2
114.221
ar-Turmeric
PTGS2
91.0625
CHRM2
88.3765
Methyl Linoleate
RXRA
138.583
PTGS2
130.838
Palmitic acid
PTGS2
114.349
Demethoxycurcumin
PTGS2
135.479
Naringin
TOP2A
141.197
PTGS2
164.305
PTGS2: Posttranscriptional silencing; TOP2A: Topoisomerase II alpha; CHRM2: Muscarinic acetylcholine receptor M2; RXRA: Retinoic X receptor alpha.
Figure 5
The virtual docking of bioactive ingredients from Ruyi Jinhuang powder for diabetic foot ulcer targets. A and D: The virtual docking of naringin with post-transcriptionalgenesilencing (PTGS2) and Topoisomerase II alpha (TOP2A); B and G: The virtual docking of methyl linoleate with Retinoid X Receptor alpha and PTGS2; C: The docking of Luteolin and TOP2A; E: The docking of Methyl palmitate and PTGS2; F and H: The virtual docking of ar-Turmeric with Muscarinic Acetylcholine receptor M2 and PTGS2; I: The virtual docking of palmitic acid and PTGS2. J represents the virtual docking of demethoxycurcumin and PTGS2.
The virtual docking of bioactive ingredients from Ruyi Jinhuang powder for diabetic foot ulcer targets. A and D: The virtual docking of naringin with post-transcriptionalgenesilencing (PTGS2) and Topoisomerase II alpha (TOP2A); B and G: The virtual docking of methyl linoleate with Retinoid X Receptor alpha and PTGS2; C: The docking of Luteolin and TOP2A; E: The docking of Methyl palmitate and PTGS2; F and H: The virtual docking of ar-Turmeric with Muscarinic Acetylcholine receptor M2 and PTGS2; I: The virtual docking of palmitic acid and PTGS2. J represents the virtual docking of demethoxycurcumin and PTGS2.Docking parameters of chemical components and target protein moleculesPTGS2: Posttranscriptional silencing; TOP2A: Topoisomerase II alpha; CHRM2: Muscarinic acetylcholine receptor M2; RXRA: Retinoic X receptor alpha.
Analysis of KEGG pathway enrichment results
The interaction between each target protein was preliminarily obtained by entering the target protein into the STRING database, as shown in Figure 6. Among the 36 targets, except for CHRM2, CHRM4, PCYT1A and Sodium voltage-gated channel alpha subunit 5, the remaining 32 target proteins are closely related to each other and may participate in multiple pathways. The target list was uploaded to the DAVID database to establish the "pathway-target" link. A total of 27 pathways related to DFU treatment were retrieved. The number of targets contained in each pathway is shown in Figure 7A. Cytoscape was used to construct the "pathway-target" network relationship diagram, which is shown in Figure 7B; the "pathway-target" relationship is shown in Figure 7C. Correspondence between channel categories is shown in Figure 7D. Most of the 27 pathway diagrams belong to the immune system and participate in the processes of inflammatory and infectious diseases, which demonstrate signal transduction effects. This finding shows that RHP can affect the body’s immune system, regulate inflammation through signal transduction, and exert curative effects in DFU. Further analysis of the literature shows that RHP is mainly related to the phosphoinositide 3-kinase (PI3K)-protein kinase B (Akt) signaling pathway, neuroactive ligand–receptor interaction and the forkhead box O (FoxO) signaling pathway in the treatment of DFU.
Figure 6
Interaction diagram of target proteins. The edges with different colors represent different association relationships. Among the 36 targets, except for Muscarinic Acetylcholine receptor M2 (CHRM2), CHRM4, Choline-phosphate cytidylyltransferase A and Sodium channel protein type 5 subunit alpha, the remaining 32 target proteins are closely related to each other and may participate in multiple pathways.
Figure 7
Target-pathway diagram of Ruyi Jinhuang powder-diabetic foot ulcer. A: A bar graph of Kyoto Encyclopedia of Genes and Genomes (KEGG) signal pathway enrichment. The left side is the pathway name, and the right side is the number of targets contained; B: The target-path network relationship diagram. The yellow circle represents 27 pathways, and the blue square represents 36 targets; C: An enriched bubble diagram of the KEGG signaling pathway. The left side is the pathway name, and the abscissa is the fold enrichment. The bubble size corresponds to the number of targets contained, and the color category corresponds to the P value; D: The correspondence diagram of pathway categories. The first column is the ownership of the pathway, the second column is the pathway, and the third column is the target.
Interaction diagram of target proteins. The edges with different colors represent different association relationships. Among the 36 targets, except for Muscarinic Acetylcholine receptor M2 (CHRM2), CHRM4, Choline-phosphate cytidylyltransferase A and Sodium channel protein type 5 subunit alpha, the remaining 32 target proteins are closely related to each other and may participate in multiple pathways.Target-pathway diagram of Ruyi Jinhuang powder-diabetic foot ulcer. A: A bar graph of Kyoto Encyclopedia of Genes and Genomes (KEGG) signal pathway enrichment. The left side is the pathway name, and the right side is the number of targets contained; B: The target-path network relationship diagram. The yellow circle represents 27 pathways, and the blue square represents 36 targets; C: An enriched bubble diagram of the KEGG signaling pathway. The left side is the pathway name, and the abscissa is the fold enrichment. The bubble size corresponds to the number of targets contained, and the color category corresponds to the P value; D: The correspondence diagram of pathway categories. The first column is the ownership of the pathway, the second column is the pathway, and the third column is the target.We established the network relationship between the three pathways of the PI3K-Akt signaling pathway, neuroactive ligand–receptor interaction and FoxO signaling pathway and the target and chemical components and analyzed the topological parameters of each node. The network diagram (Figure 8) showed that CHRM2 and RXRA are the common targets of the three pathways.
Figure 8
Critical pathway-target-component network relationship diagram of Ruyi Jinhuang powder. Pink represents the pathway, yellow represents the target and blue represents the components. It was shown that Muscarinic Acetylcholine receptor M2 and Retinoid X Receptor alpha are common targets of the three pathways. Luteolin, atractylenone and arylodone are the key components in the three pathways.
Critical pathway-target-component network relationship diagram of Ruyi Jinhuang powder. Pink represents the pathway, yellow represents the target and blue represents the components. It was shown that Muscarinic Acetylcholine receptor M2 and Retinoid X Receptor alpha are common targets of the three pathways. Luteolin, atractylenone and arylodone are the key components in the three pathways.
The prediction of anti-DFU mechanism
According to the literature, wound healing mainly involves the three processes of hemostasis and inflammation, proliferation, and remodeling. The target pathways involved are diverse and complex, as shown in Figure 9. RHP can act on leukocytes by regulating IL-6 and TNF, regulating the inflammatory response, and then participating in hemostasis to participate in inflammation and remodeling stages. RHP can also affect platelets, macrophages and fibroblasts by acting on the PI3K-Akt pathway to promote angiogenesis participation in the proliferation stage and facilitate wound healing in DFUs.
Wound healing mechanism diagram. TNF: Tumor necrosis factor; IL: Interleukin; VEGF: Vascular endothelial growth factor; NF-kB: Nuclear factor-kappaB.Further analysis of the effects of the three pathways in the human body found that if the PI3K-Akt signaling pathway is inhibited in diabetic neurons, neuronal apoptosis is increased, and diabetic neuropathic pain is induced; thus, activating this pathway can alleviate painful diabetic peripheral neuropathy. The abnormal expression or loss of genes related to the neuroactive ligand–receptor interaction pathway can cause neuropathy, such as with the downregulation of glial gene expression, which impairs nerve impulse conduction and increases the potential nerve involvement of systemic diseases. The link between the FoxO signaling pathway and diabetes is that type II diabetes mellitus causes abnormal tissue signaling due to hyperglycemia or insulin resistance. The identified important targets, such as IL-6 and TNF, participate in the abovementioned pathways and play a role in hemostasis and tissue remodeling so that the DFU wound can heal. The mechanism of action is shown in Figure 10.
The results of the three groups of experiments showed indicated that RHP has a good therapeutic effect on the healing of DFU wounds in rats. The HE staining electron microscopy results showed neutrophil infiltration and granulation tissue loosening in the model group and blank group. The fibroblasts of the RHP group had good function, the granulation tissue was mature, and nerve cells were restored. The thickness of the stratum corneum and the integrity of the epidermis were good, showing a good state of recovery, and intact lymphocytes could be observed, but there were a few mast cells. Thus, RHP promoted the healing of DFUs in rats by affecting fibroblasts and nerve cells.In this study, GC–MS and UPLC–MS were used to collect and identify the chemical components of RHP. Eighty-nine compounds were detected and analyzed, including flavonoids, terpenes, and coumarins. The components obtained by GC–MS are mainly volatile substances, and the components obtained by UPLC–MS are mostly high boiling point, nonvolatile, high molecular weight substances. The network pharmacology analysis of these compounds shows that RHP plays a role in the treatment of DFUs through multiple targets and channels. However, the manner in which specific components are combined with target proteins needs to be further explored.DFU is one of the main complications of diabetes. Current studies have found that its main causes are neuropathy, vascular disease and impaired immune function[13,14]. The specific pathogenesis is that when persistent hyperglycemia occurs, nerve cells accumulate a large number of toxic components, damage to endothelial cell function occurs, vascular tension is decreased, and nerve ischemia occurs; thus, nerve and vascular lesions develop. In addition, persistent hyperglycemia affects the inflammatory response of the wound by acting on cells and the immune system, and the wound is persistently infected and difficult to heal[15,16].The active ingredients were extracted from the ten herbs of RHP, and network pharmacological analysis of "herb-component-target pathway-disease" was performed to identify three key pathways: The PI3K-Akt signaling pathway; neuroactivated-receptor interactions; and the FoxO signaling pathway. These pathways are important targets for the treatment of DFU. Studies have found that RXRA has a negative regulatory effect on glucose-stimulated insulin secretion[17], while CHRM2 mainly plays a role in the central nervous system and peripheral nervous system and can activate endothelial cell NO lyase to relax vascular smooth muscle[18,19], which is a major cause of foot ulcers caused by diabetic vascular disease and is a key target for the treatment of vascular disease[20]. The mechanism of the PI3K-Akt signaling pathway in the treatment of DFU is reflected in the following two aspects. On the one hand, insulin initiates this signaling pathway by upregulating PI3K expression with Akt molecules and inhibits FoxO1 expression, thereby allowing normal glucose uptake by cells and improving abnormal lipid metabolism[21-23]. On the other hand, if the PI3K-Akt signaling pathway is inhibited in diabetic neurons, neuronal apoptosis is increased, and diabetic neuropathic pain is induced; thus, activation of this pathway can relieve painful diabetic peripheral neuropathy[24]. In addition, diabetic patients have difficulty healing repeated wound infections due to a prolonged wound inflammatory response due to metabolic disorders, and the PI3K-Akt signaling pathway can promote the expression of hypoxia-inducible Factor 1, an important factor in wound healing, to accelerate wound cell proliferation to promote healing[25]. The neuroligand-interreceptor interaction pathway is a collection of intracellular and extracellular related receptor ligands, and abnormal expression or deletion of related genes in this pathway can cause neuropathic lesions; for example, downregulation of Glral gene expression impairs nerve impulse conduction and increases the potential incidence of neurological diseases[26,27]. The link between the FoxO signaling pathway and diabetes is that type 2 diabetes induces skeletal muscle atrophy due to abnormal tissue signaling and protein disorders caused by hyperglycemia or insulin resistance, and insulin phosphorylates mediators of the FoxO signaling pathway and inhibits autophagy-related gene expression. Luteolin can reduce neuropathic pain by decreasing FoxO1 acetylation, regulating this pathway, and inhibiting the expression of inflammation- and pain-related genes[28]. In addition to these three pathways, the PI3K-Akt signaling pathway, neuroactivated-receptor interactions, and the FoxO signaling pathway, the network topological properties of the human T-cell leukemia virus-infection signaling pathway were also found to be prominent in the RHP network pharmacology analysis; however, there is no specific research result on how this pathway affects DFU and needs further study.
CONCLUSION
RHP, a traditional Chinese medicine formula, may play a role in the treatment of DFU through these target pathways by affecting insulin resistance, altering the nervous system and immune system, participating in inflammatory responses, and regulating cell proliferation, differentiation and apoptosis through other specific mechanisms. The key active components were luteolin, methyl palmitate, ar-turmeric, methyl linoleate, palmitic acid, demethoxycurcumin, and naringin.
ARTICLE HIGHLIGHTS
Research background
Diabetic foot ulcer (DFU) seriously affects the quality of life of patients. Traditional Chinese medicine has a unique effect in the treatment of skin ulcerative diseases. Ruyi Jinhuang powder (RHP) is one of the classic prescriptions in traditional Chinese medicine and is widely used in clinical practice.
Research motivation
Although there have been studies suggesting that RHP has a therapeutic effect on DFU, there is a lack of research that further verify the mechanism of RHP to promote wound healing.
Research objectives
The effective components of RHP were extracted and identified by chromatography-mass spectrometry, and the obtained chemical components were analyzed by network pharmacology methods to predict its therapeutic mechanism. Gas chromatography–mass spectrometry (MS) and ultra-performance liquid chromatography-MS were used to separately identify the chemical constituents.
Research methods
Sprague Dawley rats were injected with streptozotocin to establish the DFU model. Hematoxylin-eosin staining was used to observe the wound tissue under an electron microscope. Medicine Systems Pharmacology database to obtain the target information, and the molecular docking of important components and key targets was performed in Discovery Studio software. Cytoscape software was used to visualize and analyze the relationship between the chemical composition, targets and Traditional Chinese Medicine network.
Research results
RHP was shown to promote the healing of diabetic foot ulcers in rats by affecting fibroblasts and nerve cells. A total of 89 chemical components were obtained by chromatography-mass spectrometry. Network pharmacological analysis revealed that RHP was associated with 36 targets and 27 pathways in the treatment of DFU.
Research conclusions
Our results indicated that RHP may play a role in the treatment of DFU through these target pathways by affecting insulin resistance, altering the nervous system and immune system, participating in inflammatory responses and regulating cell proliferation, differentiation and apoptosis and through other specific mechanisms.
Research perspectives
We found that RHP plays a role in the treatment of diabetic foot ulcers through multiple targets and channels. However, the way in which specific components are combined with target proteins needs to be further explored.
Authors: Allan Peter Davis; Cynthia J Grondin; Robin J Johnson; Daniela Sciaky; Jolene Wiegers; Thomas C Wiegers; Carolyn J Mattingly Journal: Nucleic Acids Res Date: 2021-01-08 Impact factor: 16.971