Bistability in glycolysis pathway as a physiological switch in energy metabolism

A comprehensive prognostic and immunological implications of PFKP in pan-cancer

BACKGROUND

Phosphofructokinase P (PFKP) is a key rate-limiting enzyme in glycolysis, playing a crucial role in various pathophysiological processes. However, its specific function in tumors remains unclear. This study aims to evaluate the expression and specific role of PFKP across multiple tumor types (Pan-cancer) and to explore its potential clinical significance as a therapeutic target in cancer treatment.

METHODS

We analyzed the expression of PFKP, immune cell infiltration, and patient prognosis across various cancers using data from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Additionally, we conducted a series of experiments in lung cancer cells, including Western blot, CCK-8 assay, colony formation assay, transwell migration assay, scratch wound healing assay, LDH release assay, and flow cytometry, to evaluate the impact of PFKP on tumor cells.

RESULTS

PFKP was found to be highly expressed in most cancers and identified as a prognostic risk factor. Elevated PFKP expression is associated with poorer clinical outcomes, particularly in lung adenocarcinoma (LUAD). Receiver operating characteristic (ROC) curve analysis indicated that PFKP can effectively differentiate between cancerous and normal tissues. The expression of PFKP in most tumors showed significant correlations with tumor mutational burden (TMB), microsatellite instability (MSI), immune score, and immune cell infiltration. In vitro experiments demonstrated that PFKP overexpression promotes lung cancer cell proliferation and migration while inhibiting apoptosis, whereas PFKP deficiency results in the opposite effects.

CONCLUSION

PFKP acts as an oncogene involved in tumorigenesis and may influence the immune microenvironment within the tumor. Our findings suggest that PFKP could serve as a potential biomarker for predicting prognosis and the efficacy of immunotherapy in tumors.

TPT disrupts early embryonic development and glucose metabolism of marine medaka in different salinites

Triphenyltin (TPT) is an organotin compound frequently detected in coastal estuaries, yet studies on TPT's effects in regions with significant salinity fluctuations, such as coastal estuaries, are currently limited. To investigate the toxic effects of TPT under different salinity conditions, this study focused on marine medaka (Oryzias melastigma) embryos. Through early morphological observations, RNA-seq analysis, biochemical marker assays, and qPCR detection, we explored the impact of TPT exposure on the early embryonic development of marine medaka under varying salinities. The study found that TPT exposure significantly increased embryo mortality at salinities of 0 ppt and 30 ppt. RNA-seq analysis revealed that TPT primarily affects glucose metabolism and glycogen synthesis processes in embryos. Under high salinity conditions, TPT may inhibit glucose metabolism by suppressing glycolysis and promoting gluconeogenesis. Furthermore, TPT exposure under different salinities led to the downregulation of genes associated with the insulin signaling pathway (ins, insra, irs2b, pik3ca, pdk1b, akt1, foxo1a), which may be linked to suppressed glucose metabolism and increased embryonic mortality. In summary, TPT exposure under different salinities affects the early development of marine medaka embryos and inhibits glucose metabolism. This study provides additional data to support research on organotin compounds in coastal estuaries.

Biological and genetic determinants of glycolysis: Phosphofructokinase isoforms boost energy status of stored red blood cells and transfusion outcomes

Mature red blood cells (RBCs) lack mitochondria and thus exclusively rely on glycolysis to generate adenosine triphosphate (ATP) during aging in vivo or storage in blood banks. Here, we leveraged 13,029 volunteers from the Recipient Epidemiology and Donor Evaluation Study to identify associations between end-of-storage levels of glycolytic metabolites and donor age, sex, and ancestry-specific genetic polymorphisms in regions encoding phosphofructokinase 1, platelet (detected in mature RBCs); hexokinase 1 (HK1); and ADP-ribosyl cyclase 1 and 2 (CD38/BST1). Gene-metabolite associations were validated in fresh and stored RBCs from 525 Diversity Outbred mice and via multi-omics characterization of 1,929 samples from 643 human RBC units during storage. ATP and hypoxanthine (HYPX) levels-and the genetic traits linked to them-were associated with hemolysis in vitro and in vivo, both in healthy autologous transfusion recipients and in 5,816 critically ill patients receiving heterologous transfusions, suggesting their potential as markers to improve transfusion outcomes.

Multi-cell-line learning for the data-driven construction of mechanistic metabolic models

Mammalian cells are commonly used as hosts in cell culture for biologics production in the pharmaceutical industry. Structured mechanistic models of metabolism have been used to capture complex cellular mechanisms that contribute to varying metabolic shifts in different cell lines. However, little research has focused on the impact of temporal changes in enzyme abundance and activity on the modeling of cell metabolism. In this work, we present a framework for constructing mechanistic models of metabolism that integrate growth-signaling control of enzyme activity and transcript dynamics. The proposed approach is applied to build models for three Chinese hamster ovary (CHO) cell lines using fed-batch culture data and time-series transcript profiles. Leveraging information from the transcriptome data, we develop a parameter estimation approach based on multi-cell-line (MCL) learning, which combines data sets from different cell lines and trains the individual cell-line models jointly to improve model accuracy. The computational results demonstrate the important role of growth signaling and transcript variability in metabolic models as well as the virtue of the MCL approach for constructing cell-line models with a limited amount of data. The resulting models exhibit a high level of accuracy in predicting distinct metabolic behaviors in the different cell lines; these models can potentially be used to accelerate the process and cell-line development for the biomanufacturing of new protein therapeutics.

Homeostasis in networks with multiple inputs

Homeostasis, also known as adaptation, refers to the ability of a system to counteract persistent external disturbances and tightly control the output of a key observable. Existing studies on homeostasis in network dynamics have mainly focused on 'perfect adaptation' in deterministic single-input single-output networks where the disturbances are scalar and affect the network dynamics via a pre-specified input node. In this paper we provide a full classification of all possible network topologies capable of generating infinitesimal homeostasis in arbitrarily large and complex multiple inputs networks. Working in the framework of 'infinitesimal homeostasis' allows us to make no assumption about how the components are interconnected and the functional form of the associated differential equations, apart from being compatible with the network architecture. Remarkably, we show that there are just three distinct 'mechanisms' that generate infinitesimal homeostasis. Each of these three mechanisms generates a rich class of well-defined network topologies-called homeostasis subnetworks. More importantly, we show that these classes of homeostasis subnetworks provides a topological basis for the classification of 'homeostasis types': the full set of all possible multiple inputs networks can be uniquely decomposed into these special homeostasis subnetworks. We illustrate our results with some simple abstract examples and a biologically realistic model for the co-regulation of calcium ( Ca ) and phosphate ( PO 4 ) in the rat. Furthermore, we identify a new phenomenon that occurs in the multiple input setting, that we call homeostasis mode interaction, in analogy with the well-known characteristic of multiparameter bifurcation theory.

Biological and Genetic Determinants of Glycolysis: Phosphofructokinase Isoforms Boost Energy Status of Stored Red Blood Cells and Transfusion Outcomes

Mature red blood cells (RBCs) lack mitochondria, and thus exclusively rely on glycolysis to generate adenosine triphosphate (ATP) during aging in vivo or storage in the blood bank. Here we leveraged 13,029 volunteers from the Recipient Epidemiology and Donor Evaluation Study to identify an association between end-of-storage levels of glycolytic metabolites and donor age, sex, and ancestry-specific genetic polymorphisms in regions encoding phosphofructokinase 1, platelet (detected in mature RBCs), hexokinase 1, ADP-ribosyl cyclase 1 and 2 (CD38/BST1). Gene-metabolite associations were validated in fresh and stored RBCs from 525 Diversity Outbred mice, and via multi-omics characterization of 1,929 samples from 643 human RBC units during storage. ATP and hypoxanthine levels - and the genetic traits linked to them - were associated with hemolysis in vitro and in vivo, both in healthy autologous transfusion recipients and in 5,816 critically ill patients receiving heterologous transfusions, suggesting their potential as markers to improve transfusion outcomes.

eTOC and Highlights

Highlights

Blood donor age and sex affect glycolysis in stored RBCs from 13,029 volunteers;Ancestry, genetic polymorphisms in PFKP, HK1, CD38/BST1 influence RBC glycolysis;Modeled PFKP effects relate to preventing loss of the total AXP pool in stored RBCs;ATP and hypoxanthine are biomarkers of hemolysis in vitro and in vivo.

Glycolysis in Chronic Liver Diseases: Mechanistic Insights and Therapeutic Opportunities

Chronic liver diseases (CLDs) cover a spectrum of liver diseases, ranging from nonalcoholic fatty liver disease to liver cancer, representing a growing epidemic worldwide with high unmet medical needs. Glycolysis is a conservative and rigorous process that converts glucose into pyruvate and sustains cells with the energy and intermediate products required for diverse biological activities. However, abnormalities in glycolytic flux during CLD development accelerate the disease progression. Aerobic glycolysis is a hallmark of liver cancer and is responsible for a broad range of oncogenic functions including proliferation, invasion, metastasis, angiogenesis, immune escape, and drug resistance. Recently, the non-neoplastic role of aerobic glycolysis in immune activation and inflammatory disorders, especially CLD, has attracted increasing attention. Several key mediators of aerobic glycolysis, including HIF-1α and pyruvate kinase M2 (PKM2), are upregulated during steatohepatitis and liver fibrosis. The pharmacological inhibition or ablation of PKM2 effectively attenuates hepatic inflammation and CLD progression. In this review, we particularly focused on the glycolytic and non-glycolytic roles of PKM2 in the progression of CLD, highlighting the translational potential of a glycolysis-centric therapeutic approach in combating CLD.

Chemical Memory with Discrete Turing Patterns Appearing in the Glycolytic Reaction

Memory is an essential element in information processing devices. We investigated a network formed by just three interacting nodes representing continuously stirred tank reactors (CSTRs) in which the glycolytic reaction proceeds as a potential realization of a chemical memory unit. Our study is based on the 2-variable computational model of the reaction. The model parameters were selected such that the system has a stable limit cycle and several distinct, discrete Turing patterns characterized by stationary concentrations at the nodes. In our interpretation, oscillations represent a blank memory unit, and Turing patterns code information. The considered memory can preserve information on one of six different symbols. The time evolution of the nodes was individually controlled by the inflow of ATP. We demonstrate that information can be written with a simple and short perturbation of the inflow. The perturbation applies to only one or two nodes, and it is symbol specific. The memory can be erased with identical inflow perturbation applied to all nodes. The presented idea of pattern-coded memory applies to other reaction networks that allow for discrete Turing patterns. Moreover, it hints at the experimental realization of memory in a simple system with the glycolytic reaction.

Generation of reference cell lines, media, and a process platform for CHO cell biomanufacturing

Due to the favorable attributes of Chinese hamster ovary (CHO) cells for therapeutic proteins and antibodies biomanufacturing, companies generate proprietary cells with desirable phenotypes. One key attribute is the ability to stably express multi-gram per liter titers in chemically defined media. Cell, media, and feed diversity has limited community efforts to translate knowledge. Moreover, academic, and nonprofit researchers generally cannot study "industrially relevant" CHO cells due to limited public availability, and the time and knowledge required to generate such cells. To address these issues, a university-industrial consortium (Advanced Mammalian Biomanufacturing Innovation Center, AMBIC) has acquired two CHO "reference cell lines" from different lineages that express monoclonal antibodies. These reference cell lines have relevant production titers, key performance outcomes confirmed by multiple laboratories, and a detailed technology transfer protocol. In commercial media, titers over 2 g/L are reached. Fed-batch cultivation data from shake flask and scaled-down bioreactors is presented. Using productivity as the primary attribute, two academic sites aligned with tight reproducibility at each site. Further, a chemically defined media formulation was developed and evaluated in parallel to the commercial media. The goal of this work is to provide a universal, industrially relevant CHO culture platform to accelerate biomanufacturing innovation.

Kinetic-model-based pathway optimization with application to reverse glycolysis in mammalian cells

Over the last two decades, model-based metabolic pathway optimization tools have been developed for the design of microorganisms to produce desired metabolites. However, few have considered more complex cellular systems such as mammalian cells, which requires the use of nonlinear kinetic models to capture the effects of concentration changes and cross-regulatory interactions. In this study, we develop a new two-stage pathway optimization framework based on kinetic models that incorporate detailed kinetics and regulation information. In Stage 1, a set of optimization problems are solved to identify and rank the enzymes that contribute the most to achieving the metabolic objective. Stage 2 then determines the optimal enzyme interventions for specified desired numbers of enzyme adjustments. It also incorporates multi-scenario optimization, which allows the simultaneous consideration of multiple physiological conditions. We apply the proposed framework to find enzyme adjustments that enable a reverse glucose flow in cultured mammalian cells, thereby eliminating the need for glucose feed in the late culture stage and enhancing process robustness. The computational results demonstrate the efficacy of the proposed approach; it not only captures the important regulations and key enzymes for reverse glycolysis but also identifies differences and commonalities in the metabolic requirements for different carbon sources.

Attenuation of glutamine synthetase selection marker improves product titer and reduces glutamine overflow in Chinese hamster ovary cells

The glutamine synthetase (GS) expression system is commonly used to ensure stable transgene integration and amplification in CHO host lines. Transfected cell populations are typically grown in the presence of the GS inhibitor, methionine sulfoximine (MSX), to further select for increased transgene copy number. However, high levels of GS activity produce excess glutamine. We hypothesized that attenuating the GS promoter while keeping the strong IgG promoter on the GS-IgG expression vector would result in a more efficient cellular metabolic phenotype. Herein, we characterized CHO cell lines expressing GS from either an attenuated promoter or an SV40 promoter and selected with/without MSX. CHO cells with the attenuated GS promoter had higher IgG specific productivity and lower glutamine production compared to cells with SV40-driven GS expression. Selection with MSX increased both specific productivity and glutamine production, regardless of GS promoter strength. ¹³ C metabolic flux analysis (MFA) was performed to further assess metabolic differences between these cell lines. Interestingly, central carbon metabolism was unaltered by the attenuated GS promoter while the fate of glutamate and glutamine varied depending on promoter strength and selection conditions. This study highlights the ability to optimize the GS expression system to improve IgG production and reduce wasteful glutamine overflow, without significantly altering central metabolism. Additionally, a detailed supplementary analysis of two "lactate runaway" reactors provides insight into the poorly understood phenomenon of excess lactate production by some CHO cell cultures. This article is protected by copyright. All rights reserved.

Classification of infinitesimal homeostasis in four-node input-output networks

An input-output network has an input node [Formula: see text], an output node o, and regulatory nodes [Formula: see text]. Such a network is a core network if each [Formula: see text] is downstream from [Formula: see text] and upstream from o. Wang et al. (J Math Biol 82:62, 2021. https://doi.org/10.1007/s00285-021-01614-1 ) show that infinitesimal homeostasis can be classified in biochemical networks through infinitesimal homeostasis in core subnetworks. Golubitsky and Wang (J Math Biol 10:1-23, 2020) show that there are three types of 3-node core networks and three types of infinitesimal homeostasis in 3-node core networks. This paper uses the theory developed in Wang et al. (2021) to show that there are twenty types of 4-node core networks (Theorem 1.3) and seventeen types of infinitesimal homeostasis in 4-node core networks (Theorem 1.7). Biological contexts illustrate the classification theorems and show that the theory can be an aid when calculating homeostasis in specific biochemical networks.

Multisite Enzymes as a Mechanism for Bistability in Reaction Networks

Here, we focus on a common class of enzymes that have multiple substrate binding sites (multisite enzymes) and analyze their capacity to generate bistable dynamics in the reaction networks that they are embedded in. These networks include both substrate-product-substrate cycles and substrate-to-product conversion with subsequent product consumption. Using mathematical techniques, we show that the inherent binding and catalysis reactions arising from multiple substrate-enzyme complexes create a potential for bistable dynamics in such reaction networks. We construct a generic model of an enzyme with n-substrate binding sites and derive an analytical solution for the steady-state concentration of all enzyme-substrate complexes. By studying these expressions, we obtain a mechanistic understanding of bistability, derive parameter combinations that guarantee bistability, and show how changing specific enzyme kinetic parameters and enzyme levels can lead to bistability in reaction networks involving multisite enzymes. Thus, the presented findings provide a biochemical and mathematical basis for predicting and engineering bistability in multisite enzymes.

A mechanistic modeling framework reveals the key principles underlying tumor metabolism

While aerobic glycolysis, or the Warburg effect, has for a long time been considered a hallmark of tumor metabolism, recent studies have revealed a far more complex picture. Tumor cells exhibit widespread metabolic heterogeneity, not only in their presentation of the Warburg effect but also in the nutrients and the metabolic pathways they are dependent on. Moreover, tumor cells can switch between different metabolic phenotypes in response to environmental cues and therapeutic interventions. A framework to analyze the observed metabolic heterogeneity and plasticity is, however, lacking. Using a mechanistic model that includes the key metabolic pathways active in tumor cells, we show that the inhibition of phosphofructokinase by excess ATP in the cytoplasm can drive a preference for aerobic glycolysis in fast-proliferating tumor cells. The differing rates of ATP utilization by tumor cells can therefore drive heterogeneity with respect to the presentation of the Warburg effect. Building upon this idea, we couple the metabolic phenotype of tumor cells to their migratory phenotype, and show that our model predictions are in agreement with previous experiments. Next, we report that the reliance of proliferating cells on different anaplerotic pathways depends on the relative availability of glucose and glutamine, and can further drive metabolic heterogeneity. Finally, using treatment of melanoma cells with a BRAF inhibitor as an example, we show that our model can be used to predict the metabolic and gene expression changes in cancer cells in response to drug treatment. By making predictions that are far more generalizable and interpretable as compared to previous tumor metabolism modeling approaches, our framework identifies key principles that govern tumor cell metabolism, and the reported heterogeneity and plasticity. These principles could be key to targeting the metabolic vulnerabilities of cancer.

Sex differences in metabolic pathways are regulated by Pfkfb3 and Pdk4 expression in rodent muscle

Skeletal muscles display sexually dimorphic features. Biochemically, glycolysis and fatty acid β-oxidation occur preferentially in the muscles of males and females, respectively. However, the mechanisms of the selective utilization of these fuels remains elusive. Here, we obtain transcriptomes from quadriceps type IIB fibers of untreated, gonadectomized, and sex steroid-treated mice of both sexes. Analyses of the transcriptomes unveil two genes, Pfkfb3 (phosphofructokinase-2) and Pdk4 (pyruvate dehydrogenase kinase 4), that may function as switches between the two sexually dimorphic metabolic pathways. Interestingly, Pfkfb3 and Pdk4 show male-enriched and estradiol-enhanced expression, respectively. Moreover, the contribution of these genes to sexually dimorphic metabolism is demonstrated by knockdown studies with cultured type IIB muscle fibers. Considering that skeletal muscles as a whole are the largest energy-consuming organs, our results provide insights into energy metabolism in the two sexes, during the estrus cycle in women, and under pathological conditions involving skeletal muscles.

Baba et al. analyzed the transcriptomes from quadriceps type IIB fibers of untreated, gonadectomized, and sex steroid-treated mice of both sexes and identified Pfkfb3 and Pdk4 as differentially regulated genes between males and diestrus females. The authors found that Pfkfb3 and Pdk4 may act as metabolic switches, showed male-enriched and estradiol-enhanced expression, respectively and contributed to sexually dimorphic metabolism.

Dietary chromium modulates glucose homeostasis and induces oxidative stress in Pacific white shrimp (Litopenaeus vannamei)

While chromium (Cr) has been recognized as an essential nutrient for all animals, and dietary supplementation can be beneficial, it can also be toxic. The present study aimed to investigate the contrasting effects of dietary chromium in Pacific white shrimp Litopenaeus vannamei. Five experimental diets were formulated to contain Cr at levels of 0.82 (Cr0.82, unsupplemented diet), 1.01 (Cr1.01), 1.22 (Cu1.22), 1.43 (Cr1.43) and 1.63 (Cr1.63) mg/kg and were fed to shrimp for 8 weeks. Highest weight gain was recorded in shrimp fed the diet containing 1.22 mg/kg Cr. Shrimp fed the diet containing the highest level of Cr (1.63 mg/kg) showed the lowest weight gain and clear signs of oxidative stress and apoptosis as evidenced by higher levels of H₂O₂, malondialdehyde and 8-hydroxydeoxyguanosine, and expression of caspase 2, 3, 5, and lower contents of total and oxidized glutathione, and expression of Cu/Zn sod, cat, gpx, mt, bcl2. Chromium supplementation promoted glycolysis and inhibited gluconeogenesis as shown by increased activities of hexokinase, phosphofructokinase and pyruvate kinase, and reduced activity of phosphoenolpyruvate carboxykinase in shrimp fed the diet containing 1.43 mg/kg Cr. Shrimp fed the diet with 1.63 mg/kg Cr had lowest contents of crustacean hyperglycemic hormone and insulin like peptide in hemolymph. Expression of genes involved in insulin signaling pathway and glycose metabolism including insr, irs1, pik3ca, pdpk1, akt, acc1, gys, glut1, pk, hk were up-regulated, and foxO1, gsk-3β, g6pc, pepck were down-regulated in shrimp fed the diets supplemented with Cr. This study demonstrated that optimum dietary supplementation of Cr had beneficial effects on glucose homeostasis and growth, whereas excess caused oxidative damage and impaired growth. The results contribute to our understanding of the biological functions of chromium in shrimp.

Bistability in fatty-acid oxidation resulting from substrate inhibition

In this study we demonstrated through analytic considerations and numerical studies that the mitochondrial fatty-acid β-oxidation can exhibit bistable-hysteresis behavior. In an experimentally validated computational model we identified a specific region in the parameter space in which two distinct stable and one unstable steady state could be attained with different fluxes. The two stable states were referred to as low-flux (disease) and high-flux (healthy) state. By a modular kinetic approach we traced the origin and causes of the bistability back to the distributive kinetics and the conservation of CoA, in particular in the last rounds of the β-oxidation. We then extended the model to investigate various interventions that may confer health benefits by activating the pathway, including (i) activation of the last enzyme MCKAT via its endogenous regulator p46-SHC protein, (ii) addition of a thioesterase (an acyl-CoA hydrolysing enzyme) as a safety valve, and (iii) concomitant activation of a number of upstream and downstream enzymes by short-chain fatty-acids (SCFA), metabolites that are produced from nutritional fibers in the gut. A high concentration of SCFAs, thioesterase activity, and inhibition of the p46Shc protein led to a disappearance of the bistability, leaving only the high-flux state. A better understanding of the switch behavior of the mitochondrial fatty-acid oxidation process between a low- and a high-flux state may lead to dietary and pharmacological intervention in the treatment or prevention of obesity and or non-alcoholic fatty-liver disease.

Exometabolome profiling reveals activation of the carnitine buffering pathway in fed-batch cultures of CHO cells co-fed with glucose and lactic acid

Adjustments to CHO cell physiology were recently observed during implementation of a Raman spectroscopy-based glucose and lactate control strategy. To further understand how these cells, under monoclonal antibody (mAb) production conditions, utilized the extra lactic acid fed, we performed a comprehensive semi-quantitative and time-dependent analysis of the exometabolome. This study focused on the CHO cell's metabolic shift from the fifth day of culture. We compared relative levels of extracellular metabolites in the absence or presence of a 2 g/L lactic acid setpoint while glucose was kept at 4 g/L. Our hypothesis is that extra lactic acid would supply more pyruvate, favoring oxidative phosphorylation. We subsequentially uncovered several carnitine derivatives as biomarkers of the simultaneous activation of TCA anaplerotic pathways as well as a carbon-buffering pathway. CHO cells exhibited a balance between intermediates from (i) amino acid catabolism, (ii) fatty acid β-oxidation, and (iii) pyruvate from glycolysis and lactic acid; and the secretion of their intermediate carnitine derivatives. In addition, 3-hydroxy-methyl-glutaric acid (HMG) and mevalonate syntheses were found as biomarkers of alternative acyl group removal. Together, under a limited capacity to assimilate the surplus of acyl-CoA groups as well as an ability to maintain the acyl-CoA: free CoA ratio for proper and continuous functioning of the TCA cycle, CHO cells activate the carnitine-buffering system, HMG, and mevalonate pathways.

Flux-augmented bifurcation analysis in chemical reaction network systems

The dynamics of biochemical reaction networks are considered to be responsible for biological functions in living systems. Since real networks are immense and complicated, it is difficult to determine which reactions can cause a significant change of dynamical behaviors, namely, bifurcations. Also to what extent numerical results of network systems depend on the chosen kinetic rate parameters is not known. In this paper, an analytical setting that splits the information of the dynamics into the network structure and reaction kinetics is introduced. This setting possesses a factorization structure for some class of network systems which allows one to determine which subnetworks are responsible for the occurrence of a bifurcation. Subsequently, the bifurcation criteria are reformulated in a manner that allows the efficient determination of relevant reactions for bifurcations.

The role of miR-320 in glucose and lipid metabolism disorder-associated diseases

Glucose and lipids are important nutrients that provide the majority of energy for each organ to maintain homeostasis of the body. With the continuous improvement in living standards, the incidence of metabolic disorder-associated diseases, such as diabetes, hyperlipidemia, and atherosclerosis, is increasing worldwide. Among them, diabetes, which could be induced by both glucose and lipid metabolic disorders, is one of the five diseases with the highest incidence and mortality worldwide. However, the detailed molecular mechanisms underlying glucose and lipid metabolism disorders and target-organ damage are still not fully defined. MicroRNAs (miRNAs) are small, non-coding, single-stranded RNAs, which usually affect their target mRNAs in the cytoplasm by post-transcriptional regulation. Previously, we have found that miR-320 contributed to glucose and lipid metabolism via different signaling pathways. Most importantly, we identified that nuclear miR-320 mediated diabetes-induced cardiac dysfunction by activating the transcription of fatty acid metabolic genes to cause lipotoxicity in the heart. Here, we reviewed the roles of miR-320 in glucose and lipid metabolism and target-organ damage.

Emerging Roles and Therapeutic Interventions of Aerobic Glycolysis in Glioma

Glioma is the most common type of intracranial malignant tumor, with a great recurrence rate due to its infiltrative growth, treatment resistance, intra- and intertumoral genetic heterogeneity. Recently, accumulating studies have illustrated that activated aerobic glycolysis participated in various cellular and clinical activities of glioma, thus influencing the efficacy of radiotherapy and chemotherapy. However, the glycolytic process is too complicated and ambiguous to serve as a novel therapy for glioma. In this review, we generalized the implication of key enzymes, glucose transporters (GLUTs), signalings and transcription factors in the glycolytic process of glioma. In addition, we summarized therapeutic interventions via the above aspects and discussed promising clinical applications for glioma.

Stochastic Simulation of Cellular Metabolism

Increased technological methods have enabled the investigation of biology at nanoscale levels. Such systems require the use of computational methods to comprehend the complex interactions that occur. The dynamics of metabolic systems have been traditionally described utilizing differential equations without fully capturing the heterogeneity of biological systems. Stochastic modeling approaches have recently emerged with the capacity to incorporate the statistical properties of such systems. However, the processing of stochastic algorithms is a computationally intensive task with intrinsic limitations. Alternatively, the queueing theory approach, historically used in the evaluation of telecommunication networks, can significantly reduce the computational power required to generate simulated results while simultaneously reducing the expansion of errors. We present here the application of queueing theory to simulate stochastic metabolic networks with high efficiency. With the use of glycolysis as a well understood biological model, we demonstrate the power of the proposed modeling methods discussed herein. Furthermore, we describe the simulation and pharmacological inhibition of glycolysis to provide an example of modeling capabilities.

Computational Analysis of Insulin-Glucagon Signalling Network: Implications of Bistability to Metabolic Homeostasis and Disease states

Insulin and glucagon control plasma macronutrient homeostasis through their signalling network composed of multiple feedback and crosstalk interactions. To understand how these interactions contribute to metabolic homeostasis and disease states, we analysed the steady state response of metabolic regulation (catabolic or anabolic) with respect to structural and input perturbations in the integrated signalling network, for varying levels of plasma glucose. Structural perturbations revealed: the positive feedback of AKT on IRS is responsible for the bistability in anabolic zone (glucose >5.5 mmol); the positive feedback of calcium on cAMP is responsible for ensuring ultrasensitive response in catabolic zone (glucose <4.5 mmol); the crosstalk between AKT and PDE3 is responsible for efficient catabolic response under low glucose condition; the crosstalk between DAG and PKC regulates the span of anabolic bistable region with respect to plasma glucose levels. The macronutrient perturbations revealed: varying plasma amino acids and fatty acids from normal to high levels gradually shifted the bistable response towards higher glucose range, eventually making the response catabolic or unresponsive to increasing glucose levels. The analysis reveals that certain macronutrient composition may be more conducive to homeostasis than others. The network perturbations that may contribute to disease states such as diabetes, obesity and cancer are discussed.

Kinetic model optimization and its application to mitigating the Warburg effect through multiple enzyme alterations

Pathway engineering is a powerful tool in biotechnological and clinical applications. However, many phenomena cannot be rewired with a single enzyme change, and in a complex network like energy metabolism, the selection of combinations of targets to engineer is a daunting task. To facilitate this process, we have developed an optimization framework and applied it to a mechanistic kinetic model of energy metabolism. We then identified combinations of enzyme alternations that led to the elimination of the Warburg effect seen in the metabolism of cancer cells and cell lines, a phenomenon coupling rapid proliferation to lactate production. Typically, optimization approaches use integer variables to achieve the desired flux redistribution with a minimum number of altered genes. This framework uses convex penalty terms to replace these integer variables and improve computational tractability. Optimal solutions are identified which substantially reduce or eliminate lactate production while maintaining the requirements for cellular proliferation using three or more enzymes.

The IDO Metabolic Trap Hypothesis for the Etiology of ME/CFS

Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) is a debilitating noncommunicable disease brandishing an enormous worldwide disease burden with some evidence of inherited genetic risk. Absence of measurable changes in patients' standard blood work has necessitated ad hoc symptom-driven therapies and a dearth of mechanistic hypotheses regarding its etiology and possible cure. A new hypothesis, the indolamine-2,3-dioxygenase (IDO) metabolic trap, was developed and formulated as a mathematical model. The historical occurrence of ME/CFS outbreaks is a singular feature of the disease and implies that any predisposing genetic mutation must be common. A database search for common damaging mutations in human enzymes produces 208 hits, including IDO2 with four such mutations. Non-functional IDO2, combined with well-established substrate inhibition of IDO1 and kinetic asymmetry of the large neutral amino acid transporter, LAT1, yielded a mathematical model of tryptophan metabolism that displays both physiological and pathological steady-states. Escape from the pathological one requires an exogenous perturbation. This model also identifies a critical point in cytosolic tryptophan abundance beyond which descent into the pathological steady-state is inevitable. If, however, means can be discovered to return cytosolic tryptophan below the critical point, return to the normal physiological steady-state is assured. Testing this hypothesis for any cell type requires only labelled tryptophan, a means to measure cytosolic tryptophan and kynurenine, and the standard tools of tracer kinetics.

Dissecting N-Glycosylation Dynamics in Chinese Hamster Ovary Cells Fed-batch Cultures using Time Course Omics Analyses

N-linked glycosylation affects the potency, safety, immunogenicity, and pharmacokinetic clearance of several therapeutic proteins including monoclonal antibodies. A robust control strategy is needed to dial in appropriate glycosylation profile during the course of cell culture processes accurately. However, N-glycosylation dynamics remains insufficiently understood owing to the lack of integrative analyses of factors that influence the dynamics, including sugar nucleotide donors, glycosyltransferases, and glycosidases. Here, an integrative approach involving multi-dimensional omics analyses was employed to dissect the temporal dynamics of glycoforms produced during fed-batch cultures of CHO cells. Several pathways including glycolysis, tricarboxylic citric acid cycle, and nucleotide biosynthesis exhibited temporal dynamics over the cell culture period. The steps involving galactose and sialic acid addition were determined as temporal bottlenecks. Our results show that galactose, and not manganese, is able to mitigate the temporal bottleneck, despite both being known effectors of galactosylation. Furthermore, sialylation is limited by the galactosylated precursors and autoregulation of cytidine monophosphate-sialic acid biosynthesis.

•

Major glycosylated species exhibit temporal dynamics during fed-batch processes

•

Key metabolic pathways linked to N-glycosylation exhibit significant temporal dynamics

•

Dynamics in nucleotide sugar donors (NSDs) directly influences glycoform heterogeneity

•

Glycoform heterogeneity can be mitigated by supplementing NSD biosynthetic precursors

PFKP Signaling at a Glance: An Emerging Mediator of Cancer Cell Metabolism

Phosphofructokinase-1 (PFK-1), a rate-determining enzyme of glycolysis, is an allosteric enzyme that regulates the oxidation of glucose in cellular respiration. Glycolysis phosphofructokinase platelet (PFKP) is the platelet isoform and works as an important mediator of cell metabolism. Considering that PFKP is a crucial player in many steps of cancer initiation and metastasis, we reviewed the specificities and complexities of PFKP and its biological roles in human diseases, especially malignant tumors. The possible use of PFKP as a diagnostic marker or a drug target in the prevention or treatment of cancer is also discussed.

A Mathematical Model for Enzyme Clustering in Glucose Metabolism

We have recently demonstrated that the rate-limiting enzymes in human glucose metabolism organize into cytoplasmic clusters to form a multienzyme complex, the glucosome, in at least three different sizes. Quantitative high-content imaging data support a hypothesis that the glucosome clusters regulate the direction of glucose flux between energy metabolism and building block biosynthesis in a cluster size-dependent manner. However, direct measurement of their functional contributions to cellular metabolism at subcellular levels has remained challenging. In this work, we develop a mathematical model using a system of ordinary differential equations, in which the association of the rate-limiting enzymes into multienzyme complexes is included as an essential element. We then demonstrate that our mathematical model provides a quantitative principle to simulate glucose flux at both subcellular and population levels in human cancer cells. Lastly, we use the model to simulate 2-deoxyglucose-mediated alteration of glucose flux in a population level based on subcellular high-content imaging data. Collectively, we introduce a new mathematical model for human glucose metabolism, which promotes our understanding of functional roles of differently sized multienzyme complexes in both single-cell and population levels.

New concepts in feedback regulation of glucose metabolism

Glycolysis, the breakdown of glucose, is one of the most conserved and extensively studied biochemical pathways. Designing principles from chemistry and thermodynamics allow for energy production, biosynthesis and cellular communication. However, the kinetics or metabolic flux through the pathway also determines its function. Recently, there have been numerous developments that establish new allosteric interactions of glycolytic enzymes with small molecule metabolites and other mechanisms that may cooperate to allow for addition complex regulation of glycolysis. This review surveys these newfound sources of glycolysis regulation and discusses their possible roles in establishing kinetic design principles of glycolysis.

RNA SEQ Analysis Indicates that the AE3 Cl-/HCO3- Exchanger Contributes to Active Transport-Mediated CO2 Disposal in Heart

Loss of the AE3 Cl⁻/HCO₃⁻ exchanger (Slc4a3) in mice causes an impaired cardiac force-frequency response and heart failure under some conditions but the mechanisms are not known. To better understand the functions of AE3, we performed RNA Seq analysis of AE3-null and wild-type mouse hearts and evaluated the data with respect to three hypotheses (CO₂ disposal, facilitation of Na⁺-loading, and recovery from an alkaline load) that have been proposed for its physiological functions. Gene Ontology and PubMatrix analyses of differentially expressed genes revealed a hypoxia response and changes in vasodilation and angiogenesis genes that strongly support the CO₂ disposal hypothesis. Differential expression of energy metabolism genes, which indicated increased glucose utilization and decreased fatty acid utilization, were consistent with adaptive responses to perturbations of O₂/CO₂ balance in AE3-null myocytes. Given that the myocardium is an obligate aerobic tissue and consumes large amounts of O₂, the data suggest that loss of AE3, which has the potential to extrude CO₂ in the form of HCO₃⁻, impairs O₂/CO₂ balance in cardiac myocytes. These results support a model in which the AE3 Cl⁻/HCO₃⁻ exchanger, coupled with parallel Cl⁻ and H⁺-extrusion mechanisms and extracellular carbonic anhydrase, is responsible for active transport-mediated disposal of CO₂.

Low glucose concentrations within typical industrial operating conditions have minimal effect on the transcriptome of recombinant CHO cells

Typically, mammalian cell culture medium contains high glucose concentrations that are analogous to diabetic levels in humans, suggesting that mammalian cells are cultivated in excessive glucose. Using RNA-Seq, this study characterized the Chinese hamster ovary (CHO) cell transcriptome under two glucose concentrations to assess the genetic effects associated with metabolic pathways, in addition to other global responses. The initial extracellular glucose concentrations used represented high (30 mM) and low (10 mM) glucose conditions, where at the time the transcriptomes were compared, the glucose concentrations were approximately 24 and 4.4 mM for the mid-exponential cultures, where 4.4 mM represents a common target concentration in the biopharmaceutical industry for controlled fed-batch cultures. A recombinant CHO cell line producing a monoclonal antibody was used, such that the impact on glycosylation genes could be evaluated. Relatively few genes were identified as being significantly different (FDR ≤ 0.01) between the high and low glucose conditions, for example, only 575 genes, and only 40 of these genes had 2-fold or greater differences. Gene expression differences for glycolysis, TCA cycle, and glycosylation-related reactions were minimal and unlikely to have biological significance. This transcriptome study indicates that low glucose concentrations in the culture medium are unlikely to cause any biologically significant or detrimental changes to CHO cells at the transcriptome level. Furthermore, it is well-known that maintaining low glucose concentrations in fed-batch cultures can reduce lactate production, which in turn improves process outcomes. Taken together, the transcriptome data supports the continued development of low glucose-based processes to control lactate. © 2017 American Institute of Chemical Engineers Biotechnol. Prog., 33:771-785, 2017.

The Eukaryote-Like Serine/Threonine Kinase STK Regulates the Growth and Metabolism of Zoonotic Streptococcus suis

Like eukaryotes, bacteria express one or more serine/threonine kinases (STKs) that initiate diverse signaling networks. The STK from Streptococcus suis is encoded by a single-copy stk gene, which is crucial in stress response and virulence. To further understand the regulatory mechanism of STK in S. suis, a stk deletion strain (Δstk) and its complementary strain (CΔstk) were constructed to systematically decode STK characteristics by applying whole transcriptome RNA sequencing (RNA-Seq) and phosphoproteomic analysis. Numerous genes were differentially expressed in Δstk compared with the wild-type parental strain SC-19, including 320 up-regulated and 219 down-regulated genes. Particularly, 32 virulence-associated genes (VAGs) were significantly down-regulated in Δstk. Seven metabolic pathways relevant to bacterial central metabolism and translation are significantly repressed in Δstk. Phosphoproteomic analysis further identified 12 phosphoproteins that exhibit differential phosphorylation in Δstk. These proteins are associated with cell growth and division, glycolysis, and translation. Consistently, phenotypic assays confirmed that the Δstk strain displayed deficient growth and attenuated pathogenicity. Thus, STK is a central regulator that plays an important role in cell growth and division, as well as S. suis metabolism.

Metabolic Control in Mammalian Fed-Batch Cell Cultures for Reduced Lactic Acid Accumulation and Improved Process Robustness

Biomass and cell-specific metabolic rates usually change dynamically over time, making the "feed according to need" strategy difficult to realize in a commercial fed-batch process. We here demonstrate a novel feeding strategy which is designed to hold a particular metabolic state in a fed-batch process by adaptive feeding in real time. The feed rate is calculated with a transferable biomass model based on capacitance, which changes the nutrient flow stoichiometrically in real time. A limited glucose environment was used to confine the cell in a particular metabolic state. In order to cope with uncertainty, two strategies were tested to change the adaptive feed rate and prevent starvation while in limitation: (i) inline pH and online glucose concentration measurement or (ii) inline pH alone, which was shown to be sufficient for the problem statement. In this contribution, we achieved metabolic control within a defined target range. The direct benefit was two-fold: the lactic acid profile was improved and pH could be kept stable. Multivariate Data Analysis (MVDA) has shown that pH influenced lactic acid production or consumption in historical data sets. We demonstrate that a low pH (around 6.8) is not required for our strategy, as glucose availability is already limiting the flux. On the contrary, we boosted glycolytic flux in glucose limitation by setting the pH to 7.4. This new approach led to a yield of lactic acid/glucose (Y L/G) around zero for the whole process time and high titers in our labs. We hypothesize that a higher carbon flux, resulting from a higher pH, may lead to more cells which produce more product. The relevance of this work aims at feeding mammalian cell cultures safely in limitation with a desired metabolic flux range. This resulted in extremely stable, low glucose levels, very robust pH profiles without acid/base interventions and a metabolic state in which lactic acid was consumed instead of being produced from day 1. With this contribution, we wish to extend the basic repertoire of available process control strategies, which will open up new avenues in automation technology and radically improve process robustness in both process development and manufacturing.

Differential Proteomic Analysis Using iTRAQ Reveals Alterations in Hull Development in Rice (Oryza sativa L.)

Rice hull, the outer cover of the rice grain, determines grain shape and size. Changes in the rice hull proteome in different growth stages may reflect the underlying mechanisms involved in grain development. To better understand these changes, isobaric tags for relative and absolute quantitative (iTRAQ) MS/MS was used to detect statistically significant changes in the rice hull proteome in the booting, flowering, and milk-ripe growth stages. Differentially expressed proteins were analyzed to predict their potential functions during development. Gene ontology (GO) terms and pathways were used to evaluate the biological mechanisms involved in rice hull at the three growth stages. In total, 5,268 proteins were detected and characterized, of which 563 were differentially expressed across the development stages. The results showed that the flowering and milk-ripe stage proteomes were more similar to each other (r=0.61) than either was to the booting stage proteome. A GO enrichment analysis of the differentially expressed proteins was used to predict their roles during rice hull development. The potential functions of 25 significantly differentially expressed proteins were used to evaluate their possible roles at various growth stages. Among these proteins, an unannotated protein (Q7X8A1) was found to be overexpressed especially in the flowering stage, while a putative uncharacterized protein (B8BF94) and an aldehyde dehydrogenase (Q9FPK6) were overexpressed only in the milk-ripe stage. Pathways regulated by differentially expressed proteins were also analyzed. Magnesium-protoporphyrin IX monomethyl ester [oxidative] cyclase (Q9SDJ2), and two magnesium-chelatase subunits, ChlD (Q6ATS0), and ChlI (Q53RM0), were associated with chlorophyll biosynthesis at different developmental stages. The expression of Q9SDJ2 in the flowering and milk-ripe stages was validated by qRT-PCR. The 25 candidate proteins may be pivotal markers for controlling rice hull development at various growth stages and chlorophyll biosynthesis pathway related proteins, especially magnesium-protoporphyrin IX monomethyl ester [oxidative] cyclase (Q9SDJ2), may provide new insights into the molecular mechanisms of rice hull development and chlorophyll associated regulation.

Multiplicity of steady states in glycolysis and shift of metabolic state in cultured mammalian cells

Cultured mammalian cells exhibit elevated glycolysis flux and high lactate production. In the industrial bioprocesses for biotherapeutic protein production, glucose is supplemented to the culture medium to sustain continued cell growth resulting in the accumulation of lactate to high levels. In such fed-batch cultures, sometimes a metabolic shift from a state of high glycolysis flux and high lactate production to a state of low glycolysis flux and low lactate production or even lactate consumption is observed. While in other cases with very similar culture conditions, the same cell line and medium, cells continue to produce lactate. A metabolic shift to lactate consumption has been correlated to the productivity of the process. Cultures that exhibited the metabolic shift to lactate consumption had higher titers than those which didn't. However, the cues that trigger the metabolic shift to lactate consumption state (or low lactate production state) are yet to be identified. Metabolic control of cells is tightly linked to growth control through signaling pathways such as the AKT pathway. We have previously shown that the glycolysis of proliferating cells can exhibit bistability with well-segregated high flux and low flux states. Low lactate production (or lactate consumption) is possible only at a low glycolysis flux state. In this study, we use mathematical modeling to demonstrate that lactate inhibition together with AKT regulation on glycolysis enzymes can profoundly influence the bistable behavior, resulting in a complex steady-state topology. The transition from the high flux state to the low flux state can only occur in certain regions of the steady state topology, and therefore the metabolic fate of the cells depends on their metabolic trajectory encountering the region that allows such a metabolic state switch. Insights from such switch behavior present us with new means to control the metabolism of mammalian cells in fed-batch cultures.

Global insights into the Chinese hamster and CHO cell transcriptomes

Transcriptomics is increasingly being used on Chinese hamster ovary (CHO) cells to unveil physiological insights related to their performance during production processes. The rich transcriptome data can be exploited to provide impetus for systems investigation such as modeling the central carbon metabolism or glycosylation pathways, or even building genome-scale models. To harness the power of transcriptome assays, we assembled and annotated a set of RNA-Seq data from multiple CHO cell lines and Chinese hamster tissues, and constructed a DNA microarray. The identity of genes involved in major functional pathways and their transcript levels generated in this study will serve as a reference for future studies employing kinetic models. In particular, the data on glycolysis and glycosylation pathways indicate that the variability of gene expression level among different cell lines and tissues may contribute to their differences in metabolism and glycosylation patterns. Thereby, these insights can potentially lead to opportunities for cell engineering. This repertoire of transcriptome data also enables the identification of potential sequence variants in cell lines and allows tracing of cell lineages. Overall the study is an illustration of the potential benefit of RNA-Seq that is yet to be exploited.

Mechanism for multiplicity of steady states with distinct cell concentration in continuous culture of mammalian cells

Continuous culture for the production of biopharmaceutical proteins offers the possibility of steady state operations and thus more consistent product quality and increased productivity. Under some conditions, multiplicity of steady states has been observed in continuous cultures of mammalian cells, wherein with the same dilution rate and feed nutrient composition, steady states with very different cell and product concentrations may be reached. At those different steady states, cells may exhibit a high glycolysis flux with high lactate production and low cell concentration, or a low glycolysis flux with low lactate and high cell concentration. These different steady states, with different cell concentration, also have different productivity. Developing a mechanistic understanding of the occurrence of steady state multiplicity and devising a strategy to steer the culture toward the desired steady state is critical. We establish a multi-scale kinetic model that integrates a mechanistic intracellular metabolic model and cell growth model in a continuous bioreactor. We show that steady state multiplicity exists in a range of dilution rate in continuous culture as a result of the bistable behavior in glycolysis. The insights from the model were used to devise strategies to guide the culture to the desired steady state in the multiple steady state region. The model provides a guideline principle in the design of continuous culture processes of mammalian cells.

Heterogeneity of glycolysis in cancers and therapeutic opportunities

Upregulated glycolysis, both in normoxic and hypoxic environments, is a nearly universal trait of cancer cells. The enormous difference in glucose metabolism offers a target for therapeutic intervention with a potentially low toxicity profile. The past decade has seen a steep rise in the development and clinical assessment of small molecules that target glycolysis. The enzymes in glycolysis have a highly heterogeneous nature that allows for the different bioenergetic, biosynthetic, and signaling demands needed for various tissue functions. In cancers, these properties enable them to respond to the variable requirements of cell survival, proliferation and adaptation to nutrient availability. Heterogeneity in glycolysis occurs through the expression of different isoforms, posttranslational modifications that affect the kinetic and regulatory properties of the enzyme. In this review, we will explore this vast heterogeneity of glycolysis and discuss how this information might be exploited to better target glucose metabolism and offer possibilities for biomarker development.