Abstract
Introduction:
Huangqin decoction (HQD), a traditional Chinese medicine prescription, is used to treat gastrointestinal diseases, including ulcerative colitis (UC). However, systematic research on the components of HQD remains insufficient. Therefore, we aimed to perform chemical profiling, network pharmacology, and bioinformatics analyses of HQD to identify candidate constituents potentially associated with UC and to establish a quantitative method for their determination in HQD.
Methods:
Qualitative chemical profiling was performed to identify 51 compounds in HQD, and their potential targets were predicted. UC-related target genes were identified by combining results from public and Gene Expression Omnibus (GEO) databases. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed to investigate the biological processes and signaling pathways associated with UC. Moreover, protein–protein interaction (PPI) analysis was performed. Based on these results, 15 putative active constituents of HQD were selected and quantified. Molecular docking analysis was then performed to evaluate the binding interactions between these compounds and key target proteins.
Results:
A total of 947 HQD component-related, 1,868 UC-related, and 2,930 GEO database-related target genes were intersected to obtain 109 common target genes for HQD and UC. GO and KEGG enrichment analyses indicated that these targets were mainly associated with inflammatory and immune-related biological processes and signaling pathways. Among the identified targets, NOS2, AHR, MMP3, MMP9, and PRKCQ were highlighted as potential key targets. The analysis of three batches of HQD samples revealed that baicalin had the highest content. Molecular docking results indicated favorable predicted interactions between the putative active compounds and core target proteins, with several compound–target pairs exhibiting docking scores below −11.0 kcal/mol.
Discussion:
This study not only provides a comprehensive chemical profile of HQD but also presents a new approach for evaluating and managing quality based on putative active constituents. These findings may serve as a scientific basis for further pharmacological and clinical studies.
1 Introduction
Ulcerative colitis (UC) is a chronic inflammatory disorder characterized by inflammation of the colonic mucosal surface, with progression from mild to extensive inflammation and a clinical course of relapse and remission (; ). Patients with UC typically present with clinical symptoms of bloody stool and abdominal cramping and are at an increased risk of developing malignant tumors (). UC significantly affects the quality of life of patients and imposes an economic burden (). Although various chemical-based drugs and targeted therapies are available for clinical treatment, the global incidence of UC continues to rise, and these therapies are associated with potential side effects ().
Traditional Chinese medicine (TCM) has been widely used in the treatment of UC because of its unique advantages, such as high efficacy and safety. TCMs with various targets and synergistic effects owing to their complex composition have been prescribed for UC treatment (; ). Huangqin decoction (HQD), a TCM prescription in Shang Han Lun, consists of four herbal medicines: Scutellaria baicalensis Georgi, Glycyrrhiza uralensis Fischer, Paeonia lactiflora Pallas, and Zizyphus jujuba Miller var. Inermis Rehder (). HQD has been used to treat gastrointestinal disorders for thousands of years and is effective in treating UC (). A recent study investigated the clinical efficacy of HQD and western medicine combination therapy in treating pediatric patients with UC, reporting a low incidence of side effects and good therapeutic effects (). In addition, studies on the chemical constituents of HQD have investigated its pharmacological activity and quality (; ). TCM prescriptions, such as HQD, comprise various herbal medicines containing abundant chemical components; therefore, identifying their components is very important for clinical application and to ensure safety and consistency (). However, most studies have focused on verifying the efficacy of HQD in UC using animal experiments (; ; ), and systematic research on their components is insufficient.
Network pharmacology, which is based on a systems biology approach, is a core methodology in modern TCM research that systematically elucidates the pharmacological actions and safety of TCMs, enabling the transition from experience-based to evidence-based medicine (). Moreover, it effectively investigates TCM mechanisms in disease treatment and is widely used to explore disease and drug mechanisms and active ingredients (; ). Previous studies have reported the therapeutic effects of some TCM prescriptions, including HQD, on UC using network pharmacology (; ; ; ); however, most have explored the potential mechanism of treatment or identified the mechanism using animal models. Moreover, most studies aimed at identifying active compounds for the treatment of UC have relied on database-based predictions, whereas studies based on the actual chemical constituents present in the extracts are insufficient. Therefore, a more accurate and practical component analysis of the TCM formulations is required.
The analysis of TCM preparations is not only difficult from an analytical perspective because of its complexity, but also requires safety, making liquid chromatography-mass spectrometry-based analytical methods increasingly important (). Ultra-high-performance liquid chromatography/quadrupole Orbitrap mass spectrometry (UHPLC-Q-Orbitrap-MS) is an effective, rapid, and excellent analytical method for identifying chemical components in complex herbal formulations owing to its high sensitivity and specificity (). In addition, ultra-performance liquid chromatography coupled with triple-quadrupole tandem mass spectrometry (UPLC-TQ-MS/MS) can simultaneously detect various components within a short analysis time, providing systematic information on compounds in complex matrices (). Therefore, in this study, we conducted a qualitative analysis of HQD extracts using UHPLC-Q-Orbitrap-MS, followed by network pharmacology analysis based on the identified compounds to select putative active constituents against UC treatment, which were then quantitatively analyzed using UPLC-TQ-MS/MS.
Bioinformatic analysis is used to elucidate and explore the mechanisms of various diseases and may help discover new biomarkers for improving the diagnosis and prognosis of UC (). Additionally, collecting gene expression data using microarray technology and identifying differentially expressed genes (DEGs) through bioinformatics analysis can provide important insights into the pathogenesis of UC and contribute to the derivation of stable and reliable results (; ). Therefore, in this study, we combined the results from multiple disease databases with those obtained from DEGs derived from the Gene Expression Omnibus (GEO) database to identify key targets for UC.
Moreover, it is important to develop a more comprehensive quality control (QC) method to ensure the quality and safety of herbal medicines prescribed for disease treatment (). However, systematic and comprehensive studies on the chemical components of HQD are lacking. Therefore, in this study, we identified UC disease targets and selected putative active constituents by combining network pharmacology and bioinformatics analysis based on the chemical profiling of HQD. In addition, we quantified the putative active constituents using an optimized UPLC-TQ-MS/MS analysis method, through which we established an integrated system that enables the quality evaluation and management of HQD. Furthermore, molecular docking analysis was performed to evaluate the potential interactions between the putative active constituents and UC-related target proteins and to provide computational support for the predicted compound–target relationships. In the future, this study is expected to provide a reliable QC method for HQD, verify its therapeutic effect on UC, and help elucidate its therapeutic mechanism.
2 Materials and methods
2.1 Chemicals and sample preparation
2.1.1 Chemicals and reagents
All standards used for the qualitative and quantitative analyses were purchased from ChemFaces (Wuhan, China) and TargetMol (Wellesley Hill, MA, United States). Warfarin (internal standard (IS)) was purchased from Sigma-Aldrich (St. Louis, MO, United States), and the purity of all standards was >98%. Methanol, acetonitrile, and formic acid were all MS grade and were purchased from Thermo Fisher Scientific (Loughborough, United Kingdom).
2.1.2 Preparation of HQD extract
Scutellaria baicalensis (66.5 g), Z. jujuba (66.5 g), P. lactiflora (50.0 g), and G. uralensis (50.0 g) were weighed and mixed. Thereafter, 10 times the total weight of water was added, refluxed at 100 °C for 3 h, and filtered. The filtrate was concentrated under reduced pressure and freeze-dried to obtain the HQD extract (yield: 38.79%). The four herbal medicines in HQD, S. baicalensis, Z. jujuba, P. lactiflora, and G. uralensis, were purchased from Kwangmyungdang Pharmaceutical (Ulsan, Republic of Korea), and voucher specimens for each were deposited in the KM Convergence Research Division, Korea Institute of Oriental Medicine (Specimen No. FHM-K-02, FHM-K-09–11).
2.1.3 Sample and standard solution
Standard stock solutions of 15 compounds were accurately weighed and individually dissolved in methanol. Each stock solution was mixed and diluted with methanol to obtain standard solutions. A mixed standard solution consisting of 15 compounds was diluted with methanol to create a working solution with an appropriate concentration for use as a calibration curve and QC sample. QC samples were prepared at low, medium, and high concentrations for the recovery efficiency and method validation. HQD-1, HQD-2, and HQD-3 represent three independently prepared sample solutions obtained from the same HQD extract and were analyzed to evaluate analytical reproducibility. The sample stock solution was prepared by dissolving 20 mg of HQD extract powder in 1 mL of methanol, followed by sonication for 30 min and centrifugation at 12,500 rpm for 15 min. Prior to UPLC-TQ-MS/MS analysis, the stock solution was diluted 20,000-fold with methanol to obtain a final analytical concentration of 0.001 mg/mL. The same dilution factor was applied to all analytes, and all analyte responses were confirmed to fall within the validated calibration ranges. Warfarin was selected as the IS because it exhibited stable chromatographic retention, reproducible multiple reaction monitoring (MRM) responses, efficient ionization in both positive and negative ionization modes, and no detectable interference with the target analytes under the optimized UPLC-TQ-MS/MS conditions. The internal standard was added at a constant concentration (5 ng/mL) to the standard and sample solutions used for analysis.
2.2 Chemical profiling and quantitative analysis
2.2.1 UHPLC-Q-Orbitrap-MS
The constituent compounds of HQD were identified using a Dionex UltiMate 3000 system equipped with a Thermo Q-Exactive mass spectrometer (Thermo Fisher Scientific, Waltham, MA, United States). All analytes were separated on an Acquity BEH C18 column (100 × 2.1 mm, 1.7 µm, Waters, Milford, MA, United States), and the column temperature was maintained at 40 °C. The mobile phase consisted of A (0.1% formic acid in water) and B (acetonitrile), with a gradient elution of 3% B for 0–1 min, 3%–15% B for 1–2 min, 15%–50% B for 2–13 min, 50%–100% B for 13–20 min, 100% B for 20–23 min, 100%–3% B for 23–23.5 min, and 3% B for 23.5–27.5 min. The mobile phase was set at a flow rate of 0.25 mL/min, and the injection volume of all analytes was 3.0 µL. Mass data were acquired in both positive and negative ion modes using electrospray ionization in the full MS-ddMS2 scan mode. The full scan range was 100–1,500 m/z, and the scan resolution was set to 70,000 for MS and 17,500 for MS/MS. The normalized collision energy required to obtain fragment ions was set to 25 eV. Data acquisition and processing were performed using Xcalibur v.3.0 and Tracefinder v.3.2 software (Thermo Fisher Scientific, Waltham, MA, United States).
2.2.2 UPLC-TQ-MS/MS
An Agilent 1290 Infinity II system coupled to an Agilent 6495C triple-quadrupole mass spectrometer with a jet-stream electrospray ionization source (Agilent Technologies, Santa Clara, CA, United States) was used for quantitative analysis. The separation conditions for chromatographic analysis were the same as those used for the qualitative analysis. Fifteen compounds were detected by MRM in both positive and negative ion modes. The MRM parameters for each compound were optimized, and the selected values are summarized in Table 1. Agilent MassHunter software (v.10.1) was used for the acquisition and processing of all experimental data.
TABLE 1
| No. | Compound name | Rt (min) | Ion mode | MRM transitions (m/z) | Collision energy (V) | Cell accelerator (V) | Fragmentor (V) |
|---|---|---|---|---|---|---|---|
| Q1 → Q3 | |||||||
| 1 | Magnoflorine | 4.63 | Positive | 342.2 → 58.2 | 30 | 4 | 166 |
| 2 | Albiflorin | 4.94 | Positive | 481.2 → 105.0 | 22 | 4 | 166 |
| 3 | Paeoniflorin | 5.23 | Negative | 525.2 → 449.1 | 10 | 4 | 166 |
| 4 | Liquiritin apioside | 5.67 | Negative | 549.2 → 255.1 | 34 | 4 | 166 |
| 5 | Scutellarin | 5.78 | Positive | 463.0 → 287.0 | 22 | 4 | 166 |
| 6 | Acteoside | 5.80 | Negative | 623.2 → 161.0 | 40 | 4 | 166 |
| 7 | Liquiritin | 5.81 | Negative | 417.2 → 255.0 | 18 | 4 | 166 |
| 8 | Baicalin | 7.79 | Positive | 447.1 → 271.0 | 22 | 4 | 166 |
| 9 | Ononin | 7.89 | Positive | 431.1 → 269.1 | 18 | 4 | 166 |
| 10 | Wogonoside | 9.32 | Positive | 461.1 → 285.1 | 22 | 4 | 166 |
| 11 | Benzoylpaeoniflorin | 9.77 | Negative | 629.1 → 553.2 | 10 | 4 | 166 |
| 12 | Baicalein | 10.89 | Positive | 271.0 → 123.0 | 34 | 4 | 166 |
| 13 | Glycyrrhizic acid | 12.01 | Negative | 821.4 → 351.0 | 40 | 4 | 166 |
| 14 | Wogonin | 13.02 | Positive | 285.1 → 270.0 | 26 | 4 | 166 |
| 15 | Oroxylin A | 13.58 | Positive | 285.0 → 270.0 | 26 | 4 | 166 |
| IS1 | Warfarin | 14.00 | Negative | 307.0 → 250.0 | 22 | 4 | 166 |
| IS2 | Warfarin | 14.00 | Positive | 309.0 → 163.0 | 14 | 4 | 166 |
Optimized MRM parameters for analytes and IS in UPLC-TQ-MS/MS analysis.
MRM, multiple reaction monitoring; UPLC-TQ-MS/MS, ultra-performance liquid chromatography coupled with triple-quadrupole tandem mass spectrometry; IS, internal standards.
2.2.3 Method validation
The quantitative method was validated for the selected analytes in terms of linearity, lower limit of quantification (LLOQ), accuracy, precision, and recovery. The validation strategy followed a fit-for-purpose approach for comparative quantification of major constituents in HQD samples. System suitability was confirmed by five consecutive injections of the mixed standard solution before sample analysis to ensure stable chromatographic and mass spectrometric performance. Calibration curves were constructed using multiple concentration levels spanning the expected analytical ranges of the target compounds. The mixed stock solution containing 15 analytes was serially diluted stepwise with methanol to prepare calibration working solutions at eight concentration levels, each of which was analyzed five times. Calibration curves were constructed using the peak area ratio of each analyte to the IS over the validated concentration ranges. Linearity was evaluated using weighted (1/x2) linear regression analysis to generate the calibration equations. The lowest measured concentration that allowed accuracy and precision was selected as the lower limit of quantification. To evaluate the precision and accuracy of the analytical method, QC samples were analyzed on one and three consecutive days. The relative standard deviation (RSD%) and accuracy (%) of the QC samples were used to evaluate the intra- and inter-day precision and accuracy of the method. The acceptance criteria were set as RSD < 15% for intra- and inter-day precision, with RSD < 20% allowed at the LLOQ, and accuracy within ±15% of the nominal concentration (). Precision and accuracy at the LLOQ were evaluated together with the QC samples. Recovery tests were performed by spiking HQD samples with 15 analytes at low, medium, and high concentrations, and the added and detected amounts of each analyte were compared to calculate recovery. Blank samples were analyzed before and during sample analysis to monitor carryover, and no significant carryover was observed.
2.3 Network pharmacology analysis
2.3.1 Prediction of target genes for chemical components
For the chemical components of the HQD identified in the qualitative analysis, SMILES structures were confirmed in PubChem (https://PubChem.ncbi.nlm.nih.gov/), and potential targets corresponding to each component were predicted using the Swiss Target Prediction (http://www.swisstargetprediction.ch/, accessed on 25 November 2025) database. Compounds that were not available in SwissTargetPrediction were further searched in the STITCH database (http://stitch.embl.de/, accessed on 25 November 2025). The species was restricted to Homo sapiens, and all predicted targets, including those with a prediction probability of zero, were retained for subsequent analysis. Target names were standardized to official gene symbols, and duplicate targets were removed prior to network construction. The complete compound–target association table is provided in Supplementary Table S1.
2.3.2 Acquisition of UC-related target genes
UC-related target genes were searched on four databases using the keyword “ulcerative colitis”: OMIM database (https://omim.org/, accessed on 6 August 2025), GeneCards database (https://www.genecards.org/, accessed on 6 August 2025), DisGeNET database (https://www.disgenet.org/, accessed on 6 August 2025), and TTD database (https://idrblab.net/ttd/, accessed on 6 August 2025). In GeneCards, only genes with a relevance score ≥ 5 were retained. The UC-related target genes from the databases were merged, and duplicates were removed. Data on the remaining target genes were collected.
2.3.3 Data acquisition and analysis of DEGs
Gene expression profiles associated with UC were obtained from the GEO, a Public Functional Genomics Data Repository. Two GEO datasets (GSE92415 and GSE107499) were selected, and gene expression profiles in tissues of patients with UC and healthy controls were obtained. The GSE92415 dataset (GPL13158 platform) contained 162 and 21 colonic mucosal biopsy samples from patients with UC and healthy controls, respectively. The GSE107499 dataset (GPL15207 platform) contained 75 inflammatory colonic tissues from patients with active UC and 44 normal colonic tissues. DEGs between UC and normal control samples were identified using the GEO2R online tool. Adjusted P-values were calculated using the Benjamini–Hochberg false discovery rate correction implemented in GEO2R. Statistical significance was defined as |Log2FC| > 0.5 and an adjusted P-value < 0.05, and the DEGs were visualized using volcano plots. Probe identifiers were mapped to gene symbols according to the corresponding GEO platform annotation files. When multiple probes corresponded to the same gene symbol, duplicate entries were removed prior to subsequent analyses. DEGs common to both datasets were identified using the online tool VENNY 2.1 based on gene-symbol intersections, and only genes showing consistent directions of differential expression in both datasets were retained as UC-related target genes.
2.3.4 Gene Ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) pathway analysis
The target genes of the compounds and diseases, and the results of the GEO analysis, were combined, and duplicated target genes were derived using Venn diagrams. The selected target genes were subjected to GO and KEGG enrichment analysis using the online tool Enrichr (https://amp.pharm.mssm.edu/Enrichr/), and the top 20 items were selected based on the P-value for each item. GO and KEGG pathway enrichment analysis was performed using Enrichr’s GO Biological Process 2025 and KEGG 2021 Human libraries (accessed on 17 March 2026).
2.3.5 Protein–protein interaction (PPI) network analysis
Based on the GO and KEGG enrichment results, 12 UC-related terms associated with inflammation, immune regulation, and signaling pathways were selected. A total of 61 genes associated with these selected terms were used to construct the PPI network using the Search Tool for the Retrieval of Interacting Genes (STRING) database (http://string-db.org/). Homo sapiens was selected as the species, and the minimum required interaction score was set to 0.400 (medium confidence). Active interaction sources included textmining, experiments, databases, co-expression, neighborhood, gene fusion, and co-occurrence. No additional interactors beyond the query proteins were included, and disconnected nodes were retained in the network. Subsequently, the interaction network between the selected genes and HQD components was constructed and visualized using Cytoscape 3.10.2, and putative active constituents involved in the treatment of UC were selected.
2.4 Molecular docking
We evaluated binding interactions between putative active compounds and target proteins through molecular docking analysis. Target proteins associated with the 15 quantitatively analyzed compounds were collected, and only proteins with available three-dimensional structures in the Protein Data Bank (PDB) were selected for molecular docking analysis. The 3D structures of the target proteins were obtained from the RCSB Protein Data Bank (http://www.rcsb.org/), and the 3D structures of the putative active compounds were obtained from the PubChem Compound Database (https://pubchem.ncbi.nlm.nih.gov/). Detailed information on the protein structures used for docking is provided in Supplementary Table S2. Protein structures were prepared by removing water molecules and adding hydrogen, whereas ligand structures underwent an energy minimization process prior to docking. Molecular docking was performed using AutoDock Vina implemented in PyRx (version 0.8). Blind docking was performed using grid boxes covering the entire protein structures, and the binding energies generated by the interactions were calculated. The docking results were evaluated using the lowest binding-energy pose generated by AutoDock Vina. The docking results of the putative active compounds and target proteins were visualized as heatmaps, and for structures with some low binding energies, the docking results were visualized using PyMOL 3.1.8 software.
3 Results
3.1 Identification of chemical components and target genes in HQD
The chemical components of the HQD extract were identified using UHPLC-Q-Orbitrap-MS, and the base peak chromatograms analyzed in positive and negative ion modes are shown in Figure 1. A total of 51 compounds were identified as HQD components, and the information related to each compound is summarized in Table 2. Among them, 35 compounds were unambiguously identified by comparison with the corresponding reference compounds, and 16 compounds were tentatively identified by referring to the literature. The identified compounds included 31 flavonoids, four monoterpenoids, four alkaloids, four triterpenoids, two phenolic compounds, two phenylpropanoids, one saponin, and three other compounds. In addition, based on previously reported studies, the identified compounds were assigned to the constituent herbal medicines of HQD (; ). Among the 51 identified compounds, 15 were associated with S. baicalensis, eight with Z. jujuba, eight with P. lactiflora, and 20 with G. uralensis. A total of 947 compound-related targets were identified from 51 compounds, and detailed data on compound-target associations are provided in Supplementary Table S1.
FIGURE 1
TABLE 2
| No. | Rt (min) | Identification | Formula | Adduct | Expected (m/z) | Measured (m/z) | Error (ppm) | MS2 fragment (m/z) | Source |
|---|---|---|---|---|---|---|---|---|---|
| 1 | 1.24 | Maltose* | C12H22O11 | [M − H]− | 341.1089 | 341.1098 | 2.4676 | 89 | ZJ |
| 2 | 1.35 | Proline* | C5H9NO2 | [M + H]+ | 116.0706 | 116.0712 | 4.7987 | 70 | GU |
| 3 | 1.73 | Citric acid | C6H8O7 | [M − H]− | 191.0197 | 191.0195 | −1.4399 | 111 | ZJ |
| 4 | 1.77 | 8-Debenzoylpaeoniflorin | C16H24O10 | [M + HCO2]− | 421.1352 | 421.1362 | 2.5656 | 345, 165 | PL |
| 5 | 1.98 | Adenosine* | C10H13N5O4 | [M + H]+ | 268.1040 | 268.1046 | 2.0654 | 136 | ZJ, GU, PL |
| 6 | 2.37 | Gallic acid* | C7H6O5 | [M − H]− | 169.0142 | 169.0137 | −3.3303 | 169, 125 | ZJ, PL |
| 7 | 3.61 | Phenylalanine* | C9H11NO2 | [M + H]+ | 166.0863 | 166.0866 | 2.1006 | 120 | GU |
| 8 | 5.12 | Oxypaeoniflorin | C23H28O12 | [M − H]− | 495.1508 | 495.1518 | 1.9467 | 137 | PL |
| 9 | 5.17 | Catechin* | C15H14O6 | [M − H]− | 289.0718 | 289.0725 | 2.6993 | 289, 245 | ZJ, PL |
| 10 | 5.30 | Glucoliquiritin apioside | C32H40O18 | [M − H]− | 711.2142 | 711.2158 | 2.2923 | 549, 255 | GU |
| 11 | 5.37 | Licoagroside B | C18H24O12 | [M + H]+ | 433.1340 | 433.1350 | 2.1454 | 127 | GU |
| 12 | 5.49 | Vicenin-2* | C27H30O15 | [M − H]− | 593.1512 | 593.1530 | 3.0770 | 593, 473, 353 | GU |
| 13 | 5.63 | Epicatechin* | C15H14O6 | [M − H]− | 289.0718 | 289.0728 | 3.5438 | 289, 245 | ZJ |
| 14 | 5.79 | Magnoflorine* | C20H24NO4 | [M]+ | 342.1700 | 342.1707 | 1.9706 | 342, 297 | ZJ |
| 15 | 5.90 | Schaftoside* | C26H28O14 | [M − H]− | 563.1406 | 563.1422 | 2.7007 | 563, 353 | GU, PL, SB |
| 16 | 6.02 | Albiflorin* | C23H28O11 | [M + H]+ | 481.1704 | 481.1714 | 2.0944 | 105 | PL |
| 17 | 6.22 | Ganhuangemin | C15H12O7 | [M − H]− | 303.0510 | 303.0519 | 2.8044 | 177, 125 | SB |
| 18 | 6.31 | Paeoniflorin* | C23H28O11 | [M + HCO2]− | 525.1614 | 525.1627 | 2.5892 | 121 | PL |
| 19 | 6.34 | Chavicol | C9H10O | [M + H]+ | 135.0804 | 135.0807 | 1.9524 | 105, 79 | ZJ |
| 20 | 6.49 | Rutin* | C27H30O16 | [M − H]− | 609.1461 | 609.1475 | 2.3180 | 300, 223, 131 | ZJ |
| 21 | 6.72 | Liquiritin apioside* | C26H30O13 | [M − H]− | 549.1614 | 549.1624 | 1.9203 | 549, 255 | GU |
| 22 | 6.84 | Scutellarin* | C21H18O12 | [M + H]+ | 463.0871 | 463.0877 | 1.2460 | 287 | SB |
| 23 | 6.86 | Liquiritin* | C21H22O9 | [M − H]− | 417.1191 | 417.1200 | 2.1220 | 255, 135 | GU |
| 24 | 6.95 | Acteoside* | C29H36O15 | [M − H]− | 623.1981 | 623.1995 | 1.9250 | 623, 461, 161 | SB |
| 25 | 7.31 | Isoacteoside* | C29H36O15 | [M − H]− | 623.1981 | 623.1992 | 1.7310 | 623, 461, 161 | SB |
| 26 | 7.81 | Sophoraflavone B | C21H20O9 | [M − H]− | 415.1035 | 415.1044 | 2.1925 | 295 | GU |
| 27 | 8.07 | Iridin* | C24H26O13 | [M + H]+ | 523.1446 | 523.1459 | 2.3971 | 361 | SB |
| 28 | 8.32 | Isoliquiritin apioside | C26H30O13 | [M − H]− | 549.1614 | 549.1630 | 2.9206 | 255 | GU |
| 29 | 8.53 | Scutellarein* | C15H10O6 | [M + H]+ | 287.0550 | 287.0554 | 1.3217 | 287 | SB |
| 30 | 8.62 | Isoliquiritin* | C21H22O9 | [M + H]+ | 419.1337 | 419.1341 | 1.1088 | 257 | GU |
| 31 | 8.75 | Viscidulin III | C17H14O8 | [M + H]+ | 347.0761 | 347.0765 | 1.1470 | 347 | SB |
| 32 | 8.84 | Baicalin* | C21H18O11 | [M + H]+ | 447.0922 | 447.0927 | 1.1684 | 271 | SB |
| 33 | 8.94 | Ononin* | C22H22O9 | [M + H]+ | 431.1337 | 431.1343 | 1.5735 | 269 | GU |
| 34 | 9.23 | Liquiritigenin* | C15H12O4 | [M − H]− | 255.0663 | 255.0668 | 2.2298 | 255, 135, 119 | GU |
| 35 | 9.78 | Licorice glycoside E | C35H35NO14 | [M − H]− | 692.1985 | 692.2003 | 2.6544 | 549, 255, 160 | GU |
| 36 | 10.36 | Wogonoside* | C22H20O11 | [M + H]+ | 461.1078 | 461.1084 | 1.2773 | 285 | SB |
| 37 | 10.84 | Benzoylpaeoniflorin* | C30H32O12 | [M + HCO2]− | 629.1876 | 629.1891 | 2.3950 | 121 | PL |
| 38 | 11.46 | Glabrolide | C30H44O4 | [M + H]+ | 469.3312 | 469.3321 | 1.8160 | 469, 175 | GU |
| 39 | 11.89 | Baicalein* | C15H10O5 | [M − H]− | 269.0456 | 269.0464 | 3.1100 | 269 | SB |
| 40 | 12.03 | Isoliquiritigenin* | C15H12O4 | [M − H]− | 255.0663 | 255.0668 | 2.2298 | 255, 135, 119 | GU |
| 41 | 12.47 | Formononetin* | C16H12O4 | [M + H]+ | 269.0808 | 269.0813 | 1.7052 | 269 | GU |
| 42 | 13.10 | Glycyrrhizic acid* | C42H62O16 | [M − H]− | 821.3965 | 821.3988 | 2.7925 | 351, 193, 113 | GU |
| 43 | 13.12 | Zizyberenalic acid | C30H44O3 | [M + H]+ | 453.3363 | 453.3370 | 1.5797 | 453 | ZJ |
| 44 | 13.70 | Amentoflavone | C30H18O10 | [M − H]− | 537.0827 | 537.0842 | 2.6951 | 537, 391, 245 | GU |
| 45 | 13.84 | Licorice saponin B2 | C42H64O15 | [M − H]− | 807.4172 | 807.4192 | 2.4144 | 351 | GU |
| 46 | 13.91 | Wogonin* | C16H12O5 | [M + H]+ | 285.0758 | 285.0761 | 1.2658 | 285, 270 | SB |
| 47 | 14.09 | Chrysin* | C15H10O4 | [M + H]+ | 255.0652 | 255.0655 | 1.1790 | 255 | SB |
| 48 | 14.37 | Skullcapflavone II | C19H18O8 | [M + H]+ | 375.1074 | 375.1079 | 1.3347 | 375 | SB |
| 49 | 14.46 | Oroxylin A* | C16H12O5 | [M + H]+ | 285.0758 | 285.0761 | 1.1587 | 285, 270 | SB |
| 50 | 16.97 | Glabridin* | C20H20O4 | [M − H]− | 323.1289 | 323.1298 | 2.9146 | 323, 201, 135 | GU |
| 51 | 18.20 | Hederagenin* | C30H48O4 | [M − H]− | 471.3480 | 471.3491 | 2.2914 | 471 | PL |
Identification of chemical constituents in HQD using UHPLC-Q-Orbitrap-MS.
Compared with the reference standard. ZJ: Zizyphus jujuba Miller var. Inermis Rehder; GU: glycyrrhiza uralensis fischer; PL: paeonia lactiflora pallas; SB: scutellaria baicalensis georgi.
HQD, huangqin decoction; UHPLC-Q-Orbitrap-MS, ultra-high-performance liquid chromatography/quadrupole Orbitrap mass spectrometry.
3.2 Network pharmacology analysis
3.2.1 Prediction of UC-related target genes
UC-related target genes were searched for in four databases: OMIM with 186 genes, GeneCards with 841 genes, DisGeNET with 1,458 genes, and TTD with 46 genes. By merging the results from each database and removing duplicates, 1,868 UC-related target genes were identified. DEGs were identified between patients with UC and healthy controls in two datasets from the GEO database. From the GSE92415 and GSE107499 datasets, 4,172 and 4,441 genes were extracted, respectively, and the DEGs analysis results for each dataset were plotted as a volcano map (Figure 2A). GSE92415 identified 1,670 upregulated and 2,502 downregulated genes (Supplementary Table S3), whereas GSE107499 identified 1,914 upregulated and 2,527 downregulated genes (Supplementary Table S4). A total of 2,930 common DEGs showing consistent directions of differential expression across the two GEO datasets were identified (Figure 2B). Finally, 947 HQD component-related targets, 1,868 UC-related targets, and 2,930 GEO database-related targets were intersected to obtain 109 common HQD and UC target genes (Figure 2C; Supplementary Table S5).
FIGURE 2
3.2.2 GO and KEGG pathway analysis
GO enrichment analyses were performed on 109 duplicate target genes, and the top 20 BP terms were extracted based on statistical significance. Among them, we selected eight representative terms related to inflammation and signaling pathways, including inflammatory response, peptidyl-tyrosine phosphorylation, positive regulation of cell migration, positive regulation of the ERK1/ERK2 cascade, positive regulation of the MAPK cascade, regulation of the ERK1/ERK2 cascade, regulation of inflammatory response, and regulation of the MAPK cascade (Figure 3A). KEGG pathway analysis demonstrated that target genes appear significantly abundant in various biological pathways. Among these, pathways previously reported to be associated with inflammation, immune regulation, and the pathogenesis of ulcerative colitis were selected for focused analysis (; ). These include HIF-1, IL-17, and PI3K-Akt signaling pathways, as well as Th17 cell differentiation. (Figure 3B). The complete GO and KEGG enrichment results generated by Enrichr are provided in Supplementary Tables S6, S7. A target–pathway network was constructed to more clearly demonstrate the relationship between target genes and the GO and KEGG pathways (Figure 3C). As shown in Figure 3C, several target genes are simultaneously associated with various pathways, including HIF-1, IL-17, and PI3K–Akt signaling pathways, as well as major biological processes such as inflammatory responses and MAPK/ERK signaling.
FIGURE 3
3.2.3 PPI network analysis and putative active constituents prediction
To perform a detailed analysis of overlapping target genes between HQD and UC, we first selected 12 UC-related terms from the top 20 terms identified in the GO and KEGG pathway analysis. We analyzed the PPI network of 61 genes (Supplementary Table S8) associated with the selected terms using the STRING database. Subsequently, we constructed the network relationships between the selected genes and HQD components using Cytoscape software. As shown in Figure 4, the network consists of 128 nodes and 1,309 edges. Purple square nodes represent the four herbal medicines constituting HQD, blue circle nodes represent the 51 compounds identified in HQD, yellow circle nodes represent 61 related target genes, and orange diamond nodes represent the 12 GO and KEGG pathways. Edges indicate interactions between compounds and target genes, or associations between target genes and pathways. The results showed that there is a complex correlation between compounds and targets, and the darker the color of the circle-shaped node representing the compound, the more interactions there are between the compounds and the targets. For example, albiflorin, baicalein, and magnoflorin were associated with 13 targets, acteoside and wogonin with 12 targets, and oroxylin A with 11 targets. These results suggested that specific compounds are associated with multiple targets and are likely to be the main putative active constituents of HQD involved in the treatment of UC. Based on the network analysis results, 15 compounds (magnoflorine, albiflorin, paeoniflorin, liquiritin apioside, scutellarin, acteoside, liquiritin, baicalin, ononin, wogonoside, benzoylpaeoniflorin, baicalein, glycyrrhizic acid, wogonin, and oroxylin A) were selected as putative active constituents of the HQD among the compounds linked to more than seven targets based on the qualitative analysis results.
FIGURE 4
3.3 Quantitative analysis of putative active constituents of HQD
3.3.1 Quantitative analysis
To analyze the putative active constituents of HQD, 15 compounds were simultaneously quantified using UPLC-TQ-MS/MS. All compounds were detected within approximately 15 min. The total chromatographic run time was approximately 27.5 min, including the column washing and re-equilibration steps described in the gradient program. The detected compounds were distributed in each herbal medicine constituting the HQD as follows: scutellarin, acteoside, baicalin, wogonoside, baicalein, wogonin, and oroxylin A from S. baicalensis; magnoflorine from Z. jujuba; albiflorin, paeoniflorin, and benzoylpaeoniflorin from P. lactiflora; and liquiritin apioside, liquiritin, ononin, and glycyrrhizic acid from G. uralensis. The compounds identified in this study have also been reported as major constituents of HQD in previous chemical profiling studies (; ). Dynamic MRM was performed for accurate and sensitive quantitative analysis. The precursor ions, product ions, and collision energy were optimized in the MRM mode by considering the characteristic parameters of each compound, and the detailed analysis conditions are summarized in Table 1. All compounds were analyzed individually in positive and negative ion modes according to detection sensitivity, with magnoflorine, albiflorin, scutellarin, baicalin, ononin, wogonoside, baicalein, wogonin, and oroxylin A showing higher detection intensities in the positive ion mode than in the negative ion mode. The remaining compounds were selected in negative ion mode, which showed higher detection intensities than those in positive ion mode. Representative MRM chromatograms of the reference standard solution and HQD extract used for the quantitative analysis of the 15 compounds are shown in Figure 5.
FIGURE 5
3.3.2 Method validation
To validate the UPLC-TQ-MS/MS method for analyte quantification, linearity, precision, accuracy, and recovery were investigated. A calibration curve was generated using the peak area (Y) and concentration (X) of each analyte and the IS. The linear range was selected based on the concentration of each analyte in the sample, and good linearity was observed, with a correlation coefficient (R2) ranging from 0.9990 to 0.9994. The lower limit of quantification of 15 analytes was in the range of 0.02–1.56 ng/mL, and the detailed results are presented in Table 3. Precision was evaluated using the RSD; the RSD of the intra-day precision was in the range of 0.27%–6.28%, and the RSD of the inter-day precision was in the range of 0.25%–8.09%. Additionally, the accuracy evaluation results were in the range of 90.14%–109.79% for intra-day and 91.17%–104.54% for inter-day. The average recovery rate of 15 analytes was in the range of 91.61%–107.81%, and the RSD was in the range of 0.36%–3.87%. As presented in Table 4, the developed quantitative analysis method showed good precision, accuracy, and recovery, indicating its suitability for accurate quantification of 15 putative active compounds.
TABLE 3
| No. | Compound name | Calibration curves | Linear range (ng/mL) | R2 | LLOQ (ng/mL) |
|---|---|---|---|---|---|
| 1 | Magnoflorine | y = 0.431028x − 0.002095 | 0.05–12.50 | 0.9990 | 0.05 |
| 2 | Albiflorin | y = 0.011666x − 0.000133 | 0.05–12.50 | 0.9994 | 0.05 |
| 3 | Paeoniflorin | y = 0.058974x − 0.004921 | 0.39–100.00 | 0.9993 | 0.39 |
| 4 | Liquiritin apioside | y = 0.094094x − 0.001184 | 0.20–50.00 | 0.9993 | 0.20 |
| 5 | Scutellarin | y = 0.009589x − 0.000028 | 0.02–6.25 | 0.9990 | 0.02 |
| 6 | Acteoside | y = 0.039244x − 0.000634 | 0.05–12.50 | 0.9992 | 0.05 |
| 7 | Liquiritin | y = 0.161083x − 0.000885 | 0.10–25.00 | 0.9992 | 0.10 |
| 8 | Baicalin | y = 0.019967x + 0.000524 | 1.56–400.00 | 0.9994 | 1.56 |
| 9 | Ononin | y = 0.116863x − 0.001513 | 0.05–12.50 | 0.9992 | 0.05 |
| 10 | Wogonoside | y = 0.015843x − 0.006822 | 1.56–400.00 | 0.9992 | 1.56 |
| 11 | Benzoylpaeoniflorin | y = 0.051781x − 0.000331 | 0.02–6.25 | 0.9992 | 0.02 |
| 12 | Baicalein | y = 0.014169x + 0.000632 | 0.10–25.00 | 0.9993 | 0.10 |
| 13 | Glycyrrhizic acid | y = 0.026484x − 0.001693 | 0.39–100.00 | 0.9994 | 0.39 |
| 14 | Wogonin | y = 0.118617x − 0.003641 | 0.10–25.00 | 0.9991 | 0.10 |
| 15 | Oroxylin A | y = 0.127588x − 0.001964 | 0.05–12.50 | 0.9991 | 0.05 |
Linear regression and LLOQ values of 15 analytes in HQD for method validation.
LLOQ, lower limit of quantification; HQD, huangqin decoction.
TABLE 4
| No. | Compound name | Concentration (ng/mL) | Precision (RSD, %) | Accuracy (%) | Recovery (%) | |||
|---|---|---|---|---|---|---|---|---|
| Intra-day | Inter-day | Intra-day | Inter-day | Mean | RSD | |||
| 1 | Magnoflorine | 0.52 | 0.90 | 3.78 | 97.10 | 94.87 | 92.18 | 1.11 |
| 2.08 | 0.94 | 0.52 | 93.60 | 93.32 | 92.07 | 0.84 | ||
| 8.33 | 0.45 | 2.06 | 95.59 | 94.22 | 99.20 | 2.42 | ||
| 2 | Albiflorin | 0.52 | 1.65 | 5.46 | 109.78 | 103.39 | 102.31 | 3.47 |
| 2.08 | 1.50 | 5.86 | 107.83 | 104.54 | 106.86 | 0.63 | ||
| 8.33 | 0.34 | 3.50 | 101.30 | 99.10 | 106.46 | 1.11 | ||
| 3 | Paeoniflorin | 4.17 | 1.47 | 5.91 | 106.06 | 99.30 | 93.08 | 0.86 |
| 16.67 | 1.35 | 7.80 | 104.80 | 101.83 | 100.26 | 0.37 | ||
| 66.67 | 0.39 | 3.60 | 99.95 | 96.76 | 99.42 | 1.27 | ||
| 4 | Liquiritin apioside | 2.08 | 0.88 | 3.74 | 104.01 | 99.93 | 93.91 | 1.31 |
| 8.33 | 0.46 | 6.57 | 108.47 | 101.13 | 95.17 | 0.63 | ||
| 33.33 | 0.95 | 3.31 | 99.91 | 97.08 | 103.88 | 1.78 | ||
| 5 | Scutellarin | 0.26 | 6.28 | 5.67 | 95.23 | 100.60 | 94.99 | 0.92 |
| 1.04 | 2.74 | 1.82 | 90.14 | 91.64 | 91.61 | 2.03 | ||
| 4.17 | 0.47 | 2.28 | 96.79 | 94.53 | 106.08 | 2.44 | ||
| 6 | Acteoside | 0.52 | 1.92 | 3.53 | 100.68 | 97.41 | 92.22 | 1.63 |
| 2.08 | 1.17 | 3.85 | 99.33 | 95.15 | 92.16 | 3.89 | ||
| 8.33 | 1.98 | 1.30 | 96.05 | 95.29 | 94.21 | 2.78 | ||
| 7 | Liquiritin | 1.04 | 1.02 | 4.17 | 106.37 | 102.17 | 92.18 | 0.65 |
| 4.17 | 1.96 | 7.51 | 109.79 | 101.58 | 99.15 | 1.11 | ||
| 16.67 | 1.19 | 4.05 | 98.88 | 96.81 | 104.28 | 1.30 | ||
| 8 | Baicalin | 16.67 | 0.27 | 6.67 | 97.92 | 101.04 | 93.92 | 0.76 |
| 66.67 | 0.84 | 0.25 | 91.40 | 91.17 | 91.88 | 2.39 | ||
| 266.67 | 1.06 | 2.06 | 92.70 | 94.72 | 105.51 | 0.68 | ||
| 9 | Ononin | 0.52 | 1.24 | 1.40 | 98.10 | 99.05 | 97.51 | 0.79 |
| 2.08 | 1.10 | 4.81 | 100.78 | 99.30 | 101.63 | 1.24 | ||
| 8.33 | 0.72 | 0.99 | 99.80 | 98.67 | 104.97 | 0.39 | ||
| 10 | Wogonoside | 16.67 | 0.34 | 2.94 | 99.58 | 97.07 | 93.89 | 1.69 |
| 66.67 | 1.78 | 2.92 | 92.76 | 93.88 | 91.97 | 1.54 | ||
| 266.67 | 0.98 | 3.95 | 92.12 | 96.04 | 107.36 | 0.98 | ||
| 11 | Benzoylpaeoniflorin | 0.26 | 3.70 | 1.69 | 101.74 | 99.82 | 96.91 | 1.72 |
| 1.04 | 1.76 | 7.85 | 109.14 | 102.01 | 101.18 | 0.97 | ||
| 4.17 | 2.30 | 4.19 | 98.73 | 95.68 | 100.24 | 0.47 | ||
| 12 | Baicalein | 1.04 | 1.20 | 8.09 | 93.45 | 102.30 | 106.37 | 3.18 |
| 4.17 | 1.65 | 1.72 | 91.13 | 92.98 | 95.33 | 2.98 | ||
| 16.67 | 0.53 | 5.35 | 99.44 | 97.75 | 105.43 | 1.47 | ||
| 13 | Glycyrrhizic acid | 4.17 | 1.35 | 3.65 | 98.43 | 94.62 | 96.21 | 1.61 |
| 16.67 | 0.83 | 1.21 | 100.79 | 99.43 | 98.10 | 0.61 | ||
| 66.67 | 1.30 | 3.12 | 94.34 | 94.80 | 97.73 | 0.36 | ||
| 14 | Wogonin | 1.04 | 1.68 | 4.97 | 92.19 | 97.16 | 92.67 | 1.16 |
| 4.17 | 0.72 | 2.66 | 92.18 | 93.80 | 101.51 | 2.12 | ||
| 16.67 | 1.10 | 3.28 | 93.22 | 96.86 | 107.81 | 1.71 | ||
| 15 | Oroxylin A | 0.52 | 1.42 | 2.87 | 93.20 | 96.25 | 92.83 | 1.50 |
| 2.08 | 0.47 | 0.97 | 92.30 | 92.29 | 102.42 | 0.75 | ||
| 8.33 | 1.15 | 3.48 | 91.54 | 94.81 | 105.88 | 1.54 | ||
Precision, accuracy, and recovery of 15 analytes in HQD for method validation.
HQD, huangqin decoction; RSD, relative standard deviation.
3.3.3 Sample analysis
The contents of 15 putative active constituents were analyzed in three batches of HQD samples using a validated quantitative analytical method. The detailed results are presented in Table 5. Baicalin and wogonoside showed the highest content with a range of 61.87–64.45 mg/g, followed by liquiritin apioside and glycyrrhizic acid with a range of 8.09–8.62 mg/g.
TABLE 5
| No. | Compound | Content (Mean ± SD, mg/g) | ||
|---|---|---|---|---|
| HQD-1 (n = 6)* | HQD-2 (n = 6) | HQD-3 (n = 6) | ||
| 1 | Magnoflorine | 2.399 ± 0.050 | 2.533 ± 0.028 | 2.487 ± 0.046 |
| 2 | Albiflorin | 1.449 ± 0.033 | 1.449 ± 0.031 | 1.506 ± 0.045 |
| 3 | Paeoniflorin | 14.359 ± 0.087 | 14.485 ± 0.084 | 14.528 ± 0.069 |
| 4 | Liquiritin apioside | 8.553 ± 0.066 | 8.621 ± 0.075 | 8.537 ± 0.044 |
| 5 | Scutellarin | 0.923 ± 0.009 | 0.966 ± 0.024 | 0.967 ± 0.021 |
| 6 | Acteoside | 1.431 ± 0.038 | 1.817 ± 0.074 | 1.821 ± 0.083 |
| 7 | Liquiritin | 2.295 ± 0.042 | 2.311 ± 0.056 | 2.230 ± 0.084 |
| 8 | Baicalin | 63.140 ± 0.823 | 63.937 ± 0.858 | 64.451 ± 0.923 |
| 9 | Ononin | 1.367 ± 0.010 | 1.345 ± 0.010 | 1.362 ± 0.030 |
| 10 | Wogonoside | 62.065 ± 0.706 | 61.873 ± 0.682 | 63.162 ± 0.809 |
| 11 | Benzoylpaeoniflorin | 0.419 ± 0.005 | 0.428 ± 0.005 | 0.417 ± 0.005 |
| 12 | Baicalein | 3.681 ± 0.064 | 3.763 ± 0.098 | 3.896 ± 0.069 |
| 13 | Glycyrrhizic acid | 8.086 ± 0.054 | 8.199 ± 0.067 | 8.238 ± 0.051 |
| 14 | Wogonin | 2.476 ± 0.040 | 2.625 ± 0.034 | 2.561 ± 0.042 |
| 15 | Oroxylin A | 1.101 ± 0.014 | 1.197 ± 0.029 | 1.166 ± 0.042 |
Quantitative analysis of 15 analytes in three batches of HQD.
Replicate injections. HQD, huangqin decoction; SD, standard deviation.
3.4 Molecular docking analysis between putative active constituents and targets
Molecular docking was performed as a complementary computational approach to assess the binding potential of the putative active constituents of HQD toward UC-related target proteins identified through network pharmacology analysis. Based on the results, molecular docking was conducted between 15 candidate compounds and 39 target proteins, and docking scores were calculated for each compound-protein pair. The calculated docking scores were used to estimate the potential interactions between ligands and target proteins, with lower docking scores indicating more favorable predicted binding under the selected docking conditions. As shown in Figure 6, all compound–protein pairs yielded docking scores below −5.0 kcal/mol, suggesting potentially favorable predicted binding affinities between the selected compounds and target proteins. Among the analyzed target proteins, NOS2 showed relatively favorable docking scores with most of the selected candidate compounds. In particular, baicalin, benzoylpaeoniflorin, scutellarin, and wogonoside exhibited docking scores below −11.0 kcal/mol against NOS2. The interaction between the ligand and the receptor was primarily stabilized by several key amino acid residues, including phenylalanine (PHE369), tryptophan (TRP372), and glutamic acid (GLU377), which were commonly observed in all four compounds (Figure 7; Supplementary Figure S1). In addition, AHR–baicalin (Figure 8A), AHR–scutellarin (Figure 8B), MMP3–ononin (Figure 8C), MMP9–liquiritin (Figure 8D), and PRKCQ–glycyrrhizic acid (Figure 8E) also exhibited relatively favorable docking scores below −11.0 kcal/mol, supporting the potential interactions predicted by the network pharmacology analysis (Supplementary Figure S2).
FIGURE 6
FIGURE 7
FIGURE 8
4 Discussion
Although many studies have demonstrated that HQD is effective in treating UC, most have been limited to in vitro and animal experiments, and systematic and comprehensive studies on the components of HQD remain lacking (). In this study, we utilized network pharmacology to select putative active constituents of HQD for UC treatment and established a quantitative analysis method, thereby increasing our understanding of the effective ingredients of HQD. First, 51 chemical compounds present in HQD were identified using UHPLC-Q-Orbitrap-MS. Comprehensive identification of chemical components provides a solid foundation for performing network pharmacology, and based on this, it enables the exploration of putative active constituents through integrated analysis with network pharmacology (; ). In addition, we identified DEGs in patients with UC using the GEO database to analyze biological information and derived genes related to UC. These results have enabled the discovery of potential targets for the diagnosis and treatment of UC, as biomarkers for an accurate diagnosis of UC are still lacking (; ). Finally, 947 targets corresponding to compounds identified in HQD, 1,868 UC-related targets, and 2,930 GEO database-related targets were intersected to obtain 109 common target genes associated with for HQD and UC.
Based on these 109 common target genes, GO BP enrichment analysis revealed that these genes are primarily associated with inflammatory responses and MAPK/ERK signaling-related processes. These findings suggest that the therapeutic effects of HQD on UC may involve the regulation of intestinal inflammatory responses. Indeed, inflammatory responses have been reported to play a significant role in the pathogenesis of UC (). The MAPK/ERK signaling pathway plays a crucial role in regulating inflammatory responses by generating inflammatory cytokines and modulating immune cell activation (). Previous studies have reported that inhibition of the MAPK/ERK signaling pathway can significantly alleviate colonic inflammation and mucosal damage in UC mouse models (). Additionally, the MAPK signaling pathway is known to be an important regulator of the inflammatory process in inflammatory diseases such as UC (). Therefore, the regulation of the MAPK/ERK signaling pathway may be an important mechanism supporting the therapeutic effects of HQD on UC. KEGG pathway analysis revealed that common targets are involved in immune and inflammation-related pathways, particularly the IL-17 signaling pathway and Th17 cell differentiation. Previous studies have shown that Th17 cell numbers and IL-17 levels are significantly elevated in patients with UC and positively correlate with disease severity (). Furthermore, it has been reported that intestinal inflammation in UC mouse models can be alleviated through the regulation of Th17 cell differentiation and the IL-17 signaling pathway (). Moreover, the PI3K–Akt and HIF-1 signaling pathways have been associated with common targets and reported to alleviate symptoms of intestinal diseases through the regulation of inflammatory processes (; ). These results suggest that HQD can exert therapeutic effects on UC by regulating immune responses and major intracellular signaling pathways.
In this study, we selected terms related to UC from the GO and KEGG pathway enrichment analysis results and constructed a PPI network for target genes extracted from these terms. We then performed a network analysis of the compounds identified in HQD and the target genes of UC, and finally selected 15 putative active constituents of HQD by considering their association with multiple targets. The bioactive effects of the selected compounds on UC are supported by several previous studies. Magnoflorine, albiflorin, scutellarin, baicalin, and ononin effectively alleviate UC by reducing the production of inflammatory cytokines and regulating signal transduction pathways related to inflammation in dextran sulfate sodium-induced UC mice (; ; ; ; ). In addition, paeoniflorin, acteoside, liquiritin, wogonoside, and baicalein have protective effects on the intestinal epithelial barrier function, suggesting that they can ameliorate colitis (; ; ; ; ). Liquiritin apioside ameliorates gastrointestinal diseases by regulating intestinal metabolites, whereas glycyrrhizic acid ameliorates colonic inflammation and inhibits the development of colitis-related colorectal cancer (; ). In addition, benzoylpaeoniflorin and wogonin are useful in the treatment of inflammatory diseases by inhibiting the expression of inflammatory mediators (; ). Finally, HQD quality was assessed by measuring the content of 15 putative active constituents, and a UPLC-TQ-MS/MS method was established for their quantitative analysis. The putative active ingredient contents in the three batches of HQD were measured, and the results showed that the contents of baicalin and wogonoside were the highest. In a previously reported clinical study on patients with UC, baicalin was shown to regulate immune balance and alleviate inflammatory responses caused by UC by promoting the proliferation of CD4+CD29+ cells ().
The molecular docking results provided computational support for the potential interactions between the candidate bioactive constituents and target proteins predicted by the network pharmacology analysis. All selected compounds yielded docking scores below −5.0 kcal/mol with their corresponding targets, suggesting potentially favorable interactions under the selected docking conditions. These findings support the potential compound–target associations identified through network analysis and indicate that these interactions may contribute to the therapeutic effects of HQD. Among the identified targets, NOS2 showed relatively favorable docking scores with several compounds. In particular, baicalin and wogonoside, which were present at relatively high levels in the quantitative analysis of HQD, exhibited docking scores below −11.0 kcal/mol against NOS2. NOS2 is an important inflammatory mediator involved in the pathophysiology of UC. Previous studies have shown that NOS2 expression is significantly increased in UC models and that NOS2 inhibition can alleviate the progression of UC by regulating inflammatory signaling pathways (; ). Therefore, the results of this study suggest that major compounds of HQD, such as baicalin and wogonoside, may play an important role in the treatment of UC by regulating NOS2 activity. In addition to NOS2, several other target proteins exhibited strong binding affinities with specific putative active compounds. For example, AHR showed strong interactions with baicalin and scutellarin, and the regulation of AHR is associated with anti-inflammatory and immunomodulatory effects (). MMP3 and MMP9 showed strong binding affinities with ononin and liquiritin, respectively, and MMPs have been reported to play important roles in mucosal degradation and inflammation in inflammatory bowel disease (). Additionally, PRKCQ exhibited strong binding affinity with glycyrrhizic acid, and PKCθ is known to play an important role in regulating T cell-mediated colitis and intestinal inflammation (). Collectively, these findings suggest that the candidate bioactive constituents of HQD may exhibit therapeutic effects on UC by targeting key regulators of inflammation and immune responses, including NOS2, AHR, MMPs, and PRKCQ. However, as molecular docking is a predictive method based on static structure, it is insufficient to assess these interactions and biological efficacy. Furthermore, the predicted targets were derived from public databases, and the bioavailability of the identified compounds was not evaluated in this study. Therefore, additional experimental studies, including in vitro and in vivo studies, are required to confirm the pharmacological effects and mechanisms of action of HQD.
These results demonstrate that the proposed method can be extended beyond simple quality assessments based on the quantitative analysis of compounds within a prescription. In this study, we integrated network pharmacology and molecular docking to evaluate the quality of HQD based on putative active constituents and to explore the potential compound–target relationships associated with UC. In particular, the favorable docking scores further supported the predicted compound–target relationships. These results suggest that the identified putative active constituents may contribute to the therapeutic effects of HQD on UC and that the established analytical method is reliable for both quality assessment and the characterization of putative active constituents. However, this study is limited to the prediction of putative active constituents and potential targets. Therefore, further experimental verification using cell and animal models is necessary to elucidate the mechanism of action of HQD on UC and confirm the involvement of the proposed pathways and targets.
5 Conclusion
In this study, we identified 51 compounds in HQD using UHPLC-Q-Orbitrap-MS profiling. Based on the identified composition, we combined network pharmacology and bioinformatic analyses to identify 15 putative active constituents closely related to UC-related target genes. We established a quantitative analysis method for these putative active constituents and measured their contents in HQD extracts. In addition, molecular docking analysis predicted favorable binding interactions between these putative active compounds and key target proteins, providing computational support for the predicted compound–target interactions. These results provide a quantitative analytical approach for the characterization of HQD constituents and offer preliminary insights into the chemical basis of HQD associated with its potential therapeutic effects against UC. Further studies are needed to elucidate the mechanism of action of HQD in UC treatment and experimentally validate its therapeutic effects.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
SJ: Data curation, Writing – original draft, Writing – review and editing. YJK: Data curation, Writing – review and editing. Y-HH: Conceptualization, Funding acquisition, Writing – review and editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This research was funded by the Korea Institute of Oriental Medicine (grant number: KSN2225021).
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2026.1897628/full#supplementary-material
Abbreviations
HQD, Huangqin decoction; UC, Ulcerative colitis; GEO, Gene Expression Omnibus; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; PPI, Protein–protein interaction; TCM, Traditional Chinese medicine; UHPLC-Q-Orbitrap-MS, Ultra-high-performance liquid chromatography/quadrupole Orbitrap mass spectrometry; UPLC-TQ-MS/MS, Ultra-performance liquid chromatography coupled with a triple-quadrupole tandem mass spectrometry; DEGs, Differentially expressed genes; IS, Internal standard; QC, Quality control; MRM, Multiple reaction monitoring; RSD, Relative standard deviation; BP, Biological process.
References
1
BaumgartD. C.SandbornW. J. (2007). Inflammatory bowel disease: clinical aspects and established and evolving therapies. Lancet369 (9573), 1641–1657. 10.1016/S0140-6736(07)60751-X
2
ChenY. L.XuB.PanZ. F.CaiY. P.YangC. Y.CaoS. L.et al (2025). Glycyrrhizic acid reduces neutrophil extracellular trap formation to ameliorate colitis-associated colorectal cancer by inhibiting peptidylarginine deiminase 4. J. Ethnopharmacol.341, 119337. 10.1016/j.jep.2025.119337
3
ChengF.LiQ.WangJ.ZengF.WangK.ZhangY. (2020). Identification of differential intestinal mucosa transcriptomic biomarkers for ulcerative colitis by bioinformatics analysis. Dis. Markers2020, 8876565. 10.1155/2020/8876565
4
Di MattiaM.SalleseM.NeriM.LopetusoL. R. (2024). Hypoxic functional regulation pathways in the GI tract: focus on the HIF-1alpha and microbiota's crosstalk. Inflamm. Bowel Dis.30 (8), 1406–1418. 10.1093/ibd/izae046
5
DuanZ. L.WangY. J.LuZ. H.TianL.XiaZ. Q.WangK. L.et al (2023). Wumei wan attenuates angiogenesis and inflammation by modulating RAGE signaling pathway in IBD: network pharmacology analysis and experimental evidence. Phytomedicine111, 154658. 10.1016/j.phymed.2023.154658
6
GanzeraM.SturmS. (2018). Recent advances on HPLC/MS in medicinal plant analysis-An update covering 2011-2016. J. Pharm. Biomed. Anal.147, 211–233. 10.1016/j.jpba.2017.07.038
7
GaoW.WangC.YuL.ShengT.WuZ.WangX.et al (2019). Chlorogenic acid attenuates dextran sodium sulfate-induced ulcerative colitis in mice through MAPK/ERK/JNK pathway. Biomed. Res. Int.2019, 6769789. 10.1155/2019/6769789
8
GaoW. Y.SiN.LiM. L.GuX. R.ZhangY.ZhouY. Y.et al (2021). The integrated study on the chemical profiling and in vivo course to explore the bioactive constituents and potential targets of Chinese classical formula qingxin lianzi yin decoction by UHPLC-MS and network pharmacology approaches. J. Ethnopharmacol.272, 113917. 10.1016/j.jep.2021.113917
9
GuoM.WangX. (2023). Pathological mechanism and targeted drugs of ulcerative colitis: a review. Med. Baltim.102 (37), e35020. 10.1097/MD.0000000000035020
10
HuangS.FuY.XuB.LiuC.WangQ.LuoS.et al (2020). Wogonoside alleviates colitis by improving intestinal epithelial barrier function via the MLCK/pMLC2 pathway. Phytomedicine68, 153179. 10.1016/j.phymed.2020.153179
11
HuangJ.ZhengY.MaJ.MaJ.LuM.MaX.et al (2021). Exploration of the potential mechanisms of wumei pill for the treatment of ulcerative colitis by network pharmacology. Gastroenterol. Res. Pract.2021, 4227668. 10.1155/2021/4227668
12
KawalecP. (2016). Indirect costs of inflammatory bowel diseases: crohn's disease and ulcerative colitis. A systematic review. Arch. Med. Sci.12 (2), 295–302. 10.5114/aoms.2016.59254
13
KimC.SimH.BaeJ. S. (2022). Benzoylpaeoniflorin activates anti-inflammatory mechanisms to mitigate sepsis in cell-culture and mouse sepsis models. Int. J. Mol. Sci.23 (21), 13130. 10.3390/ijms232113130
14
LiT.ZhuangS.WangY.WangY.WangW.ZhangH.et al (2016). Flavonoid profiling of a traditional Chinese medicine formula of huangqin tang using high performance liquid chromatography. Acta Pharm. Sin. B6 (2), 148–157. 10.1016/j.apsb.2016.01.001
15
LiM. Y.LuoH. J.WuX.LiuY. H.GanY. X.XuN.et al (2019). Anti-inflammatory effects of huangqin decoction on dextran sulfate sodium-induced ulcerative colitis in mice through regulation of the gut microbiota and suppression of the Ras-PI3K-Akt-HIF-1alpha and NF-kappaB pathways. Front. Pharmacol.10, 1552. 10.3389/fphar.2019.01552
16
LiM. Y.LiM. X.XuN.LiZ. H.ZhangY. M.GanY. X.et al (2021). Effects of huangqin decoction on ulcerative colitis by targeting estrogen receptor alpha and ameliorating endothelial dysfunction based on system pharmacology. J. Ethnopharmacol.271, 113886. 10.1016/j.jep.2021.113886
17
LiM. Y.WuY. Z.QiuJ. G.LeiJ. X.LiM. X.XuN.et al (2023). Huangqin decoction ameliorates ulcerative colitis by regulating fatty acid metabolism to mediate macrophage polarization via activating FFAR4-AMPK-PPARalpha pathway. J. Ethnopharmacol.311, 116430. 10.1016/j.jep.2023.116430
18
LiY.YanM.ZhangM.ZhangB.XuB.DingX.et al (2024). Scutellarin alleviated ulcerative colitis through gut microbiota-mediated cAMP/PKA/NF-kappaB pathway. Biochem. Biophys. Res. Commun.735, 150837. 10.1016/j.bbrc.2024.150837
19
LiM.WuJ.LiD.LeiJ.HuangS.MoP.et al (2025). Baicalin ameliorates dextran sulfate solidum-induced colitis by modulating Th17-macrophage immune network via JAK2/STAT3/IL-17/NF-kappaB pathway. Phytomedicine148, 157414. 10.1016/j.phymed.2025.157414
20
LiM. X.LiM. Y.LeiJ. X.WuY. Z.LiZ. H.ChenL. M.et al (2022). Huangqin decoction ameliorates DSS-Induced ulcerative colitis: role of gut microbiota and amino acid metabolism, mTOR pathway and intestinal epithelial barrier. Phytomedicine100, 154052. 10.1016/j.phymed.2022.154052
21
LiY. Y.WangX. J.SuY. L.WangQ.HuangS. W.PanZ. F.et al (2022). Baicalein ameliorates ulcerative colitis by improving intestinal epithelial barrier via AhR/IL-22 pathway in ILC3s. Acta Pharmacol. Sin.43 (6), 1495–1507. 10.1038/s41401-021-00781-7
22
LiuS.LiQ.LiuF.CaoH.LiuJ.ShanJ.et al (2021). Uncovering the mechanism of curcuma in the treatment of ulcerative colitis based on network pharmacology, molecular docking technology, and experiment verification. Evid. Based Complement. Altern. Med.2021, 6629761. 10.1155/2021/6629761
23
LiuJ.SunL.ChenD.HuoX.TianX.LiJ.et al (2022). Prdx6-induced inhibition of ferroptosis in epithelial cells contributes to liquiritin-exerted alleviation of colitis. Food Funct.13 (18), 9470–9480. 10.1039/d2fo00945e
24
LiuY.LiB. G.SuY. H.ZhaoR. X.SongP.LiH.et al (2022). Potential activity of traditional Chinese medicine against Ulcerative colitis: a review. J. Ethnopharmacol.289, 115084. 10.1016/j.jep.2022.115084
25
LuoJ.ChenG.LiuD.WangY.QiQ.HuH.et al (2019). Study on the material basis of Houpo Wenzhong decoction by HPLC fingerprint, UHPLC-ESI-LTQ-Orbitrap-MS, and network pharmacology. Molecules24 (14), 2561. 10.3390/molecules24142561
26
MaY.LangX.YangQ.HanY.KangX.LongR.et al (2023). Paeoniflorin promotes intestinal stem cell-mediated epithelial regeneration and repair via PI3K-AKT-mTOR signalling in ulcerative colitis. Int. Immunopharmacol.119, 110247. 10.1016/j.intimp.2023.110247
27
MaronekM.MarafiniI.GardlikR.LinkR.TronconeE.MonteleoneG. (2021). Metalloproteinases in inflammatory bowel diseases. J. Inflamm. Res.14, 1029–1041. 10.2147/JIR.S288280
28
NagahamaK.OgawaA.ShiraneK.ShimomuraY.SugimotoK.MizoguchiA. (2008). Protein kinase C theta plays a fundamental role in different types of chronic colitis. Gastroenterology134 (2), 459–469. 10.1053/j.gastro.2007.11.005
29
NiuS.ZhangY. (2025). Targeting the PI3K/Akt pathway for the treatment of ulcerative colitis: integrative regulatory features of traditional Chinese medicine. Front. Pharmacol.16, 1620138. 10.3389/fphar.2025.1620138
30
OrdasI.EckmannL.TalaminiM.BaumgartD. C.SandbornW. J. (2012). Ulcerative colitis. Lancet380 (9853), 1606–1619. 10.1016/S0140-6736(12)60150-0
31
PanZ.LinH.FuY.ZengF.GuF.NiuG.et al (2023). Identification of gene signatures associated with ulcerative colitis and the association with immune infiltrates in colon cancer. Front. Immunol.14, 1086898. 10.3389/fimmu.2023.1086898
32
SagginiR. P. R. (2024). MAPK is implicated in sepsis, immunity, and inflammation. Int. J. Infect.8 (3), 100–104.
33
ShahV. P.MidhaK. K.FindlayJ. W.HillH. M.HulseJ. D.McGilverayI. J.et al (2000). Bioanalytical method validation--a revisit with a decade of progress. Pharm. Res.17 (12), 1551–1557. 10.1023/a:1007669411738
34
ShenC.WangY.ZhangH.LiW.ChenW.KuangM.et al (2023). Exploring the active components and potential mechanisms of Rosa roxburghii Tratt in treating type 2 diabetes mellitus based on UPLC-Q-exactive Orbitrap/MS and network pharmacology. Chin. Med.18 (1), 12. 10.1186/s13020-023-00713-z
35
ShiL.HanX.LiJ. X.LiaoY. T.KouF. S.WangZ. B.et al (2020). Identification of differentially expressed genes in ulcerative colitis and verification in a colitis mouse model by bioinformatics analyses. World J. Gastroenterol.26 (39), 5983–5996. 10.3748/wjg.v26.i39.5983
36
StockingerB.DiazO. E.WincentE. (2024). The influence of AHR on immune and tissue biology. EMBO Mol. Med.16 (10), 2290–2298. 10.1038/s44321-024-00135-w
37
StroberW.FussI. J. (2011). Proinflammatory cytokines in the pathogenesis of inflammatory bowel diseases. Gastroenterology140 (6), 1756–1767. 10.1053/j.gastro.2011.02.016
38
TangX.CuiY.FengB. (2024). The chemical constituents and metabolite profiles of huangqin decoction in normal and ulcerative colitis rats by UHPLC-Q-TOF/MS analysis. J. Pharm. Biomed. Anal.237, 115763. 10.1016/j.jpba.2023.115763
39
UngaroR.ColombelJ. F.LissoosT.Peyrin-BirouletL. (2019). A treat-to-target update in ulcerative colitis: a systematic review. Am. J. Gastroenterol.114 (6), 874–883. 10.14309/ajg.0000000000000183
40
WangC.TangX.ZhangL. (2017). Huangqin-Tang and ingredients in modulating the pathogenesis of ulcerative colitis. Evid. Based Complement. Altern. Med.2017, 7016468. 10.1155/2017/7016468
41
WangX.ZhouW.WangQ.ZhangY.LingY.ZhaoT.et al (2021). A novel and comprehensive strategy for quality control in complex Chinese medicine formula using UHPLC-Q-Orbitrap HRMS and UHPLC-MS/MS combined with network pharmacology analysis: take Tangshen formula as an example. J. Chromatogr. B Anal. Technol. Biomed. Life Sci.1183, 122889. 10.1016/j.jchromb.2021.122889
42
WangX.SuL.TanJ.DingT.YueY. (2023). Albiflorin alleviates DSS-induced ulcerative colitis in mice by reducing inflammation and oxidative stress. Iran. J. Basic Med. Sci.26 (1), 48–56. 10.22038/IJBMS.2022.66678.14624
43
Wedebye SchmidtE. G.LarsenH. L.KristensenN. N.PoulsenS. S.Lynge PedersenA. M.ClaessonM. H.et al (2013). TH17 cell induction and effects of IL-17A and IL-17F blockade in experimental colitis. Inflamm. Bowel Dis.19 (8), 1567–1576. 10.1097/MIB.0b013e318286fa1c
44
WenJ.YangY.LiL.XieJ.YangJ.ZhangF.et al (2024). Magnoflorine alleviates dextran sulfate sodium-induced ulcerative colitis via inhibiting JAK2/STAT3 signaling pathway. Phytother. Res.38 (9), 4592–4613. 10.1002/ptr.8271
45
WuY.LiuX.LiG. (2022). Integrated bioinformatics and network pharmacology to identify the therapeutic target and molecular mechanisms of huangqin decoction on ulcerative Colitis. Sci. Rep.12 (1), 159. 10.1038/s41598-021-03980-8
46
WuH.JiangF.YuanW.ZhaoY.LiuN.MiaoX. (2024). Rauwolfia polysaccharide can inhibit the progress of ulcerative colitis through NOS2-mediated JAK2/STAT3 pathway. PLoS One19 (4), e0301660. 10.1371/journal.pone.0301660
47
XiaX.ZhangY.ZhuL.YingY.HaoW.WangL.et al (2023). Liquiritin apioside alleviates colonic inflammation and accompanying depression-like symptoms in colitis by gut metabolites and the balance of Th17/Treg. Phytomedicine120, 155039. 10.1016/j.phymed.2023.155039
48
XuL.ZhangJ.WangY.ZhangZ.WangF.TangX. (2021). Uncovering the mechanism of ge-gen-qin-lian decoction for treating ulcerative colitis based on network pharmacology and molecular docking verification. Biosci. Rep.41 (2), BSR20203565. 10.1042/BSR20203565
49
YinjuanY.QiuG.MianyongW. (2024). Efficacy of Huangqin decoction combined with western medicine in the treatment of pediatric ulcerative colitis. Chin. J. Prim. Med. Pharm.31 (2), 267–272. 10.3760/cma.j.cn341190-20230711-00607
50
YuF. Y.HuangS. G.ZhangH. Y.YeH.ChiH. G.ZouY.et al (2014). Effects of baicalin in CD4 + CD29 + T cell subsets of ulcerative colitis patients. World J. Gastroenterol.20 (41), 15299–15309. 10.3748/wjg.v20.i41.15299
51
YuT.LuX.LiangY.YangL.YinY.ChenH. (2023). Ononin alleviates DSS-Induced colitis through inhibiting NLRP3 inflammasome via triggering mitophagy. Immun. Inflamm. Dis.11 (2), e776. 10.1002/iid3.776
52
ZhangC. L.ZhangS.HeW. X.LuJ. L.XuY. J.YangJ. Y.et al (2017). Baicalin May alleviate inflammatory infiltration in dextran sodium sulfate-induced chronic ulcerative colitis via inhibiting IL-33 expression. Life Sci.186, 125–132. 10.1016/j.lfs.2017.08.010
53
ZhangR.ZhuX.BaiH.NingK. (2019). Network pharmacology databases for traditional Chinese medicine: review and assessment. Front. Pharmacol.10, 123. 10.3389/fphar.2019.00123
54
ZhangH.ZhangY.LiJ.FanC.GuJ.JinY.et al (2023). UHPLC Q-Orbitrap mass spectrometry-based molecular networking for identification of chemical constituents in the multi-herb formula runyan Mixture. ACS Omega8 (7), 6515–6522. 10.1021/acsomega.2c06885
55
ZhangJ. L.ZhangM. N.WangH. G.YangX. Z.YuC. G. (2022). Jatrorrhizine alleviates ulcerative colitis via regulating gut microbiota and NOS2 expression. Gut Pathog.14 (1), 41. 10.1186/s13099-022-00514-z
56
ZhangX.ZhangL.ChanJ. C. P.WangX.ZhaoC.XuY.et al (2022). Chinese herbal medicines in the treatment of ulcerative colitis: a review. Chin. Med.17 (1), 43. 10.1186/s13020-022-00591-x
57
ZhangY.ShenJ.MaX.YaoM.ZhangY.CaoD. (2022). Anti-infammatory and antioxidant activities of acteoside isolated from Acanthus ilicifolius var. xiamenensis. Appl. Biol. Chem.65 (1), 29. 10.1186/s13765-022-00695-w
58
ZhaoK.LiJ.AsaiY.ZhangL.HuZ.GaoH. (2026). Vitamin D3 mediates amelioration of ulcerative colitis via the TRPV1-MAPK signaling pathway. Front. Immunol.17, 1727932. 10.3389/fimmu.2026.1727932
59
ZhengS.XueT.WangB.GuoH.LiuQ. (2022). Application of network pharmacology in the study of mechanism of Chinese medicine in the treatment of ulcerative colitis: a review. Front. Bioinform2, 928116. 10.3389/fbinf.2022.928116
60
ZhouY.DouF.SongH.LiuT. (2022). Anti-ulcerative effects of wogonin on ulcerative colitis induced by dextran sulfate sodium via Nrf2/TLR4/NF-kappaB signaling pathway in BALB/c mice. Environ. Toxicol.37 (4), 954–963. 10.1002/tox.23457
61
ZhuJ.WangZ.ChenF.LiuC. (2019). Identification of genes and functional coexpression modules closely related to ulcerative colitis by gene datasets analysis. PeerJ7, e8061. 10.7717/peerj.8061
62
ZhuJ.JinR.SuM.PeiJ.ChangY.JiangM. (2024). Quantitative analysis of active components in Rhodiola species based on disease module-guided network pharmacology. Arab. J. Chem.17 (2), 105570. 10.1016/j.arabjc.2023.105570
Summary
Keywords
bioinformatics analysis, chemical constituents, Huangqin decoction, network pharmacology, UHPLC-Q-Orbitrap-MS, UPLC-TQ-MS/MS
Citation
Jang S, Kim YJ and Hwang Y-H (2026) Identification and quantification of putative active constituents of huangqin decoction for ulcerative colitis treatment through the integration of chemical profiling, network pharmacology, and bioinformatics. Front. Chem. 14:1897628. doi: 10.3389/fchem.2026.1897628
Received
02 June 2026
Revised
03 July 2026
Accepted
09 July 2026
Published
31 July 2026
Volume
14 - 2026
Edited by
Zhisheng Wu, Beijing University of Chinese Medicine, China
Reviewed by
Mu Kai-lang, Guizhou University of Traditional Chinese Medicine, China
Lili Tang, Gansu Medical College, China
Updates
Copyright
© 2026 Jang, Kim and Hwang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Youn-Hwan Hwang, hyhhwang@kiom.re.kr
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.