Abstract

Qufeng Zhitong capsule (QZC), a Chinese patent medicine officially approved in China for the treatment of rheumatoid arthritis (RA) and other diseases, possesses the primary effects of dispelling wind, relieving pain, and promoting blood circulation, whose clinical applications have been confined owing to the incomplete elucidation of its chemical compositions and the underlying molecular mechanism for the treatment of RA. In this study, 61 compounds including 16 phenylpropanoids, 15 organic acids, 13 alkaloids, seven flavonoids, six iridoids, one saccharide, two aldehydes, and one saponin in QZC were simultaneously identified and traced to their herbal origins by ultra-high performance liquid chromatography tandem Q-Exactive Orbitrap high-resolution mass spectrometry (UHPLC/Q-Orbitrap-MS), where 31 of them were unambiguously identified by reference compounds, and the other 30 were tentatively characterized. Besides, all these compounds were proven to have potential pharmacological activity in the treatment of RA based on network pharmacology analysis. In conclusion, this study first investigated the chemical composition and potential pharmacological effects of the main chemical compounds in QZC, which will contribute to the revelation of bioactive compounds in QZC and provide evidence for clinical application.

1. Introduction

Rheumatoid arthritis (RA) is a chronic inflammatory joint disease of an autoimmune nature, featuring persistent inflammatory synovitis, whose progression is influenced by multiple factors including environmental, genetic, and infection factors [13]. With the accelerated pace of life and increased work pressure, the number of people suffering from RA in China is increasing annually, which is difficult to heal and prone to recurrence, seriously affecting their physical and mental health, and quality of life. Besides, nonsteroidal anti-inflammatory drugs and slow-acting antirheumatic drugs, mainly adopted by western medicine, usually lead to gastrointestinal side effects and psychological disorders. RA belongs to the category of “obstinate rheumatism” in traditional Chinese medicine (TCM), which is usually treated by removing wind and cold, dispelling dampness and clearing heat, eliminating blood stasis, and promoting blood circulation with remarkable efficacy [4], where Qufeng Zhitong capsule (QZC) is one of the effective drugs listed in the national essential medicine for the treatment of RA [5].

QZC, an effective patent medicine for treating RA, is composed of seven herbs that play different roles. Herba Geranii (Laoguancao, HG) and Radix Aconiti Kusnezoffii Cocta (Zhicaowu, RAKC) are monarch herbs that mainly exert the effect of dispelling wind and dampness, relaxing tendons, and relieving pain. Flos Carthami (Honghua, FC) and Radix Angelicae Pubescentis (Duhuo, RAP) serve as minister drugs with the function of dispelling wind and dredging collaterals, activating blood circulation and stimulating meridians, and reducing swelling and relieving pain. Radix Dipsaci (Xuduan, RD), Radix et Rhizoma Clematidis (Weilingxian, ReRC), and Herba Visci (Hujisheng, HV) all act as assistant drugs to possess the effects of dispelling wind and dampness, dredging collaterals and relieving pain, invigorating the liver and kidney, and strengthening tendons and bones [6]. With the integrated effects of herbs, QZC can be widely used to alleviate the symptoms of RA, such as limb numbness, waist and knee pain, owing to the synergism of these herbal medicines [7].

At present, most studies of QZC focus on its clinical applications with few analytical studies on its chemical compounds, causing poorly clarified chemical substances. In addition, little attention has been paid to the identification and quality evaluation of the overall medical product. QZC, as a TCM prescription, also has the characteristics of multicomponents, multitargets, and multipathways, and its comprehensive action mechanism on RA is still unclear. Li et al. reported that QZC has a positive effect on cortical bone and can prevent bone erosion in patients with RA [8]. With the rapid development of bioinformatics, network pharmacology, characterized by the holistic idea of TCM, has become a powerful tool and has been successfully applied to forecast the complicated action mechanism of TCM in the remedy of multiple illnesses by constructing and analyzing biological networks. Using existing databases, the potential mechanisms of Chinese herbs can be predicted, which can provide a theoretical reference for in-depth studies on the underlying mechanisms of TCM prescriptions [9].

In this study, a fast and effective method was developed for the chemical characterization of QZC using ultra-high performance liquid chromatography tandem Q-Exactive Orbitrap high-resolution mass spectrometry (UHPLC/Q-Orbitrap-MS). The systematic pharmacological mechanism of QZC in the treatment of RA was predicted via network pharmacology. This is the first study focusing on the analysis of the chemical types, herb origins, and potential pharmacological effects of the main chemical components in QZC, whose results can provide a reference for future research on the chemical compositions, pharmacological effects, and clinical applications of QZC.

2. Materials and Methods

2.1. Reagents and Materials

Methanol and formic acid were purchased from Fisher Scientific (Fair Lawn, NJ, USA). Dimethyl sulfoxide (DMSO) was acquired from Meridian Medical Technologies (MREDA, New York, NY, USA). Water for UHPLC/Q-Orbitrap-MS analysis was purified by the Milli-Q water purification system (Millipore, Billerica, MA, USA). QZC samples were provided by Shaanxi Buchang Pharmaceutical Co., Ltd. (Shaanxi, China). Reference compounds were obtained from Shanghai Yuanye Bio-Technology Co., Ltd. (Shanghai, China), including gallic acid, protocatechuic acid, neochlorogenic acid, chlorogenic acid, cryptochlorogenic acid, loganic acid, corilagin, loganin, isochlorogenic acids A−C, angelol A, columbianetin acetate, osthole, columbianadin, akebia saponin D, 5-hydroxymethyl furfural, 4-hydroxybenzoic acid, caffeic acid, hydroxysafflor yellow A, syringin, ferulic acid, secoxyloganin, ellagic acid, columbianetin, angelol G, and isoimperatorin. Benzoylmesaconine, benzoylaconine, and benzoylhypaconine were obtained from Shanghai Acmec Biochemical Co., Ltd. (Shanghai, China). Kaempferol-3-O-rutinoside was obtained from Chengdu Pufei De Biotech Co., Ltd. (Chengdu, China). The purities of these reference compounds were all determined to be above 98% by UPLC analysis.

2.2. Standard Solution Preparation

Thirty-one standards were accurately weighed and dissolved using methanol in volumetric flasks to obtain individual stock solutions. A certain amount of each stock solution was placed in a 10 mL volumetric flask and diluted to volume with 10% methanol aqueous solution for preparing a mixed reference solution at a final concentration of 0.0603 mg·mL−1 gallic acid, 0.0323 mg·mL−1 protocatechuic acid, 0.0714 mg·mL−1 neochlorogenic acid, 0.0451 mg·mL−1 chlorogenic acid, 0.0506 mg ·mL−1 cryptochlorogenic acid, 0.0603 mg·mL−1 loganic acid, 0.0537 mg·mL−1 corilagin, 0.0539 mg·mL−1 loganin, 0.0498 mg·mL−1 isochlorogenic acid B, 0.0499 mg·mL−1 isochlorogenic acid A, 0.0399 mg·mL−1 isochlorogenic acid C, 0.0543 mg·mL−1 angelol A, 0.0502 mg·mL−1 columbianetin acetate, 0.0494 mg·mL−1 osthole, 0.0320 mg·mL−1 columbianadin, 0.0467 mg·mL−1 akebia saponin D, 0.0500 mg·mL−1 5-hydroxymethyl furfural, 0.0500 mg·mL−1 4-hydroxybenzoic acid, 0.0501 mg·mL−1 caffeic acid, 0.0501 mg·mL−1 hydroxysafflor yellow A, 0.0500 mg·mL−1 syringin, 0.0505 mg·mL−1 ferulic acid, 0.0498 mg·mL−1 secoxyloganin, 0.0503 mg·mL−1 ellagic acid, 0.0506 mg·mL−1 benzoylmesaconine, 0.0500 mg·mL−1 kaempferol-3-O-rutinoside, 0.0505 mg·mL−1 benzoylaconine, 0.0504 mg mL−1 benzoylhypaconine, 0.0502 mg ·mL−1 columbianetin, 0.0498 mg·mL−1 angelol G, and 0.0498 mg·mL−1 isoimperatorin, respectively. All solutions were stored at 4°C when not in use.

2.2.1. Sample Solution Preparation

The accurately weighed QZC powder (0.5 g) was transferred into a 25 mL volumetric flask and ultrasonically extracted with a certain amount of 75% methanol aqueous solution at 60°C for 30 min, which was subsequently diluted to scale by 75% methanol aqueous solution after cooling down to room temperature.

Weilingxian and Duhuo were pulverized into homogeneous powder, respectively. Then, 0.5 g was respectively weighed and ultrasonically extracted with 75% methanol aqueous solution (25 mL) at 60°C for 30 min. Honghua (10 g), Laoguancao (10 g), Hujisheng (10 g), Zhicaowu (10 g), and Xuduan (10 g) were respectively added to 80 mL of water and extracted with reflux extraction twice (3 h each time).

QZC, Weilingxian, and Duhuo sample solutions were diluted 25 times with water, while Honghua, Laoguancao, Hujisheng, Zhicaowu, and Xuduan sample solutions were diluted 25 times with 10% methanol aqueous solution for UHPLC/Q-Orbitrap-MS analysis, respectively.

2.3. UHPLC/Q-Orbitrap-MS Analysis

The qualitative analysis was performed on a UHPLC/Q-Orbitrap-MS (Thermo Fisher Scientific, San Jose, CA, USA) with chromatographic separation achieved using an ACQUITY UPLC® BEH C18 (2.1 × 100 mm, 1.7 µm, Waters, Milford, MA, USA) at 40°C. The mobile phase system was composed of 0.1% formic acid aqueous solution (v/v) (A) and methanol (B) at 0.3 mL min−1, which was performed according to following optimized gradient program: 0–5 min, 3%–9% B; 5–9 min, 9%–15% B; 9–11 min, 15%–17% B; 11–15 min, 17%–27% B; 15–18 min, 27%–28% B; 18–26 min, 28%–50% B; 26–31 min, 50%–66% B; 31–35 min, 66%–74% B; and 35–36 min, 74%–3% B. The mass spectrometer was carried out in both positive and negative ion modes with centroided MS and MS2 spectra recorded from 100 to 1500 m/z in full MS and dd-MS2 (TopN) modes at a resolution of 70000 and 17500, respectively. The maximum injection time was set at 50 and 100 ms for MS1 and MS2, respectively. The automatic gain control (AGC) target of MS1 and MS2 was 3e6 and 1e5, respectively. The optimal MS parameters were as follows: spray voltage, −2.5 kV in negative ion mode and 3.5 kV in positive ion mode; sheath gas flow rate, 35 Arb; aux gas flow rate, 10 Arb; capillary temperature, 350°C; aux gas heater temperature, 350°C; collision energy at 20, 40, and 60 eV. All the data were acquired and processed by the Thermo Scientific™ Xcalibur™ system.

2.4. Network Pharmacology Analysis

The structural formulae and IUPAC-NIST chemical identifiers (InChIs) of the compounds were obtained from the PubChem database (https://pubchem.ncbi.nlm.nih.gov/), which were uploaded to the SwissTargetPrediction (https://www.swisstargetprediction.ch/) and BATMAN-TCM database (https://bionet.ncpsb.org.cn/batman-tcm/) to acquire their targets. All the targets were de-duplicated and then imported to Ingenuity Pathway Analysis (IPA) for core analysis so as to identify canonical pathways and networks of compounds, targets, and diseases.

3. Results and Discussion

3.1. Characterization of Chemical Constituents in Qufeng Zhitong Capsule by UHPLC/Q-Orbitrap-MS

The chemical compounds of QZC were identified by UHPLC/Q-Orbitrap-MS analysis in both positive and negative ion modes. As a result, 61 compounds were characterized by comparison with reported literatures, including 16 phenylpropanoids, 15 organic acids, 13 alkaloids, seven flavonoids, six iridoids, one saccharide, two aldehydes, and one saponin, among which 31 compounds were confirmed by comparison with the reference substances. The fragment ions of identified compounds are shown in Table 1, and the MS spectra of QZC sample and reference compounds are displayed in Figure 1, respectively. The MS spectra of individual herbs are shown in Figure 2. The chemical structures of the compounds identified by comparison with the reference substances in QZC are shown in Figure S1.

Coumarins, which are widely distributed in the plant families of Rutaceae and Umbelliferae, have been reported to have various therapeutic properties, such as antibacterial, antitumor, and antidepression activities [22]. Coumarins in QZC were found to mainly come from Duhuo. Ions at m/z 187.03911 [M+H]+ of compound 40 was observed in positive ion mode, whose product ions at m/z 159.04410 [M + H−CO]+, 131.04938 [M + H−CO−CO]+, 143.04921 [M + H−CO2]+, and 115.05457 [M + H−CO2−CO]+ were also observed. Compound 40 was identified as psoralen following the cleavage rule of coumarin, which is consistent with the reported literature [13]. Additionally, compounds 8, 37, 45, 47, 49, 50, 51, 53, 54, 55, 56, 59, 60, and 61 were identified as umbelliferone, nodakenetin, columbianetin, angelol D, angelol A, angelol B, angelol K, angelol G, columbianetin acetate, isoangenomalin, angenomalin, osthole, isoimperatorin, and columbianadin, respectively.

With the continuous development of natural product chemistry and the constant improvement of analytical techniques, organic acid components in TCM, such as chlorogenic acid and gallic acid, have been discovered to have strong activity, drawing increasing attention globally. Organic acids refer to acidic compounds containing the carboxylic group (−COOH), mostly distributed in herbs with a sour taste, and have been demonstrated to possess various pharmacological activities, including anti-inflammatory, antioxidation, and inhibition of platelet aggregation, which have important clinical application value in the prevention and treatment of cardiovascular diseases [23]. In negative ion mode, organic acids exhibited similar cleavage behaviours by the loss of CO2, H2O, or CO. Compounds 1, 2, 5, 7, 10, 15, 16, 19, 25, 28, 33, 35, 36, 42, and 57 were identified as quinic acid, gallic acid, protocatechuic acid, neochlorogenic acid, 4-hydroxybenzoic acid, caffeic acid, chlorogenic acid, cryptochlorogenic acid, corilagin, ferulic acid, ellagic acid, isochlorogenic acid B, isochlorogenic acid A, isochlorogenic acid C, and 9,12,13-trihydroxy-10-octadecenoic acid, respectively. Quinic acid is mainly present in Hujisheng. The mass spectrum of compound 1 in negative ion mode showed a quasimolecular ion at m/z 191.05511 [M−H] with its molecular formula speculated as C7H12O6, where its fragment ions were detected at m/z 173.00810 [M−H−H2O], 127.03873 [M−H−2H2O−CO], and 109.02816 [M−H−3H2O−CO]. Compared with the cleavage pattern reported in the literature, compound 1 was confirmed as quinic acid [10]. Isochlorogenic acids A−C, neochlorogenic acid, chlorogenic acid, and cryptochlorogenic acid have been shown to contain quinic acid. Compound 33 was identified as ellagic acid, a common dimer of gallic acid, which was characterized by the presence of a quasimolecular ions at m/z 300.99890 [M−H] with fragment ion at m/z 257.00940 [M−H−CO2] and 229.01354 [M−H−CO2−CO] [19].

Iridoids, a type of monoterpene derivative in TCM, are widely distributed and have a variety of biological activities, such as antitumor, antiviral, and anti-inflammatory. Iridoids feature complex structures, mostly containing hemiacetal hydroxyl groups and cyclopentane rings, which are prone to the formation of iridoid glycosides with saccharides. The loss of monosaccharides or monosaccharide residues and neutral molecules such as CO2 (44), H2O (18), C2H2 (26), and CH2O (30), usually occurs in the mass spectra of iridoids [18]. The iridoids in QZC are mostly derived from XuDuan. In negative ion mode, compound 30 exhibited an ion with m/z value of 403.12555 [M−H], and produced fragment ions at m/z 223.06062 [M−H−Glc], 165.05470 [M−H−Glc−CH2COO], 121.02822 [M−H−Glc−CH2COO−CO2], and 95.04887 [M−H−Glc−CH2COO−CO2−C2H2], which are consistent with the mass spectrum fragmentation of secoxyloganin. According to the quasimolecular ions and the secondary fragment ions information, the fragmentation patterns of iridoid glycosides were analyzed. Compounds 20, 27, 29, 46, and 52 were identified as loganic acid, sweroside, loganin, cantleyoside, and triplostoside A.

Aconitine-type alkaloids are the primary active ingredients of Zhicaowu in QZC, which are both effective and toxic components with analgesic, anti-inflammatory, and circulation-promoting effects [24]. Zhicaowu has the effects of dispelling wind, removing dampness, warming menstruation, and relieving pain, which is related to its alkaloids, such as mesaconitine, hypaconitine, and karakoline. In total, 13 alkaloids from QZC were only detected in positive ion mode. Alkaloids with acetyl and benzoyl groups at C-8 and C-14, respectively, are called diester diterpenoid alkaloids (DDAs). Given that C-8 is the first active site of DDAs, DDAs tend to lose acetyl side chains so as to form stable fragment ions. AcOH (60) is usually removed from C-8 when the diester alkaloids are cleaved. Compound 48 showed quasimolecular ion at m/z 616.31207 [M+H]+ in positive ion mode with the product ions at m/z 556.29083 [M + H−AcOH]+, 524.26385 [M + H−AcOH−CH3OH]+, and 496.27423 [M + H−AcOH−CH3OH−CO]+, which was preliminary identified as hypaconitine [16]. Compounds 6, 12, 13, 14, 17, 18, 23, 24, 26, 39, 43, and 44 were identified as karakoline, songorine, mesaconine, isotalatizidine, napelline, aconine, hypaconine, neoline, talatisamine, benzoylmesaconine, benzoylaconine, and benzoylhypaconine, respectively.

Flavonoids are an essential class of active ingredients in TCM with a wide range of biological activities including vasodilation, antihepatic injury, antispasmodic, antibacterial, antiviral, antitumor, and antioxidant, which can also be used as food additives, natural antioxidants, natural pigments, and so on. The Retro-Diels-Alder (RDA) cleavage reaction involves the loss and rearrangement of flavonoid aglycone C ring in different ways in the positive ion mode, where [M+H]+ ions are subjected to collision-induced dissociation to lose neutral radicals or molecules, such as H2O, CO, and C2H2O2 [25], while flavonoid glycosides are prone to lose glycosyl during multistage cleavage. Flavonoids in QZC are mainly derived from Honghua. Compound 41 was tentatively identified as kaempferol-3-O-rutinoside due to the observation of quasimolecular ion at m/z 593.15094 [M–H]in negative ion mode and the product ions at m/z 285.04004 [M−H−rutinose residue], 163.00266 [M−H−rutinose residue–C7H6O2], and 151.00261 [M−H−rutinose residue−C8H6O2]. Compounds 9, 21, 31, 32, 34, and 38 were identified as safflomin A, hydroxysafflor yellow A, 6-hydroxykaempferol-3-O-glucoside, 6-hydroxykaempferol, carthamone, and homoeriodictyol-7-O-β-D-apiosyl-(1 ⟶ 2)-β-D-glucoside, respectively.

To further confirm the structures of the compounds, a total of 31 compounds were verified by comparison with reference standards, where compounds 2, 4, 5, 7, 10, 15, 16, 19−22, 25, 28, 29, 30, 33, 35, 36, 39, 41−45, 49, 53, 54, and 58−61 were unambiguously assigned as gallic acid, 5-hydroxymethylfurfural, protocatechuic acid, neochlorogenic acid, 4-hydroxybenzoic acid, caffeic acid, chlorogenic acid, cryptochlorogenic acid, loganic acid, hydroxysafflor yellow A, syringin, corilagin, ferulic acid, loganin, secoxyloganin, ellagic acid, isochlorogenic acid B, isochlorogenic acid A, benzoylmesaconine, kaempferol-3-O-rutinoside, isochlorogenic acid C, benzoylaconine, benzoylhypaconine, columbianetin, angelol A, angelol G, columbianetin acetate, akebia saponin D, osthole, isoimperatorin, and columbianadin, respectively (Table S1).

As an important ingredient in QZC, RAKC was reported to have alkaloids, such as aconine, mesaconine, and hypaconine, which were responsible for anti-inflammation and analgesia functions. Notably, toxicity would happen when RAKC was misused, which must be concerned about. The pharmacopoeia of the People’s Republic of China stipulates that the total amount of aconitine, mesaconitine, and hypaconine must be no more than 0.04%, and the total amount of benzoylaconitine, benzoylhypaconine, and benzoylmesaconine must range from 0.02% to 0.07% [7]. Importantly, the dosage and the course of treatment of QZC must be paid more attention to in clinical use by taking the detected risky alkaloids from RAKC into account.

HG, RAP, RD, and FC are effective ingredients in QZC. From HG, corilagin was reported to exert antirheumatoid arthritis via downregulation of NF-κB and MAPK signaling pathways [26]. The osthole in RAP can reduce proinflammatory cytokines, such as TNF-α, IL-1β, and IL-6, showing promising analgesic effects [27]. Hydroxysafflor yellow A, which derived from FC, can reduce expression of inflammatory mediators in synovial tissue, thus reducing inflammation and synovitis [28]. In addition, asperosaponin VI, one of the saponins in RD, can significantly inhibit RANKL-induced osteoclast formation to relieve arthritis in vivo [29]. The abovementioned compounds play an important role in the treatment of RA.

3.2. Network Pharmacology

As an open-source database, the SwissTargetPrediction is a web-based tool which can predict the targets of any bioactive small molecule as ligands [30]. IPA is a graphical interface bioinformatics software based on cloud computing, which can search for some important information, such as genes, proteins, diseases, and classical pathways [31].

Based on the 61 compounds identified by UHPLC/Q-Orbitrap-MS, a network pharmacological analysis was carried out to preliminarily predict their potential pharmacological activities. The chemical structures and InChIs for 61 compounds from QZC were obtained from the PubChem database and literatures, which were uploaded to the SwissTargetPrediction database to obtain the targets of HOMO except for hypaconine, angelol A, angelol B, angelol D, angelol G, angelol K, columbianetin, mesaconine, and nodakenetin. The BATMAN-TCM database was used to predict the targets of the aforementioned nine compounds, with a cutoff score of 3.351. A total of 1002 targets were acquired after removing duplicates, which were imported to IPA software for core analysis with 926 targets detected. The 1002 targets corresponding to the 61 index compounds in QZC are shown in Table S2. RA was selected for analysis in the related diseases, resulting in 173 related targets, which are shown in Table S3. The network of “QZC-components-targets-RA” is shown in Figure 3(a). “Pathogen-influenced signaling” was selected for the reason that pathogen infection is an important problem in the research of RA [32]. The canonical pathway analysis for pathogen-influenced signaling of the predicted QZC targets is shown in Figure 3(b).

The pathogenesis of rheumatoid arthritis is still unclear. Current studies suggest that inflammation caused by an imbalance of Th1/Th2 cytokines is an important pathological mechanism [33]. Th1 cells mainly secrete proinflammatory cytokines, such as interferon-γ (IFN-γ) and interleukin-1 (IL-1), to promote the production of inflammatory cytokines. Th2 cells mainly secrete anti-inflammatory cytokines such as interleukin-4 (IL-4) and transforming growth factor-β (TGF-β), and inhibit the production of inflammatory cytokines [34]. Th1 and Th2 cells are balanced in healthy people, while Th1 cells increase and Th2 cells decrease in patients with RA [35]. NF-κB, ICAM-1, JAK, IL-6, IL-23, CCR1, and CCR3 are closely associated with the Th1/Th2 signaling pathway, which may be regulated by the chemical compounds from QZC, suggesting that the chemical compounds from QZC may participate in the balance of the Th1/Th2 signaling pathway, thus, playing a therapeutic role in RA.

The results indicated that the identified compounds detected by UHPLC/Q-Orbitrap-MS have potential pharmacological activity for the treatment of RA, which sheds light on the modern approaches to investigate QZC and promoting the development of QZC.

4. Conclusions

In this study, 61 compounds in QZC were identified by UHPLC/Q-Orbitrap-MS, including 16 phenylpropanoids, 15 organic acids, 13 alkaloids, seven flavonoids, six iridoids, one saccharide, two aldehydes, and one saponin, with their herbal sources traced. Among them, 31 compounds were unambiguously identified by reference compounds, and another 30 compounds were tentatively characterized based on their fragmentation pathways. Network pharmacology analysis revealed that all detected compounds exhibited potential pharmacological activity in the treatment of RA via different signaling pathways. In conclusion, this is the first study focusing on the analysis of the overall chemical basis of QZC, which is conducive to revealing the bioactive compounds and potential application value of QZC.

Data Availability

The data used to support the findings of this study are included within the article and the supplementary information file.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Mengjie Xue and Lihua Zhang contributed equally to this work.

Acknowledgments

This work was supported by grants from the National Key Research and Development Program of China (2021YFC1712905), Science and Technology Project of Haihe Laboratory of Mordern Chinese Medicine (22HHZYJC00003), and Innovation Group of Component-based Chinese Medicine and Intelligent Manufacturing with multicrossed disciplines (No. ZYYCXTD-D-202002).

Supplementary Materials

Figure S1: chemical structures of the compounds identified by comparison with reference substances. Table S1: characterization of chemical constituents from QZC by comparison with the reference substances. Table S2: 1002 targets corresponding to the 61 index compounds in QZC. Table S3: 173 targets related to rheumatoid arthritis. (Supplementary Materials)