Coordination of adaptive metabolism through cellular signaling networks and metabolic response is essential for balanced flow of energy and homeostasis. Post-translational modifications such as phosphorylation offer a rapid, efficient, and dynamic mechanism to regulate metabolic networks. Although numerous phosphorylation sites have been identified on metabolic enzymes, much remains unknown about their contribution to enzyme function and systemic metabolism. In this study, we stratify phosphorylation sites on metabolic enzymes based on their location with respect to functional and dimerization domains. Our analysis reveals that the majority of published phosphosites are on oxidoreductases, with particular enrichment of phosphotyrosine (pY) sites in proximity to binding domains for substrates, cofactors, active sites, or dimer interfaces. We identify phosphosites altered in obesity using a high fat diet (HFD) induced obesity model coupled to multiomics, and interrogate the functional impact of pY on hepatic metabolism. HFD induced dysregulation of redox homeostasis and reductive metabolism at the phosphoproteome and metabolome level in a sex-specific manner, which was reversed by supplementing with the antioxidant butylated hydroxyanisole (BHA). Partial least squares regression (PLSR) analysis identified pY sites that predict HFD or BHA induced changes of redox metabolites. We characterize predictive pY sites on glutathione S-transferase pi 1 (GSTP1), isocitrate dehydrogenase 1 (IDH1), and uridine monophosphate synthase (UMPS) using CRISPRi-rescue and stable isotope tracing. Our analysis revealed that sites on GSTP1 and UMPS inhibit enzyme activity while the pY site on IDH1 induces activity to promote reductive carboxylation. Overall, our approach provides insight into the convergence points where cellular signaling fine-tunes metabolism.
Keywords: GSTP1; IDH1; UMPS; computational modeling; metabolism; metabolomics; oxidative stress response; phosphoproteomics; reductive metabolism.