Abstract
Introduction
Rheumatoid arthritis (RA) is an autoimmune disease characterized by the presence of aggressive arthritis. The basic pathological changes in RA involve chronic inflammation of joint synovia, formation of pannus, and destruction of cartilage caused by inflammatory cytokines and mediators. 1 Clinically, RA patients often suffer from joint stiffness, swelling, pain, and other symptoms, which may lead to joint deformity and loss of function in severe cases. 2 At present, RA is mainly treated with non-steroidal anti-inflammatory drugs (NSAIDs), corticosteroids, disease-modifying anti-rheumatic drugs, and biological agents in the field of Western medicine. Although these drugs can alleviate the clinical symptoms of RA to a certain extent, they are prone to adverse reactions when taken in large quantities for a long time, resulting in intolerance and unsatisfactory long-term effect.3,4
Previous studies have proved that there are pathological changes in blood coagulation and fibrinolysis indexes of RA patients, which are clinically manifested as a hypercoagulable state.5‐7 Coagulation and fibrinolysis indexes mainly include D-dimer (D-D), fibrinogen (FBG), and platelet (PLT) count. The coagulation and fibrinolysis indexes FBG and D-D in RA patients were significantly higher than those in the healthy control group.8,9 The abnormal changes in coagulation and fibrinolysis indexes in RA patients are also closely related to their immune and inflammatory indexes.10,11 Therefore, the abnormality of the coagulation and fibrinolysis system may be closely related to the occurrence and development of inflammation in RA patients.
Xinfeng capsule (XFC; Patent No.: Z20050062), a traditional Chinese medicine (TMC) compound, is composed of Huangqi (Radix Astragali Mongolici), Yiyiren (Semen Coicis), Wugong (Scolopendra), and Leigongteng (Radix et Rhizoma Tripterygii). It has been applied in clinical practice for 30 years with remarkable curative effects. 12 In the early stages, our team conducted systematic studies on the fingerprint and pharmacokinetics of XFC. Relatively complete chromatographic information was obtained, and the relative peak area ratio of each substance was within a certain range, which was the unique rule of XFC. In terms of the fingerprint, it could not only reflect the type and quantity of chemical components in XFC, but also describe and evaluate XFC to control its quality in an all-round way. Therefore, it can be used as one of the quality control methods of XFC. 13
Apriori algorithm is a popular approach used in data mining and association rule analysis. 14 Network pharmacology has been applied to reveal the network characteristics of drugs, indicating that multiple drugs and targets may have better clinical efficacy and less toxic side effects compared to single-target drugs. Based on Hospital Information System (HIS) and data mining, this study focused on the changes of immune-inflammation indexes and coagulation indexes of RA patients who received TCM combined with XFC treatment, and explored the potential targets, main active ingredients and core targets of XFC in the treatment of RA through network pharmacology and molecular docking technology. This study opens up new ideas for the follow-up experimental study of XFC in the treatment of RA and provides a scientific basis for its clinical application.
Materials and Methods
Materials
Hospitalization data were collated from RA inpatients with blood stasis syndrome before October 2021 in the Department of Rheumatology and Immunology of The First Affiliated Hospital of Anhui University of Chinese Medicine. The dataset contained records of the use of TCM, formulations of XFC, and disease-related laboratory indicators, including immune-inflammation indexes, such as C-reactive protein (CRP), erythrocyte sedimentation rate (ESR), rheumatoid factor (RF), complement component 3 (C3), and C4, and coagulation indexes, such as PLT, FBG, thrombin time (TT), prothrombin time (PT), and activated partial thromboplastin time (APTT). This study was approved by the Ethics Research Committee of the First Affiliated Hospital of Anhui University of Traditional Chinese Medicine. A total of 524 patients were enrolled and treated with XFC.
Methods
Cluster Analysis
The designation for the use of Chinese herbal medicine was 1 while that for nonuse was 0. Systematic clustering in SPSS v. 21.0 (IBM Corp.) was applied to investigate the compatibility of the Chinese herbal medicine. In the clustering analysis algorithm, each herb was regarded as a cluster. According to the similarity between objects, N clusters were combined to form a new class. Euclidean metric was used to calculate the similarity between herbs
14
:
Association Rules
Apriori algorithm: The designation for the use of Chinese herbal medicine was 1 while that for nonuse was 0. The Apriori module in SPSS Clementine v. 11.1 (IBM Corp.) was employed to identify the correlations among Chinese herbal medicines. The minimum support and confidence were set to 80% and the degree of improvement to >1. The Apriori algorithm was implemented to establish the relationships among items within a dataset. It is also known as a shopping blue analysis. In this dataset, each drug was treated as a variable. The formulae applied were as follows
15
:
Random Walking
The random walking model is based on the therapeutic effect of each patient. When the scale of patients reaches a certain level, random walking routes can be formed to reflect the improvement of the overall efficacy of the treatment. Oracle Developer Suite 10 g was utilized to evaluate the random walking model of international BioMed immune-inflammation indexes, and the improvement of drug compatibility of laboratory indicators was observed. The formula used was as follows
16
:
Screening of Active Ingredients of XFC and Acquisition of Targets
The pharmacology database of the TCM system was used to search the chemical composition and pharmacodynamic parameters of each herb. XFC mainly consists of Huangqi (Radix Astragali Mongolici), Yiyiren (Semen Coicis), Wugong (Scolopendra), and Leigongteng (Radix et Rhizoma Tripterygii). The Absorption, Distribution, Metabolism, and Excretion (ADME) model was utilized to screen the active ingredients of XFC more efficiently and economically. Virtual screening of potential active compounds of each herb in XFC was conducted using 3 parameters, oral bioavailability (OB), druglikeness (DL), and half-life (HL), in which the OB screening threshold was set to ≥30% and the DL ≥ 0.18.
Acquisition of XFC Active Components—RA Common Targets
With “rheumatoid arthritis” as the keyword, we searched GeneCards, Online Mendelian Inheritance in Man (OMIM), DrugBank, and other databases to obtain all known targets related to RA. In addition, relevant target names were standardized through UniProt database. FunRich software was used to conduct Venn analysis on the main targets of XFC and RA-related targets to obtain the potential targets of XFC in the treatment of RA.
Construction of Protein-Protein Interaction Network
A network diagram of XFC active component—potential target—disease was constructed using Cytoscape3.7.2 software. The common targets of XFC and RA were imported into STRING online platform, with “multiple proteins” chosen as the analysis model and the species defined as “Homo sapiens”, to obtain the information of protein–protein interaction (PPI) and construct a PPI interaction network. The downloaded PPI network files were imported into Cytoscape, and 5 core targets were obtained using Maximal Clique Centrality (MCC) algorithm in the cytoHubba plug-in.
Gene Ontology and Kyoto Encyclopedia of Genes and Genomes Enrichment Analyses
Gene Ontology (GO) function enrichment analysis was performed on the common targets of XFC-RA through the Database for Annotation, Visualization and Integrated Discovery (DAVID), mainly including biological process (BP), cellular component (CC), and molecular function (MF). The results were sorted by -log10 (P value). Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was carried out to examine the key role of XFC in RA treatment. According to the size of the P value, the top 10 signal pathways were selected for visual analysis.
Statistical Processing
All data were analyzed in SPSS v. 21.0 (IBM Corp.). No-parametric tests were performed on 2 related samples before and after the treatment for the control group and the experimental group. Mann–Whitney U test was used to identify the differences between groups before and after the treatment. The difference was considered statistically significant at P < .05.
Results
Changes of Immune-Inflammatory Indexes before and after the Treatment
The immune inflammatory reaction in RA patients can activate the coagulation and fibrinolysis system, increase the levels of FBG and PLT, and lead to a hypercoagulable state. Our results demonstrated that the immune-inflammatory indexes of RA patients, including high-sensitivity (Hs)-CRP, ESR, RF, C3 and C4, and coagulation indexes, including FBG, PLT, TT, PT, and APPT, were decreased significantly after treatment compared with those before treatment (P < .01) (Table 1), indicating that the inflammatory response and hypercoagulable state of RA patients were effectively improved.
Changes of Immune-Inflammatory Indexes and Coagulation Indexes Before and After the Treatment.
Abbreviations: APTT, activated partial thromboplastin time; CRP, C-reactive protein; C3, complete component 3; ESR, erythrocyte sedimentation rate; FBG, fibrinogen; Hs, high sensitivity; PT, prothrombin time; PLT, platelet; RF, rheumatoid factor; TT, thrombin time.
Correlation of XFC in Immune-Inflammation and Coagulation Indexes
The latter item was set as the indicator and the former as XFC. The minimum support was set to 60% and the minimum confidence to 70%. The following results were obtained in the Apriori module analysis of the association rules between XFC and improvement of RF, Hs-CRP, ESR, C3, C4, FBG, and PLT: Support > 60%, Confidence > 70%, and Lift > 1, indicating that XFC had a strong correlation with immune inflammation and coagulation indicators (Table 2).
Correlation of Xinfeng Capsule (XFC) in Immune-Inflammatory Indicators.
Abbreviations: CRP, C-reactive protein; C3, complete component 3; ESR, erythrocyte sedimentation rate; FBG, fibrinogen; Hs, high sensitivity; LHS, left-hand side; PLT, platelet; RF, rheumatoid factor; RHS, right-hand side.
Evaluation of Immune-Inflammation and Coagulation Indicators by Random Walking Model
After XFC treatment, 675 comprehensive evaluations of RF were recorded, and its improvement coefficient was 0.330, while its clinical significance was that 5.130 steps were required for each improvement in the comprehensive index. C3 had 533 records in the comprehensive evaluation, and its improvement coefficient was 0.215, while its clinical significance was that 9.050 steps were required for each improvement in the patient comprehensive index. C4 had 533 records in the comprehensive evaluation, and its improvement coefficient was 0.268, while its clinical significance indicated that 7.280 steps were required for each improvement of the comprehensive index. HS-CRP had 862 records in the comprehensive evaluation, and its improvement coefficient was 0.354, while the clinical significance suggested that 4.540 steps were required for each improvement in the patient comprehensive index. ESR had 794 records in the comprehensive evaluation, and its improvement coefficient was 0.353, while the clinical significance suggested that 4.710 steps were required for each improvement in the patient comprehensive index. The comprehensive evaluation records of FBG were 151, and its improvement coefficient was 0.133, while the clinical significance was that 33.75 steps were required for each improvement in the patient comprehensive index. PLT showed 775 records in terms of the comprehensive evaluation and 0.088 in the improvement coefficient, while its clinical significance was that 19.100 steps were required for each improvement in the patient comprehensive index. TT had 152 records in the comprehensive evaluation and showed 0.138 in the improvement coefficient, while its clinical significance was that 32.190 steps were required for each improvement in the patient comprehensive index. PT had 152 records in the comprehensive evaluation, and its improvement coefficient was 0.118, while its clinical significance suggested that 37.560 steps were required for each improvement in the patient comprehensive index. APTT had 151 records in the comprehensive evaluation and 0.205 in the improvement coefficient, and its clinical significance was that 21.77 steps were required for each improvement in the patient comprehensive index (Table 3) (Figure 1).

Random walking model of immune-inflammatory indices in rheumatoid arthritis (RA) patients.
Random Walking Model of Immune-Inflammatory and Coagulation Indicators.
Abbreviations: APTT,#activated partial thromboplastin time; CRP, C-reactive protein; C3, complete component 3; ESR, erythrocyte sedimentation rate; FBG, fibrinogen; Hs, high sensitivity; PT, prothrombin time; PLT, platelet; RF, rheumatoid factor; TT, thrombin time.
Active Components of XFC
By searching the traditional chinese medicine systems pharmacology database and analysis platform (TCMSP) database and reviewing relevant literature, 92 chemically active components in XFC formula were obtained. Among them, there were 20 effective components of Huangqi (Radix Astragali Mongolici), 12 of Yiyiren (Semen Coicis), 51 of Leigongteng (Radix et Rhizoma Tripterygii), and 9 of Wugong (Scolopendra).
Action Targets of XFC and RA and Their Common Targets
A total of 197 targets corresponding to XFC active ingredients were retrieved from the TCMSP database. Moreover, 2810, 27, 571, and 169 RA-related targets were collected, respectively, by searching GeneCards, OMIM, DrugBank, and TTD databases. A total of 152 RA-related targets were eventually obtained after the removal of duplicate genes. Venn analysis was performed on 197 targets corresponding to XFC active ingredients and 152 RA-related targets using RData software. Thirty intersected genes were found which were the potential targets of XFC in the treatment of RA (Figure 2).

Common targets of Xinfeng Capsule (XFC) and rheumatoid arthritis (RA).
Target Prediction of Active Ingredients of XFC
Cytoscape3.7.2 software was used to construct a network diagram of “Chinese Medicine – Active Ingredient – Potential Target”, in which pink represents the drugs in XFC, yellow represents an active ingredient of Chinese medicine, and purple represents the potential target (Figure 3).

Active components of Xinfeng Capsule (XFC) targets network diagram.
Construction of PPI Network
The intersection genes were imported into the STRING database to draw the PPI network diagram, as shown in Figure 4. The downloaded PPI network files were imported into Cytoscape, and 5 core targets were obtained using the MCC algorithm in the cytoHubba plug-in (Figure 5).

Protein–protein interaction (PPI) network of Xinfeng Capsule (XFC) in the treatment of rheumatoid arthritis (RA).

Core hub targets (top 5).
GO Functional Analysis and KEGG Pathway Enrichment Analysis
In this study, the enrichment analysis of the abovementioned PPI protein interaction network targets was completed through the DAVID database, suggesting that the overall treatment process had multiple signaling pathways that took part in the regulation of various BPs. The results of GO functional enrichment analysis showed that the key targets of XFC in the treatment of RA were mainly enriched in biological reactions, such as inflammation, apoptosis, I-kappaB phosphorylation, and immune reaction. KEGG signaling pathway enrichment analysis manifested that the related targets were significantly enriched in interleukin-17 (IL-17), tumor necrosis factor (TNF), nuclear factor-κB (NF-κB), and other signaling pathways (Figures 6 and 7).

Gene Ontology (GO) enrichment analysis.

Kyoto Encyclopedia of Genes and Genomes (KEGG) signal pathway bubble map (top10).
Molecular Docking of key Targets
To conduct an in-depth investigation into the key molecular mechanism of XFC in the treatment of RA, we selected 10 active components as ligands and 5 core targets as receptors. The docking results are shown in Figure 8. The abscissa represents the active component of XFC, and the ordinate represents the core target. The color depth indicates the binding energy. A binding energy <−7.0 kcal·mol−1 indicated a strong binding activity and a negative correlation between the docking effect and binding energy. The results unveiled that the minimum binding energy obtained by β-sitosterol and stigmasterol docking with TNF and prostaglandin-endoperoxide synthase 2 (PTGS2) was less than −8.0 kJ mol−1, andβ-sitosterol had a better binding ability to TNF. In addition, PyMOL was used to visualize the docking results ofβ-sitosterol and TNF, as shown in Figure 9. Therefore, these 2 compounds and 2 proteins may be the key compounds and core targets in the treatment of RA by XFC, which may also provide a molecular biological basis for further experimental verification.

Heat map of the minimum binding energy of docking of active compounds with core targets.

Molecular docking diagram of β-sitosterol and tumor necrosis factor (TNF).
Discussion
Synovial vasculitis is regarded as the major pathological change in RA patients. In vivo inflammation can activate the body’s coagulation and fibrinolysis system and increase the levels of FBG and PLT, leading to hypercoagulability. Therefore, it is of great value to clarify the hypercoagulability of patients for the treatment of RA. FBG is a protein synthesized by the liver, which can result in PLT aggregation in the body and induce red blood cell adhesion and thrombosis when it plays a role. The increase of this index indicates that the human body is in a hypercoagulable state. 17 According to the study carried out by Lu Xuedan, 18 RA patients had high serum D-D activity and their FBG levels were higher compared to the normal group due to the high activity group joint inflammation. Joint inflammation can induce prothrombin and FBG, such as plasma components, into the joint cavity, while the joint cavity inflammation factor could lead to the activation of blood coagulation responses. Moreover, the higher the RA activity is, the more severe is the inflammatory response. Meanwhile, the inflammatory response can stimulate the synthesis of large amounts of FBG in the liver, which further participates in the formation of a microthrombus and activation of the coagulation-fibrinolytic system, resulting in a vicious cycle. As one of the tangible components in blood, PLTs play an important role in physiological and pathological processes, such as coagulation, inflammatory reaction, and thrombosis. 19 Modern studies have shown that an elevated PLT count is one of the main characteristics of RA, and RA patients have a significantly higher PLT count than that of normal subjects.9,20
In this study, laboratory indicators of clinical RA patients were observed before and after the treatment based on the HIS system. According to our results, after the treatment, the inflammatory indicators, including RF, ESR, C3, C4, and Hs-CRP, improved significantly, while the coagulation indicators, including PLT, FBG, TT, PT, and APTT, improved notably. Association rule analysis results showed that XFC was strongly correlated with the improvement of RF, Hs-CRP, ESR, C3, C4, FBG, and PLT indexes, with Support > 70%, Confidence > 70%, and Lift > 1, as shown in Table 2. According to the data mining results, XFC had the effect of improving immune-inflammatory and coagulation indicators in RA patients and regulate the immune-inflammatory response.
XFC is a TCM compound. Our team's previous research results show that XFC accelerates the apoptosis of synovial fibroblasts by regulating long noncoding RNA MAPKAPK5-AS1, thus repressing the response of RA. 21 XFC can improve the symptoms of RA patients, improve cardiopulmonary function, and reduce the parameters of anemia and PLT factor. 22 XFC can reduce the immune inflammatory response to improve the extra-articular lesions in RA patients.23‐27 XFC can regulate the level of B cell activating factor and its receptor in active RA and increase the level of immunoglobulin in RA patients. 28 XFC improves blood stasis in patients with RA through the Actl/NF-κB signaling pathway. 29 XFC has played a good role in the treatment of RA. However, due to the large number of chemical active ingredients in the TCM and the lack of molecular-level research, the mechanism of its treatment process is not clear. Based on network pharmacology and molecular docking technology, this study further analyzed the active ingredients, potential targets, and pathways of XFC in the treatment of RA, and studied its molecular mechanism of action.
According to the results of network pharmacology analysis, the main target proteins involved in the treatment of RA by XFC, such as PTGS2, CASP3, TNF, AKT1, and JUN, were closely related to inflammation, tumor, and immunity. Among them, PTGS2 and TNF were the active target proteins with optimal binding in molecular docking. Moreover, AKT1 is involved in the regulation of cell proliferation, invasion, and differentiation, as well as tissue fibrosis. The transforming growth factor-β1 (TGFβ1) can bind to the fibroblast membrane receptors to activate the PI3K/Akt pathway, leading to phosphorylation and activation of the AKT factor. 30 As an important cytokine, TNF has multiple biological functions, which can be involved in inducing inflammatory and immune responses, cell survival and apoptosis. 31
Through screening the effective compounds of XFC in the treatment of RA in the early stage and verifying by molecular docking technology, our study found that β-sitosterol and stigmasterol, the active ingredients of XFC, had the strongest binding with TNF and PTGS2 in molecular docking, suggesting that XFC may play a critical role in the treatment of RA. Research shows thatβ-sitosterol has neuroprotective, antimutagenic, anti-diabetes, and antifungal effects. 32 Its effect is especially related to the reduction of oxidative stress in cardiovascular diseases, aging, and cancer. 33 Stigmasterol can also remarkably reduce the mRNA expressions of cyclooxygenase-2 (COX-2) and nitric oxide synthase (iNOS), which are inflammatory factors stimulated by lipopolysaccharide (LPS), and has an anti-inflammatory effect by reducing PGE2 and NO release. 34
Conclusion
This study elucidated that the immune-inflammatory indicators and coagulation indicators in RA patients were effectively improved after XFC treatment, suggesting the critical role of XFC in the improvement of immune inflammatory response. In addition, the molecular mechanism of XFC in the treatment of RA was preliminarily discussed by adopting the network pharmacology method, which demonstrated that the key active components of XFC, such as β-sitosterol and stigmasterol, mainly acted on TNF, PTGS2, CASP3, AKT1, JUN, and other targets. Through IL-17, TNF, NF-κB, and other signaling pathways, the drug participated in the activity of transcription factors and immune response, thus contributing to anti-PLT aggregation, inhibition of inflammatory response, and anti-myocardial ischemia. Overall, the research results can comprehensively reflect the synergistic characteristics of Chinese medicinal compounds with multiple components, targets, and pathways. However, due to the limitations of network pharmacology analysis, this study only studied the mechanism of action of XFC in treating RA from the theoretical level. Therefore, relevant pharmacological experiments are needed for verification so as to further clarify the pharmacological mechanism of XFC in treating RA and provide new research ideas for the prevention and treatment of RA in clinical practice.
Footnotes
Acknowledgments
This research was supported by the following projects: National Nature Fund Program (81973655), the University Synergy Innovation Program of Anhui Province (GXXT-2020-025), Anhui Famous Traditional Chinese Medicine Liu Jian Studio Construction Project (Traditional Chinese Medicine Development Secret [2018] No. 11), The 12th Batch of “115” Innovation Team in Anhui Province (Anhui Talent Office [2019] No. 1), Anhui Natural Science Youth Foundation (2008085QH386); The Key Laboratory of Xin’an Ministry of Medical Education (No.2020xayx10), National Key Research and Development Project (No. 2018YFC1705203), and Collaborative research project of Chinese and Western medicine for major and difficult diseases in Anhui Province (TCM Development Secretary [2021] No. 70).
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Nature Fund Program (grant number 81973655).
