Coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) poses a mortal threat to human health. The elucidation of the relationship between peripheral immune cells and the development of inflammation is essential for revealing the pathogenic mechanism of COVID-19 and developing related antiviral drugs. The immune cell metabolism-targeting therapies exhibit a desirable anti-inflammatory effect in some treatment cases. In this study, based on differentially expressed gene (DEG) analysis, a genome-scale metabolic model (GSMM) was reconstructed by integrating transcriptome data to characterize the adaptive metabolic changes in peripheral blood mononuclear cells (PBMCs) in severe COVID-19 patients. Differential flux analysis revealed that metabolic changes such as enhanced aerobic glycolysis, impaired oxidative phosphorylation, fluctuating biogenesis of lipids, vitamins (folate and retinol), and nucleotides played important roles in the inflammation adaptation of PBMCs. Moreover, the main metabolic enzymes such as the solute carrier (SLC) family 2 member 3 (SLC2A3) and fatty acid synthase (FASN), responsible for the reactions with large differential fluxes, were identified as potential therapeutic targets. Our results revealed the inflammation regulation potentials of partial metabolic reactions with differential fluxes and their metabolites. This study provides a reference for developing potential PBMC metabolism-targeting therapy strategies against COVID-19.

The network module-based method has been used for drug repositioning. The traditional drug repositioning method only uses the gene characteristics of the drug but ignores the drug-triggered metabolic changes. The metabolic network systematically characterizes the connection between genes, proteins, and metabolic reactions. The differential metabolic flux distribution, as drug metabolism characteristics, was employed to cluster the agents with similar MoAs (mechanism of action). In this study, agents with the same pharmacology were clustered into one group, and a total of 1309 agents from the CMap database were clustered into 98 groups based on differential metabolic flux distribution. Transcription factor (TF) enrichment analysis revealed the agents in the same group (such as group 7 and group 26) were confirmed to have similar MoAs. Through this agent clustering strategy, the candidate drugs which can inhibit (Japanese encephalitis virus) JEV infection were identified. This study provides new insights into drug repositioning and their MoAs.