Tao Huang1†
Ziwan Ning1†
Dongdong Hu1,2,3†Man Zhang1,4†
Ling Zhao1
Chengyuan Lin1,5
Linda L. D. Zhong1,6Zhijun Yang1Hongxi Xu2,3
Zhaoxiang Bian1,6* for MZRW Research Group- 1Lab of Brain and Gut Research, School of Chinese Medicine, Hong Kong Baptist University, Kowloon Tong, Hong Kong
- 2School of Pharmacy, Shanghai University of Traditional Chinese Medicine, Shanghai, China
- 3Engineering Research Center of Shanghai Colleges for TCM New Drug Discovery, Shanghai, China
- 4Guangzhou Research Institute of Snake Venom, Guangzhou Medical University, Guangzhou, China
- 5YMU-HKBU Joint Laboratory of Traditional Natural Medicine, Yunnan Minzu University, Kunming, China
- 6Hong Kong Chinese Medicine Clinical Study Centre, Hong Kong Baptist University, Kowloon Tong, Hong Kong
MaZiRenWan (MZRW, also known as Hemp Seed Pill) is a Chinese Herbal Medicine which has been demonstrated to safely and effectively alleviate functional constipation (FC) in a randomized, placebo-controlled clinical study with 120 subjects. However, the underlying pharmacological actions of MZRW for FC, are still largely unknown. We systematically analyzed the bioactive compounds of MZRW and mechanism-of-action biological targets through a novel approach called “focused network pharmacology.” Among the 97 compounds identified by UPLC-QTOF-MS/MS in MZRW extract, 34 were found in rat plasma, while 10 were found in rat feces. Hierarchical clustering analysis suggest that these compounds can be classified into component groups, in which compounds are highly similar to each other and most of them are from the same herb. Emodin, amygdalin, albiflorin, honokiol, and naringin were selected as representative compounds of corresponding component groups. All of them were shown to induce spontaneous contractions of rat colonic smooth muscle in vitro. Network analysis revealed that biological targets in acetylcholine-, estrogen-, prostaglandin-, cannabinoid-, and purine signaling pathways are able to explain the prokinetic effects of representative compounds and corresponding component groups. In conclusion, MZRW active components enhance colonic motility, possibly by acting on multiple targets and pathways.
Introduction
Functional constipation (FC) is a common gastrointestinal (GI) disorder characterized by slow bowel movement and defecation difficulties. It influences over 14.0% adults worldwide (Suares and Ford, 2011). Female gender, old age, and low socioeconomic status are risk factors for people to have FC (Suares and Ford, 2011). FC cause worse life quality and significant financial burden of whole society (Liem et al., 2009; Peery et al., 2012, 2015). People with mild or moderate FC can be treated with high-fiber or laxatives (Rao et al., 2015; Christodoulides et al., 2016), while patients with severe FC need special cares and aggressive therapies. Several pharmacological therapeutics have been approved for FC, including diphenyl mechanics or derivatives, anthraquinone, 5-hydroxytryptamine receptor 4 (5-HT4) agonist, chloride channel type 2 activator, guanylate cyclase C receptor agonist, apical sodium bile acid inhibitors (Nelson et al., 2016). More than 50% FC patients are not completely satisfied with current therapies (Nelson et al., 2016), and alternative therapies for FC are required.
MaZiRenWan (MZRW, also known as Hemp Seed Pill) is an herbal formula for constipation from Traditional Chinese Medicine (TCM). It was firstly recorded in a TCM classic book, Discussion of Cold-Induced Disorders (Shang Han Lun) (Jing, 1997; Mitchell et al., 1998), in about 2,000 years ago. MZRW is comprised of six herbs, including Fructus cannabis (Huo Ma Ren, HMR), Radix et rhizoma rhei (Da Huang, DH), Semen Armeniacae Amarum (Ku Xing Ren, KXR), Radix paeoniae Alba (Bai Shao, BS), Cortex magnolia officinalis (Hou Pu, HP), and Fructus aurantii immaturus (Zhi Shi, ZS) (Cheng et al., 2011).
A systematic review of the published literature showed that MZRW is the most frequently used TCM formula for constipation (Zhong et al., 2016), but there is little strict clinical evidence to prove its efficacy. In the randomized, placebo-controlled clinical study with 120 FC patients, we demonstrated that MZRW is significantly better than placebo in improvement of bowel movement during the drug treatment period, and such effect is more sustainable than placebo during the 8 weeks follow-up period (Cheng et al., 2011). We also identified 10 major compounds from MZRW in rat plasma by UPLC-MS/MS (Hu et al., 2015) to facilitate the pharmacokinetic study of MZRW in healthy volunteers (Hu et al., 2017). However, the active components and the mechanism-of-actions by which MZRW utilized to alleviate FC, are still unclear.
In this study, we investigated the pharmacological actions of MZRW for FC with a novel idea: focused network pharmacology. In the classic network pharmacology studies for TCM, hundreds of herbal compounds were retrieved from database and/or literatures, and hundreds/thousands of targets were mapped (Zeng et al., 2017; Chen et al., 2018). The network comprised of huge number of compounds and targets could explain the holistic and complex effects of TCM, but they will also result in numerous mechanistic hypotheses to be tested. To confer such defects, we only included the compounds validated by quantitative methods, then further focused on the representative compounds with chemical space analysis.
To precisely define the effective compounds, we used UPLC-QTOF-MS/MS to identify compounds in MZRW extract (after preparation), and compounds in biological samples of rats after oral administration of MZRW. Followed-by computational analysis, we found that the identified compounds can be clustered as “component groups” based on pairwise chemical similarity. Five representative compounds for the component groups were selected and their effects on rat intestinal tract contractions were tested. The relevant disease targets and signaling pathways were investigated.
Materials and Methods
Compounds Identification
Reagents and Materials
All the six herbal medicines (HMR: RM130121-08; DH: 130524; BS: RM130128-11; KXR: RM121026-04; HP: RM121026-05; ZS: RM121210-11) contained in MZRW were supported and authenticated by PuraPharm Company, Nanning department (Nanning, China). The reference standards of gallic acid, amygdalin, paeoniflorin, hesperidin, neohesperidin, naringin, magnolol, honokiol, aloe emodin, physcion, chrysophanol were purchased from Shanghai Yuan-Ye Bio-technology company (Shanghai, China). Emodin, rhein, albiflorin were purchased from the National Institute for the Control of Pharmaceutical and Biological Products (Beijing, China). The purity of all these reference standards was ≥98.0% (HPLC). HPLC grade of acetonitrile, methanol and formic acid were purchased from Merck (Darmstadt, Germany). Deionized water was purified by the Millipore water purification system (Millipore, Milford, MA, United States). All other reagents used were of analytical grade. Twelve male Sprague-Dawley rats (220–250 g) were kindly provided by the Experimental Animal Center of Chinese University of Hong Kong (Hong Kong). The rats were bred in an environmentally controlled room (22 ± 2°C, relative humidity 50 ± 20%) with a natural light-dark cycle for 7 days before the experiment carried out. The animal study was carried out in accordance with the Guideline for Animal Experimentation of Hong Kong.
Instrumentation and Analytical Conditions
Liquid Chromatographic analysis was performed on an Agilent 1290 UPLC system, consisting of a 1290 binary pump solvent management system, an 1290 TCC, and an 1290 auto-sampler. A waters Acquity BEH C18 column (100 mm × 2.1 mm, 1.7 μm) was employed and the column temperature was maintained at 40°C. The mobile phase was composed of A (0.1% formic acid in water) and B (0.1% formic acid in acetonitrile) using a gradient elution of 2–5% B at 0–2.5 min, 5–35% B at 2.5–10 min, 35–75% B at 10–14 min, 75–100% B at 14–16 min, 100% B at 16–20 min, 100–2% B at 20–20.1 min and 2% B at 20.1–24 min with a flow rate set at 0.40 mL/min. The auto-sampler was conditioned at 4°C and the injection volume was 3 μL.
Mass detections were performed using an Agilent 6540 Quadruple-Time of flight mass spectrometer (6540 Q-TOF-MS, Agilent Technologies, Santa Clara, CA, United States) equipped with an AJS electrospray ionization source (ESI). The ESI source was set in negative ionization mode. The parameters in the source were set as follows: Gas temperature, 300°C; gas flow, 8.0 L/min; nebulizer, 45 psi; sheath gas temperature, 350°C; sheath gas flow, 8.0 L/min; capillary voltage, 3000 V; nozzle voltage, 500 V; fragmentor voltage, 150 V; Oct RFV, 600 V; and skimmer voltage 65.0 V. charging, 3000 V; The operations, acquisition, and data analysis were operated under MassHunter Acquisition Software Version B.06.00 (Agilent Technologies). Accurate mass of all mass peaks were measured and recorded with a full scan mode set at a mass range from m/z 100 to 1700. The molecular masses of the precursor and product ions were accurately determined. Reference masses 966.00072500 and 112.98558700 in negative ionization mode.
To optimize signals and obtain maximal structural information, the Q-TOF MS/MS scan was also conducted. The energies for collision-induced dissociation (CID) were set at two fixed collision energies of 15 and 35 V, other Q-TOF fragmentation parameters were the same as those for Q-TOF MS, and data were also acquired over a mass range from m/z 100 to 1700.
Preparation of MZRW Herbal Extract for Administration
Six pulverized herbal ingredients, Fructus Cannabis (892.9 g), Semen Armeniacae Amarum (446.4 g), Radix et Rhizoma Rhei (446.4 g), Fructus Aurantii Immaturus (223.2 g), Cortex Magnoliae Officinalis (267.9 g) and Radix Paeoniae Alba (223.2 g) were mixed together and decocted with 8-flood volumes of water (1:8, w/v) for 2.5 h, then the decoction was filtered. The filtered decoction was concentrated by rotary evaporation under vacuum at 55°C. The residue was then freeze dried. Finally, the dried residue was dissolved in water to obtain MZRW oral solution with a concentration of 0.5 g/mL. Raw materials were stored at room temperature. MZRW preparations were stored at -20°C.
In Vivo Study
Twelve male Sprague-Dawley (SD) rats, weighing 250–300 g, were fasted for 12 h with free access to water before experiments. The rats were randomly divided into two groups. One group was for plasma collection, and the other group for feces sample collection. Blood and feces samples before administration of each rat were collected for self-control study. Single dosage of MZRW oral solution was orally administered to each rat at a dose of 1g/100 g bodyweight. Then blood samples of about 500 μL were collected at 15 min, 30 min, and 3 h after administration by cutting tail. Feces samples were collected in 12 h after administration by metabolic cages. The blood samples were then centrifuged at 5000 rcf for 5 min to obtain plasma samples, which were frozen immediately and stored at -20°C until analysis; meanwhile, the feces sample were directly stored at -20°C until analysis.
Sample Preparation
MaZiRenWan dried extract was accurately weighted and dissolved in 70% methanol-water (30 mg/mL, w/v) following with ultrasonic extraction at room temperature for 30 min. Then the extract solution was centrifuged at 13000 rcf for 10 min. The supernatant was injected for components identification of MZRW.
Reference compounds of gallic acid, amygdalin, paeoniflorin, hesperidin, neohesperidin, naringin, magnolol, honokiol, aloe emodin, physcion, chrysophanol, emodin, albiflorin, and rhein were accurately weighted and dissolved in methanol to give a stock solution of 200 μg/ml. Then mixed standard solution of this 14 compounds at 5 μg/ml were prepared for the confirmation of compounds in MZRW.
Rat plasma sample (300 μL) was added to a 1.5 mL tube (Eppendorf). Then extraction was performed by adding 700 μL methanol to precipitate protein and extract components from the plasma. The sample was vortexed for 2 min, and then centrifuged for 10 min at 13000 rcf. The supernatant was transferred into another tube and dried under a flow of nitrogen gas. The residue was re-constituted in 100 μL methanol solutions, and centrifuged (13000 rcf for 10 min). The supernatant was transferred to an auto-sampler vial, 3 μL was injected into the UPLC-MS/MS system for analysis.
Feces samples were accurately weighted and mixed with normal saline as 1g/2ml. Then ultrasonic extraction (30 min at room temperature) was conducted to extract the components in feces. The mixed solution was then centrifuged. Three microliter supernatant was injected for LC-MS/MS analysis.
Chemical Space Analysis
The Venn’s diagram was drawn by online tool (Venny 2.11) (Oliveros, 2007–2015). Fingerprint similarity between compounds was calculated with ECFP-4-like Morgan (radius = 2) (Rogers and Hahn, 2010) in RDKit2. The clustering map was made with seaborn (version 0.8.0) package in Python. The PCA analysis was performed with scikit-learn (Pedregosa et al., 2011) package in Python.
Organ Bath
Adult Sprague-Dawley male rats (250–300 g) were obtained from Guangdong Medical Laboratory Animal Centre. Rats were kept on a 12 h light/dark cycle at 23 ± 2°C with free access to food and water. Experimental protocols were approved by the Committee on the Use of Human and Animal Subjects in Teaching and Research, Hong Kong Baptist University (HASC/16-17/0331). Emodin, Naringin, Amygdalin, Albiflorin, and Honokiol were purchased from Shenzhen ChemStrong Scientific Co., Ltd. The purity of the drugs were equal or above 99%. The drugs were dissolved in DMSO and the final concentration of DMSO in solution were limited in 1‰.
The rats were euthanized with CO2, and a 4-cm long segment of distal colon was dissected out. The content of the colonic segment was flushed with Krebs solution and then cut into 4 segments, each 1 cm in length. Each segment was mounted longitudinally and the contraction of longitudinal smooth muscle was recorded. Briefly, the lower end of the colonic segment was tied with a holder at the bottom of the organ bath, while the upper end was connected to an isometric force transducer with a silk thread (Braided silk wax, US 5/0, Pearsalls Ltd., United Kingdom) for recording mechanical activity. Colonic contractions were recorded using the PowerLab system and Chart5 software (AD instrument Ltd., Australia). The segment was bathed in an oxygenated (95%O2–5% CO2) Krebs solution at 37°C. Krebs solution contained (in mM) 119 NaCl, 4.5 KCl, 1.2 MgCl2, 25 NaHCO3, 1.2 KH2PO4, 2.5 CaCl2, and 11.1 Glucose. The colonic segments were allowed to equilibrate for 1 h before the experiment started. During the equilibration period, the segment was rinsed every 20 min with Krebs solution and the basal tension was maintained (Lin et al., 2015).
The amplitude (in g) of contractions was measured and expressed as active tension (force/area, g/mm2), using the following equation: force/area = grams tension/ [gram wet wt/ (1.05 × L0)], with 1.05 as the density of smooth muscle (Bossone et al., 2001). The optimal length (Lo) was obtained using the similar method described before.
Network Analysis of Representative Compounds
Referenced Targets Analysis
For each of the representative compounds of MZRW, the corresponding biological targets were searched from BindingDB3 (Chen et al., 2001; Liu et al., 2007; Gilson et al., 2016). These reported targets were manually checked and curated by searching relevant references from PubMed. Full-text of publications were retrieved and read by one expert. If no relevant publications were founds with that target, it was removed from the referenced targets list. If relevant publications were found, the associated bioactivity data were also gathered.
Predicted Targets Analysis
The targets of MZRW representative compounds were predicted by the MOST method (Huang et al., 2017). Based on the bioactivity type, the pKi, pIC50, and pEC50 datasets were generated from the CHEMBL20 database (Gaulton et al., 2012) by using procedures described in previous work (Huang et al., 2017). The machine learning model was training with dataset from CHEMBL19 and the performance of prediction were summarized in Supplementary Table S8. For each compound, the SMILES string was used for dataset searching. The results were filtered by the following criterions: Tc ≥ 0.45 and p < 0.05.
Target-Constipation Link Analysis
For all the referenced and predicted targets, their links with constipation were validated by searching PubMed4. The following key words were used: “[target gene symbol or full name] AND constipation.” Only targets with constipation-relevant publications were remained.
Results
Identification of Compounds in MZRW Extract, and in Rat Plasma and Feces After Oral Administration of MZRW
We set out to identify compounds in the MZRW extract, as well as in plasma and feces samples of rats orally administrated with MZRW extract. In total, 97 compounds were identified by UPLC-QTOF-MS/MS in MZRW extract (Supplementary Table S1). Among them, 7 were from HMR, 48 were from DH, 3 were from KXR, 14 were from BS, 10 were from HP, and 15 were from ZS (Figure 1A and Supplementary Table S1). We also identified the compounds in rat plasma and feces samples after oral administration: 34 compounds can be found in rat plasma, and 10 can be found in rat feces (Figure 1B). Nine compounds were commonly identified in rat plasma and feces, while only 1 (naringenin) was found in feces (Figure 1B).
FIGURE 1. Identification of compounds of MZRW in extract and biological samples and selection of representative compounds. (A) Histogram of number of identified compounds by herbs. (B) Venn diagram of number of identified compounds by extract and biological samples. (C) Clustering analysis of identified compounds based on chemical similarity. Component groups were highlighted by yellow rectangle. (D) PCA analysis of identified compounds according to chemical similarity. (E) Representative compounds of each component group.
Next we analyzed the chemical spaces of the identified 97 compounds. By calculating the fingerprint similarity and analyzing the clusters, these compounds were classified into hierarchical groups (Figure 1C). In general, they are fatty acids and non-fatty acid compounds. The fatty acids, including linoleic acid, oleic acid, and sativic acid, are oil components from HMR. The non-fatty acid compounds can be further categorized into compounds with glucoside (e.g., albiflorin), and without glucoside (e.g., emodin) (Figures 1C,D). Very interestingly, these compounds were clustered into small groups based on chemical similarity. Within such group, the compounds are highly similar to each other in terms of chemical structure, and most of them are from the same herb (Supplementary Tables S2–S6). It is well-known that similar compounds are likely to have similar bioactivity profiles (Bender and Glen, 2004), the clustered compounds suggest a way to abstract the complex chemical space of MZRW. We use the concept of “component group” to define the clustered compounds. Firstly, a component group contain compounds with highly similar structures, and most of the compounds are from the same herb. Secondly, a component group can be treated as a whole in investigation of pharmacological actions. Thirdly, the biological function of a component group can be demonstrated by a representative compound, which is abundant in herb, in formula extract, and in biological samples.
With these definitions, we were able to locate five component groups, which are from DH, KXR, BS, HP, and ZS (Figure 1C), respectively (Supplementary Tables S2–S6). Five representative compounds were selected because of their high abundance in MZRW (Hu et al., 2015), commercial availability, and frequent use in pharmacology study: emodin for DH, amygdalin for KXR, albiflorin for BS, honokiol for HP, and naringin for ZS (Figure 1E). These representative compounds were used with organ bath experiments to elucidate the pharmacological actions of MZRW. The component group from HMR were excluded, because they are fatty acids which are hard to be tested in organ bath.
Effects of the Representative Compounds of MZRW on the Contractions of Rat Colonic Segments
A major part of FC patients are with symptom of slow bowel movement, which can be alleviated by MZRW treatment. Thus, we tested if the representative compounds could induce contractions of the rat colonic segments in organ bath experiments.
Emodin (1 μM to 0.3 mM) concentration-dependently increased colonic smooth muscle contractions (Figure 2A). Amygdalin (1 nM to 1 μM) concentration-dependently increased colonic smooth muscle contractions (Figure 2B). Albiflorin (1–30 μM) concentration-dependently increased colonic smooth muscle contractions while decreased the contraction at the dose range from 100 μM to 1 mM (Figure 2C). Honokiol (1–30 nM) concentration-dependently increased colonic smooth muscle contractions while decreased the contraction at the dose range from 300 nM to 100 μM (Figure 2D). Naringin (1–100 nM) concentration-dependently increased colonic smooth muscle contractions while decreased the contraction at the dose range from 1 μM to 1 mM (Figure 2E).
FIGURE 2. Effects of MZRW representative compounds on spontaneous contraction of colonic smooth muscle in rats in vitro. Effects of emodin (A), amygdalin (B), albiflorin (C), honokiol (D), and naringin (E) on spontaneous contraction of colonic smooth muscle. Data are mean ± SEM (n = 4–6). Significant difference is indicated by ∗P < 0.05 and ∗∗P < 0.01 compared with control.
These results suggested that all the five representative compounds have prokinetic effects in vitro. Next we asked, what’s the biological targets mediating such prokinetic effect?
Representative Compound-Target-Disorder Network of MZRW for FC
We used two complementary methods to investigate the constipation-relevant targets of representative compounds.
Firstly, the targets of representative compounds were retrieved from the BindingDB (see footnote 3) (Chen et al., 2001; Liu et al., 2007; Gilson et al., 2016), and curated manually, resulting in the “referenced targets” list (Supplementary Table S7). No referenced targets can be found with amygdalin and albiflorin, while 7 targets for emodin [ESR1 (Matsuda et al., 2001), ESR2 (Matsuda et al., 2001), CSNK2A1 (Sarno et al., 2003), EPHX2 (Lee et al., 2015), PTP4A3 (Han et al., 2012), LCK (Chang and Geahlen, 1992), and ELANE (Zembower et al., 1992)], 7 targets for honokiol [PTGS1 (Schuhly et al., 2009), PTGS2 (Schuhly et al., 2009), ALOX5 (Schuhly et al., 2009), CNR1 (Rempel et al., 2013), CNR2 (Rempel et al., 2013), GABAA (Taferner et al., 2011), and RXR (Kotani et al., 2010)], and 1 target for naringin [CYP19A1 (Endringer et al., 2008)] were found (Table 1).
TABLE 1. Referenced and predicted targets of MZRW representative compounds.
Secondly, a computational method, namely MOST (Huang et al., 2017), was used to predict the possible targets of representative compounds. The machine learning models of MOST were trained by datasets from CHEMBL19 database (Gaulton et al., 2012). The prediction accuracies for pIC50, pEC50, and pKi datasets are 73.5, 80.4, and 75.5%, respectively in temporal validation, where pIC50/pEC50/pKi ≥ 6 are defined as “active” (Supplementary Table S8). Then, the representative compounds were searched against the datasets from CHEMBL20 (Supplementary Table S9) and corresponding predictions were made. In general, 14 predicted targets for emodin, 10 predicted targets for amygdalin, 7 predicted targets for albiflorin, 7 predicted targets for honokiol, and 11 predicted targets for naringin were found (Table 1).
Next, we checked which biological targets (from both referenced and predicted targets) of MZRW representative compounds have direct links with the symptom of constipation. By searching with PubMed, seven targets have been reported to be associated with constipation (Table 2).
TABLE 2. Referenced and predicted targets correlated with constipation.
Finally, we generated the component group-target-disorder network of MZRW for FC (Figure 3). Albiflorin was not connected into this network since no targets of albiflorin have been reported to be linked with constipation. For the remaining four compounds, emodin, amygdalin, honokiol, and naringin, they were linked with constipation through seven targets, which can be categorized into five signaling pathways: acetylcholine (ACh), estrogen, prostaglandin, cannabinoid, and purine. All of these signaling pathways are related with constipation or gastrointestinal (GI) motility.
FIGURE 3. The drug-target-disorder network of MZRW representative compounds for constipation. The prokinetic effects of MZRW representative compounds can be explained by several targets. The compound is linked with biological target if it has reported (solid line) or predicted (dashed line) to be interacted with that target. All the biological targets have been reported to have links with constipation. Overall, these biological targets are belonging to the ACh signaling (red), estrogen signaling (orange), prostaglandin signaling (blue), cannabinoid signaling (green), and purine signaling (purple) pathways.
Discussion
To investigate the pharmacology of TCM formula is a daunting task because of the complicated constituents and mechanism-of-actions. Network pharmacology was thought to be an important step to TCM modernization and a new research paradigm for the translation of TCM into evidence-based-medicine system (Liu and Sun, 2012; Li and Zhang, 2013). To date, more than 100 articles have been published in the cross-section of Chinese Medicine and network pharmacology. However, with this classic approach, typically hundreds of compounds are mapped to hundreds/thousands of targets, which finally would generate too many hypotheses to be tested. In our work studying the pharmacology of MZRW for FC, we used an improved strategy, namely, focused network pharmacology. Instead of looking for all the herbal compounds from public database and literatures, only compounds in the extracts that can be validated by quantitative method (UPLC-QTOF-MS/MS), were included. Such improvement provide a ‘real chemical space’ for subsequent analysis. We also proposed a hierarchical representation, herb-component group-representative compound, to abstract the chemical space. Based on the hierarchical clustering analysis, MZRW were seen as combination of several component groups, while each component group was represented by a major compound: that is the reason we named it “focused.” Compared with the traditional network pharmacology approach, the focused network pharmacology approach generate feasible number of mechanistic hypotheses, which are much easier to be tested experimentally. Very recently, several groups also used similar strategy to investigate the mechanism-of-actions of TCM. For example, Gao et al. (2018) proposed a new idea, “main active compound-based network pharmacology” to investigate the anti-cancer mechanism of Compound Kushen Injection. Four main active compounds were confirmed by UPLC-MS and cell proliferation assay, and the potential targets and pathways related with anti-HCC effects were predicted by network pharmacology approach. In another study performed by Zhao and He (2018), the main active compounds of Ganoderma lucidum extract were identified by HPLC, EI-MS, and NMR, and related targets were predicted by reverse-docking. Taken together, the focused network pharmacology approach could suggest a new direction in TCM pharmacology study: from explainable science to testable science.
In TCM theory, MZRW can drain heat, unblock the bowel, promote the movement of Qi, and moisten the intestines (Cheng et al., 2011). With modern pharmacology study, it was found that administration of MZRW increase the fecal pellet number and weight in mice. MZRW also enhance the contractions of intestinal segments of rabbit and guinea pig (ref?). Our work identified MZRW active components enhancing the colonic motility and analyzed the potential targets and pathways. It is an important step to understand the mechanisms of MZRW for FC at the molecular levels. It will also promote the quality control and modernization of MZRW.
Functional constipation is a highly heterogeneous disorders: psychological and psychiatric problems, personality traits, neuroendocrine and neurological changes, have been suggested to contribute to the development of slow transit constipation (STC), a major subtype of FC (Surrenti et al., 1995; Velio and Bassotti, 1996; El-Salhy, 2003; Frattini and Nogueras, 2008; Shahid et al., 2012). We found that several component groups (DH and ZS) could improve the constipation symptoms through modulating the estrogen signaling pathway. Naringin from ZS and Emodin from DH, through acting on estrogen production (CYP19A1) and estrogen-induced transcriptional activation (ESR2), cooperatively inhibit the estrogen signaling pathway (Figure 3). It has been found that FC is about two-times prevalent in women [17.4% (95%CI: 13.4–21.8%)] than in men [9.2% (95%CI: 6.5–12.2%)] (Suares and Ford, 2011). Consistently, progesterone and estrogen are implicated in the development of constipation (Chen et al., 1995; Knowles et al., 2003; Xiao et al., 2005, 2009; Oh et al., 2013). Since none of estrogen therapy has been approved for FC, the lesson from MZRW provide the clues to develop such novel therapies.
Our work also have several limitations. Firstly, due to the technical problems, fatty acids (mainly from HMR) were excluded from organ bath test and network analysis. Traditionally, the fatty acids are thought to have lubricant effects in treating constipation. However, there are increasing evidences suggest that the fatty acids may also have pharmacological effects in increasing bowel movement. For example, ricinoleic acid, major compound from castor oil, was found to induce laxation and uterus contraction by activating the EP3 receptors (Tunaru et al., 2012). Thus, it will be interested to investigate the pharmacological actions of fatty acid compounds of MZRW in the future. Secondly, the representative compounds of MZRW were identified but the combinational effect of them, have not been investigate yet. It will be of great research interest, to test if they have synergistic, additive, or antagonistic effect in future studies. Thirdly, several signaling pathways have been predicted by the focused network pharmacology approach, but the contribution of each pathway in enhancing colonic motility have not been measured. For individual animal subject, the contribution of single pathway could be quantitatively measured with highly specific pathway inhibition during the MZRW treatment. For the FC population, the capability of targeting multiple signaling pathway may offer MZRW the advantage over single-compound-agent in treating FC patients with different pathogenic causes.
Conclusion
We developed a novel idea, namely “focused network pharmacology,” to investigate the pharmacological actions of MZRW for FC. Representative compounds of component groups in MZRW, emodin from DH, amygdalin from KXR, albiflorin from BS, honokiol from HP, and naringin from ZS induce spontaneous contractions of rat colonic smooth muscle in vitro. Biological targets of the representative compounds are within ACh-, estrogen-, prostaglandin-, cannabinoid-, and purine signaling pathways, which explain the prokinetic effects of corresponding component groups and herbs. MZRW active components enhance colonic motility, possibly by acting on multiple targets and pathways.
Author Contributions
ZB and TH designed the whole study. TH carried out the chemical space analysis and target analysis. DH and ZN performed compound identification experiments. MZ and ZN did animal testing with the help from LZ and CL. ZB and TH analyzed the data and wrote the manuscript. LLDZ, HX, and ZY made contributions to the manuscript preparation.
Funding
This study was supported by Health and Health Services Fund, Hong Kong SAR, China (project no. 09101501) and Shenzhen Science and Technology Innovation Committee (No. JCYJ20140419130444178).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors thank all the members of MZRW research group, including Shu-hai Lin, Yan-hong Li, and Chung-Wah Cheng, for their great help and support.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2018.00270/full#supplementary-material
Footnotes
- ^ http://bioinfogp.cnb.csic.es/tools/venny/
- ^ RDKit: Open-source cheminformatics: http://www.rdkit.org
- ^ https://www.bindingdb.org
- ^ https://www.ncbi.nlm.nih.gov/pubmed/
References
Aviello, G., Romano, B., and Izzo, A. A. (2008). Cannabinoids and gastrointestinal motility: animal and human studies. Eur. Rev. Med. Pharmacol. Sci. 12(Suppl. 1), 81–93.
Bassotti, G., Villanacci, V., Bellomi, A., Fante, R., Cadei, M., Vicenzi, L., et al. (2012). An assessment of enteric nervous system and estroprogestinic receptors in obstructed defecation associated with rectal intussusception. Neurogastroenterol. Motil. 24, e155–e161. doi: 10.1111/j.1365-2982.2011.01850.x
Bender, A., and Glen, R. C. (2004). Molecular similarity: a key technique in molecular informatics. Org. Biomol. Chem. 2, 3204–3218. doi: 10.1039/b409813g
Bossone, C., Hosseini, J. M., Piñeiro-Carrero, V., and Shea-Donohue, T. (2001). Alterations in spontaneous contractions in vitro after repeated inflammation of rat distal colon. Am. J. Physiol. Gastrointest. Liver Physiol. 280, G949–G957. doi: 10.1152/ajpgi.2001.280.5.G949
Chang, C. J., and Geahlen, R. L. (1992). Protein-tyrosine kinase inhibition: mechanism-based discovery of antitumor agents. J. Nat. Prod. 55, 1529–1560. doi: 10.1021/np50089a001
Chen, L., Cao, Y., Zhang, H., Lv, D., Zhao, Y., Liu, Y., et al. (2018). Network pharmacology-based strategy for predicting active ingredients and potential targets of Yangxinshi tablet for treating heart failure. J. Ethnopharmacol. (in press). doi: 10.1016/j.jep.2017.12.011
Chen, T. S., Doong, M. L., Chang, F. Y., Lee, S. D., and Wang, P. S. (1995). Effects of sex steroid hormones on gastric emptying and gastrointestinal transit in rats. Am. J. Physiol. 268, G171–G176. doi: 10.1152/ajpgi.1995.268.1.G171
Chen, X., Liu, M., and Gilson, M. K. (2001). BindingDB: a web-accessible molecular recognition database. Comb. Chem. High Throughput Screen. 4, 719–725. doi: 10.2174/1386207013330670
Cheng, C. W., Bian, Z. X., Zhu, L. X., Wu, J. C., and Sung, J. J. (2011). Efficacy of a Chinese herbal proprietary medicine (Hemp Seed Pill) for functional constipation. Am. J. Gastroenterol. 106, 120–129. doi: 10.1038/ajg.2010.305
Choijookhuu, N., Sato, Y., Nishino, T., Endo, D., Hishikawa, Y., and Koji, T. (2012). Estrogen-dependent regulation of sodium/hydrogen exchanger-3 (NHE3) expression via estrogen receptor β in proximal colon of pregnant mice. Histochem. Cell Biol. 137, 575–587. doi: 10.1007/s00418-012-0935-2
Christodoulides, S., Dimidi, E., Fragkos, K. C., Farmer, A. D., Whelan, K., and Scott, S. M. (2016). Systematic review with meta-analysis: effect of fibre supplementation on chronic idiopathic constipation in adults. Aliment. Pharmacol. Ther. 44, 103–116. doi: 10.1111/apt.13662
Cong, P., Pricolo, V., Biancani, P., and Behar, J. (2007). Abnormalities of prostaglandins and cyclooxygenase enzymes in female patients with slow-transit constipation. Gastroenterology 133, 445–453. doi: 10.1053/j.gastro.2007.05.021
Coupar, I. M., and Tran, B. L. (2002). Effects of adenosine agonists on consumptive behaviour and body temperature. J. Pharm. Pharmacol. 54, 289–294. doi: 10.1211/0022357021778330
Czajka-Oraniec, I., and Simpson, E. R. (2010). Aromatase research and its clinical significance. Endokrynol. Pol. 61, 126–134.
El-Salhy, M. (2003). Chronic idiopathic slow transit constipation: pathophysiology and management. Colorectal Dis. 5, 288–296. doi: 10.1046/j.1463-1318.2003.00498.x
Endringer, D. C., Guimaraes, K. G., Kondratyuk, T. P., Pezzuto, J. M., and Braga, F. C. (2008). Selective inhibition of aromatase by a dihydroisocoumarin from Xyris pterygoblephara. J. Nat. Prod. 71, 1082–1084. doi: 10.1021/np800098f
Fichna, J., Sibaev, A., Salaga, M., Sobczak, M., and Storr, M. (2013). The cannabinoid-1 receptor inverse agonist taranabant reduces abdominal pain and increases intestinal transit in mice. Neurogastroenterol. Motil. 25, e550–e559. doi: 10.1111/nmo.12158
Frattini, J. C., and Nogueras, J. J. (2008). Slow transit constipation: a review of a colonic functional disorder. Clin. Colon Rectal Surg. 21, 146–152. doi: 10.1055/s-2008-1075864
Gao, L., Wang, K. X., Zhou, Y. Z., Fang, J. S., Qin, X. M., and Du, G. H. (2018). Uncovering the anticancer mechanism of Compound Kushen Injection against HCC by integrating quantitative analysis, network analysis and experimental validation. Sci. Rep. 8:624. doi: 10.1038/s41598-017-18325-7
Gaulton, A., Bellis, L. J., Bento, A. P., Chambers, J., Davies, M., Hersey, A., et al. (2012). ChEMBL: a large-scale bioactivity database for drug discovery. Nucleic Acids Res. 40, D1100–D1107. doi: 10.1093/nar/gkr777
Gilson, M. K., Liu, T., Baitaluk, M., Nicola, G., Hwang, L., and Chong, J. (2016). BindingDB in 2015: a public database for medicinal chemistry, computational chemistry and systems pharmacology. Nucleic Acids Res. 44, D1045–D1053. doi: 10.1093/nar/gkv1072
Han, Y. M., Lee, S. K., Jeong, D. G., Ryu, S. E., Han, D. C., Kim, D. K., et al. (2012). Emodin inhibits migration and invasion of DLD-1 (PRL-3) cells via inhibition of PRL-3 phosphatase activity. Bioorg. Med. Chem. Lett. 22, 323–326. doi: 10.1016/j.bmcl.2011.11.008
Hu, D.-D., Han, Q.-B., Zhong, L. L.-D., Li, Y.-H., Lin, C.-Y., Ho, H.-M., et al. (2015). Simultaneous determination of ten compounds in rat plasma by UPLC-MS/MS: application in the pharmacokinetic study of Ma-Zi-Ren-Wan. J. Chromatogr. B 1000, 136–146. doi: 10.1016/j.jchromb.2015.07.003
Hu, D. D., Zhong, L. L. D., Kun, W., Mu, H. X., Lin, C. Y., Zhao, L., et al. (2017). Evaluation of the pharmacokinetics and renal excretion of Ma-Zi-Ren-Wan in health subjects. World J. Tradit. Chin. Med. 3, 8–15. doi: 10.15806/j.issn.2311-8571.2016.0050
Huang, T., Mi, H., Lin, C. Y., Zhao, L., Zhong, L. L., Liu, F. B., et al. (2017). MOST: most-similar ligand based approach to target prediction. BMC Bioinformatics 18:165. doi: 10.1186/s12859-017-1586-z
Knowles, C. H., Scott, S. M., Rayner, C., Glia, A., Lindberg, G., Kamm, M. A., et al. (2003). Idiopathic slow-transit constipation: an almost exclusively female disorder. Dis. Colon Rectum 46, 1716–1717. doi: 10.1007/BF02660783
Kobayashi, H., Li, Z., Yamataka, A., Lane, G. J., Yokota, H., Watanabe, A., et al. (2002). Acetylcholinesterase distribution and refractory constipation - a new criterion for diagnosis and management. Pediatr. Surg. Int. 18, 349–353. doi: 10.1007/s00383-002-0776-5
Kotani, H., Tanabe, H., Mizukami, H., Makishima, M., and Inoue, M. (2010). Identification of a naturally occurring rexinoid, honokiol, that activates the retinoid X receptor. J. Nat. Prod. 73, 1332–1336. doi: 10.1021/np100120c
Lee, G. Y., Kim, J. H., Choi, S. K., and Kim, Y. H. (2015). Constituents of the seeds of Cassia tora with inhibitory activity on soluble expoxide hydrolease. Bioorg. Med. Chem. Lett. 25, 5097–5101. doi: 10.1016/j.bmcl.2015.10.014
Li, S., and Zhang, B. (2013). Traditional Chinese medicine network pharmacology: theory, methodology and application. Chin. J. Nat. Med. 11, 110–120. doi: 10.3724/SP.J.1009.2013.00110
Liem, O., Harman, J., Benninga, M., Kelleher, K., Mousa, H., and Di Lorenzo, C. (2009). Health utilization and cost impact of childhood constipation in the United States. J. Pediatr. 154, 258–262. doi: 10.1016/j.jpeds.2008.07.060
Lin, C. Y., Zhang, M., Huang, T., Yang, L. L., Fu, H. B., Zhao, L., et al. (2015). Spexin enhances bowel movement through activating L-type voltage-dependent calcium channel via galanin receptor 2 in mice. Sci. Rep. 5:12095. doi: 10.1038/srep12095
Liu, T., Lin, Y., Wen, X., Jorissen, R. N., and Gilson, M. K. (2007). BindingDB: a web-accessible database of experimentally determined protein-ligand binding affinities. Nucleic Acids Res. 35, D198–D201. doi: 10.1093/nar/gkl999
Liu, Z. H., and Sun, X. B. (2012). Network pharmacology: new opportunity for the modernization of traditional Chinese medicine. Yao Xue Xue Bao 47, 696–703.
Matsuda, H., Shimoda, H., Morikawa, T., and Yoshikawa, M. (2001). Phytoestrogens from the roots of Polygonum cuspidatum (Polygonaceae): structure-requirement of hydroxyanthraquinones for estrogenic activity. Bioorg. Med. Chem. Lett. 11, 1839–1842. doi: 10.1016/S0960-894X(01)00318-3
Mitchell, C., Chang, C.-C., Ye, F., Wiseman, N., and Zhang, Z. (1998). Shang Han Lun: On Cold Damage, Translation & Commentaries. Brookline, MA: Paradigm Publications.
Nagahama, K., Matsunaga, Y., Kawachi, M., Ito, K., Tanaka, T., Hori, Y., et al. (2012). Acotiamide, a new orally active acetylcholinesterase inhibitor, stimulates gastrointestinal motor activity in conscious dogs. Neurogastroenterol. Motil. 24, 566–574, e256. doi: 10.1111/j.1365-2982.2012.01912.x
Nelson, A. D., Camilleri, M., Chirapongsathorn, S., Vijayvargiya, P., Valentin, N., Shin, A., et al. (2016). Comparison of efficacy of pharmacological treatments for chronic idiopathic constipation: a systematic review and network meta-analysis. Gut 66, 1611–1622. doi: 10.1136/gutjnl-2016-311835
Nolan, M. L., and Scott, L. J. (2013). Acotiamide: first global approval. Drugs 73, 1377–1383. doi: 10.1007/s40265-013-0100-9
Oh, J. E., Kim, Y. W., Park, S. Y., and Kim, J. Y. (2013). Estrogen rather than progesterone cause constipation in both female and male mice. Korean J. Physiol. Pharmacol. 17, 423–426. doi: 10.4196/kjpp.2013.17.5.423
Oliveros, J. C. (2007–2015). VENNY. In Interactive Tool for Comparing Lists with Venn’s Diagrams. http://bioinfogp.cnb.csic.es/tools/venny/index.html
Pedregosa, F., Varoquaux, G., Gramfort, A., Michel, V., Thirion, B., Grisel, O., et al. (2011). Scikit-learn: machine learning in python. J. Mach. Learn. Res. 12, 2825–2830.
Peery, A. F., Crockett, S. D., Barritt, A. S., Dellon, E. S., Eluri, S., Gangarosa, L. M., et al. (2015). Burden of gastrointestinal, liver, and pancreatic diseases in the United States. Gastroenterology 149, 1731.e3–1741.e3. doi: 10.1053/j.gastro.2015.08.045
Peery, A. F., Dellon, E. S., Lund, J., Crockett, S. D., Mcgowan, C. E., Bulsiewicz, W. J., et al. (2012). Burden of gastrointestinal disease in the United States: 2012 update. Gastroenterology 143, 1179.e3–1187.e3. doi: 10.1053/j.gastro.2012.08.002
Rao, S. S., Yu, S., and Fedewa, A. (2015). Systematic review: dietary fibre and FODMAP-restricted diet in the management of constipation and irritable bowel syndrome. Aliment. Pharmacol. Ther. 41, 1256–1270. doi: 10.1111/apt.13167
Rempel, V., Fuchs, A., Hinz, S., Karcz, T., Lehr, M., Koetter, U., et al. (2013). Magnolia extract, magnolol, and metabolites: activation of cannabinoid CB2 receptors and blockade of the related GPR55. ACS Med. Chem. Lett. 4, 41–45. doi: 10.1021/ml300235q
Rogers, D., and Hahn, M. (2010). Extended-connectivity fingerprints. J. Chem. Inf. Model. 50, 742–754. doi: 10.1021/ci100050t
Sarno, S., De Moliner, E., Ruzzene, M., Pagano, M. A., Battistutta, R., Bain, J., et al. (2003). Biochemical and three-dimensional-structural study of the specific inhibition of protein kinase CK2 by [5-oxo-5,6-dihydroindolo-(1,2-a)quinazolin-7-yl]acetic acid (IQA). Biochem. J. 374, 639–646. doi: 10.1042/bj20030674
Schuhly, W., Hufner, A., Pferschy-Wenzig, E. M., Prettner, E., Adams, M., Bodensieck, A., et al. (2009). Design and synthesis of ten biphenyl-neolignan derivatives and their in vitro inhibitory potency against cyclooxygenase-1/2 activity and 5-lipoxygenase-mediated LTB4-formation. Bioorg. Med. Chem. 17, 4459–4465. doi: 10.1016/j.bmc.2009.05.018
Shahid, S., Ramzan, Z., Maurer, A. H., Parkman, H. P., and Fisher, R. S. (2012). Chronic idiopathic constipation more than A simple colonic transit disorder. J. Clin. Gastroenterol. 46, 150–154. doi: 10.1097/MCG.0b013e318231fc64
Suares, N. C., and Ford, A. C. (2011). Prevalence of, and risk factors for, chronic idiopathic constipation in the community: systematic review and meta-analysis. Am. J. Gastroenterol. 106, 1582–1591; quiz 1581, 1592. doi: 10.1038/ajg.2011.164
Surrenti, E., Rath, D. M., Pemberton, J. H., and Camilleri, M. (1995). Audit of constipation in a tertiary referral gastroenterology practice. Am. J. Gastroenterol. 90, 1471–1475.
Taferner, B., Schuehly, W., Huefner, A., Baburin, I., Wiesner, K., Ecker, G. F., et al. (2011). Modulation of GABAA-receptors by honokiol and derivatives: subtype selectivity and structure-activity relationship. J. Med. Chem. 54, 5349–5361. doi: 10.1021/jm200186n
Tunaru, S., Althoff, T. F., Nusing, R. M., Diener, M., and Offermanns, S. (2012). Castor oil induces laxation and uterus contraction via ricinoleic acid activating prostaglandin EP3 receptors. Proc. Natl. Acad. Sci. U.S.A. 109, 9179–9184. doi: 10.1073/pnas.1201627109
Velio, P., and Bassotti, G. (1996). Chronic idiopathic constipation: pathophysiology and treatment. J. Clin. Gastroenterol. 22, 190–196. doi: 10.1097/00004836-199604000-00008
Xiao, Z. L., Biancani, P., and Behar, J. (2009). Effects of progesterone on motility and prostaglandin levels in the distal guinea pig colon. Am. J. Physiol. Gastrointest. Liver Physiol. 297, G886–G893. doi: 10.1152/ajpgi.00184.2009
Xiao, Z. L., Pricolo, V., Biancani, P., and Behar, J. (2005). Role of progesterone signaling in the regulation of G-protein levels in female chronic constipation. Gastroenterology 128, 667–675. doi: 10.1053/j.gastro.2004.12.001
Zembower, D. E., Kam, C. M., Powers, J. C., and Zalkow, L. H. (1992). Novel anthraquinone inhibitors of human leukocyte elastase and cathepsin G. J. Med. Chem. 35, 1597–1605. doi: 10.1021/jm00087a014
Zeng, L., Yang, K., and Ge, J. (2017). Uncovering the pharmacological mechanism of astragalus salvia compound on pregnancy-induced hypertension syndrome by a network pharmacology approach. Sci. Rep. 7:16849. doi: 10.1038/s41598-017-17139-x
Zhao, R. L., and He, Y. M. (2018). Network pharmacology analysis of the anti-cancer pharmacological mechanisms of Ganoderma lucidum extract with experimental support using Hepa1-6-bearing C57 BL/6 mice. J. Ethnopharmacol. 210, 287–295. doi: 10.1016/j.jep.2017.08.041
Keywords: MaZiRenWan, component group, functional constipation, focused network pharmacology, representative compound, Chinese Herbal Medicine
Citation: Huang T, Ning Z, Hu D, Zhang M, Zhao L, Lin C, Zhong LLD, Yang Z, Xu H and Bian Z (2018) Uncovering the Mechanisms of Chinese Herbal Medicine (MaZiRenWan) for Functional Constipation by Focused Network Pharmacology Approach. Front. Pharmacol. 9:270. doi: 10.3389/fphar.2018.00270
Received: 09 September 2017; Accepted: 09 March 2018;
Published: 26 March 2018.
Edited by:
Ruixin Zhu, Tongji University, ChinaReviewed by:
Yibin Feng, The University of Hong Kong, Hong KongYong Wang, Beijing University of Chinese Medicine, China
Copyright © 2018 Huang, Ning, Hu, Zhang, Zhao, Lin, Zhong, Yang, Xu and Bian. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhaoxiang Bian, bianzxiang@gmail.com
†These authors have contributed equally to this work.