Next Article in Journal
Antimelanogenic Effects of Polygonum tinctorium Flower Extract from Traditional Jeju Fermentation via Upregulation of Extracellular Signal-Regulated Kinase and Protein Kinase B Activation
Previous Article in Journal
Cellular Mechanisms Responsible for Success and Failure of Bone Substitute Materials
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 19
Issue 10
10.3390/ijms19102894
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Inhibitory Effect of Methotrexate on Rheumatoid Arthritis Inflammation and Comprehensive Metabolomics Analysis Using Ultra-Performance Liquid Chromatography-Quadrupole Time of Flight-Mass Spectrometry (UPLC-Q/TOF-MS)

by
Zhiqiang Pang
1,
Guoqiang Wang
1,
Nan Ran
1,
Hongqiang Lin
2,
Ziyan Wang
1,
Xuewa Guan
1,
Yuze Yuan
1,
Keyong Fang
3,
Jinping Liu
2 and
Fang Wang
1,*
1
Department of Pathogen Biology, College of Basic Medical Sciences, Jilin University, Changchun 130021, China
2
Research Center of Natural Drug, School of Pharmaceutical Sciences, Jilin University, Changchun 130012, China
3
Department of Pharmacology, College of Basic Medical Sciences, Jilin University, Changchun 130012, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2018, 19(10), 2894; https://doi.org/10.3390/ijms19102894
Submission received: 19 August 2018 / Revised: 20 September 2018 / Accepted: 21 September 2018 / Published: 23 September 2018
(This article belongs to the Section Molecular Pathology, Diagnostics, and Therapeutics)
Browse Figures
Versions Notes

Abstract

:
Rheumatoid arthritis (RA) is a common autoimmune disease. The inflammation in joint tissue and system endanger the human health seriously. Methotrexate have exhibited a satisfactory therapeutic effect in clinical practice. The aim of this research was to establish the pharmacological mechanism of methotrexate on RA therapy. Collagen induced arthritic rats were used to identify how methotrexate alleviates inflammation in vivo. Lipopolysaccharide-induced inflammatory proliferation in macrophages was also be detected in vitro. The activation level of Nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) and Nucleotide binding domain and leucine-rich repeat pyrin 3 domain (NLRP3)/Caspase-1 and related cytokines were examined by real-time PCR and western blotting or quantified with the enzyme-linked immunosorbent assay. Comprehensive metabolomics analysis was performed to identify the alteration of metabolites. Results showed that treating with methotrexate could alleviate the inflammatory condition, downregulate the activation of NF-κB and NLRP3/Caspase-1 inflammatory pathways and reduce the level of related cytokines. Docking interaction between methotrexate and caspase-1 was visualized as six H-bonds indicating a potential inhibitory effect. Metabolomics analysis reported three perturbed metabolic inflammation related pathways including arachidonic acid, linoleic acid and sphingolipid metabolism. These findings indicated that methotrexate could inhibit the onset of inflammation in joint tissue by suppressing the activation of NF-κB and NLRP3/Caspase-1 pathways and regulating the inflammation related metabolic networks.

Graphical Abstract

1. Introduction

Rheumatoid arthritis (RA) is a chronic inflammatory joint disease that is characterized by synovial proliferation, joint destruction, and systemic inflammation [1,2]. It affects about 1–2% adult population globally [3]. RA begins with an autoimmune phase and is followed by an inflammation, which is mediated by the citrullinated antigen [4,5,6]. Despite the autoimmune reaction remains clinically silent, inflammation has the potential to cause severe pain and swelling [5].
The pathogenesis of RA may be associated with the involvement of environmental, infectious, and genetic factors [4]. Nuclear factor (NF)-κB activation had been demonstrated to be a crucial central role in the development and progression of RA because it not only participates in the production of many pro-inflammatory mediators but also regulates osteoclast genesis [7]. Besides, there also exists a high transcriptional and activated level of Nucleotide binding domain and leucine-rich repeat pyrin 3 domain (NLRP3) inflammasome in patients with RA [8]. The caspase-1 activation induced by NLRP3 inflammasome has been reported to be involved in RA [9]. Many inflammatory mediators contribute together to the pathogenesis of RA, despite the mechanism being still incompletely elucidated.
Metabolomics, as an emerging high-throughput technique for comprehensive analysis of endogenous small metabolites (<1 kDa) in biological system has been widely used to identify novel biomarkers and explore the potential pharmacological targets [10]. It usually focuses on the metabolic profile in vivo and directly shows us the metabolic changes that are associated with the pathological process of disease. Previous metabolomics studies have shown the potential of metabolic profiling to diagnose RA [11]. Recently, it has been applied extensively to investigate the pathogenesis of RA in patients and experimental arthritis models [11,12]. In addition, the pharmacological principles of many multi-targets medicine have also been clarified with non-targeted metabolomics techniques.
Conventional strategy of medical therapy for RA involves the Nonsteroidal anti-inflammatory drugs (NSAIDs), such as diclofenac [13]. But the side-effect and resistance of these pharmaceuticals usually cause a chronic pain and poor life quality and require a larger amount of resources for further management [14]. Methotrexate as an anti-folic acid metabolism medicine is widely used to inhibit the growth of carcinomas such as acute leukemia by blocking dihydrofolate reductase. As a result, humoral immunity is inhabited after the treatment [15]. Despite the patients with variant genotypes showed variable responsiveness towards methotrexate, it is still considered as a standard first-line medicine for RA in clinic practice and it has achieved a satisfactory therapeutic effect in over 60 percent RA cases [16]. However, the detailed mechanism behind this therapeutic effect has not been clarified completely yet.
In the present study, we established RA murine models and assessed the therapeutic effect of methotrexate in vivo and in vitro. At the same time, the inflammatory responses, mainly the activation of the typical NLRP3/Caspase-1 and NF-κB pathways in the whole joint tissue and macrophages were investigated. We tested the hypothesis that the treatment with methotrexate will alleviate the inflammatory condition by down-regulating the NF-κB and NLRP3/Caspase-1 pathways. Besides, the global view of metabolic profiles after the treatment of methotrexate was also identified. The molecular and protein interaction that is based on molecular dynamics virtual docking will contribute to the further pharmacological exploration.

2. Results

How methotrexate alleviate inflammation was studied in collagen induced arthritic (CIA) murine model and Lipopolysaccharide (LPS)-induced macrophages in present study.

2.1. Inhibitory Effect of Methotrexate on Murine CIA Paw Swelling

Swelling of CIA paws showed a statistical significance at 11th day according to the visual scores after the booster injection. Paw swelling volume was examined for all experimental animals then. As shown in Figure 1A, the experimental arthritis groups that were injected with collagen emulsion showed a significant increase in paw swelling volume when compared with the control (Control) group (p < 0.05). Drugs administration began at day 11 and continued until the visual scores of the CIA model (Model) group came to a platform stage. As expected, the endpoint of the administration came at day 25. The low dosage (Low) group showed a trend of inhibitory effect on paw swelling when compared with the model group (p = 0.14), as shown in Figure 1B,C. However, middle dosage (Middle), high dosage (High), and positive drug diclofenac sodium (DCS, Positive) group showed a significantly reduced swelling volume and visual scores, as indicted by Figure 1B–D.

2.2. Downregulation of Methotrexate on NF-κB/NLRP3 Pathways in Murine Model

NF-κB and NLRP3/Caspase-1 pathways mRNA expression was remarkably higher in CIA ankle joint tissue in Model compared with Control indicating a significant activation (Figure 2). In detail, NF-κB p65 and TNF-α gene expression in all the drug administrated groups decreased significantly (p < 0.05). Level of NLRP3 gene expression in Middle and High group is different from Model group significantly (p < 0.01). The mRNA level of Caspase-1 in drug-administrated groups showed an obvious decrease compared with the model animals (p < 0.01) in a dose-independent manner. However, IL-1β and IL-18 gene expression only decreased significantly in Middle, High, and Positive groups (p < 0.01). The IL-1β and IL-18 gene expression in the Low group simply showed a trend of decrease.
NLRP3 protein in Model rats increased significantly (p < 0.01). So, as the protein expression level of caspase-1 and NF-κB p65 in Model (p < 0.05), as displayed in Figure 3. All of the drug administrated groups showed a downregulated expression of NLRP3, caspase-1, and NF-κB p65 protein when compared with Model. Especially, they were downregulated remarkably in Middle, High and Positive group compared with Model (p < 0.05 or p < 0.01). Moreover, NF-κB p65 protein in Low group also showed a statistical significance (p < 0.05).
Besides, the quantification assay results of IL-1β, IL-18, and TNF-α in plasma were represented in Figure 4. The NF-κB and NLRP3/Caspase-1 pathways related cytokines in plasma of Model increased significantly (p < 0.01 vs. Control) and decreased in all the drug-administrated groups (p < 0.05 for Low vs. Model and p < 0.01 for Middle, High, and Positive groups vs. Model, respectively).

2.3. Methotrexate Suppressed the Inflammatory Proliferation of Macrophages

The cytotoxicity of both methotrexate and DCS on RAW 264.7 was monitored with Real Time Cellular Analysis (RTCA) system and displayed as Figure 5. After a continuous incubation with different drug concentration, the cytotoxicity was evaluated at 72 h for methotrexate and at 96 h for DCS in case of proliferation recession. As shown in Figure 5A1,A2, methotrexate significantly inhibited the cell proliferation over the range of 1 μM to 1 mM indicating a cytotoxic effect (p < 0.01). DCS did not exhibit any suppression on the proliferation with all drug concentrations (Figure 5B, p > 0.05). The non-cytotoxic concentrations were set as 50, 20, 10, 5, 2 nM for methotrexate and 1 mM for DCS to test their anti-inflammatory effects. As a result, compared to the Model, the Control cells exhibited a similar inflammatory proliferative status (p > 0.05). DCS and methotrexate in 50, 20 and 10 nM significantly reduced the inflammatory proliferation (p < 0.01). However, there did not exist any difference on proliferative characteristics compared to the Model after treating with methotrexate of 5 and 2 nM (p > 0.05).

2.4. Methotrexate Inhibited NF-κB and NLRP3/Caspase-1 Pathways in Macrophages

The mRNA expression level of all related genes in RAW 264.7 inflammatory cellular model after the treatment was shown in Figure 6. All genes including were remarkably upregulated in inflammatory cellular models (p < 0.01). After the treatment with methotrexate, the expression of all investigated genes was significantly inhibited in dose-dependent manner (p < 0.01 vs. Model). However, methotrexate with 2 and 5 nM only displayed a potential to downregulate the expression of NF-κB p65 and IL-18 genes (p > 0.05 vs. Model).
Consistent with the results in CIA murine model, proteins including NLRP3, Caspase-1 and NF-κB p65 in inflammation cellular Model rose sharply compared with Control (p < 0.01). Their levels decreased significantly after the therapy with methotrexate and DCS in a dose-independent manner for NLRP3 and Caspase-1 or in a dose-dependent manner for NF-κB p65 (p < 0.01 vs. Model) as shown in Figure 7. However, the pathways related cytokines including IL-1β, IL-18 and TNF-α increased significantly in LPS induced cellular inflammation Model (p < 0.01 vs. Control) and decreased in a dose-dependent manner (p < 0.01 vs. Model, Figure 8).

2.5. Multivariate Analysis and Differentiated Metabolic Patterns

Metabolomics analysis on the metabolic perturbation after the treatment was performed with Ultra-Performance Liquid Chromatography-Quadrupole Time of Flight-Mass Spectrometry (UPLC-Q/TOF-MS) systems. The stability and suitability of spectrum and chromatogram was evaluated based on the 10 ions in both ESI+ and ESI− modes from six quality control (QC) sample injections respectively. All relative standard deviation (RSD) values including repeatability and intermediate precision were calculated as less than 6.0% (Table S1), which indicated a good confidence for subsequent analysis.
Principle component analysis (PCA) as an unsupervised lowering-dimension pattern recognition model revealed the inherent similarity among the different samples in ESI+ and ESI− modes visually (Figure 9A1,B1). Sample from different groups displayed as intrinsic clusters indicating a lager similarity within them under both ESI+ (R2 = 64.05%, Q2 = 47.80%) and ESI− (R2 = 54.68%, Q2 = 33.90%) modes.
4 orthogonal projections to latent structures discriminant analysis (OPLS-DA) models were also established with a supervise based on the experimental treatment groups. As shown in OPLS-DA plots of Figure 9, all samples in different models were appreciably separated from each other indicating there exists no extremely abnormal sample. As a result, OPLS-DA models established can be used to identify the differentiated metabolites and potential biomarkers (R2 = 99.99%, Q2 = 78.32% for C vs. M in ESI+, R2 = 99.32%, Q2 = 82.46% for M vs. T in ESI+; R2 = 99.99%, Q2 = 45.02% for C vs. M in ESI−, R2 = 99.56%, Q2 = 97.10% for M vs. T in ESI−). Besides, the validation of all OPLS-DA models for differentiated metabolites analysis have been confirmed by cross validation ANONA analysis (p < 0.0001).

2.6. Global Profiles of Distinct Metabolites and Perturbed Metabolic Pathways

Distinct metabolites were selected according to the Variable Importance for the Projection (VIP) indexes that were calculated based on the OPLS-DA models above. Metabolites with VIP values more than 1.0 and p values were less than 0.05 were chosen for tandem mass spectrometry (MS/MS) tandem identification and metabolic pathway analysis. A total of 24 remarkably changed metabolites in plasma were determined with the details listed in Table 1. Comparison details were provided in Figures S1–S24.
All of the identified metabolites ions were also displayed in S-plots (Figure 10), which combined correlation and covariation values of the OPLS-DA models. The edging distribution of all the distinct metabolites in the loading plot revealed enough significance and confidence to illustrate the group difference within the OPLS-DA models. All of the potential biomarkers were evaluated with receiver operating characteristic curve (ROC), as shown in Figure S25. Metabolites with AUC > 0.9 were defined as high diagnostic accuracy. As a result, 19 differentiated metabolites were identified as potential biomarkers with high sensitivity and specialty.
Besides, the relative levels of the distinct metabolites were displayed with a heatmap after the mathematic LOG transformation (Figure 11A). Samples were naturally clustered into three major branches, which is consistent with the biological groups. It is noted that the drug-administrated samples share a large similarity on their metabolic patterns.
Eight metabolic pathways were identified to be involved (Table S3, Figure 11B) according the results from MetaboAnalyst 4.0 platform. Three of them, including arachidonic acid, linoleic acid, and sphingolipid metabolisms were significantly perturbed (Impacts > 0.1 and −log(P) > 5). The comprehensive interaction between inflammatory pathways and metabolic pathways were summarized and displayed in Figure 10C.

2.7. Potential Docking Mode of Methotrexate within Caspase-1 Protein

Molecular docking performed to investigate the potential interaction between methotrexate and the critical inflammatory cytokines catalyzation related enzyme, caspase-1 revealed an accurate docking site and the appropriate conformation. The precise modes of interaction and relative binding energy were measured as Table 2. Moreover, the number of H-bonds and the docking scores was calculated with the Schrödinger demo. As shown in Figure 12A,B, both the molecular conformation of methotrexate and the polar contacts (seven contacts, including six H-bonds and one other polar contacts) were displayed. The residue ARG-341 has the most contacts with the ligand.

3. Discussion

Collagen-induced arthritis (CIA) in rats is the most widely used model for RA because it shares many similar immunological and pathological features with RA in patients [17]. Innate immune components, especially macrophage is a prime source of pro-inflammatory cytokines and critically contribute to RA pathology [3,18]. In present study, we established the CIA murine model and the inflammatory macrophage cellular model for pharmacological research on the anti-rheumatic effect of methotrexate. We discovered that methotrexate could alleviate inflammatory swelling in CIA rats and inhibit the inflammatory proliferation of macrophages by suppressing NF-κB and NLRP3/Caspase-1 pathway and downregulate the level of pathway related cytokines. Besides, methotrexate has the potential to inhibit catalytic activity of caspase-1 directly. Three inflammation related metabolic pathways including arachidonic acid, linoleic acid and sphingolipid metabolism also participated in the pathogenesis of RA and were regulated significantly after the treatment with methotrexate.
NF-κB is responsible for the regulation of inflammation and immune responses and it is involved in the pathogenesis of RA. The damage-associated molecular patterns binding to Toll-like receptors, as the priming signal activate the upregulation of the NF-κB nucleus transfer [3,19]. The translocated free NF-κB then activates its target genes to transcribe various kinds of inflammatory mediators, such as proIL-1β, proIL-18, and NLRP3 [19,20]. NLRP3 protein was constructed as NLRP3 inflammasome assembly with other components including caspase recruitment domain containing protein 8, pro-caspase-1 and so on together. The activation of caspase1, which is mediated by the inflammasome further induces maturation of IL-1β and IL-18 in RA [21]. The similar activation of NF-κB and NLRP3/Caspase-1 pathways were also observed in CIA ankle joint tissues and LPS-induced macrophages in present study. A great deal of evidence has been accumulated and revealed the beneficial effect of direct inhibition of NF-κB on inflammatory disease [22,23]. Herein, methotrexate inhibited the activation of NF-κB and NLRP3/Caspase-1 pathways in joint tissue and decreased the level of systematic cytokines in a dose dependent manner. The activity-inhibiting binding of methotrexate within caspase-1 may have the potential to partly explain the therapeutic effect, but the detailed mechanism is not clear.
Metabolomics techniques as a high-throughput approach have been extensively utilized to investigate the pathogenesis of RA and therapeutic targets of anti-RA drugs [11,24]. In the present study, we discovered that arachidonic acid, linoleic acid, and sphingolipid metabolism perturbed and they were significantly regulated after administrating with methotrexate.
Arachidonic acid, a component that is released from cellular membrane would be metabolized as a variety of metabolites, including EETs, HETEs, prostaglandins, thromboxane, and leukotrienes by cyclooxygenases or lipoxygenases and promote the nuclear translocation of NF-κB, which is at least partly induced by eicosanoids [25,26]. At the same time, the NF-κB pathway is the most critical activator for NLRP3 inflammasome and the massive production of lipid mediators [27]. EETs are a series of lipid mediators with critical physiological functions that include anti-inflammation, anti-hypertension, and organ protection effects [28]. The anti-inflammation effect is associated with the regulation of NF-κB translocation, as described previously [29]. Methotrexate also increased the production of EETs and inhibited the activation of NF-κB in present study. Prostaglandins and leukotrienes, which was believed to cause pain and activate the formation of cytokine were reduced after the treatment of methotrexate. So did the level of thromboxane A2. HETEs as a class of lipids with distinct inflammation bioactive roles were also been influenced by methotrexate. In brief, methotrexate modulate the NF-κB and NLRP3/caspase-1 inflammasome-mediated arachidonic acid metabolic lipid profiles, inhibit the release of pro-inflammatory lipids and promote the production of anti-inflammatory mediators in CIA murine rats.
The precursor of arachidonic acid, linoleic acid can be converted to arachidonic acid and epoxides of linoleic acid (including 9,10-EpOME and 12,13-EpOME) [30]. The EpOMEs that are produced by neutrophils and macrophages are considered as leukotoxins and they contribute to the pathogenesis of inflammatory disease, such as acute respiratory distress syndrome [31]. We observed a significant increase in CIA rats and a decrease of EpOMEs after the treatment of methotrexate. Linoleic acid itself seems to have potent anti-inflammatory effect, which is associated with the diminished NF-κB binding activity and delayed translocation of NF-κB [32,33]. A lower level of linoleic acid reported before in RA serum is consistent with our results in CIA [34]. Linoleic acid also has the potential to ameliorate the high expression level of NLRP3 genes and reduced the release of mature IL-1β and other cytokines [35,36]. Despite that methotrexate could damage the dietary absorption of linoleic acid [37], it had already been reported to increase linoleic acid in carcinoma cell and was demonstrated in present study in CIA rats to play its anti-inflammation role [38].
Sphingolipids as a class of bioactive lipids, including ceramides, sphingomyelins, and lactosylceramides play a key modulating role in cell cycle, apoptosis, and inflammatory responses [39]. They could also participate in the osteoblast crosstalk on bone homeostasis [40]. Sphingomyelin, a critical structural component of biological membranes is tightly associated with the activation of NF-κB [39,41]. The hydrolyzed metabolites of sphingomyelin, ceramides have already been linked to the inflammation process and act as a mediator inducing the activation of NF-κB and the sensitivity of NLRP3 inflammasome [42]. Another pro-inflammatory key branching point of sphingolipid biosynthesis, lactosylceramides were also reported to involve the activation of Ras/NF-κB and inflammatory gene expression [43]. Similar to our results, the level of these sphingolipids increased significantly in RA synovial fluid [44]. Therapeutic effect of methotrexate on CIA murine sphingolipid metabolism was reported for the first time in the present study. We demonstrated that methotrexate might have downregulated the NF-κB and NLRP3/caspase-1 pathways at least partly by intervening the sphingolipid metabolism.
RA was accompanied with a differentiated energy metabolism and expenditure [45]. Trehalose was reported to inhibit the inflammation by regulating NF-κB [46]. Herein, we observed a mild alteration of energy metabolism in CIA rats after the administration of methotrexate, indicating a potent mechanism that is unestablished. Glycerophospholipids are a critical series compounds for participating in many biological regulatory processes, including the pathogenesis of RA [47]. In present study, we identified large series phosphatidylcholines, most of which had long polyunsaturated side chains. The detailed role of every phosphatidylcholines metabolites was unknown, but their potential predictive effects for diagnosis and therapy will benefit for the development and investigation of biomarkers in clinical practice.
Although the subjects with RA mainly manifests a pathological change in synovial joint and cartilage, recent studies reported that the autoimmune response in patients with RA was initiated outside the joint, affected different parts of the body, and resulted in multi-organ disorders inevitably [48]. The inflammation related pathways and cytokines that we detected in the joint ankle tissue could reflect the pathophysiological changes more sufficiently. Given that many immune components such as innate B lymphocyte and adaptive T lymphocyte have been reported to participate the pathogenesis of RA [49], whether methotrexate could potentially regulate the inflammatory mechanism of these cells need further elucidation. Untargeted metabolomics in present study was utilized to analyze all the measurable molecules in a sample, including chemical unknowns. But, the extremely limited ability to quantify with the untargeted metabolomics require us to perform the targeted lipidomic research on the RA patients treated with methotrexate. Small size for metabolomics study is the other limitation of this study. Besides, as a classical chemical therapeutic drug, the unknown side effects after inhibiting the potential target also need our further evaluation. Nevertheless, this pilot study will be a pioneer for further determination on the discovered metabolism pathways.

4. Materials and Methods

4.1. Experimental Animal and Ethical Approval

Female Wistar rats, aged six to eight weeks and weighing 160–180 g were provided by Yi-Si Animal Corporation (Changchun, China, Certificate No. SCXK (JI)-2016-0003). All of the rats were housed individually in a 12 h’ light and dark circular room. The environment of the room was under control (temperature 23 ± 2 °C and humidity 40–60%). The animals could obtain food and water ad libitum. All the animal experiments protocols were conducted based on the guidelines for the administration of laboratory animals (Directive 86/609/EEC in the Protection of Animals Used for Experimental and Other Scientific Purposes, 1986) and were approved by the Institutional Animal Care and Use Committee of Jilin University (No. SCXK-2013-0001, 27th May 2017).

4.2. RA induction and Drug Administration

Animals were randomly divided into six groups: Control Group (n = 12), Model Group (n = 12), Low Concentration Group (n = 12), Middle Concentration Group (n = 12), High Concentration Group (n = 12), and Positive Medicine Group (n = 12). After adaptive feeding for five days, fresh emulsion was prepared by emulsifying the Incomplete Freund’s Adjuvant (Chondrex, Inc.2607 151st Place NE Redmond, WA 98052, USA. Cat. 7002) with the bovine collagen solution (Chondrex, Cat. 20022) in the ice water. The induction of CIA was performed according to the manufacturer’s instructions. Briefly, the animals were induced twice at 0th and 7th day, respectively. After the booster injection, animals were administrated intragastrically with methotrexate (Tong Hua Mao Xiang Pharmaceutical Co., Ltd., Tong Hua City, Jilin Province, China, Approval Number: H22022674, 15th July 2011) twice a week (0.103 mg/mL for Low, 0.206 mg/mL for Middle, 0.412 mg/mL for High, 5 mL/kg Body weight) or diclofenac sodium (Warrant Pharmaceutical, Hunan Province, China, Approval Number: H20067776, 23rd March 2016) everyday (2.08 mg/mL for Positive, 5 mL/kg Body weight) once the swelling of paw appeared in the Model. Model and Control were administrated with distilled water.

4.3. Assessment on RA Swelling

The experimenter recorded the visual scores according to the recommended scoring system from manufacturer of Chondrex every two days. Paw swelling volume was also recorded before and after drug administration with paw volume measuring instrument (PV-200, Chengdu Tai Meng Software Co., Ltd., Chengdu City, Sichuan Province, China). The end-line of the paw volume measurement was defined as 0.3 cm above the ankle. Experimenter continued to score the swelling level until the animals were sacrificed. The observer was blind to the biological groups of experimental animals.

4.4. Sacrifice and Sample Collection

All of the animals were sacrificed when the average visual score of Model measured nearly and stably up to 4. The animals were anesthetized with sodium pentobarbital and sacrificed by collecting whole blood from abdominal aorta directly into the heparin sodium treated vacuum tube. Plasma samples was prepared by centrifuging the blood at 3000 rpm for 15 min at 4 °C and stored at −80 °C. The preparation of ankle joints tissue was performed as described previously [50]. All of the ankle joints tissues were stored at liquid nitrogen for RNA and protein extraction.

4.5. Cell Viability Assay with Real Time Cellular Analysis

Proliferation of RAW 264.7 cells, as obtained from American Type Culture Collection (ATCC TIB-71) was monitored with Real Time Cellular Analysis system (xCELLigence RTCA S16, ACEA Biosciences, San Diego, USA) in Dulbecco’s modified Eagle’s medium (DMEM, Caisson, Smithfield, UT, USA) containing 10% fetal bovine serum (FBS, Genedirex, Las Vegas, NV, USA) with 1% penicillin/streptomycin in a humidified atmosphere with 5% CO2 in a 37 °C incubator. The toxicity assays of both methotrexate and diclofenac sodium were performed with RTCA system at a series of geometric concentrations (1 mM, 100 μM, 10 μM, 1 μM, 100 nM, 10 nM, 1 nM, 0 fM). The non-toxic dosages of methotrexate and diclofenac sodium were selected to investigate their anti-inflammation effect. 3500 cells in every well were cultured in the E-plates of the RTCA system and stimulated with lipopolysaccharide (LPS, 1 μg/mL, Sigma, Darmstadt, Germany. Cat. L2630) at 12th h for 4 h. After stimulating with LPS, the medium of Model group was continued to be under the stimulation of LPS. The anti-inflammation groups were treated with methotrexate in different concentrations or diclofenac sodium respectively. Control was cultured with 10% DMEM+PBS after the four hours’ stimulation with LPS.

4.6. Cell Culture and Induction of Inflammation

RAW 264.7 cells were seeded in the six-well plates with 1 × 105 cell/well and cultured in the medium described above. The stimulation and treatment steps are consistent with the RTCA. All cells were harvested at the endpoint which was determined according to the highest peak of RTCA monitor for RNA and protein extraction. The supernatants were collected after removing cell debris by centrifugation prior to use for enzyme-linked immunosorbent assay (ELISA).

4.7. qPCR Assay

About 100–150 mg ankle tissues or 2 × 106 cells were homogenized in 1.0 mL or 200 μL TRIzol reagent (Invitrogen, Waltham, MA, USA. Cat. 15596-026). The dissected tissue was homogenized with magnetic beads homogenizer (Gering Scientific Instruments (Beijing) Ltd., Beijing, China) at 4 °C. Total RNA was extracted to be reversely transcribed and amplified according to the manufactural instructions (RevertAcid First Strand cDNA Synthesis Kit, Thermo Scientific, Waltham, MA, USA. Cat. 1662). All of the primers are listed in Table S2. Gene expression level of NLRP3, Caspase-1, IL-1β, IL-18, TNF-α, and NF-κB p65, β-actin was analyzed after 40 cycles while using Applied Biosystems StepOne Real-Time PCR System. The qPCR amplification mix was prepared as the recommended protocol of FastStart Universal SYBR Green Master (ROX, Roche, Basel, Switzerland. Cat. 04913914001,). The expressive level was calculated with the 2−ΔΔCt method.

4.8. Western Blotting and ELISA Assay

About 100–150 mg arthritic ankle joint tissue or 2 × 106 cells were homogenized as described above in 1.0 mL or 150 μL reagent of tissue protein lysate (RIPA: PMSF = 100:1) for protein extraction. The concentration of protein was determined with BCA kit (Beyotime Biotechnology Shanghai City, China, Cat. P0010). Briefly, all of the denatured protein samples were subjected by 12% SDS-PAGE gel and transferred to a polyvinylidene difluoride (PVDF) membrane (Millipore, Biosharp, Hefei City, China.) followed by blocking for 1 h and incubating with antibodies against NLRP3 (abcam, Cambridge, USA, Cat. ab214185), Caspase-1 (abcam, ab54932), NF-κB p65 (abcam, Cat. ab16502), and β-actin (Bioss, Beijing City, China. Cat.bs-0061R), respectively, overnight at 4 °C. After being incubated with secondary antibodies (Beyotime, Cat. A0208 and Cat. A0216) labeled with horseradish peroxidase (HRP), all the membranes were detected with the chemiluminescence (ECL) system and analyzed with Image-Pro Plus 6.0 software.
The plasma and cell supernatants were thawed on ice for the quantification of inflammatory cytokines. IL-1β, IL-18, TNF-α were determined by ELISA (Muti Sciences, Hangzhou City, China, Cat. EK301B5/2, Cusabio, Wuhan, China, Cat. CSB-E04610r, and Muti Sciences, Cat. EK3825/2 for murine samples and Invitrogen, Cat. BMS6002, BMS618-3 and BMS607-2INST for cellular supernant samples respectively) according to the manufactural instructions.

4.9. Metabolomic Procedure and Data Processing

All of the plasma samples, 200 μL for each sample, were added methanol (600 μL, Methanol (HPLC), Fisher Chemical, Cat. A452-1) to make free of protein. Then the supernatant (500 μL) was lyophilized at −60 °C and 10.0 pa. The residue was re-dissolved in 100 μL of methanol-water (4:1, v/v). An aliquot of 2 μL was injected for UPLC-MS/MS analysis with Waters ACQUITY UPLC system (Waters Corporation, Milford, MA, USA), which had been equipped with a BEH C18 column (2.1 mm × 100 mm, 1.7 mm, Waters Corporation). A Waters Xevo G2-S Quadrupole Time-Of-Flight (QTOF) mass spectrometer (Waters Corporation) coupled to the UPLC system was utilized to carry out the mass spectrometry with an electrospray ionization in both and positive (ESI+) and negative (ESI−) ion modes. The chromatographic conditions and spectrometric parameters were set and optimized, as described previously [51].
To ensure the stability and suitability of MS analysis, a QC sample was prepared by pooling the same volume (20 μL) from every plasma samples. Ten chromatographic peaks of ions from the QC sequencing datasheet were selected to evaluate the validation of systematic method. Repeatability and intermediate precision on spectrum and chromatogram were also estimated with six successive QC injections or 6 replicates of a plasma in both ESI+ and ESI− modes before the injection of the plasma samples. Moreover, four QC injections were performed randomly through the whole worklist.
The pre-processing of raw data was finished with MarkerLynx XS V4.1 software for alignment, deconvolution, and data reduction, so as to pair the mass and retention time with the corresponding intensities of all detected peaks. The main parameters were set similarly as before [51]. The processed files from ESI+ and ESI− modes were exported for multivariate analysis.

4.10. Molecular Docking of Methotrexate within Caspase-1

Molecular docking of methotrexate within caspase-1 was performed with Schrödinger software packages. Maestro (version 2015-2 demo, Schrödinger) was used for protein and ligand preparation, receptor grid generation and docking. The X-ray crystal structure of inhibited interleukin-1β converting enzyme (Protein Data Bank (PDB) code: 1IBC) was retrieved and downloaded from the PDB database (available online: http://www.rcsb.org/pdb). Caspase-1 protein was prepared in the Protein Preparation Wizard to optimize the structure with assigning bond orders and water orientations, removing water molecules, adding hydrogens, creating zero-order bonds to metals, and disulfide bonds. Energy of the protein was minimized with the default parameters. Receptor grid of the protein ligand was generated with the limited size of 20 Å at the active site. The crystal structure of ligand, methotrexate, was downloaded from National Center for Biotechnology Information Open Chemistry Database (NCBI, PubChem CID: 126941), and optimized in the Ligand Preparation Wizard prior to the molecular docking. Docking was performed on the Workspace using GLIDE v.6.7. At least 10 conformations were set for methotrexate. The docking result was plotted with Pymol 1.8 software, which was based on python 2.7.

4.11. Statistic and Bioinformatic Analysis

Multivariate analysis for metabolomics data pre-possessed was finished with SIMCA-P software (v14.1, Umetric, Umeå, Sweden). PCA and OPLS-DA models were established. Variable importance of project (VIP) values, which indicated a significant difference between the groups were also estimated statistically. RSD also was calculated for the pooled QC injections to assess the quality and stability of MS data. All the distinct metabolites were identified by matching accurate mass (ppm < 10) to the Human Metabolome Database (HMDB Version 4.0) [52], with confirmation determined by comparing characteristic tandem mass spectrometry (MS/MS) fragmentation patterns with METLIN [53] and HMDB Database [52] or demonstrated by referring the chemical standards (HPLC grade, Sigma). All distinct metabolites were analyzed with MetaboAnalyst 4.0 for metabolic pathways [54]. The comprehensive metabolic network was constructed with Cytoscape Software (v3.6.1) [55] based on the data from Kyoto Encyclopedia of Genes and Genomes (KEGG) database [56].
Statistical significance of the difference was estimated using the F-test. Kolmogorov-Smirnov test was used to ensure the normality of the data. Student’s t test for data with homogeneity of variance or Welch’s t-test was performed for pairwise two-group analysis. Moreover, multiple comparisons among the groups were performed by one-way analyses of variance (ANOVA). Mann-Whitney-Wilcoxon test was performed for the dataset, which does not follow the normality. All statistical significance was accepted at p < 0.05. The statistical analysis was completed with R (v3.3.3) basic statistical packages. All of the statistical bar charts were finished with R’s Basic Graphic and ggplot2 (v2.2.1) packages. The heatmap of all the differentiated metabolites were prepared with ComplexHeatmap (v1.19.1) package. The specificity and sensitivity of biomarkers was displayed by ROC estimated with pROC (v1.12.1) package.

5. Conclusions

We demonstrated that methotrexate play an inhibitory effect on rheumatoid arthritis swelling and inflammation by downregulating NF-κB and NLRP3/Caspase-1 pathways and potentially affecting the activity of caspase-1. The critical perturbed metabolic pathways, including arachidonic acid, linoleic acid, and sphingolipid interacting with NF-κB and NLRP3/Caspase-1 could contribute to the pathogenesis of RA and were also involved in the pharmacological regulation of methotrexate. Various differentiated mediators provide us with series potential clinical biomarkers for future investigation. The elucidating pharmacological mechanism of methotrexate will benefit its rational administration in clinical practice.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/19/10/2894/s1.

Author Contributions

Conceptualization, Z.P. and J.L.; Methodology, Z.P., H.L. and G.W.; Software, Z.P. and N.R.; Validation, Z.P. and Z.W., G.W. and K.F.; Formal Analysis, Z.P.; Data Curation, X.G.; Writing-Original Draft Preparation, Z.P. and Y.Y.; Writing-Review & Editing, J.L. and F.W; Visualization, H.L.; Supervision, F.W.; Project Administration, F.W. and J.L.; Funding Acquisition, F.W.

Funding

This research is supported by Graduate Innovation Fund of Jilin University [Grant No. 2017139], funded by Fundamental Research Funds for the General Universities [Grant No. 451170306058], funded by Key Laboratory of Health and Family Planning Commission of Jilin Province of China [Grant No. 3D5172303426] and Provincial School Co-construction Industrialization Demonstration Project of Jilin Province [Grant No. SXGJSF2017-1-1(01)].

Acknowledgments

The authors would like to thank Cuizhu Wang from Research Center of Natural Drug, School of Pharmaceutical Sciences, Jilin University for providing her advice and help on the metabolomics researching work.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ANOVAAnalysis of Variance
ARGArginine
ATCCAmerican Type Culture Collection
AUCArea Under Curve
CIACollagen induced arthritic
DCSDiclofenac sodium
DMEMDulbecco’s Modified Eagle’s Medium
ECLChemiluminescence
EETEpoxyoctadecenoic acid
EpOMEEpoxide of linoleic acid
ESIElectrospray Ionization
FBSFetal Bovine Serum
HETEHydroxyeicosatetraenoic acid
HMDBHuman Metabolome Database
HRPHorseradish Peroxidase
IODIntegrated optical density
KEGGKyoto Encyclopedia of Genes and Genomes
LPSLipopolysaccharide
LTC4Leukotrienes C4
MSMass Spectrometry
MS/MStandem mass spectrometry
NF-κBNuclear factor kappa-light-chain-enhancer of activated B cells
NLRP3Nucleotide binding domain and leucine-rich repeat pyrin 3 domain
NSAIDsNonsteroidal anti-inflammatory drugs
OPLS-DAOrthogonal Projections to Latent Structures Discriminant Analysis
PCPrinciple Component
PCAPrinciple Component Analysis
PcPhosphatidylcholine
PDBProtein Data Bank
PDG2Prostaglandin D2
PPMMass error ppm
PVDFPolyvinylidene Difluoride
Q/TOFQuadrupole time of flight mass spectrometry
QCQuality Control
RARheumatoid arthritis
ROCReceiver Operating Characteristic Curve
RSDRelative Standard Deviation
RTCAReal Time Cellular Analysis
sPBSStimulated with LPS followed by PBS
SMSm(d18:0/18:1)
TXA2Thromboxane A2
UPLCUltra-Performance Liquid Chromatography
VIPVariable Importance for the Projection

References

  1. Taylor, P.C.; Keystone, E.C.; van der Heijde, D.; Weinblatt, M.E.; Del Carmen Morales, L.; Reyes Gonzaga, J.; Yakushin, S.; Ishii, T.; Emoto, K.; Beattie, S.; et al. Baricitinib versus Placebo or Adalimumab in Rheumatoid Arthritis. N. Engl. J. Med. 2017, 376, 652–662. [Google Scholar] [CrossRef] [PubMed]
  2. Yousri, N.A.; Bayoumy, K.; Elhaq, W.G.; Mohney, R.P.; Emadi, S.A.; Hammoudeh, M.; Halabi, H.; Masri, B.; Badsha, H.; Uthman, I.; et al. Large Scale Metabolic Profiling identifies Novel Steroids linked to Rheumatoid Arthritis. Sci. Rep. 2017, 7, 9137. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Vande Walle, L.; Van Opdenbosch, N.; Jacques, P.; Fossoul, A.; Verheugen, E.; Vogel, P.; Beyaert, R.; Elewaut, D.; Kanneganti, T.D.; van Loo, G.; et al. Negative regulation of the NLRP3 inflammasome by A20 protects against arthritis. Nature 2014, 512, 69–73. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. McInnes, I.B.; Schett, G. The pathogenesis of rheumatoid arthritis. N. Engl. J. Med. 2011, 365, 2205–2219. [Google Scholar] [CrossRef] [PubMed]
  5. Engdahl, C.; Bondt, A.; Harre, U.; Raufer, J.; Pfeifle, R.; Camponeschi, A.; Wuhrer, M.; Seeling, M.; Martensson, I.L.; Nimmerjahn, F.; et al. Estrogen induces St6gal1 expression and increases IgG sialylation in mice and patients with rheumatoid arthritis: A potential explanation for the increased risk of rheumatoid arthritis in postmenopausal women. Arthritis Res. Ther. 2018, 20, 84. [Google Scholar] [CrossRef] [PubMed]
  6. Sohn, D.H.; Rhodes, C.; Onuma, K.; Zhao, X.; Sharpe, O.; Gazitt, T.; Shiao, R.; Fert-Bober, J.; Cheng, D.; Lahey, L.J.; et al. Local Joint inflammation and histone citrullination in a murine model of the transition from preclinical autoimmunity to inflammatory arthritis. Arthritis Rheumatol. 2015, 67, 2877–2887. [Google Scholar] [CrossRef] [PubMed]
  7. Van Loo, G.; Beyaert, R. Negative regulation of NF-kappaB and its involvement in rheumatoid arthritis. Arthritis Res. Ther. 2011, 13, 221. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Yang, Z.; Cao, J.; Yu, C.; Yang, Q.; Zhang, Y.; Han, L. Caspase-1 mediated interleukin-18 activation in neutrophils promotes the activity of rheumatoid arthritis in a NLRP3 inflammasome independent manner. Joint Bone Spine 2016, 83, 282–289. [Google Scholar] [CrossRef] [PubMed]
  9. Choulaki, C.; Papadaki, G.; Repa, A.; Kampouraki, E.; Kambas, K.; Ritis, K.; Bertsias, G.; Boumpas, D.T.; Sidiropoulos, P. Enhanced activity of NLRP3 inflammasome in peripheral blood cells of patients with active rheumatoid arthritis. Arthritis Res. Ther. 2015, 17, 257. [Google Scholar] [CrossRef] [PubMed]
  10. Wang, Z.; Liu, J.Q.; Xu, J.D.; Zhu, H.; Kong, M.; Zhang, G.H.; Duan, S.M.; Li, X.Y.; Li, G.F.; Liu, L.F.; et al. UPLC/ESI-QTOF-MS-based metabolomics survey on the toxicity of triptolide and detoxication of licorice. Chin. J. Nat. Med. 2017, 15, 474–480. [Google Scholar] [CrossRef]
  11. Li, R.; Guo, L.X.; Li, Y.; Chang, W.Q.; Liu, J.Q.; Liu, L.F.; Xin, G.Z. Dose-response characteristics of Clematis triterpenoid saponins and clematichinenoside AR in rheumatoid arthritis rats by liquid chromatography/mass spectrometry-based serum and urine metabolomics. J. Pharm. Biomed. Anal. 2017, 136, 81–91. [Google Scholar] [CrossRef] [PubMed]
  12. Li, J.; Che, N.; Xu, L.; Zhang, Q.; Wang, Q.; Tan, W.; Zhang, M. LC-MS-based serum metabolomics reveals a distinctive signature in patients with rheumatoid arthritis. Clin. Rheumatol. 2018, 37, 1493–1502. [Google Scholar] [CrossRef] [PubMed]
  13. Van Walsem, A.; Pandhi, S.; Nixon, R.M.; Guyot, P.; Karabis, A.; Moore, R.A. Relative benefit-risk comparing diclofenac to other traditional non-steroidal anti-inflammatory drugs and cyclooxygenase-2 inhibitors in patients with osteoarthritis or rheumatoid arthritis: A network meta-analysis. Arthritis Res. Ther. 2015, 17, 66. [Google Scholar] [CrossRef] [PubMed]
  14. Huffman, K.M.; Jessee, R.; Andonian, B.; Davis, B.N.; Narowski, R.; Huebner, J.L.; Kraus, V.B.; McCracken, J.; Gilmore, B.F.; Tune, K.N.; et al. Molecular alterations in skeletal muscle in rheumatoid arthritis are related to disease activity, physical inactivity, and disability. Arthritis Res. Ther. 2017, 19, 12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Huang, Z.; Tong, H.F.; Li, Y.; Qian, J.C.; Wang, J.X.; Wang, Z.; Ruan, J.C. Effect of the Polymorphism of Folylpolyglutamate Synthetase on Treatment of High-Dose Methotrexate in Pediatric Patients with Acute Lymphocytic Leukemia. Med. Sci. Monit. 2016, 22, 4967–4973. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Peres, R.S.; Liew, F.Y.; Talbot, J.; Carregaro, V.; Oliveira, R.D.; Almeida, S.L.; Franca, R.F.; Donate, P.B.; Pinto, L.G.; Ferreira, F.I.; et al. Low expression of CD39 on regulatory T cells as a biomarker for resistance to methotrexate therapy in rheumatoid arthritis. Proc. Natl. Acad. Sci. USA 2015, 112, 2509–2514. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Cheng, B.; Zheng, H.; Wu, F.; Wu, J.; Liu, X.; Tang, C.; Lu, S.; Chen, Z.; Song, F.; Ruan, J.; et al. Metabolomics analysis of Danggui Sini decoction on treatment of collagen-induced arthritis in rats. J. Chromatogr. B Analyt. Technol. Biomed. Life Sci. 2017, 1061–1062, 282–291. [Google Scholar] [CrossRef] [PubMed]
  18. Schett, G. Review: Immune cells and mediators of inflammatory arthritis. Autoimmunity 2008, 41, 224–229. [Google Scholar] [CrossRef] [PubMed]
  19. Bauernfeind, F.G.; Horvath, G.; Stutz, A.; Alnemri, E.S.; MacDonald, K.; Speert, D.; Fernandes-Alnemri, T.; Wu, J.; Monks, B.G.; Fitzgerald, K.A.; et al. Cutting edge: NF-kappaB activating pattern recognition and cytokine receptors license NLRP3 inflammasome activation by regulating NLRP3 expression. J. Immunol. 2009, 183, 787–791. [Google Scholar] [CrossRef] [PubMed]
  20. Ghosh, S.; Hayden, M.S. New regulators of NF-kappaB in inflammation. Nat. Rev. Immunol. 2008, 8, 837–848. [Google Scholar] [CrossRef] [PubMed]
  21. Mathews, R.J.; Robinson, J.I.; Battellino, M.; Wong, C.; Taylor, J.C.; Eyre, S.; Churchman, S.M.; Wilson, A.G.; Isaacs, J.D.; Hyrich, K.; et al. Evidence of NLRP3-inflammasome activation in rheumatoid arthritis (RA); genetic variants within the NLRP3-inflammasome complex in relation to susceptibility to RA and response to anti-TNF treatment. Ann. Rheum. Dis. 2014, 73, 1202–1210. [Google Scholar] [CrossRef] [PubMed]
  22. Liu, X.; Cai, F.; Zhang, Y.; Yang, A.; Liu, L. Celastrol, an NF-kappaB inhibitor, ameliorates hypercalciuria and articular cartilage lesions in a mouse model of secondary osteoporosis. J. Pharmacol. Sci. 2016, 130, 204–211. [Google Scholar] [CrossRef] [PubMed]
  23. Sun, W.X.; Liu, Y.; Zhou, W.; Li, H.W.; Yang, J.; Chen, Z.B. Shikonin inhibits TNF-alpha production through suppressing PKC-NF-kappaB-dependent decrease of IL-10 in rheumatoid arthritis-like cell model. J. Nat. Med. 2017, 71, 349–356. [Google Scholar] [CrossRef] [PubMed]
  24. Tao, H.X.; Xiong, W.; Zhao, G.D.; Peng, Y.; Zhong, Z.F.; Xu, L.; Duan, R.; Tsim, K.W.K.; Yu, H.; Wang, Y.T. Discrimination of three Siegesbeckiae Herba species using UPLC-QTOF/MS-based metabolomics approach. Food Chem. Toxicol. 2018, 119, 400–406. [Google Scholar] [CrossRef] [PubMed]
  25. Kim, Y.; Gromovsky, A.D.; Brown, J.M.; Chung, S. Gamma-tocotrienol attenuates the aberrant lipid mediator production in NLRP3 inflammasome-stimulated macrophages. J. Nutr. Biochem. 2018, 58, 169–177. [Google Scholar] [CrossRef] [PubMed]
  26. Camandola, S.; Leonarduzzi, G.; Musso, T.; Varesio, L.; Carini, R.; Scavazza, A.; Chiarpotto, E.; Baeuerle, P.A.; Poli, G. Nuclear factor kB is activated by arachidonic acid but not by eicosapentaenoic acid. Biochem. Biophys. Res. Commun. 1996, 229, 643–647. [Google Scholar] [CrossRef] [PubMed]
  27. Wen, H.; Gris, D.; Lei, Y.; Jha, S.; Zhang, L.; Huang, M.T.; Brickey, W.J.; Ting, J.P. Fatty acid-induced NLRP3-ASC inflammasome activation interferes with insulin signaling. Nat. Immunol. 2011, 12, 408–415. [Google Scholar] [CrossRef] [PubMed]
  28. Morisseau, C.; Hammock, B.D. Impact of soluble epoxide hydrolase and epoxyeicosanoids on human health. Annu. Rev. Pharmacol. Toxicol. 2013, 53, 37–58. [Google Scholar] [CrossRef] [PubMed]
  29. Node, K.; Huo, Y.; Ruan, X.; Yang, B.; Spiecker, M.; Ley, K.; Zeldin, D.C.; Liao, J.K. Anti-inflammatory properties of cytochrome P450 epoxygenase-derived eicosanoids. Science 1999, 285, 1276–1279. [Google Scholar] [CrossRef] [PubMed]
  30. Tam, V.C. Lipidomic profiling of bioactive lipids by mass spectrometry during microbial infections. Seminars Immunol. 2013, 25, 240–248. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Ishizaki, T.; Shigemori, K.; Nakai, T.; Miyabo, S.; Ozawa, T.; Chang, S.W.; Voelkel, N.F. Leukotoxin, 9,10-epoxy-12-octadecenoate causes edematous lung injury via activation of vascular nitric oxide synthase. Am. J. Physiol. 1995, 269, L65–L70. [Google Scholar] [CrossRef] [PubMed]
  32. De Lima-Salgado, T.M.; Alba-Loureiro, T.C.; do Nascimento, C.S.; Nunes, M.T.; Curi, R. Molecular mechanisms by which saturated fatty acids modulate TNF-alpha expression in mouse macrophage lineage. Cell Biochem. Biophys. 2011, 59, 89–97. [Google Scholar] [CrossRef] [PubMed]
  33. Reynolds, C.M.; Roche, H.M. Conjugated linoleic acid and inflammatory cell signalling. Prostaglandins Leukot. Essent. Fatty Acids 2010, 82, 199–204. [Google Scholar] [CrossRef] [PubMed]
  34. Luczaj, W.; Gindzienska-Sieskiewicz, E.; Jarocka-Karpowicz, I.; Andrisic, L.; Sierakowski, S.; Zarkovic, N.; Waeg, G.; Skrzydlewska, E. The onset of lipid peroxidation in rheumatoid arthritis: Consequences and monitoring. Free Radic. Res. 2016, 50, 304–313. [Google Scholar] [CrossRef] [PubMed]
  35. Reynolds, C.M.; Segovia, S.A.; Zhang, X.D.; Gray, C.; Vickers, M.H. Maternal high-fat diet-induced programing of gut taste receptor and inflammatory gene expression in rat offspring is ameliorated by CLA supplementation. Physiol. Rep. 2015, 3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Aggarwal, N.; Korenbaum, E.; Mahadeva, R.; Immenschuh, S.; Grau, V.; Dinarello, C.A.; Welte, T.; Janciauskiene, S. α-Linoleic acid enhances the capacity of alpha-1 antitrypsin to inhibit lipopolysaccharide-induced IL-1beta in human blood neutrophils. Mol. Med. 2016, 22. [Google Scholar] [CrossRef] [PubMed]
  37. Fijlstra, M.; Tissing, W.J.; Stellaard, F.; Verkade, H.J.; Rings, E.H. Reduced absorption of long-chain fatty acids during methotrexate-induced gastrointestinal mucositis in the rat. Clin. Nutr. 2013, 32, 452–459. [Google Scholar] [CrossRef] [PubMed]
  38. Dai, J.; Shen, J.; Pan, W.; Shen, S.; Das, U.N. Effects of polyunsaturated fatty acids on the growth of gastric cancer cells in vitro. Lipids Health Dis. 2013, 12, 71. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Adada, M.; Luberto, C.; Canals, D. Inhibitors of the sphingomyelin cycle: Sphingomyelin synthases and sphingomyelinases. Chem. Phys. Lipids 2016, 197, 45–59. [Google Scholar] [CrossRef] [PubMed]
  40. Meshcheryakova, A.; Mechtcheriakova, D.; Pietschmann, P. Sphingosine 1-phosphate signaling in bone remodeling: Multifaceted roles and therapeutic potential. Expert Opin. Ther. Targets 2017, 21, 725–737. [Google Scholar] [CrossRef] [PubMed]
  41. Luberto, C.; Yoo, D.S.; Suidan, H.S.; Bartoli, G.M.; Hannun, Y.A. Differential effects of sphingomyelin hydrolysis and resynthesis on the activation of NF-kappa B in normal and SV40-transformed human fibroblasts. J. Biol. Chem. 2000, 275, 14760–14766. [Google Scholar] [CrossRef] [PubMed]
  42. Maceyka, M.; Spiegel, S. Sphingolipid metabolites in inflammatory disease. Nature 2014, 510, 58–67. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Won, J.S.; Singh, A.K.; Singh, I. Lactosylceramide: A lipid second messenger in neuroinflammatory disease. J. Neurochem. 2007, 103 (Suppl. 1), 180–191. [Google Scholar] [CrossRef] [PubMed]
  44. Kosinska, M.K.; Liebisch, G.; Lochnit, G.; Wilhelm, J.; Klein, H.; Kaesser, U.; Lasczkowski, G.; Rickert, M.; Schmitz, G.; Steinmeyer, J. Sphingolipids in human synovial fluid--a lipidomic study. PLoS ONE 2014, 9, e91769. [Google Scholar] [CrossRef] [PubMed]
  45. Choe, J.Y.; Park, S.H.; Lee, H.; Kwon, H.H.; Kim, S.K. Resting energy expenditure is not associated with disease activity in women with rheumatoid arthritis: Cross-sectional study. Korean J. Intern. Med. 2014, 29, 516–524. [Google Scholar] [CrossRef] [PubMed]
  46. Echigo, R.; Shimohata, N.; Karatsu, K.; Yano, F.; Kayasuga-Kariya, Y.; Fujisawa, A.; Ohto, T.; Kita, Y.; Nakamura, M.; Suzuki, S.; et al. Trehalose treatment suppresses inflammation, oxidative stress, and vasospasm induced by experimental subarachnoid hemorrhage. J. Transl. Med. 2012, 10, 80. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Fuchs, B.; Schiller, J.; Wagner, U.; Hantzschel, H.; Arnold, K. The phosphatidylcholine/ lysophosphatidylcholine ratio in human plasma is an indicator of the severity of rheumatoid arthritis: Investigations by 31P NMR and MALDI-TOF MS. Clin. Biochem. 2005, 38, 925–933. [Google Scholar] [CrossRef] [PubMed]
  48. Yun, X.; Dong, S.; Hu, Q.; Dai, Y.; Xia, Y. 1H NMR-based metabolomics approach to investigate the urine samples of collagen-induced arthritis rats and the intervention of tetrandrine. J. Pharm. Biomed. Anal. 2018, 154, 302–311. [Google Scholar] [CrossRef] [PubMed]
  49. Semerano, L.; Minichiello, E.; Bessis, N.; Boissier, M.C. Novel Immunotherapeutic Avenues for Rheumatoid Arthritis. Trends Mol. Med. 2016, 22, 214–229. [Google Scholar] [CrossRef] [PubMed]
  50. Tanwar, A.; Chawla, R.; Ansari, M.M.; Thakur, P.; Chakotiya, A.S.; Goel, R.; Ojha, H.; Asif, M.; Basu, M.; Arora, R.; et al. In vivo anti-arthritic efficacy of Camellia sinensis (L.) in collagen induced arthritis model. Biomed. Pharmacother. 2017, 87, 92–101. [Google Scholar] [CrossRef] [PubMed]
  51. Wang, C.; Zhang, N.; Wang, Z.; Qi, Z.; Zhu, H.; Zheng, B.; Li, P.; Liu, J. Nontargeted Metabolomic Analysis of Four Different Parts of Platycodon grandiflorum Grown in Northeast China. Molecules 2017, 22, 1280. [Google Scholar] [CrossRef] [PubMed]
  52. Wishart, D.S.; Feunang, Y.D.; Marcu, A.; Guo, A.C.; Liang, K.; Vazquez-Fresno, R.; Sajed, T.; Johnson, D.; Li, C.; Karu, N.; et al. HMDB 4.0: The human metabolome database for 2018. Nucleic Acids Res. 2018, 46, D608–D617. [Google Scholar] [CrossRef] [PubMed]
  53. Guijas, C.; Montenegro-Burke, J.R.; Domingo-Almenara, X.; Palermo, A.; Warth, B.; Hermann, G.; Koellensperger, G.; Huan, T.; Uritboonthai, W.; Aisporna, A.E.; et al. METLIN: A Technology Platform for Identifying Knowns and Unknowns. Anal. Chem. 2018, 90, 3156–3164. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Chong, J.; Soufan, O.; Li, C.; Caraus, I.; Li, S.; Bourque, G.; Wishart, D.S.; Xia, J. MetaboAnalyst 4.0: Towards more transparent and integrative metabolomics analysis. Nucleic Acids Res. 2018, 46, W486–W494. [Google Scholar] [CrossRef] [PubMed]
  55. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef] [PubMed]
  56. Kanehisa, M.; Furumichi, M.; Tanabe, M.; Sato, Y.; Morishima, K. KEGG: New perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017, 45, D353–D361. [Google Scholar] [CrossRef] [PubMed]
Ijms 19 02894 g001 550
Figure 1. Evaluation on paw swelling of all groups. Swelling volume on Day 11 (A) and Day 25 (B). (C) The visual scores line of all groups including the whole time points. (D) The representative swelling states of all groups on Day 25. # p < 0.05 vs. Control. ## p < 0.01 vs. Control. * p < 0.05 vs. Model. ** p < 0.01 vs. Model. The free small circles in (B) represent the mild outlier samples.
Figure 1. Evaluation on paw swelling of all groups. Swelling volume on Day 11 (A) and Day 25 (B). (C) The visual scores line of all groups including the whole time points. (D) The representative swelling states of all groups on Day 25. # p < 0.05 vs. Control. ## p < 0.01 vs. Control. * p < 0.05 vs. Model. ** p < 0.01 vs. Model. The free small circles in (B) represent the mild outlier samples.
Ijms 19 02894 g001
Ijms 19 02894 g002 550
Figure 2. Quantification results of real-time PCR in all murine groups. Gene expression level of (A) NLRP3, (B) caspase-1, (C) IL-1β, (D) IL-18, (E) NF-κB p65, (F) TNF-α. All values were normalized against β-actin. DCS, Diclofenac sodium. ## p < 0.01 vs. Control. * p < 0.05 vs. Model. ** p < 0.01 vs. Model.
Figure 2. Quantification results of real-time PCR in all murine groups. Gene expression level of (A) NLRP3, (B) caspase-1, (C) IL-1β, (D) IL-18, (E) NF-κB p65, (F) TNF-α. All values were normalized against β-actin. DCS, Diclofenac sodium. ## p < 0.01 vs. Control. * p < 0.05 vs. Model. ** p < 0.01 vs. Model.
Ijms 19 02894 g002
Ijms 19 02894 g003 550
Figure 3. Protein levels of NF-κB and NLRP3/Caspase-1 pathway in CIA rats. Results of optical density of NLRP3 (A), caspase-1 (B), NF-κB p65 (C) estimated by integrated optical density (IOD index) and normalized against β-actin in all groups. (D)The representative Chemiluminescencez panels of all protein membranes. DCS, Diclofenac sodium. ## p < 0.01 vs. Control. * p < 0.05 vs. Model. ** p < 0.01 vs. Model.
Figure 3. Protein levels of NF-κB and NLRP3/Caspase-1 pathway in CIA rats. Results of optical density of NLRP3 (A), caspase-1 (B), NF-κB p65 (C) estimated by integrated optical density (IOD index) and normalized against β-actin in all groups. (D)The representative Chemiluminescencez panels of all protein membranes. DCS, Diclofenac sodium. ## p < 0.01 vs. Control. * p < 0.05 vs. Model. ** p < 0.01 vs. Model.
Ijms 19 02894 g003
Ijms 19 02894 g004 550
Figure 4. Quantification of plasma cytokines related with NF-κB and NLRP3/Caspase-1 pathways in all CIA murine groups. (A) IL-1β, (B) IL-18, (C) TNF-α levels in all groups. DCS, Diclofenac sodium. ## p < 0.01 vs. Control. * p < 0.05 vs. Model. ** p < 0.01 vs. Model.
Figure 4. Quantification of plasma cytokines related with NF-κB and NLRP3/Caspase-1 pathways in all CIA murine groups. (A) IL-1β, (B) IL-18, (C) TNF-α levels in all groups. DCS, Diclofenac sodium. ## p < 0.01 vs. Control. * p < 0.05 vs. Model. ** p < 0.01 vs. Model.
Ijms 19 02894 g004
Ijms 19 02894 g005 550
Figure 5. Proliferation curves of RAW 264.7 monitored with RTCA system and the statistic results. Cytotoxic assay of methotrexate (A) and Diclofenac sodium (DCS) (B) with a series of geometric concentrations from 1 nM to 1 mM. The treatment of methotrexate was performed at 24th h. (A1) The proliferation curve of methotrexate with different dosages. (A2) Statistics of cell indexes in A1 at 72nd h The treatment of DCS was performed at 48th h. (B1) The proliferation curve of DCS with different dosages. (B2) Statistics of cell indexes in B1 at 96th h. (C1) Anti-inflammation curves of methotrexate with different dosages and DCS (1 nM). All cellular groups were treated with LPS at 12th h for 4 h. Cellular Control (Stimulated with LPS followed by PBS, sPBS) was free of LPS after the 4 h’ stimulation. The cell indexes were statistically evaluated at 72nd h. ** p < 0.01 vs. Model. ▲ p > 0.05 vs. 0 fM in A2 and B2, vs. Model in C2.
Figure 5. Proliferation curves of RAW 264.7 monitored with RTCA system and the statistic results. Cytotoxic assay of methotrexate (A) and Diclofenac sodium (DCS) (B) with a series of geometric concentrations from 1 nM to 1 mM. The treatment of methotrexate was performed at 24th h. (A1) The proliferation curve of methotrexate with different dosages. (A2) Statistics of cell indexes in A1 at 72nd h The treatment of DCS was performed at 48th h. (B1) The proliferation curve of DCS with different dosages. (B2) Statistics of cell indexes in B1 at 96th h. (C1) Anti-inflammation curves of methotrexate with different dosages and DCS (1 nM). All cellular groups were treated with LPS at 12th h for 4 h. Cellular Control (Stimulated with LPS followed by PBS, sPBS) was free of LPS after the 4 h’ stimulation. The cell indexes were statistically evaluated at 72nd h. ** p < 0.01 vs. Model. ▲ p > 0.05 vs. 0 fM in A2 and B2, vs. Model in C2.
Ijms 19 02894 g005
Ijms 19 02894 g006 550
Figure 6. Quantification results of real-time PCR in all cellular groups. Gene expression level of (A) NLRP3, (B) caspase-1, (C) IL-1β, (D) IL-18, (E) NF-κB p65, (F) TNF-α. All values were normalized against β-actin. DCS, Diclofenac sodium. ## p < 0.01 vs. Control. * p < 0.05 vs. Model. ** p < 0.01 vs. Model.
Figure 6. Quantification results of real-time PCR in all cellular groups. Gene expression level of (A) NLRP3, (B) caspase-1, (C) IL-1β, (D) IL-18, (E) NF-κB p65, (F) TNF-α. All values were normalized against β-actin. DCS, Diclofenac sodium. ## p < 0.01 vs. Control. * p < 0.05 vs. Model. ** p < 0.01 vs. Model.
Ijms 19 02894 g006
Ijms 19 02894 g007 550
Figure 7. Protein levels of NF-κB and NLRP3/Caspase-1 pathway in RAW 264.7 cells. Results of optical density of NLRP3 (A), caspase-1 (B), NF-κB p65 (C) estimated by integrated optical density (IOD index) and normalized against β-actin. (D)The representative Chemiluminescencez panels of all protein membranes. DCS, Diclofenac sodium. ## p < 0.01 vs. Control. ** p < 0.01 vs. Model.
Figure 7. Protein levels of NF-κB and NLRP3/Caspase-1 pathway in RAW 264.7 cells. Results of optical density of NLRP3 (A), caspase-1 (B), NF-κB p65 (C) estimated by integrated optical density (IOD index) and normalized against β-actin. (D)The representative Chemiluminescencez panels of all protein membranes. DCS, Diclofenac sodium. ## p < 0.01 vs. Control. ** p < 0.01 vs. Model.
Ijms 19 02894 g007
Ijms 19 02894 g008 550
Figure 8. Quantification of supernatants cytokines related with NF-κB and NLRP3/Caspase-1 pathways in all cellular groups at 48th h (A) IL-1β, (B) IL-18, (C) TNF-α levels in all groups. DCS, Diclofenac sodium. ## p < 0.01 vs. Control. ** p < 0.01 vs. Model.
Figure 8. Quantification of supernatants cytokines related with NF-κB and NLRP3/Caspase-1 pathways in all cellular groups at 48th h (A) IL-1β, (B) IL-18, (C) TNF-α levels in all groups. DCS, Diclofenac sodium. ## p < 0.01 vs. Control. ** p < 0.01 vs. Model.
Ijms 19 02894 g008
Ijms 19 02894 g009 550
Figure 9. Multivariate analysis results of Ultra-Performance Liquid Chromatography-Mass Spectrometry (UPLC-MS). Principle Component Analysis (PCA) and orthogonal projections to latent structures discriminant analysis (OPLS-DA) models were displayed as two-dimensional plots. PCA of all samples in electrospray ionization in positive (ESI+) mode (A1) and in electrospray ionization in negative (ESI−) mode (B1) (PC1 = 51.1%, PC2 = 8.25% for ESI+ and PC1 = 28%, PC2 = 11.6% for ESI−). The OPLS-DA discrimination score plots between different groups: (A2) Control vs. Model in ESI+, p < 0.00001; (A3) Model vs. Treatment groups in ESI+, p = 0.015; (B2) Control vs. Model in ESI−, p < 0.00001; (B3) Model vs. Treatment groups in ESI−, p < 0.0001. Control was abbreviated as C, Model as M, Treatment with low, middle, high dosage as T_L, T_M and T_H respectively. Treatment with positive DCS as T_P, QC as Q. In PCA, t [1] are the first variables summarized from the original data and t [2] are the second ones. In OPLS-DA, t [1] is the predictive X-scores and to [1] is the Orthogonal X-scores. All p values were estimated with CV-ANOVA test. * in A3 means the multiple magnification to the coordinate axes.
Figure 9. Multivariate analysis results of Ultra-Performance Liquid Chromatography-Mass Spectrometry (UPLC-MS). Principle Component Analysis (PCA) and orthogonal projections to latent structures discriminant analysis (OPLS-DA) models were displayed as two-dimensional plots. PCA of all samples in electrospray ionization in positive (ESI+) mode (A1) and in electrospray ionization in negative (ESI−) mode (B1) (PC1 = 51.1%, PC2 = 8.25% for ESI+ and PC1 = 28%, PC2 = 11.6% for ESI−). The OPLS-DA discrimination score plots between different groups: (A2) Control vs. Model in ESI+, p < 0.00001; (A3) Model vs. Treatment groups in ESI+, p = 0.015; (B2) Control vs. Model in ESI−, p < 0.00001; (B3) Model vs. Treatment groups in ESI−, p < 0.0001. Control was abbreviated as C, Model as M, Treatment with low, middle, high dosage as T_L, T_M and T_H respectively. Treatment with positive DCS as T_P, QC as Q. In PCA, t [1] are the first variables summarized from the original data and t [2] are the second ones. In OPLS-DA, t [1] is the predictive X-scores and to [1] is the Orthogonal X-scores. All p values were estimated with CV-ANOVA test. * in A3 means the multiple magnification to the coordinate axes.
Ijms 19 02894 g009
Ijms 19 02894 g010 550
Figure 10. S-plots of all differentiated metabolites based on the OPLS-DA models. C vs. M in ESI+ (A) and in ESI− (B); M vs. T in ESI+ (C) and in ESI− (D). All of the significantly changed metabolites identified were marked as red four-polar stars. Other components were green points. The p [1]-axis describes the magnitude of each variable in X. The p(corr) [1]-axis represents the reliability of each variable in X. Control was abbreviated as C, Model as M, Treatment with drugs as T.
Figure 10. S-plots of all differentiated metabolites based on the OPLS-DA models. C vs. M in ESI+ (A) and in ESI− (B); M vs. T in ESI+ (C) and in ESI− (D). All of the significantly changed metabolites identified were marked as red four-polar stars. Other components were green points. The p [1]-axis describes the magnitude of each variable in X. The p(corr) [1]-axis represents the reliability of each variable in X. Control was abbreviated as C, Model as M, Treatment with drugs as T.
Ijms 19 02894 g010
Ijms 19 02894 g011 550
Figure 11. Differentiated metabolites and perturbed metabolic pathways. (A) The heatmap all of differentiated metabolites. (B) Results of metabolic pathways produced by MetaboAnalyst 4.0. The pathway Impact value described the significance of the metabolic pathways. (C) The comprehensive analysis about the invention from methotrexate and interaction between inflammatory pathways (including NF-κB and NLRP3/Caspase-1) and metabolic pathways (arachidonic acid, linoleic acid, sphingolipid, glycerophospholipid, Galactose, and starch and sucrose metabolic pathways). EET, epoxyoctadecenoic acid; HETE, hydroxyeicosatetraenoic acid; LTC4, Leukotriene C4; SM, Sm(d18:0/18:1); TXA2, Thromboxane A2. PGD2, Prostaglandin D2.
Figure 11. Differentiated metabolites and perturbed metabolic pathways. (A) The heatmap all of differentiated metabolites. (B) Results of metabolic pathways produced by MetaboAnalyst 4.0. The pathway Impact value described the significance of the metabolic pathways. (C) The comprehensive analysis about the invention from methotrexate and interaction between inflammatory pathways (including NF-κB and NLRP3/Caspase-1) and metabolic pathways (arachidonic acid, linoleic acid, sphingolipid, glycerophospholipid, Galactose, and starch and sucrose metabolic pathways). EET, epoxyoctadecenoic acid; HETE, hydroxyeicosatetraenoic acid; LTC4, Leukotriene C4; SM, Sm(d18:0/18:1); TXA2, Thromboxane A2. PGD2, Prostaglandin D2.
Ijms 19 02894 g011
Ijms 19 02894 g012 550
Figure 12. Docking interaction modes of methotrexate within the binding pockets of caspase-1. (A) The 3-dimensional polar contacts between methotrexate and caspase-1. The light green cartoon is caspase-1 prepared by pymol. The red dotted line is the polar contacts. The thin orange sticks are the interactive residues of caspase-1. The thick sticks are the methotrexate (Yellow are carbon atoms, Blue are nitrogen atoms, Red are oxygen atoms, White are some of the hydrogen atoms). (B) two-dimensional display of the interaction between the ligand and the protein (The colorful balls at the edge are related amino acids. Polar contacts are displayed with different lines. Purple dotted lines represent the polar contacts with Oxygen anion, green line represents with Benzene, blue red line represents with ammonium).
Figure 12. Docking interaction modes of methotrexate within the binding pockets of caspase-1. (A) The 3-dimensional polar contacts between methotrexate and caspase-1. The light green cartoon is caspase-1 prepared by pymol. The red dotted line is the polar contacts. The thin orange sticks are the interactive residues of caspase-1. The thick sticks are the methotrexate (Yellow are carbon atoms, Blue are nitrogen atoms, Red are oxygen atoms, White are some of the hydrogen atoms). (B) two-dimensional display of the interaction between the ligand and the protein (The colorful balls at the edge are related amino acids. Polar contacts are displayed with different lines. Purple dotted lines represent the polar contacts with Oxygen anion, green line represents with Benzene, blue red line represents with ammonium).
Ijms 19 02894 g012
Table
Table 1. Detailed results of all differentiated compounds identified in plasma.
Table 1. Detailed results of all differentiated compounds identified in plasma.
NOtR/minMass/DaESI §VIPIdentified CompoundKEGG IDPPMMSE Fragmentation †Pathways
S1 *17.26295.2306−1.479,10- & 12,13-EpOMEC14825 ▲9134.8979, 152.9976, 167.1678, 183.0139, 195.8134, 197.2099, 223.1726, 277.2354, 295.2302Linoleic acid
S218.15303.2299+2.47Linoleic acidC015951105.1102, 146.9803, 149.0211, 165.0605, 185.0739, 221.9592, 281.1889Linoleic acid
S322.98303.2358−1.43Arachidonic acidC002199100.9381, 160.8434, 219.8489, 243.9024, 255.2361, 259.2451, 303.2365Arachidonic acid
S4 *18.26319.2309−1.02EETs # & HETEs ##C14768 ▲▲8114.9385, 127.4831, 136.9957, 163.1423, 174.9591, 219.1486, 235.1270, 255.6255, 257.6209, 319.2300, etc.Arachidonic acid
S529.05341.1083−2.26Alpha-lactoseC002432116.9300, 145.1682, 160.8460, 265.1563, 277.4372, 281.2538, 311.1749Starch and Galactose
S629.24341.1122−2.01SucroseC000892101.0466, 146.9683, 160.9446, 249.5477, 341.1111Starch and Galactose
S729.23387.1167−1.25TrehaloseC010836101.0466, 146.9683, 160.9446, 179.0878, 341.1111Starch and Galactose
S8 *19.10397.2239−1.15ProstanoidsC02198 ▲▲▲2115.1325, 134.9978, 166.8685, 183.0141, 205.5410, 207.1378, 279.2340, 303.2343, 351.1919Arachidonic acid
S918.69624.3016−1.12Leukotriene C4C021669146.1692, 167.0458, 225.0720, 508.3372, 562.3419Arachidonic acid
S10 *2.55648.6284−1.13Ceramide (d18:1/24:0)C001952248.9640, 316.9509, 452.9229, 649.8531Sphingolipid
S1124.76730.5941+1.97Sm(d18:0/18:1)C005506149.0213, 183.0713, 282.2614, 310.2769, 561.3276Sphingolipid
S1225.75850.5526−1.82Lactosylceramide (d18:1/12:0)C012901146.9697, 179.0561, 198.2738, 480.3100, 626.3899Sphingolipid
S1328.68734.5693+2.97Pc (16:0/16:0)C001570184.0714, 496.2644, 551.5518, 734.5715Linoleic acid & Glycerophospholipid
S1426.62742.5342-3.87Pc (15:0/18:2)C001577141.0705, 279.2360, 281.2520, 655.4543, 742.5328
S1527.63752.5248+2.78Pc (16:1/16:1)C001576184.0159, 237.0634, 279.2367, 281.0246, 464.3080, 474.2838
S1624.79756.5538+1.85Pc (14:0/20:3)C001570183.0151, 247.2348, 271.2355, 468.3029, 628.3584
S1723.84758.5713+2.89Pc (18:1/16:1)C001572184.0157, 255.2365, 590.5201, 655.4091, 740.5466
S1826.38782.5714+3.71Pc (20:4/16:0)C001573184.0706, 257.2607, 478.3263, 679.5030
S1927.41786.6048+18.77Pc (18:1/18:1)C001575184.0718, 283.2617, 339.3408, 604.5342, 786.6055
S2027.08798.5323−2.59Pc (20:4/14:0)C001574227.2040, 257.2359, 303.2356, 439.2255, 669.4767
S2126.90804.5745−1.42Pc (14:1/20:0)C001572181.1628, 265.2368, 307.2393, 758.5609
S2228.72810.6019+23.25Pc (20:3/18:1)C001571184.0713, 221.0816, 283.2609, 307.2764, 627.3745
S2325.79832.5876+6.78Pc (18:2/20:2)C001576184.0728, 309.2785, 627.3777, 684.4404, 752.5748
S2426.87854.5839−4.99Pc (20:4/18:0)C001579255.2352, 279.2362, 283.2639, 303.2360, 725.5103
§ ESI modes, Positive mode (+) and Negative mode (−) of electrospray ionization; † Spectrum ion mass fragments in MS/MS under MSE mode. The characteristic fragments are highlighted as bold font; * Metabolites identified against to the standard chemicals; # EET, epoxyoctadecenoic acid; ## HETE, hydroxyeicosatetraenoic acid; ▲, include not only C14825, but also C14826; ▲▲, include not only C14768, but also C14769, C14771, C04742 and C04805; ▲▲▲, include not only C02198, but also C00696.
Table
Table 2. Molecular docking values of methotrexate within caspase-1.
Table 2. Molecular docking values of methotrexate within caspase-1.
ItemsPotential Energy-OPLS3Docking ScoreXP GScoreH-Bond
Values209.403−4.683−5.9536

Share and Cite

MDPI and ACS Style

Pang, Z.; Wang, G.; Ran, N.; Lin, H.; Wang, Z.; Guan, X.; Yuan, Y.; Fang, K.; Liu, J.; Wang, F. Inhibitory Effect of Methotrexate on Rheumatoid Arthritis Inflammation and Comprehensive Metabolomics Analysis Using Ultra-Performance Liquid Chromatography-Quadrupole Time of Flight-Mass Spectrometry (UPLC-Q/TOF-MS). Int. J. Mol. Sci. 2018, 19, 2894. https://doi.org/10.3390/ijms19102894

AMA Style

Pang Z, Wang G, Ran N, Lin H, Wang Z, Guan X, Yuan Y, Fang K, Liu J, Wang F. Inhibitory Effect of Methotrexate on Rheumatoid Arthritis Inflammation and Comprehensive Metabolomics Analysis Using Ultra-Performance Liquid Chromatography-Quadrupole Time of Flight-Mass Spectrometry (UPLC-Q/TOF-MS). International Journal of Molecular Sciences. 2018; 19(10):2894. https://doi.org/10.3390/ijms19102894

Chicago/Turabian Style

Pang, Zhiqiang, Guoqiang Wang, Nan Ran, Hongqiang Lin, Ziyan Wang, Xuewa Guan, Yuze Yuan, Keyong Fang, Jinping Liu, and Fang Wang. 2018. "Inhibitory Effect of Methotrexate on Rheumatoid Arthritis Inflammation and Comprehensive Metabolomics Analysis Using Ultra-Performance Liquid Chromatography-Quadrupole Time of Flight-Mass Spectrometry (UPLC-Q/TOF-MS)" International Journal of Molecular Sciences 19, no. 10: 2894. https://doi.org/10.3390/ijms19102894

APA Style

Pang, Z., Wang, G., Ran, N., Lin, H., Wang, Z., Guan, X., Yuan, Y., Fang, K., Liu, J., & Wang, F. (2018). Inhibitory Effect of Methotrexate on Rheumatoid Arthritis Inflammation and Comprehensive Metabolomics Analysis Using Ultra-Performance Liquid Chromatography-Quadrupole Time of Flight-Mass Spectrometry (UPLC-Q/TOF-MS). International Journal of Molecular Sciences, 19(10), 2894. https://doi.org/10.3390/ijms19102894

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop