Polyunsaturated fatty acids (PUFAs) are essential FAs for human health. Cytochrome P450 oxygenates PUFAs to produce anti-inflammatory and pain-resolving epoxy fatty acids (EpFAs) and other oxylipins whose epoxide ring is opened by the soluble epoxide hydrolase (sEH/Ephx2), resulting in the formation of toxic and pro-inflammatory vicinal diols (dihydroxy-FAs). Pharmacological inhibition of sEH is a promising strategy for the treatment of pain, inflammation, cardiovascular diseases, and other conditions. We tested the efficacy of a potent, selective sEH inhibitor, 1-trifluoromethoxyphenyl-3-(1-propionylpiperidin-4-yl) urea (TPPU), in an animal model of multiple sclerosis (MS), experimental autoimmune encephalomyelitis (EAE). Prophylactic TPPU treatment significantly ameliorated EAE without affecting circulating white blood cell counts. TPPU accumulated in the spinal cords (SCs), which was correlated with plasma TPPU concentration. Targeted lipidomics in EAE SCs and plasma identified that TPPU blocked production of dihydroxy-FAs efficiently and increased some EpFA species including 12(13)-epoxy-octadecenoic acid (12(13)-EpOME) and 17(18)-epoxy-eicosatrienoic acid (17(18)-EpETE). TPPU did not alter levels of cyclooxygenase (COX-1/2) metabolites, while it increased 12-hydroxyeicosatetraenoic acid (12-HETE) and other 12/15-lipoxygenase metabolites. These analytical results are consistent with sEH inhibitors that reduce neuroinflammation and accelerate anti-inflammatory responses, providing the possibility that sEH inhibitors could be used as a disease modifying therapy, as well as for MS-associated pain relief.

Arachidonic acid (AA), a representative ω6-polyunsaturated fatty acid (PUFA), is a precursor of 2-series prostaglandins (PGs) that play important roles in inflammation, pain, fever, and related disorders including cardiovascular diseases. Eating fish or supplementation with the ω3-PUFAs such as eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) is widely assumed to be beneficial in preventing cardiovascular diseases. A proposed mechanism for a cardio-protective role of ω3-PUFAs assumes competition between AA and ω3-PUFAs for cyclooxygenases (COX), leading to reduced production of 2-series PGs. In this study, we have used a systems biology approach to integrate existing knowledge and novel high-throughput data that facilitates a quantitative understanding of the molecular mechanism of ω3- and ω6-PUFA metabolism in mammalian cells. We have developed a quantitative computational model of the competitive metabolism of AA and EPA via the COX pathway through a two-step matrix-based approach to estimate the rate constants. This model was developed by using lipidomic data sets that were experimentally obtained from EPA-supplemented ATP-stimulated RAW264.7 macrophages. The resulting model fits the experimental data well for all metabolites and demonstrates that the integrated metabolic and signaling networks and the experimental data are consistent with one another. The robustness of the model was validated through parametric sensitivity and uncertainty analysis. We also validated the model by predicting the results from other independent experiments involving AA- and DHA-supplemented ATP-stimulated RAW264.7 cells using the parameters estimated with EPA. Furthermore, we showed that the higher affinity of EPA binding to COX compared with AA was able to inhibit AA metabolism effectively. Thus, our model captures the essential features of competitive metabolism of ω3- and ω6-PUFAs.

Eicosanoids, including prostaglandins (PG) and leukotrienes, are lipid mediators derived from arachidonic acid. A quantitative and biochemical level understanding of eicosanoid metabolism would aid in understanding the mechanisms that govern inflammatory processes. Here, we present a combined experimental and computational approach to understanding the biochemical basis of eicosanoid metabolism in macrophages. Lipidomic and transcriptomic measurements and analyses reveal temporal and dynamic changes of the eicosanoid metabolic network in mouse bone marrow-derived macrophages (BMDM) upon stimulation of the Toll-like receptor 4 with Kdo2-Lipid A (KLA) and stimulation of the P2X7 purinergic receptor with adenosine 5'-triphosphate. Kinetic models were developed for the cyclooxygenase (COX) and lipoxygenase branches of arachidonic acid metabolism, and then the rate constants were estimated with a data set from ATP-stimulated BMDM, using a two-step matrix-based approach employing a constrained least-squares method followed by nonlinear optimization. The robustness of the model was validated through parametric sensitivity, uncertainty analysis, and predicting an independent dataset from KLA-primed ATP-stimulated BMDM by allowing the parameters to vary within the uncertainty range of the calculated parameters. We analyzed the functional coupling between COX isozymes and terminal enzymes by developing a PGH2-divided model. This provided evidence for the functional coupling between COX-2 and PGE2 synthase, between COX-1/COX-2 and PGD2 synthase, and also between COX-1 and thromboxane A2 synthase. Further, these functional couplings were experimentally validated using COX-1 and COX-2 selective inhibitors. The resulting fluxomics analysis demonstrates that the "multi-omics" systems biology approach can define the complex machinery of eicosanoid networks.