Global metabolic effects of glycerol kinase overexpression in rat hepatoma cells

Transcription profiling reveals tissue-specific metabolic pathways in the fat body and ovary of the diapausing mosquito Culex pipiens

The northern house mosquito, Culex pipiens, employs diapause as an essential survival strategy during winter, inducing important phenotypic changes such as enhanced stress tolerance, lipid accumulation, and extended longevity. During diapause, the cessation of reproductive development represents another distinctive phenotypic change, underlining the need for adjusted modulation of gene expressions within the ovary. Although considerable advancements in screening gene expression profiles in diapausing and non-diapausing mosquitoes, there remains a gap in tissue-specific transcriptomic profiling that could elucidate the complicated formation of diverse diapause features in Cx. pipiens. Here, we filled this gap by utilizing RNA sequencing, providing a detailed examination of gene expression patterns in the fat body and ovary during diapause compared to non-diapause conditions. Functional annotation of upregulated genes identified associations with carbohydrate metabolism, stress tolerance, immunity, and epigenetic regulation. The validation of candidate genes using quantitative real-time PCR verified the differentially expressed genes identified in diapausing mosquitoes. Our findings contribute novel insights into potential regulators during diapause in Cx. pipiens, thereby opening possible avenues for developing innovative vector control strategies.

Exploring a Novel Role of Glycerol Kinase 1 in Prostate Cancer PC-3 Cells

Clinically, prostate cancer is infamous for its histological and molecular heterogeneity, which causes great challenges to pinpoint therapy and pharmaceutical development. To overcome these difficulties, researchers are focusing on modulating tumor microenvironment and immune responses in addition to genetic alteration and epigenetic regulation. Here, we aimed to identify potential biomarkers or modulators of prostate cancer by investigating genes specifically altered in prostate cancer cells treated with established anti-cancer agents. Glycerol kinase 1 (GK1) is phosphotransferase encoded on the X chromosome, is associated with the synthesis of triglycerides and glycerophospholipids, and has been mainly studied for X-linked metabolic disorder GK deficiency (GKD). Interestingly, our DNA microarray analysis showed that several anti-cancer agents highly induced the expression of GK1, especially GK1a and GK1b isoforms, in human prostate cancer PC-3 cells. To elucidate the relationship between GK1 and cancer cell death, a human GK1b-specific expression vector was constructed and transfected into the PC-3 cells. Surprisingly, GK1b overexpression dramatically reduced cell viability and significantly accelerated apoptotic cell death. These findings suggest that GK1b may serve as a promising modulator and biomarker of cell death in prostate cancer, offering potential avenues for therapeutic intervention.

Glycerol kinase enzyme is a prognostic predictor in esophageal carcinoma and is associated with immune cell infiltration

The influence of lipid metabolism on tumorigenesis and progression has garnered significant attention. However, the role of Glycerol Kinase (GK), a key enzyme in glycerol metabolism, in Esophageal Carcinoma (ESCA) remains unclear. To further elucidate the relationship between GK and ESCA, we investigated GK expression levels using database information. Controlled studies employing immunohistochemistry were conducted on clinical ESCA tumor samples and normal specimens, confirming GK's elevated expression in ESCA. Analysis of The Cancer Genome Atlas (TCGA) data via Kaplan-Meier (KM) survival plots revealed that increased GK expression correlates with poorer ESCA patient outcomes, particularly in overall survival (OS) and disease-specific survival (DSS). Multiple regression analysis indicated that elevated GK expression is an independent risk factor affecting ESCA prognosis. Statistical analysis of prognostic data from clinical samples further corroborated this finding. Moreover, there appears to be a significant correlation between GK expression and immune infiltration, specifically involving certain T and B lymphocytes. In conclusion, elevated GK expression in ESCA is strongly linked to poor prognosis and increased immune cell infiltration, highlighting its potential as an independent prognostic biomarker and a viable therapeutic target.

Study on the acute toxicity of sodium taurocholate via zebrafish mortality, behavioral response, and NMR-metabolomics analysis

Sodium taurocholate (NaT) is a hydrophobic bile salt that exhibits varying toxicity and antimicrobial activity. The accumulation of BSs during their entero-hepatic cycle causes cytotoxicity in the liver and intestine and could also alter the intestinal microbiome leading to various diseases. In this research, the acute toxicity of sodium taurocholate in different concentrations (3000 mg/L, 1500 mg/L, 750 mg/L, 375 mg/L, and 0 mg/L) was investigated on four months old zebrafish by immersion in water for 96 h. The results were determined based on the fish mortality, behavioral response, and NMR metabolomics analysis which revealed LC50 of 1760.32 mg/L and 1050.42 mg/L after 72 and 96 h treatment, respectively. However, the non-lethal NaT concentrations of 750 mg/L and 375 mg/L at 96 h exposure significantly (p ≤ 0.05) decreased the total distance traveled and the activity duration, also caused surface respiration on the zebrafish. Orthogonal Projections to Latent Structures Discriminant Analysis (OPLS-DA) revealed that the metabolome of the fish treated with 750 mg/L was discriminated from that of the control by PC1. Major significantly downregulated metabolites by NaT-induction include valine, isoleucine, 2-hydroxyvalerate, glycine, glycerol, choline, glucose, pyruvate, anserine, threonine, carnitine and homoserine. On the contrary, taurine, creatine, lactate, acetate and 3-hydroxybutyrate were upregulated suggesting cellular consumption of lipids, glucose and amino acids for adenosine triphosphate (ATP) generation during immune and inflammatory response. whereby these metabolites were released in the process. In conclusion, the research revealed the toxic effect of NaT and its potential to trigger changes in zebrafish metabolism.

Genome wide transcriptome analysis provides bases on hepatic lipid metabolism disorder affected by increased dietary grain ratio in fattening lambs

BACKGROUND

The liver is a principal metabolic organ and has a major role in regulating lipid metabolism. With the development of rapidly fattening livestock in the modern breeding industry, the incidence of hepatic steatosis and accumulation in animals was significantly increased. However, the molecular mechanisms responsible for hepatic lipid metabolic disturbances in a high concentrate diet remain unclear. The objective of this study was to evaluate the effects of increasing concentrate level in a fattening lamb diet on biochemical indices, hepatic triglycerides (TG) concentration, and hepatic transcriptomic profiles. In the present study, 42 weaned lambs (about 3 ± 0.3 months old) were randomly assigned to the GN60 group (60% concentrate of dry matter, GN60, n = 21) or GN70 group (70% concentrate of dry matter, n = 21) for a 3-months feeding trial.

RESULTS

No difference was observed in the growth performance or plasma biochemical parameters between the GN60 group and the GN70 group. The hepatic TG concentration was higher in the GN70 group than GN60 group (P < 0.05). Hepatic transcriptomic analysis showed that there were 290 differentially expressed genes identified between GN60 and GN70 groups, with 125 genes up-regulated and 165 genes down-regulated in the GN70 group. The enriched Gene Ontology (GO) items and KEGG pathways and protein-protein interaction (PPI) network of differentially expressed genes (DEGs) revealed that the majority of enriched pathways were related to lipid metabolism. Further analysis revealed that the fatty acid synthesis was up-regulated, while fatty acid transport, oxidation, and TG degradation were down-regulated in the GN70 group when compared with the GN60 group.

CONCLUSIONS

These results indicated that GN70 induced excess lipid deposition in the liver of lambs during the fattening period, with high synthesis rates and low degradation rates of TG. The identified mechanisms may help understand hepatic metabolism in lambs with a high concentrate diet and provide insight into decreasing the risk of liver metabolism disorder in animals.

Isotope-Assisted Metabolic Flux Analysis: A Powerful Technique to Gain New Insights into the Human Metabolome in Health and Disease

Cell metabolism represents the coordinated changes in genes, proteins, and metabolites that occur in health and disease. The metabolic fluxome, which includes both intracellular and extracellular metabolic reaction rates (fluxes), therefore provides a powerful, integrated description of cellular phenotype. However, intracellular fluxes cannot be directly measured. Instead, flux quantification requires sophisticated mathematical and computational analysis of data from isotope labeling experiments. In this review, we describe isotope-assisted metabolic flux analysis (iMFA), a rigorous computational approach to fluxome quantification that integrates metabolic network models and experimental data to generate quantitative metabolic flux maps. We highlight practical considerations for implementing iMFA in mammalian models, as well as iMFA applications in in vitro and in vivo studies of physiology and disease. Finally, we identify promising new frontiers in iMFA which may enable us to fully unlock the potential of iMFA in biomedical research.

Recent findings in the regulation of G6PD and its role in diseases

Glucose-6-phosphate dehydrogenase (G6PD) is the only rate-limiting enzyme in the pentose phosphate pathway (PPP). Rapidly proliferating cells require metabolites from PPP to synthesize ribonucleotides and maintain intracellular redox homeostasis. G6PD expression can be abnormally elevated in a variety of cancers. In addition, G6PD may act as a regulator of viral replication and vascular smooth muscle function. Therefore, G6PD-mediated activation of PPP may promote tumor and non-neoplastic disease progression. Recently, studies have identified post-translational modifications (PTMs) as an important mechanism for regulating G6PD function. Here, we provide a comprehensive review of various PTMs (e.g., phosphorylation, acetylation, glycosylation, ubiquitination, and glutarylation), which are identified in the regulation of G6PD structure, expression and enzymatic activity. In addition, we review signaling pathways that regulate G6PD and evaluate the role of oncogenic signals that lead to the reprogramming of PPP in tumor and non-neoplastic diseases as well as summarize the inhibitors that target G6PD.

Hepatic Glycerol Metabolism-Related Genes in Carnivorous Rainbow Trout (Oncorhynchus mykiss): Insights Into Molecular Characteristics, Ontogenesis, and Nutritional Regulation

Glycerol metabolism in rainbow trout is poorly studied even though it is at the interface between lipid and glucose metabolism. Moreover, glycerol can be an important ingredient in new aquafeed formulation to decrease the catabolism of dietary amino acids. Thus, the present study aimed to characterize for the first time the different genes coding for key enzymes and proteins involved in hepatic glycerol metabolism. From the trout genomes, all the paralogous genes coding for glycerol transport (aqp9b), glycerol kinase (gk2a and gk5), glycerol-3-phosphate phosphatase (pgp), and glycerol-3-phosphate dehydrogenase (gpd1a, gpd1b, and gpd1c) were identified. The ontogenesis determined that the capacity to metabolize glycerol begins with the apparition of the liver during the development (stage 22) and are more expressed at the endogenous–exogenous feeding period (stage 35). The postprandial regulation of the expression of these genes in juvenile trout showed that the postprandial peak of expression is between 4 and 24 h after the last meal for many of the genes, demonstrating that glycerol metabolism could be nutritionally regulated at a molecular level. However, surprisingly, no regulation of the mRNA abundance for the glycerol metabolism-related genes by different levels of dietary glycerol (0, 2.5, and 5%) have been detected, showing that hepatic glycerol metabolism is poorly regulated at a molecular level by dietary glycerol in rainbow trout juveniles.

Glycerol kinase enhances hepatic lipid metabolism by repressing nuclear receptor subfamily 4 group A1 in the nucleus

Glycerol kinase (GYK) plays a critical role in hepatic metabolism by converting glycerol to glycerol 3-phosphate in an ATP-dependent reaction. GYK isoform b is the only glycerol kinase present in whole cells, and has a non-enzymatic moonlighting function in the nucleus. GYK isoform b acts as a co-regulator of nuclear receptor subfamily 4 group A1 (NR4A1) and participates in the regulation of hepatic glucose metabolism by protein–protein interaction with NR4A1. Herein, GYK expression was found to upregulate the expression of NR4A1-mediated lipid metabolism-related genes (SREBP1C, FASN, ACACA, and GPAM) in HEK293T and L02 cells, and in mouse in vivo studies. GYK expression increased blood levels of cholesterol, triglyceride, and high-density lipoprotein cholesterol, but not low-density lipoprotein cholesterol levels. It enhanced the transcriptional activity of Nr4a1 target genes by negatively cooperating with NR4A1 and its enzymatic activity or by other undefined moonlighting functions. This enhancement was observed in both normal and diabetic mice. We also found a feed-forward regulation loop between GYK and NR4A1, serving as part of a GYK-NR4A1 regulatory mechanism in hepatic metabolism. Thus, GYK regulates the effect of NR4A1 on hepatic lipid metabolism in normal and diabetic mice, partially through the cooperation of GYK and NR4A1.

Global liver proteomic analysis of Wistar rats chronically exposed to low-levels of bisphenol A and S

Exposure to bisphenol A (BPA) and bisphenol S (BPS) has been associated with the development of metabolic disorders, such as obesity, dyslipidemias, and nonalcoholic fatty liver disease. Nonetheless, the associated mechanisms are still not fully understood. BPS is being used with no restrictions to replace BPA, which increases the concern regarding its safety and claims for further investigation on its potential mechanisms of toxicity. The present study aims to access liver molecular disturbances which could be associated with systemic metabolic disorders following exposure to BPA or BPS. Therefore, body weight gain and serum biochemical parameters were measured in male Wistar rats chronically exposed to 50 or 500 µg/kg/day of BPA or BPS, while an extensive evaluation of liver protein expression changes was conducted after exposure to 50 µg/kg/day of both compounds. Exposure to the lowest dose of BPA led to the development of hyperglycemia and hypercholesterolemia, while the BPS lowest dose led to the development of hypertriglyceridemia. Besides, exposure to 500 µg/kg/day of BPS significantly increased body weight gain and LDL-cholesterol levels. Hepatic proteins differentially expressed in BPA and BPS-exposed groups compared to the control group were mostly related to lipid metabolism and synthesis, with upregulation of glucokinase activity-related sequence 1 (1.8-fold in BPA and 2.4-fold in BPS), which is involved in glycerol triglycerides synthesis, and hydroxymethylglutaryl-CoA synthase cytoplasmic (2-fold in BPS), an enzyme involved in mevalonate biosynthesis. Essential mitochondrial proteins of the electron transport chain were upregulated after exposure to both contaminants. Also, BPA and BPS dysregulated expression of liver antioxidant enzymes, which are involved in cellular reactive oxygen species detoxification. Altogether, the results of the present study contribute to expand the scientific understanding of how BPA and BPS lead to the development of metabolic disorders and reinforce the risks associated with exposure to these contaminants.

Oncogenic KRAS Induces NIX-Mediated Mitophagy to Promote Pancreatic Cancer

Activating KRAS mutations are found in nearly all cases of pancreatic ductal adenocarcinoma (PDAC), yet effective clinical targeting of oncogenic KRAS remains elusive. Understanding of KRAS-dependent PDAC-promoting pathways could lead to the identification of vulnerabilities and the development of new treatments. We show that oncogenic KRAS induces BNIP3L/NIX expression and a selective mitophagy program that restricts glucose flux to the mitochondria and enhances redox capacity. Loss of Nix restores functional mitochondria to cells, increasing demands for NADPH reducing power and decreasing proliferation in glucose-limited conditions. Nix deletion markedly delays progression of pancreatic cancer and improves survival in a murine (KPC) model of PDAC. Although conditional Nix ablation in vivo initially results in the accumulation of mitochondria, mitochondrial content eventually normalizes via increased mitochondrial clearance programs, and pancreatic intraepithelial neoplasia (PanIN) lesions progress to PDAC. We identify the KRAS-NIX mitophagy program as a novel driver of glycolysis, redox robustness, and disease progression in PDAC. SIGNIFICANCE: NIX-mediated mitophagy is a new oncogenic KRAS effector pathway that suppresses functional mitochondrial content to stimulate cell proliferation and augment redox homeostasis. This pathway promotes the progression of PanIN to PDAC and represents a new dependency in pancreatic cancer.This article is highlighted in the In This Issue feature, p. 1143.

Glycerol kinase interacts with nuclear receptor NR4A1 and regulates glucose metabolism in the liver

Glycerol kinase (Gyk), consisting of 4 isoforms, plays a critical role in metabolism by converting glycerol to glycerol 3-phosphate in an ATP-dependent reaction. Only Gyk isoform b is present in whole cells, but its function in the nucleus remains elusive. Previous studies have shown that nuclear orphan receptor subfamily 4 group A member (NR4A)-1 is an important regulator of hepatic glucose homeostasis and lipid metabolism in adipose tissue. We aimed to elucidate the functional interaction between nuclear Gyk and NR4A1 during hepatic gluconeogenesis in the unfed state and diabetes. We identified nuclear Gyk as a novel corepressor of NR4A1 in the liver; moreover, this recruitment was dependent on the C-terminal ligand-binding domain instead of the N-terminal activation function 1 domain, which interacts with other NR4A1 coregulators. NR4A1 transcriptional activity was inhibited by Gyk via protein-protein interaction but not enzymatic activity. Moreover, Gyk overexpression suppressed NR4A1 ability to regulate the expression of target genes involved in hepatic gluconeogenesis in vitro and in vivo as well as blood glucose regulation, which was observed in both unfed and diabetic mice. These results highlight the moonlighting function of nuclear Gyk, which was found to act as a coregulator of NR4A1, participating in the regulation of hepatic glucose homeostasis in the unfed state and diabetes.-Miao, L., Yang, Y., Liu, Y., Lai, L., Wang, L., Zhan, Y., Yin, R., Yu, M., Li, C., Yang, X., Ge, C. Glycerol kinase interacts with nuclear receptor NR4A1 and regulates glucose metabolism in the liver.

The challenge and potential of photosynthesis: unique considerations for metabolic flux measurements in photosynthetic microorganisms

Photosynthetic microorganisms have the potential for sustainable production of chemical feedstocks and products but have had limited success due to a lack of tools and deeper understanding of metabolic pathway regulation. The application of instationary metabolic flux analysis (INST-MFA) to photosynthetic microorganisms has allowed researchers to quantify fluxes and identify bottlenecks and metabolic inefficiencies to improve strain performance or gain insight into cellular physiology. Additionally, flux measurements can also highlight deviations between measured and predicted fluxes, revealing weaknesses in metabolic models and highlighting areas where a lack of understanding still exists. In this review, we outline the experimental steps necessary to successfully perform photosynthetic flux experiments and analysis. We also discuss the challenges unique to photosynthetic microorganisms and how to account for them, including: light supply, quenching, concentration, extraction, analysis, and flux calculation. We hope that this will enable a larger number of researchers to successfully apply isotope assisted metabolic flux analysis (¹³C-MFA) to their favorite photosynthetic organism.

The ATP-stimulated translocation promoter (ASTP) activity of glycerol kinase plays central role in adipogenesis

Glycerol kinase (GK) is a multifunctional enzyme located at the interface of carbohydrate and fat metabolism. It contributes to both central carbon metabolism and adipogenesis; specifically, through its role as the ATP-stimulated translocation promoter (ASTP). GK overexpression leads to increased ASTP activity and increased fat storage in H4IIE cells. We performed metabolic flux analysis in human GK-overexpressing H4IIE cells and found that overexpressing cells had significantly altered fluxes through central carbon and lipid metabolism including increased flux through the pentose phosphate pathway and increased production of lipids. We also observed an equal contribution of glycerol to carbohydrate metabolism in all cell lines, suggesting that GK's alternate functions rather than its enzymatic function are important for these processes. To further elucidate the contributions of the enzymatic (phosphorylation) and alternative (ASTP) functions of GK in adipogenesis, we performed experiments on mammalian GK and E. coli GK. We determined that the ASTP function of GK (which is absent in E. coli GK) plays a greater role than the enzymatic activity in these processes. These studies further emphasize GK's diverse functionality and provides fundamental insights into the multiple protein functions of glycerol kinase.

Concurrent isotope-assisted metabolic flux analysis and transcriptome profiling reveal responses of poplar cells to altered nitrogen and carbon supply

Reduced nitrogen is indispensable to plants. However, its limited availability in soil combined with the energetic and environmental impacts of nitrogen fertilizers motivates research into molecular mechanisms toward improving plant nitrogen use efficiency (NUE). We performed a systems-level investigation of this problem by employing multiple 'omics methodologies on cell suspensions of hybrid poplar (Populus tremula × Populus alba). Acclimation and growth of the cell suspensions in four nutrient regimes ranging from abundant to deficient supplies of carbon and nitrogen revealed that cell growth under low-nitrogen levels was associated with substantially higher NUE. To investigate the underlying metabolic and molecular mechanisms, we concurrently performed steady-state ¹³ C metabolic flux analysis with multiple isotope labels and transcriptomic profiling with cDNA microarrays. The ¹³ C flux analysis revealed that the absolute flux through the oxidative pentose phosphate pathway (oxPPP) was substantially lower (~threefold) under low-nitrogen conditions. Additionally, the flux partitioning ratio between the tricarboxylic acid cycle and anaplerotic pathways varied from 84%:16% under abundant carbon and nitrogen to 55%:45% under deficient carbon and nitrogen. Gene expression data, together with the flux results, suggested a plastidic localization of the oxPPP as well as transcriptional regulation of certain metabolic branchpoints, including those between glycolysis and the oxPPP. The transcriptome data also indicated that NUE-improving mechanisms may involve a redirection of excess carbon to aromatic metabolic pathways and extensive downregulation of potentially redundant genes (in these heterotrophic cells) that encode photosynthetic and light-harvesting proteins, suggesting the recruitment of these proteins as nitrogen sinks in nitrogen-abundant conditions.

Steady-state and instationary modeling of proteinogenic and free amino acid isotopomers for flux quantification

Metabolic flux analysis (MFA) is a powerful tool for exploring and quantifying carbon traffic in metabolic networks. Accurate flux quantification requires (1) high-quality isotopomer measurements, usually of biomass components including proteinogenic/free amino acids or central carbon metabolites, and (2) a mathematical model that relates the unknown fluxes to the measured isotopomers. Modeling requires a thorough knowledge of the structure of the underlying metabolic network, often available from many databases, as well as the ability to make reasonable assumptions that will enable simplification of the model. Here we describe a general methodology underlying computer-aided mathematical modeling of a flux-isotopomer relationship and some of the accompanying data-processing steps. One of two modeling strategies will need to be employed, depending on the type of isotope labeling experiment performed. These strategies-steady-state modeling and instationary modeling-have different experimental and computational demands. We discuss the concepts underlying these two types of modeling and demonstrate steady-state modeling in a step-by-step manner. Our methodology should be applicable to most isotope-assisted MFA applications and should serve as a general framework applicable to many realistic metabolic networks with little modification.

Flux and reflux: metabolite reflux in plant suspension cells and its implications for isotope-assisted metabolic flux analysis

Isotope-assisted metabolic flux analysis (MFA) is a powerful methodology to quantify intracellular fluxes via isotope labeling experiments (ILEs). In batch cultures, which are often convenient, inexpensive or inevitable especially for eukaryotic systems, MFA is complicated by the presence of the initially present biomass. This unlabeled biomass may either mix with the newly synthesized labeled biomass or reflux into the metabolic network, thus masking the true labeling patterns in the newly synthesized biomass. Here, we report a detailed investigation of such metabolite reflux in cell suspensions of the tree poplar. In ILEs supplying 28% or 98% U-(13)C glucose as the sole organic carbon source, biomass components exhibited lower (13)C enrichments than the supplied glucose as well as anomalous isotopomers not explainable by simple mixing of the initial and newly synthesized biomass. These anomalous labeling patterns were most prominent in a 98% U-(13)C glucose ILE. By comparing the performance of light- and dark-grown cells as well as by analyzing the isotope labeling patterns in aspartic and glutamic acids, we eliminated photosynthetic or anaplerotic fixation of extracellular (12)CO2 as explanations for the anomalous labeling patterns. We further investigated four different metabolic models for interpreting the labeling patterns and evaluating fluxes: (i) a carbon source (glucose) dilution model, (ii) an isotopomer correction model with uniform dilution for all amino acids, (iii) an isotopomer correction model with variable dilution for different amino acids, and (iv) a comprehensive metabolite reflux model. Of these, the metabolite reflux model provided a substantially better fit for the observed labeling patterns (sum of squared residues: 538) than the other three models whose sum of squared residues were (i) 4626, (ii) 4983, and (iii) 1748, respectively. We compared fluxes determined using the metabolite reflux model to those determined using an independent methodology involving an excessively long ILE to wash out initial biomass and a minimal reflux model. This comparison showed identical or similar distributions for a majority of fluxes, thus validating our comprehensive reflux model. In summary, we have demonstrated the need for quantifying interactions between initially present biomass and newly synthesized biomass in batch ILEs, especially through the use of ≈100% U-(13)C carbon sources. Our ILEs reveal a high amount of metabolite reflux in poplar cell suspensions, which is well explained by a comprehensive metabolite reflux model.

Mathematical modeling of isotope labeling experiments for metabolic flux analysis

Isotope labeling experiments (ILEs) offer a powerful methodology to perform metabolic flux analysis. However, the task of interpreting data from these experiments to evaluate flux values requires significant mathematical modeling skills. Toward this, this chapter provides background information and examples to enable the reader to (1) model metabolic networks, (2) simulate ILEs, and (3) understand the optimization and statistical methods commonly used for flux evaluation. A compartmentalized model of plant glycolysis and pentose phosphate pathway illustrates the reconstruction of a typical metabolic network, whereas a simpler example network illustrates the underlying metabolite and isotopomer balancing techniques. We also discuss the salient features of commonly used flux estimation software 13CFLUX2, Metran, NMR2Flux+, FiatFlux, and OpenFLUX. Furthermore, we briefly discuss methods to improve flux estimates. A graphical checklist at the end of the chapter provides a reader a quick reference to the mathematical modeling concepts and resources.

Experimental evidence and isotopomer analysis of mixotrophic glucose metabolism in the marine diatom Phaeodactylum tricornutum

BACKGROUND

Heterotrophic fermentation using simple sugars such as glucose is an established and cost-effective method for synthesizing bioproducts from bacteria, yeast and algae. Organisms incapable of metabolizing glucose have limited applications as cell factories, often despite many other advantageous characteristics. Therefore, there is a clear need to investigate glucose metabolism in potential cell factories. One such organism, with a unique metabolic network and a propensity to synthesize highly reduced compounds as a large fraction of its biomass, is the marine diatom Phaeodactylum tricornutum (Pt). Although Pt has been engineered to metabolize glucose, conflicting lines of evidence leave it unresolved whether Pt can natively consume glucose.

RESULTS

Isotope labeling experiments in which Pt was mixotrophically grown under light on 100% U-(13)C glucose and naturally abundant (~99% (12)C) dissolved inorganic carbon resulted in proteinogenic amino acids with an average 13C-enrichment of 88%, thus providing convincing evidence of glucose uptake and metabolism. The dissolved inorganic carbon was largely incorporated through anaplerotic rather than photosynthetic fixation. Furthermore, an isotope labeling experiment utilizing 1-(13)C glucose and subsequent metabolic pathway analysis indicated that (i) the alternative Entner-Doudoroff and Phosphoketolase glycolytic pathways are active during glucose metabolism, and (ii) during mixotrophic growth, serine and glycine are largely synthesized from glyoxylate through photorespiratory reactions rather than from 3-phosphoglycerate. We validated the latter result for mixotrophic growth on glycerol by performing a 2-(13)C glycerol isotope labeling experiment. Additionally, gene expression assays showed that known, native glucose transporters in Pt are largely insensitive to glucose or light, whereas the gene encoding cytosolic fructose bisphosphate aldolase 3, an important glycolytic enzyme, is overexpressed in light but insensitive to glucose.

CONCLUSION

We have shown that Pt can use glucose as a primary carbon source when grown in light, but cannot use glucose to sustain growth in the dark. We further analyzed the metabolic mechanisms underlying the mixotrophic metabolism of glucose and found isotopic evidence for unusual pathways active in Pt. These insights expand the envelope of Pt cultivation methods using organic substrates. We anticipate that they will guide further engineering of Pt towards sustainable production of fuels, pharmaceuticals, and platform chemicals.

Metabolic analyses elucidate non-trivial gene targets for amplifying dihydroartemisinic acid production in yeast

Synthetic biology enables metabolic engineering of industrial microbes to synthesize value-added molecules. In this, a major challenge is the efficient redirection of carbon to the desired metabolic pathways. Pinpointing strategies toward this goal requires an in-depth investigation of the metabolic landscape of the organism, particularly primary metabolism, to identify precursor and cofactor availability for the target compound. The potent antimalarial therapeutic artemisinin and its precursors are promising candidate molecules for production in microbial hosts. Recent advances have demonstrated the production of artemisinin precursors in engineered yeast strains as an alternative to extraction from plants. We report the application of in silico and in vivo metabolic pathway analyses to identify metabolic engineering targets to improve the yield of the direct artemisinin precursor dihydroartemisinic acid (DHA) in yeast. First, in silico extreme pathway (ExPa) analysis identified NADPH-malic enzyme and the oxidative pentose phosphate pathway (PPP) as mechanisms to meet NADPH demand for DHA synthesis. Next, we compared key DHA-synthesizing ExPas to the metabolic flux distributions obtained from in vivo (13)C metabolic flux analysis of a DHA-synthesizing strain. This comparison revealed that knocking out ethanol synthesis and overexpressing glucose-6-phosphate dehydrogenase in the oxidative PPP (gene YNL241C) or the NADPH-malic enzyme ME2 (YKL029C) are vital steps toward overproducing DHA. Finally, we employed in silico flux balance analysis and minimization of metabolic adjustment on a yeast genome-scale model to identify gene knockouts for improving DHA yields. The best strategy involved knockout of an oxaloacetate transporter (YKL120W) and an aspartate aminotransferase (YKL106W), and was predicted to improve DHA yields by 70-fold. Collectively, our work elucidates multiple non-trivial metabolic engineering strategies for improving DHA yield in yeast.

Isotopically nonstationary 13C flux analysis of Myc-induced metabolic reprogramming in B-cells

We assessed several methods of (13)C metabolic flux analysis (MFA) and found that isotopically nonstationary MFA achieved maximum flux resolution in cultured P493-6 B-cells, which have been engineered to provide tunable expression of the Myc oncoprotein. Comparison of metabolic flux maps obtained under oncogenic (High) and endogenous (Low) Myc expression levels revealed network-wide reprogramming in response to ectopic Myc expression. High Myc cells relied more heavily on mitochondrial oxidative metabolism than Low Myc cells and globally upregulated their consumption of amino acids relative to glucose. TCA cycle and amphibolic mitochondrial pathways exhibited 2- to 4-fold flux increases in High Myc cells, in contrast to modest increases in glucose uptake and lactate excretion. Because our MFA approach relied exclusively upon isotopic measurements of protein-bound amino acids and RNA-bound ribose, it is readily applicable to more complex tumor models that are not amenable to direct extraction and isotopic analysis of free intracellular metabolites.

Genotype-specific responses in Atlantic salmon (Salmo salar) subject to dietary fish oil replacement by vegetable oil: a liver transcriptomic analysis

Background

Expansion of aquaculture is seriously limited by reductions in fish oil (FO) supply for aquafeeds. Terrestrial alternatives such as vegetable oils (VO) have been investigated and recently a strategy combining genetic selection with changes in diet formulations has been proposed to meet growing demands for aquaculture products. This study investigates the influence of genotype on transcriptomic responses to sustainable feeds in Atlantic salmon.

Results

A microarray analysis was performed to investigate the liver transcriptome of two family groups selected according to their estimated breeding values (EBVs) for flesh lipid content, 'Lean' or 'Fat', fed diets containing either FO or a VO blend. Diet principally affected metabolism genes, mainly of lipid and carbohydrate, followed by immune response genes. Genotype had a much lower impact on metabolism-related genes and affected mostly signalling pathways. Replacement of dietary FO by VO caused an up-regulation of long-chain polyunsaturated fatty acid biosynthesis, but there was a clear genotype effect as fatty acyl elongase (elovl2) was only up-regulated and desaturases (Δ5 fad and Δ6 fad) showed a higher magnitude of response in Lean fish, which was reflected in liver fatty acid composition. Fatty acid synthase (FAS) was also up-regulated by VO and the effect was independent of genotype. Genetic background of the fish clearly affected regulation of lipid metabolism, as PPARα and PPARβ were down-regulated by the VO diet only in Lean fish, while in Fat salmon SREBP-1 expression was up-regulated by VO. In addition, all three genes had a lower expression in the Lean family group than in the Fat, when fed VO. Differences in muscle adiposity between family groups may have been caused by higher levels of hepatic fatty acid and glycerophospholipid synthesis in the Fat fish, as indicated by the expression of FAS, 1-acyl-sn-glycerol-3-phosphate acyltransferase and lipid phosphate phosphohydrolase 2.

Conclusions

This study has identified metabolic pathways and key regulators that may respond differently to alternative plant-based feeds depending on genotype. Further studies are required but data suggest that it will be possible to identify families better adapted to alternative diet formulations that might be appropriate for future genetic selection programmes.

Tracking cellular metabolomics in lipoapoptosis- and steatosis-developing liver cells

Palmitate (PA) is known to induce reactive oxygen species (ROS) formation and apoptosis in liver cells, whereas concurrent treatment of oleate (OA) with PA predominately induces steatosis without ROS in liver cells. We previously reported that PA treatment induces the decoupling of glycolysis and tricarboxylic acid cycle (TCA cycle) fluxes, but OA co-treatment restored most metabolic fluxes to their control levels. However, the mechanisms by which metabolites are linked to metabolic fluxes and subsequent lipoapoptotic or steatotic phenotypes remain unclear. To determine the link, we used GC-MS-based polar and non-polar metabolic profiling in lipoapoptosis- or steatosis-developing H4IIEC3 hepatoma cells, to examine the metabolome at different time points after treatment with either PA alone (PA cells) or both PA and OA (PA/OA cells). Metabolic profiles revealed various changes in metabolite levels for TCA cycle intermediates, pentose phosphate pathway (PPP) intermediates, and energy storage metabolites between PA and PA/OA cells. For example, adenosine was markedly increased only in PA cells, whereas gluconate was increased in PA/OA cells. To assess the interaction among these metabolites, the metabolite-to-metabolite correlations were calculated and correlation networks were visualized. These correlation networks demonstrate that a dissociation among PPP metabolites was introduced in PA-treated cells, and this dissociation was restored in PA/OA-treated cells. Thus, our data suggest that abnormal PPP fluxes, in addition to increased adenosine levels, might be related to the decoupling of glycolysis and the resulting lipoapoptotic phenotype.

Active and inactive metabolic pathways in tumor spheroids: determination by GC-MS

Active metabolic pathways in three-dimensional cancer-cell cultures are potential chemotherapeutic targets that would be effective throughout tumors. Chaotic vasculature creates cellular regions in tumors with distinct metabolic behavior that are only present in aggregate cell masses. To quantify cancer cell metabolism, transformed mouse fibroblasts were grown as spheroids and fed isotopically labeled culture medium. Metabolite uptake and production rates were measured as functions of time. Gas chromatography-mass spectrometry was used to quantify the extent of labeling on amino acids present in cytoplasmic extracts. The labeling pattern identified several active and inactive metabolic pathways: Glutaminolysis was found to be active, and malic enzyme and gluconeogenesis were inactive. Transformed cells in spheroids were also found to actively synthesize serine, cysteine, alanine, aspartate, glutamate, and proline; and not synthesize glutamine. The activities of these pathways suggest that cancer cells consume glutamine for biosynthesis and not to provide cellular energy. Determining active metabolic pathways indicates how cells direct carbon flow and may lead to the discovery of novel molecular targets for anticancer therapy.

Moonlighting function of glycerol kinase causes systems-level changes in rat hepatoma cells

Glycerol kinase (GK) is an enzyme with diverse (moonlighting) cellular functions. GK overexpression affects central metabolic fluxes substantially; therefore, to elucidate the mechanism underlying these changes, we employed a systems-level evaluation of GK overexpression in H4IIE rat hepatoma cells. Microarray analysis revealed altered expression of genes in metabolism (central carbon and lipid), which correlated with previous flux analysis, and of genes regulated by the glucocorticoid receptor (GR). Oil Red O staining showed that GK overexpression leads to increased fat storage in H4IIE cells. Network component analysis revealed that activities of peroxisome proliferator-activated receptor alpha, GR, and seven other transcription factors were altered by GK overexpression. The increased activity of GR was experimentally verified by quantitative RT-PCR of GR-responsive genes in the presence and absence of the glucocorticoid agonist, dexamethasone. This systems biology approach further emphasizes GK's essential role in central and lipid metabolism and experimentally verifies GK's alternative (moonlighting) function of affecting GR transcription factor activity.

Fenofibrate simultaneously induces hepatic fatty acid oxidation, synthesis, and elongation in mice

A growing body of evidence indicates that peroxisome proliferator-activated receptor alpha (PPARalpha) not merely serves as a transcriptional regulator of fatty acid catabolism but also exerts a much broader role in hepatic lipid metabolism. We determined adaptations in hepatic lipid metabolism and related aspects of carbohydrate metabolism upon treatment of C57Bl/6 mice with the PPARalpha agonist fenofibrate. Stable isotope procedures were applied to assess hepatic fatty acid synthesis, fatty acid elongation, and carbohydrate metabolism. Fenofibrate treatment strongly induced hepatic de novo lipogenesis and chain elongation (+/-300, 150, and 600% for C16:0, C18:0, and C18:1 synthesis, respectively) in parallel with an increased expression of lipogenic genes. The lipogenic induction in fenofibrate-treated mice was found to depend on sterol regulatory element-binding protein 1c (SREBP-1c) but not carbohydrate response element-binding protein (ChREBP). Fenofibrate treatment resulted in a reduced contribution of glycolysis to acetyl-CoA production, whereas the cycling of glucose 6-phosphate through the pentose phosphate pathway presumably was enhanced. Altogether, our data indicate that beta-oxidation and lipogenesis are induced simultaneously upon fenofibrate treatment. These observations may reflect a physiological mechanism by which PPARalpha and SREBP-1c collectively ensure proper handling of fatty acids to protect the liver against cytotoxic damage.

Evaluation of 13C isotopic tracers for metabolic flux analysis in mammalian cells

(13)C metabolic flux analysis (MFA) is the most comprehensive means of characterizing cellular metabolic states. Uniquely labeled isotopic tracers enable more focused analyses to probe specific reactions within the network. As a result, the choice of tracer largely determines the precision with which one can estimate metabolic fluxes, especially in complex mammalian systems that require multiple substrates. Here we have experimentally determined metabolic fluxes in a tumor cell line, successfully recapitulating the hallmarks of cancer cell metabolism. Using these data, we computationally evaluated specifically labeled (13)C glucose and glutamine tracers for their ability to precisely and accurately estimate fluxes in central carbon metabolism. These methods enabled us to identify the optimal tracer for analyzing individual fluxes, specific pathways, and central carbon metabolism as a whole. [1,2-(13)C(2)]glucose provided the most precise estimates for glycolysis, the pentose phosphate pathway, and the overall network. Tracers such as [2-(13)C]glucose and [3-(13)C]glucose also outperformed the more commonly used [1-(13)C]glucose. [U-(13)C(5)]glutamine emerged as the preferred isotopic tracer for the analysis of the tricarboxylic acid (TCA) cycle. These results provide valuable, quantitative information on the performance of (13)C-labeled substrates and can aid in the design of more informative MFA experiments in mammalian cell culture.

Transcriptomic and network component analysis of glycerol kinase in skeletal muscle using a mouse model of glycerol kinase deficiency

Glycerol kinase (GK) is at the interface of fat and carbohydrate metabolism and has been linked to obesity and type 2 diabetes mellitus (T2DM). The purpose of this study was to investigate the role of GK in fat metabolism and insulin signaling in skeletal muscle (an important end organ tissue in T2DM). Microarray analysis determined that there were 525 genes that were differentially expressed (1.2-fold, p value<0.05) between knockout (KO) and wild-type (WT) mice. Quantitative PCR (qPCR) confirmed the differential expression of genes including glycerol kinase (Gyk), phosphatidylinositol 3-kinase regulatory subunit, polypeptide 1 (p85 alpha) (Pik3r1), insulin-like growth factor 1 (Igf1), and growth factor receptor bound protein 2-associated protein 1 (Gab1). Network component analysis demonstrated that transcription factor activities of myogenic differentiation 1 (MYOD), myogenic regulatory factor 5 (MYF5), myogenin (MYOG), nuclear receptor subfamily 4, group A, member 1 (NUR77) are decreased in the Gyk KO whereas the activity of paired box 3 (PAX3) is increased. The activity of MYOD was confirmed using a DNA binding assay. In addition, myoblasts from Gyk KO had less ability to differentiate into myotubes compared to WT myoblasts. These findings support our previous studies in brown adipose tissue and demonstrate that the role of Gyk in muscle is due in part to its non-metabolic (moonlighting) activities.