High-Throughput Chinmedomics Strategy Discovers the Quality Markers and Mechanisms of Wutou Decoction Therapeutic for Rheumatoid Arthritis

Wutou decoction (WTD) is a traditional Chinese medicine prescription for the treatment of rheumatoid arthritis (RA), and this study systematically analyzed the metabolic mechanism and key pharmacodynamic components of WTD in RA rats by combining untargeted metabolomics and serum pharmacochemistry of traditional Chinese medicine to enrich the evidence of WTD quality markers (Q-markers) studies. WTD prevented synovial edema in RA rats and reduced tumor necrosis factor-alpha and interleukin 6 levels in rat serum, according to the results of an enzyme-linked immunosorbent examination and histopathological inspection. In model rats, pattern recognition and multivariate statistical analysis revealed 24 aberrant metabolites that disrupted linoleic acid metabolism, arachidonic acid metabolism, arginine and proline metabolism, etc. However, continued dosing of WTD for 28 days reversed 13 abnormal metabolites, which may be an important therapeutic mechanism from a metabolomic perspective. Importantly, 12 prototypical components and 16 metabolites from WTD were characterized in RA rat serum. The results of Pearson correlation analysis showed that aconitine, L-ephedrine, L-methylephedrine, quercetin, albiflorin, paeoniflorigenone, astragaline A, astragaloside II, glycyrrhetic acid, glycyrrhizic acid, licurazide, and isoliquiritigenin are the key pharmacological components that regulate the metabolism of RA rats, and they are identified as Q-markers. In sum, utilizing metabolomics and serum pharmacochemistry of traditional Chinese medicine, the metabolic mechanisms and Q-markers of WTD therapy in RA rats were revealed, providing a theoretical basis for the quality control investigation of WTD.

Introduction

Traditional Chinese medicine has a long history of application and accumulated rich clinical experience, TCM plays a pivotal role in the protection of human health, of which herbal medicine is the key medium. A large number of chemicals in herbs and prescriptions poses a difficult problem for quality control. Unlike in the past, the quality control of herbal medicines is not only concerned with the high content of the components in the herbal medicines, but the remarkable medicinal effects are more noticeable. The concept and requirements of quality markers (Q-markers) proposed by Academician Liu have promoted the development of quality control systems for TCM and herbal medicines. In general, ideal herbal medicines and prescriptions need to have five characteristics: 1) quality transfer and traceability, 2) component specificity, 3) component effectiveness, 4) component measurability, 5) formula compatibility environment (He et al., 2018; Bai et al., 2018; Zhang et al., 2018).

Chinmedomics comprehensively blends metabolomics and serum pharmacochemistry of TCM and serves as a way for the discovery of Q-markers for herbal medicines and formulations (Zhang et al., 2019). In short, the metabolic mechanisms of herbal medicines and prescriptions for the treatment of diseases are gained through metabolomics insights (Zhang et al., 2017; Zhang et al., 2020; Qiu et al., 2020), and key pharmacodynamic components and Q-markers that regulate abnormal metabolism are identified from a large number of chemical components by correlation analysis (Wang et al., 2015; Han et al., 2020). Currently, chinmedomics was used as a research method for Q-markers to analyze the active components and therapeutic mechanisms of Shengmai San, Yinchenhao Decoction, Sijunzi Decoction, and Kaixin San (Zhang et al., 2018; Sun et al., 2019; Wang et al., 2019; Zhao et al., 2020).

Wutou decoction (WTD) is a classic prescription for treatment rheumatology, derived from “Shang Han Lun” of Zhang Zhongjing in the Han Dynasty and consisting of Aconitum carmichaeli Debeaux [Ranunculaceae], Ephedra sinica Stapf [Ephedraceae], Paeonia lactiflora Pall. [Ranunculaceae], Astragalus mongholicus Bunge [Leguminosae], and Glycyrrhiza uralensis Fisch. [Leguminosae] with a mass ratio of 10 g: 15 g: 15 g: 15 g: 15 g. Pharmacological studies have shown WTD analgesic and anti-inflammatory effects (Guo et al., 2017; Guo et al., 2020). Besides, WTD mainly contains alkaloids, triamcinolone, and flavonoid compounds that can regulate the biological signal pathways associated with rheumatoid arthritis (RA) (Qi et al., 2014; Xu et al., 2017). To enhance the Q-markers and quality control evidence of WTD, this study adopted metabolomics and serum pharmacochemistry of the TCM to carefully evaluate the global metabolic profile of WTD-treated RA rats and the components of WTD that are closely related to the therapeutic effect (Figure 1).

FIGURE 1

FIGURE 1. The workflow of quality markers research for wutou decoction in the treatment of adjuvant-induced arthritis rat based on chinmedomics strategy.

Material and Methods

Reagents and Medicines

HPLC grade methanol and acetonitrile were purchased from Thermo Fisher Scientific (Massachusetts, United States), analytical grade formic acid, leucine enkephalin, and Freund’s complete adjuvant were purchased from Sigma-Aldrich (Shanghai, China), ultrapure water purchased from Watsons (Guangzhou, China), the kit of tumor necrosis factor-alpha (TNF-α) and interleukin-6 (IL-6) purchased from Jiancheng Biologics (Nanjing, China), the five botanical drugs that compose WTD were purchased from Tongrentang Pharmacy (Harbin, China) and were authenticated by Xijun Wang from the Department of Pharmacognosy of Heilongjiang University of Chinese Medicine.

Preparation of Wutou Decoction Extracts

The crude Aconitum carmichaeli Debeaux [Ranunculaceae], Ephedra sinica Stapf [Ephedraceae], Paeonia lactiflora Pall. [Ranunculaceae], Astragalus mongholicus Bunge [Leguminosae], and Glycyrrhiza uralensis Fisch. [Leguminosae] cut into small round pieces according to the mass ratio 2: 3: 3: 3: 3 and were immersed in distilled water (eight-fold volume) for 1 h and then boiled for 1.5 h. The extracted solutions were filtered through eight layers of gauze, and the filter residue was boiled with distilled water (eight-fold volume) again for another 1 h. Finally, the two parts of the filtrates were combined, and the decoction was transformed into powder by freeze-drying. The frozen powder was dissolved with distilled water and kept at 4°C until use.

Animal Grouping and Administration

Experiments were performed with eight-week-old male Sprague-Dawley rats weighing 180–220 g (Liaoning Changsheng Biotechnology Co., Ltd., License No. SCXK [Liao] 2015-0001) under controlled environmental conditions (12/12 light/dark, 23–25°C, 50 ± 5% humidity) with unlimited access to food and water. After 7 days of acclimatization, all animals were randomly divided into the following three groups (n = 8): control group (C); model group (M), and WTD treatment group (T). The M and T rats were injected with 0.1 ml of Freund’s complete adjuvant at a concentration of 2.0 mg.ml⁻¹ at the base of the tail, and the C rats were injected with the same volume of normal saline (Zhang et al., 2014). On the 14th day of model preparation, rats in the T were orally administered a 7.56 g.kg⁻¹ lyophilized powder solution of WTD once a day for 28 consecutive days (Dong et al., 2015). The C rats received the same volume of ultrapure water. The protocol was approved by the Committee on the Ethics of Animal Experiments of the College of Pharmacy of the Heilongjiang University of Chinese Medicine.

Collect and Prepare Biological Samples

Throughout the experimental period, 12 h urine from 8 p.m. to 8 a.m. the next day was collected in the metabolic cage, and the blood samples were obtained by the final stage of administration: fresh urine centrifugation for 15 min (4°C, 13,000 r·min⁻¹), serum centrifugation for 15 min (4°C, 3500 r·min⁻¹). Urine samples were diluted 10 times with ultrapure water, serum samples were precipitated with four times the volume of methanol and the supernatant was taken. All samples were filtered through a 0.22 μm filter before injection into the ultra-performance liquid chromatography quadrupole time of flight coupled with mass spectrometry (UPLC-Q/TOF-MS) instrument. Separately, 100 μl of urine and serum were obtained from each animal for the preparation of quality control samples in examining and optimizing the metabolomics analysis parameters.

Inflammatory Factors and Histopathological Examination

The concentrations of IL-6 and TNF- in serum were determined using the kit’s instructions. X-rays of the rat feet were taken to evaluate the effects of WTD. Fresh ankles were immersed in the formalin solution, and pathological changes were observed by the affiliated hospital of the Heilongjiang University of Chinese Medicine under proper microscopic circumstances.

Metabolomic Conditions

Chromatography

The metabolomic analysis was performed using an ACQUITY UPLC^® system coupled with a time-of-flight mass spectrometer (Waters Corp.). The HSS T₃ column (2.1 mm × 100 mm id, 1.8 µm; Waters Corp.) separates urine samples, whereas the BEH C₁₈ column (2.1 mm × 100 mm id, 1.7 µm, Waters Corp.) separates serum samples. The column temperature is 35°C. The gradient mobile phase consisted of solvent A (0.1% FA-H₂O) and solvent B (0.1% FA-ACN). Optimized urine gradient elution status: 0–4.5 min, 1–21 A%; 4.5–7 min, 21–40% A; 7–7.5 min, 40–50% A; 7.5–9.5 min, 50–99% A. Optimized serum gradient elution status: 0–1.2 min, 1–39% A; 1.2–4.3 min, 39–68% A; 4.3–4.5 min, 68–71% A; 4.5–7.0 min, 71–78% A; 7.0–8.2 min, 78–90% A; 8.2–9.5 min, 90–99% A. The injection volume is 4 μl and the flow rate is 0.5 ml.min⁻¹ in both modes.

Mass Spectrometric

The Waters G2-Si mass spectrometer system (Waters Corp.) is equipped with an electrospray ion source to detect all metabolites, and the ion source temperature is 110°C. The capillary voltage in the positive ion mode (ESI⁺) is 3.0 kV; in the negative ion mode (ESI⁻), it is 2.7 kV. The cone voltage is 20 V, and the gas flow rate is 50 L.h⁻¹. Dry nitrogen is desolventizing at a flow rate of 800 L.h⁻¹ and a temperature of 350°C. Use the full scan mode of m/z 50–1200 Da to collect MS data. Leucine enkephalin with a concentration of 1 ng/μl and a flow rate of 10 μl.min⁻¹ in ESI⁺ [M+H]⁺ = 556.2771 and ESI⁻ [M-H]⁻ = 554.2615 was used as a reference compound to obtain an accurate mass.

Constituent Analysis of Wutou Decoction

After lyophilized powder of WTD decoction was dissolved in ultrapure water and was centrifuged 15 min (4°C, 13,000 r·min⁻¹), filtered with a 0.22 μm filter membrane, and injected into UPLC-MS for analysis chemical profile of WTD. The serum containing components of WTD was added with a quarter volume of 20% phosphoric acid, utilizing an OASIS MCX solid-phase extraction column (30 mg, 30 μm, Waters Corp.) to enrich the WTD components in vivo, dissolved in methanol for UPLC-MS data collection. Injection volume is 2 μl, the gradient mobile phase consisted of solvent A (0.1% FA-H₂O) and solvent B (0.1% FA-ACN), optimized gradient elution conditions: 0–3.0 min, 1–10% A; 3.0–7.0 min, 10–12% A; 12.0–20.0 min, 20–40% A; 20.0–22.0 min, 40–100% A. The MS^E data scanning mode of UPLC-Q/TOF-MS was combined with powerful UNIFI software (Waters Corp.) for efficient and accurate characterization of the chemical profile of the WTD and the components absorbed in vivo. The adduct ion is [M+H]⁺ in ESI⁺ and [M-H]^- in ESI⁻, respectively. To reduce false positives, select the embedded TCM database and set the tolerance error of components and fragments to 5 mDa. The study compared the serum between the M and the T and preferred the results of the high response value only in the serum of the T.

Multivariate Statistics and Data Analysis

The raw data is imported into Progenesis QI software (Waters Corp.) for noise reduction, peak alignment, and normalization to get an information matrix that includes m/z, retention time, and peak relative intensity. The EZinfo software (Waters Corp.) performs a principal component analysis (PCA) and judges differences between groups. Orthogonal partial least squares discriminant analysis (OPLS-DA) calculates variable importance in projection value (VIP), a statistical normalization abundance of different ions. Finally, MS/MS information combined with HMDB (https://hmdb.ca/) and MassBank (http://www.massbank.jp/) databases identifies metabolites, VIP >1, and metabolites normalized abundance with intergroup significance are considered potential biomarkers of RA rats. Comparisons between the two groups of the quantitative data were completed using Graphpad software (California, United States) independent t-test and presented as mean ± standard deviation, with p < 0.05 indicating significant differences and p < 0.01 indicating extremely significant differences.

Analysis of Key Medicinal Constituents

The PCMS (plotting of correlation between marker metabolites and serum components) patented software calculates the correlation coefficient between the components absorbed into the blood of WTD and the potential biomarkers of RA rats to determine which component is most likely to have a therapeutic effect. Normalize the abundance of biomarkers and components, and then use Pearson’s statistical correlation to set the reference standard for the correlation coefficient R. In this study, 0.7 ≤│R│<0.8 is selected as the highly correlated component; 0.8 ≤│R│≤ 1 is regarded as extremely correlated, and extremely correlated components are regarded as Q-markers.

Results

Therapeutic Effects of Wutou Decoction

Compared with the C rats, 14 days after the injection of Freund’s complete adjuvant, the rats in the M reached the uttermost state of inflammation. Compared with the M, the T significantly inhibited the levels of IL-6 and TNF-α in rats and reduced the swelling of the feet (Figure 2).

FIGURE 2

FIGURE 2. Evaluation of the effect of wutou decoction on adjuvant-induced arthritis rat. (A) X-ray examination (B) Synovial histopathological examination (*200) (C) changes in paw volume before and after administration (D) Determination of inflammatory factors in serum. Results are expressed as mean ± SD, n = 8. ^*p < 0.05, ^**p < 0.01 vs control group. ^#p < 0.05, ^##p < 0.01 vs model group.

Characterize Potential Biomarkers

Multivariate statistical analysis was performed on the metabolic fingerprint information collected from serum and urine at the end of the experiment. There are differences in serum total ion chromatograms (TIC) in C, M, and T rats in ESI⁺ and ESI⁻ and these are reflected in PCA (Figure 3). Similarly, there are differences in the TIC of the urine between the three groups, suggesting a change in the metabolic profile (Figure 4). The OPLS-DA statistical model calculates the VIP of differential ions and is employed to screen potential biomarkers, the R2Y and Q2 values of the urine model (R2Y (cum) = 0.995, Q2(cum) = 0.927 in ESI⁺, and R2Y(cum) = 0.994, Q2(cum) = 0.615 in ESI⁻) and serum model (R2Y (cum) = 0.858, Q2(cum) = 0.798 in ESI⁺, and R2Y (cum) = 0.872, Q2(cum) = 0.795 in ESI⁻) indicated that the statistic model has good quality and predictability.

FIGURE 3

FIGURE 3. Serum metabolic profile collection and multivariate statistical analysis. (A–C) TIC spectra of C, M, and T rats in ESI⁺, respectively (D–F) TIC spectra of C, M, and T rats in ESI⁻, respectively (G,J) PCA plot of C, M, and T in ESI⁺ and ESI⁻, respectively (H,K) OPLS-DA plot of C and M in ESI⁺ and ESI⁻, respectively (I,L) VIP plot of C and M in ESI⁺ and ESI⁻, respectively. Note: : C rats; : M rats; : T rats.

FIGURE 4

FIGURE 4. Urine metabolic profile collection and multivariate statistical analysis. (A–C) TIC spectra of C, M, and T rats in ESI⁺, respectively (D–F) TIC spectra of C, M, and T rats in ESI⁻, respectively (G,J) PCA plot of C, M, and T in ESI⁺ and ESI⁻, respectively (H,K) OPLS-DA plot of C and M in ESI⁺ and ESI⁻, respectively (I,L) VIP plot of C and M in ESI⁺ and ESI⁻, respectively. Note: : C rats; : M rats; : T rats.

Take arachidonic acid as an example to illustrate the identification process of potential biomarkers. First, OPLS-DA calculated that the VIP of 8.09_303.2320m/z was 4.004 and the normalized abundance of this ion in the serum of the C and the M was different (p < 0.05). Secondly, chromatographic peak extraction and elemental composition analysis showed that the tolerance of C₂₀H₃₂O₂ was accepted, while m/z 304 and m/z 279 were simultaneously observed in MS extraction at 8.09 min. Finally, it was identified that m/z 304 Da is the parent ion of arachidonic acid, and the fragments of 303[M-H]⁻ and 279[M-COOH]^- were generated by high-energy bombardment in ESI⁻ (Supplementary Figure S1). Eventually, 24 biomarkers are identified (Table 1).

TABLE 1

TABLE 1. Identification results of potential biomarker groups in adjuvant-induced arthritis rats.

Enrichment Analysis of Potential Biomarkers

The above 24 potential biomarkers of serum and urine are mostly engaged in 17 metabolic pathways, which belong to lipid metabolism, carbohydrate metabolism, and amino acid metabolism (Table 2). Where the impact value is greater than 0.1, linoleic acid metabolism, arachidonic acid metabolism, glyoxylate and dicarboxylate metabolism, arginine and proline metabolism, arginine biosynthesis, steroid hormone biosynthesis, and glycerophospholipid metabolism are considered important metabolic pathways. The WTD significantly reversed potential biomarkers with abnormal levels of 13/24 (Figure 5), including 4-Hydroxynonenal, L-Arginine, Prostaglandin E3, Aracidonic acid, N-Acetyltaurine, Kynurenic acid, D-Glucose, Linoleyl carnitine, Galactose-beta-1,4-xylose, Linoleic acid, Octadecenoylcarnitine, L-Palmitoylcarnitine, Sphingosine L-phosphate. WTD may correct the disturbed Linoleic acid metabolism, Arachidonic acid metabolism, Arginine and proline metabolism, Arginine biosynthesis, Steroid hormone biosynthesis, Sphingolipid metabolism, and Glycolysis/Gluconeogenesis by reversing the levels of the above markers, thereby treating Freund’s complete adjuvant-induced RA rats (Figure 6).

TABLE 2

TABLE 2. Metabolic pathway analysis based on MetPA.

FIGURE 5

FIGURE 5. The normalized abundance of potential biomarkers before and after wutou decoction administration. (A) Normalized abundance values between different groups (B) Heat map of the differences in the potential biomarkers of the C, M and T obtained in the ESI⁺ and ESI⁻. Results are expressed as mean ± SD, n = 8. ^*p < 0.05, ^**p < 0.01 vs control group. ^#p < 0.05, ^##p < 0.01 vs model group. C represents control group, M represents model group, T represents wutou decoction group.

FIGURE 6

FIGURE 6. Enrichment analysis of metabolic pathway of potential biomarkers. (A) Pathway impact enrichment (B) The proportion of biomarkers in the overall pathway (C) Association network of biomarker metabolic pathways.

Blood Components of Wutou Decoction

The UPLC-Q/TOF-MS collects the TIC of WTD before and after entering the blood in ESI⁺ and ESI⁻ (Figure 7). The UNIFI software, with a powerful analysis module, can automatically identify the structure information obtained under high and low energy, including m/z, retention time, and MS/MS fragment response values, and deduce the cleavage relationship of compounds. Compound aconitine was used as an example to describe the process of applying UNIFI software to analyze the chemical profile of WTD. The UNIFI software matched the MS^E data collected in the continuum mode. Both ion response values between groups and all m/z information at 14.12 min were obtained, and m/z 646.32 Da with a high response value in WTD samples was scanned for small molecule ions. After elemental composition and fragment matching, m/z 646.32 Da was recognized as the hydrogenation ion of aconitine and identified to 586.30 Da [M+H-C₂H₄O₂]⁺, 568.29 Da [M+H-C₂H₆O₃]⁺, 554.27 Da [M+H-C₃H₈O₃]⁺ (Supplementary Figure S2). Finally, they characterized and identified 61 components of the WTD (Supplementary Table S1). Importantly, 19 components were absorbed into the blood, including 12 prototype components (Supplementary Table S2).

FIGURE 7

FIGURE 7. TIC spectrum of chemical profile and absorbed components from wutou decoction. (A,B) The chemical profile of wutou decoction in ESI⁺ and ESI⁻, respectively (C,D) The blood components of wutou decoction in ESI⁺ and ESI⁻, respectively.

Discover the Key Medicinal Components

PCMS software was utilized in this investigation to uncover the important pharmacodynamic components of WTD in the treatment of RA rats, as well as biological activity tracking and literature analysis to identify Q-markers. The findings revealed that 19 of the WTD components were absorbed into the serum. According to correlation analysis, 12 components were confirmed as Q-markers (Figure 8), of which aconitine was from Aconitum carmichaeli Debeaux [Ranunculaceae], L-ephedrine, L-methylephedrine, quercetin, and kaempferol-3-O-rhamnoside were from Ephedra sinica Stapf [Ephedraceae], albiflorin and paeoniflorigenone were from Paeonia lactiflora Pall. [Ranunculaceae], astragaline A, and astragaloside II were from Astragalus mongholicus Bunge [Leguminosae], glycyrrhetic acid, glycyrrhizic acid, licurazide, and isoliquiritigenin came from Glycyrrhiza uralensis Fisch. [Leguminosae].

FIGURE 8

FIGURE 8. Pearson correlation analysis of blood components of wutou decoction and potential biomarkers in adjuvant-induced arthritis rat.

Discussion

RA is an autoimmune disease that primarily affects synovial joints and is characterized by chronic synovitis, systemic inflammation, and varying degrees of erosion of bone and cartilage. While the cause of RA is unknown, it involves a combination of genetic and environmental factors, potential mechanisms involving the body’s immune system attacking joints, and researchers have focused on exploring the immune mechanisms and molecular networks of RA (Lee and Weinblatt, 2001). Nonsteroidal anti-inflammatory drugs and biologics are currently used to treat RA, but traditional medicine has attracted the attention of doctors and patients in many developing countries because these treatment methods have similar therapeutic goals: anti-inflammatory and pain suppression, slow joint injury, and prevent the development of inflammation-related complications (Moudgil and Berman, 2014; Wang et al., 2021). WTD is a Chinese prescription for rheumatoid arthritis and joint pain which was described in ancient medical writings. As a result, scientific evidence for WTD is accumulating, including clinical efficacy evaluations, characterization of complex components, and molecular mechanism analysis. The RA rat model was effectively duplicated in this work using a variety of assessment methods after injection of Freund’s complete adjuvant in a short period. Through X-ray examination, histopathological analysis, and the detection of inflammatory factors, it was discovered that the model rats’ feet, tails, and spinal joints were gradually swollen and stiff and that there was a severe inflammatory response in the joints of the rats, so the inflammatory response evolved into bone deformation and stiffness due to unregulated, which was the same pathological phenomenon as in previous similar studies (Pan et al., 2017; Dai et al., 2022). Although clinical features of patients with RA cannot be entirely simulated, the current results may explain why rat models of adjuvant-induced arthritis are commonly implemented in drug development.

Studies indicate metabolic disorders in biofluids in RA patients or animal models, covering lipid and glucose metabolism, and amino acid metabolism (Falconer et al., 2018). Untargeted metabolomics with pattern recognition was used to evaluate the metabolites in the serum and urine of RA rats, and 24 different metabolites were detected. Metabolic pathway enrichment results showed crosstalk of three types of pathways dominated by lipid metabolism, of which the important pathway associated with RA rats with an impact value greater than 0.1 was identified (Table 2). Arachidonic acid is a polyunsaturated essential fatty acid located in the body’s fats, the liver, brain, and a range of other organs. Arachidonic acid performs an important role in the structure and function of the body’s cell membranes; it can directly or indirectly synthesize prostaglandins and leukotrienes, and other biologically active substances, and it can cause inflammation in the body by a cascade of arachidonic acid (Charles-Schoeman et al., 2018; Szczuko et al., 2020). The level of arachidonic acid in the model group was abnormal when compared to the control group, which severely disrupted linoleic acid metabolism and arachidonic acid metabolism in lipid metabolism. However, WTD reversed the level of arachidonic acid and may modulate related metabolic pathways, suggesting that this mechanism is both underlying and noteworthy. Lysolecithin is formed by phosphatase A2 hydrolysis of lecithin, and changes in lysolecithin reflect cell renewal and phospholipid consumption (Charles-Schoeman et al., 2018). These compounds are critical in cell proliferation and inflammatory response. WTD has no positive feedback effect on LysoPC(20:5) in the pathway, but it does have a regulatory influence on steroid hormone biosynthesis of lipid metabolism via L-palmitoylcarnitine, and high levels of L-palmitoylcarnitine may be detrimental to energy metabolism in mitochondria (Imagawa et al., 1984). Both RA and osteoarthritis patients exhibit anomalies in arginine and proline metabolism in their synovial fluid (Carlson et al., 2018). According to KEGG, hydroxypyruvic acid interfered with glyoxylate and dicarboxylate metabolism, and it might cause cascading fluctuations by acting on arginine and proline metabolism via the downstream enzyme 5.3.1.22 (Figure 6C). Unfortunately, WTD only regulates L-arginine and cannot regulate hydroxypyruvic acid.

To discover the Q-markers of WTD, this study employed a chinmedomics strategy to analyze the correlation between the components absorbed into the blood of WTD and the potential biomarkers of RA rats, highlighting the importance of certain components in regulating RA metabolism and markers. Finally, 12 key components were determined as Q-markers for the efficacy of WTD in the treatment of RA. Aconitine, the main active component of Aconitum carmichaeli Debeaux [Ranunculaceae], has immunomodulatory properties that might be useful in treating autoimmune diseases, such as rheumatoid arthritis and systemic lupus erythematosus (Li et al., 2017). Ephedrine can mediate toll-like receptor 4 to encourage the production of interleukin 10 and inhibit the secretion of TNF-α to decrease inflammation (Zheng et al., 2012). Quercetin inhibits arthritis mice even stronger than methotrexate (Haleagrahara et al., 2017), and quercetin regulates the enteric nervous system to protect the jejunal inflammation caused by RA (Piovezana et al., 2019). Albiflorin and paeoniflorigenone could modulate the functions and activation of immune cells, decreases inflammatory medium production, and restore abnormal signal pathway (Zhang and Wei, 2020). Astragalin could attenuate synovial inflammation and joint destruction in RA at least partially by restraining the phosphorylation of mitogen-activated protein kinases and the activating of N-terminal kinase/activated protein 1 (Zhang and Wei, 2020). Moreover, kaempferol-3-O-rhamnoside, glycyrrhetic acid, and isoliquiritigenin all have the potential and activity to resist RA (Puchner et al., 2012; Zhu et al., 2012; Ahmed and Abd Elkarim, 2021). With the help of the correlation analysis of chinmedomics, we finally identified 12 Q-markers of WTD, which have a certain degree of anti-inflammatory and immune-regulating biological activity, but it is unclear whether the combination of these compounds has a better anti-RA effect, this needs further verification.

Conclusion

To summarize, this work used chinmedomics to evaluate the therapeutic effect of WTD on adjuvant-induced arthritic rats and to investigate the Q-markers of WTD. WTD helps alleviate disease states and pathological conditions in inflammatory rats, and chinmedomics as a research tool for TCM effectiveness and pharmacology promotes the discovery of WTD Q-markers and contributes to the theoretical evidence for WTD quality control.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Ethics Statement

The animal study was reviewed and approved by Animal Experiment Center, Heilongjiang University of Chinese Medicine.

Author Contributions

TL, FW, AZ, and IU in charge of experimental design, HD, HL, JM, YiH, and YaH in charge of data collection, XW in charge of funding.

Funding

This work was supported by grants from the Key Program of Natural Science Foundation of State (Grant No. 81830110, 81861168037), Scientific and Technology Development Program of Guangxi (AD18126013), the Ba Gui Scholars program of Guangxi, the Central Government Guides Local Science and Technology Development Fund Projects (ZY21195044), Natural Science Foundation of Heilongjiang Province (LH2019H056), Heilongjiang Touyan Innovation Team Program.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

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/fphar.2022.854087/full#supplementary-material

Abbreviations

ESI, electrospray ion source; IL-6, interleukin-6; OPLS-DA, orthogonal partial least squares discriminant analysis; PCA, principal component analysis; PCMS, A plotting of correlation between marker metabolites and serum constituents; Q-markers, quality markers; RA, rheumatoid arthritis; TCM, traditional Chinese medicine; TIC, total ion chromatograms; TNF-α, tumor necrosis factor-alpha; UPLC-Q/TOF-MS, ultraperformance liquid chromatography quadrupole time of flight mass spectrometry; VIP, variable importance in projection value; WTD, wutou decoction.

References

Ahmed, A. H., and Abd Elkarim, A. S. (2021). Bioactive Compounds with Significant Anti-Rheumatoid Arthritis Effect Isolated for the First Time from Leaves of Bougainvillea Spectabilis. Curr. Pharm. Biotechnol. 22, 2048–2053. doi:10.2174/1389201021666201229111825

PubMed Abstract | CrossRef Full Text | Google Scholar

Bai, G., Zhang, T., Hou, Y., Ding, G., Jiang, M., and Luo, G. (2018). From Quality Markers to Data Mining and Intelligence Assessment: A Smart Quality-Evaluation Strategy for Traditional Chinese Medicine Based on Quality Markers. Phytomedicine 44, 109–116. doi:10.1016/j.phymed.2018.01.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Carlson, A. K., Rawle, R. A., Adams, E., Greenwood, M. C., Bothner, B., and June, R. K. (2018). Application of Global Metabolomic Profiling of Synovial Fluid for Osteoarthritis Biomarkers. Biochem. Biophys. Res. Commun. 499 (2), 182–188. doi:10.1016/j.bbrc.2018.03.117

PubMed Abstract | CrossRef Full Text | Google Scholar

Charles-Schoeman, C., Meriwether, D., Lee, Y. Y., Shahbazian, A., and Reddy, S. T. (2018). High Levels of Oxidized Fatty Acids in HDL Are Associated with Impaired HDL Function in Patients with Active Rheumatoid Arthritis. Clin. Rheumatol. 37 (3), 615–622. doi:10.1007/s10067-017-3896-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, Y., Sheng, J., He, S., Wu, Q., Wang, Y., and Su, L. (2022). Dehydroevodiamine Suppresses Inflammatory Responses in Adjuvant-Induced Arthritis Rats and Human Fibroblast-like Synoviocytes. Bioengineered 13 (1), 268–279. doi:10.1080/21655979.2021.1999554

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, H., Yan, G. L., Han, Y., Sun, H., Zhang, A. H., Li, X. N., et al. (2015). UPLC-Q-TOF/MS-based Metabolomic Studies on the Toxicity Mechanisms of Traditional Chinese Medicine Chuanwu and the Detoxification Mechanisms of Gancao, Baishao, and Ganjiang. Chin. J. Nat. Med. 13 (9), 687–698. doi:10.1016/S1875-5364(15)30067-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Falconer, J., Murphy, A. N., Young, S. P., Clark, A. R., Tiziani, S., Guma, M., et al. (2018). Review: Synovial Cell Metabolism and Chronic Inflammation in Rheumatoid Arthritis. Arthritis Rheumatol. 70 (7), 984–999. doi:10.1002/art.40504

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, Q., Zheng, K., Fan, D., Zhao, Y., Li, L., Bian, Y., et al. (2017). Wu-Tou Decoction in Rheumatoid Arthritis: Integrating Network Pharmacology and In Vivo Pharmacological Evaluation. Front. Pharmacol. 8, 230. doi:10.3389/fphar.2017.00230

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, Q., Mizuno, K., Okuyama, K., Lin, N., Zhang, Y., Hayashi, H., et al. (2020). Antineuropathic Pain Actions of Wu-tou Decoction Resulted from the Increase of Neurotrophic Factor and Decrease of CCR5 Expression in Primary Rat Glial Cells. Biomed. Pharmacother. 123, 109812. doi:10.1016/j.biopha.2020.109812

PubMed Abstract | CrossRef Full Text | Google Scholar

Haleagrahara, N., Miranda-Hernandez, S., Alim, M. A., Hayes, L., Bird, G., and Ketheesan, N. (2017). Therapeutic Effect of Quercetin in Collagen-Induced Arthritis. Biomed. Pharmacother. 90, 38–46. doi:10.1016/j.biopha.2017.03.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, Y., Sun, H., Zhang, A., Yan, G., and Wang, X. J. (2020). Chinmedomics, a New Strategy for Evaluating the Therapeutic Efficacy of Herbal Medicines. Pharmacol. Ther. 216, 107680. doi:10.1016/j.pharmthera.2020.107680

PubMed Abstract | CrossRef Full Text | Google Scholar

He, J., Feng, X., Wang, K., Liu, C., and Qiu, F. (2018). Discovery and Identification of Quality Markers of Chinese Medicine Based on Pharmacokinetic Analysis. Phytomedicine 44, 182–186. doi:10.1016/j.phymed.2018.02.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Imagawa, T., Watanabe, T., and Nakamura, T. (1984). Inhibition of Palmitate Oxidation in Mitochondria by Lipid Hydroperoxides. J. Biochem. 95 (3), 771–778. doi:10.1093/oxfordjournals.jbchem.a134668

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, D. M., and Weinblatt, M. E. (2001). Rheumatoid Arthritis. Lancet 358 (9285), 903–911. doi:10.1016/S0140-6736(01)06075-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Gu, L., Yang, L., Zhang, D., and Shen, J. (2017). Aconitine: A Potential Novel Treatment for Systemic Lupus Erythematosus. J. Pharmacol. Sci. 133 (3), 115–121. doi:10.1016/j.jphs.2017.01.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Moudgil, K. D., and Berman, B. M. (2014). Traditional Chinese Medicine: Potential for Clinical Treatment of Rheumatoid Arthritis. Expert Rev. Clin. Immunol. 10 (7), 819–822. doi:10.1586/1744666X.2014.917963

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, T., Cheng, T. F., Jia, Y. R., Li, P., and Li, F. (2017). Anti-rheumatoid Arthritis Effects of Traditional Chinese Herb Couple in Adjuvant-Induced Arthritis in Rats. J. Ethnopharmacol. 205, 1–7. doi:10.1016/j.jep.2017.04.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Piovezana Bossolani, G. D., Silva, B. T., Colombo Martins Perles, J. V., Lima, M. M., Vieira Frez, F. C., Garcia de Souza, S. R., et al. (2019). Rheumatoid Arthritis Induces Enteric Neurodegeneration and Jejunal Inflammation, and Quercetin Promotes Neuroprotective and Anti-inflammatory Actions. Life Sci. 238, 116956. doi:10.1016/j.lfs.2019.116956

PubMed Abstract | CrossRef Full Text | Google Scholar

Puchner, A., Hayer, S., Niederreiter, B., Hladik, A., Blueml, S., Bonelli, M., et al. (2012). Effects of 18β-Glycyrrhetinic Acid in hTNFtg Mice - a Model of Rheumatoid Arthritis. Wien. Klin. Wochenschr. 124 (5-6), 170–176. doi:10.1007/s00508-011-0103-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Qi, Y., Li, S., Pi, Z., Song, F., Lin, N., Liu, S., et al. (2014). Chemical Profiling of Wu-tou Decoction by UPLC-Q-TOF-MS. Talanta 118, 21–29. doi:10.1016/j.talanta.2013.09.054

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiu, S., Zhang, A.-h., Guan, Y., Sun, H., Zhang, T.-l., Han, Y., et al. (2020). Functional Metabolomics Using UPLC-Q/TOF-MS Combined with Ingenuity Pathway Analysis as a Promising Strategy for Evaluating the Efficacy and Discovering Amino Acid Metabolism as a Potential Therapeutic Mechanism-Related Target for Geniposide against Alcoholic Liver Disease. RSC Adv. 10 (5), 2677–2690. doi:10.1039/c9ra09305b

CrossRef Full Text | Google Scholar

Sun, H., Zhang, A. H., Yang, L., Li, M. X., Fang, H., Xie, J., et al. (2019). High-throughput Chinmedomics Strategy for Discovering the Quality-Markers and Potential Targets for Yinchenhao Decoction. Phytomedicine 54, 328–338. doi:10.1016/j.phymed.2018.04.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Szczuko, M., Kikut, J., Komorniak, N., Bilicki, J., Celewicz, Z., and Ziętek, M. (2020). The Role of Arachidonic and Linoleic Acid Derivatives in Pathological Pregnancies and the Human Reproduction Process. Int. J. Mol. Sci. 21 (24), 9628. doi:10.3390/ijms21249628

CrossRef Full Text | Google Scholar

Wang, X., Zhang, A., Sun, H., and Yan, G. (2015). “Origin of Chinmedomics,” in Chinmedomics. Editors X. Wang, A. Zhang, and H. Sun (Boston: Academic Press), 1–15. doi:10.1016/b978-0-12-803117-9.00001-9

CrossRef Full Text | Google Scholar

Wang, X. J., Zhang, A. H., Kong, L., Yu, J. B., Gao, H. L., Liu, Z. D., et al. (2019). Rapid Discovery of Quality-Markers from Kaixin San Using Chinmedomics Analysis Approach. Phytomedicine 54, 371–381. doi:10.1016/j.phymed.2017.12.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Chen, S., Du, K., Liang, C., Wang, S., Owusu Boadi, E., et al. (2021). Traditional Herbal Medicine: Therapeutic Potential in Rheumatoid Arthritis. J. Ethnopharmacol. 279, 114368. doi:10.1016/j.jep.2021.114368

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, T., Li, S., Sun, Y., Pi, Z., Liu, S., Song, F., et al. (2017). Systematically Characterize the Absorbed Effective Substances of Wutou Decoction and Their Metabolic Pathways in Rat Plasma Using UHPLC-Q-TOF-MS Combined with a Target Network Pharmacological Analysis. J. Pharm. Biomed. Anal. 141, 95–107. doi:10.1016/j.jpba.2017.04.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang AH, A. H., Yu, J. B., Sun, H., Kong, L., Wang, X. Q., Zhang, Q. Y., et al. (2018). Identifying Quality-Markers from Shengmai San Protects against Transgenic Mouse Model of Alzheimer's Disease Using Chinmedomics Approach. Phytomedicine 45, 84–92. doi:10.1016/j.phymed.2018.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, L., and Wei, W. (2020). Anti-inflammatory and Immunoregulatory Effects of Paeoniflorin and Total Glucosides of Paeony. Pharmacol. Ther. 207, 107452. doi:10.1016/j.pharmthera.2019.107452

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, W., Gong, L., Zhou, L. L., Shan, J. J., Chen, L. T., Hui-Qin, X. U., et al. (2014). Effect of different compatibility of Daphnes Giraldii Cortex and Glycyrrhizae Radix et Rhizoma on adjuvant-induced arthritis in rats. Chin. Tradit. Herbal Drugs 45, 1418. doi:10.7501/j.issn.0253-2670.2014.10.013

CrossRef Full Text | Google Scholar

Zhang, A., Fang, H., Wang, Y., Yan, G., Sun, H., Zhou, X., et al. (2017). Discovery and Verification of the Potential Targets from Bioactive Molecules by Network Pharmacology-Based Target Prediction Combined with High-Throughput Metabolomics. RSC Adv. 7, 51069–51078. doi:10.1039/c7ra09522h

CrossRef Full Text | Google Scholar

Zhang, A.-H., Sun, H., Yan, G.-L., Han, Y., Zhao, Q.-Q., and Wang, X.-J. (2019). Chinmedomics: A Powerful Approach Integrating Metabolomics with Serum Pharmacochemistry to Evaluate the Efficacy of Traditional Chinese Medicine. Engineering 5, 60–68. doi:10.1016/j.eng.2018.11.008

CrossRef Full Text | Google Scholar

Zhang, A. H., Ma, Z. M., Kong, L., Gao, H. L., Sun, H., Wang, X. Q., et al. (2020). High-throughput Lipidomics Analysis to Discover Lipid Biomarkers and Profiles as Potential Targets for Evaluating Efficacy of Kai-Xin-San against APP/PS1 Transgenic Mice Based on UPLC-Q/TOF-MS. Biomed. Chromatogr. 34 (2), e4724. doi:10.1002/bmc.4724

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang T, T., Bai, G., Han, Y., Xu, J., Gong, S., Li, Y., et al. (2018). The Method of Quality Marker Research and Quality Evaluation of Traditional Chinese Medicine Based on Drug Properties and Effect Characteristics. Phytomedicine 44, 204–211. doi:10.1016/j.phymed.2018.02.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Q., Gao, X., Yan, G., Zhang, A., Sun, H., Han, Y., et al. (2020). Chinmedomics Facilitated Quality-Marker Discovery of Sijunzi Decoction to Treat Spleen Qi Deficiency Syndrome. Front. Med. 14 (3), 335–356. doi:10.1007/s11684-019-0705-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, Y., Guo, Z., He, W., Yang, Y., Li, Y., Zheng, A., et al. (2012). Ephedrine Hydrochloride Protects Mice from LPS challenge by Promoting IL-10 Secretion and Inhibiting Proinflammatory Cytokines. Int. Immunopharmacol. 13 (1), 46–53. doi:10.1016/j.intimp.2012.03.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, L., Wei, H., Wu, Y., Yang, S., Xiao, L., Zhang, J., et al. (2012). Licorice Isoliquiritigenin Suppresses RANKL-Induced Osteoclastogenesis In Vitro and Prevents Inflammatory Bone Loss In Vivo. Int. J. Biochem. Cel Biol. 44 (7), 1139–1152. doi:10.1016/j.biocel.2012.04.003

CrossRef Full Text | Google Scholar