The increasing availability of metabolomics data necessitates novel methods for deeper data analysis and interpretation. We present a flux balance analysis method that allows for the computation of dynamic intracellular metabolic changes at the cellular scale through integration of time-course absolute quantitative metabolomics. This approach, termed "unsteady-state flux balance analysis" (uFBA), is applied to four cellular systems: three dynamic and one steady-state as a negative control. uFBA and FBA predictions are contrasted, and uFBA is found to be more accurate in predicting dynamic metabolic flux states for red blood cells, platelets, and Saccharomyces cerevisiae. Notably, only uFBA predicts that stored red blood cells metabolize TCA intermediates to regenerate important cofactors, such as ATP, NADH, and NADPH. These pathway usage predictions were subsequently validated through ¹³C isotopic labeling and metabolic flux analysis in stored red blood cells. Utilizing time-course metabolomics data, uFBA provides an accurate method to predict metabolic physiology at the cellular scale for dynamic systems.

Background

Recent years have witnessed a rising trend in exploring microalgae for valuable carotenoid products as the demand for lutein and many other carotenoids in global markets has increased significantly. In green microalgae lutein is a major carotenoid protecting cellular components from damage incurred by reactive oxygen species under stress conditions. In this study, we investigated the effects of abiotic stressors on lutein accumulation in a strain of the marine microalga D. salina which had been selected for growth under stress conditions of combined blue and red lights by adaptive laboratory evolution.

Results

Nitrate concentration, salinity and light quality were selected as three representative influencing factors and their impact on lutein production in batch cultures of D. salina was evaluated using response surface analysis. D. salina was found to be more tolerant to hyper-osmotic stress than to hypo-osmotic stress which caused serious cell damage and death in a high proportion of cells while hyper-osmotic stress increased the average cell size of D. salina only slightly. Two models were developed to explain how lutein productivity depends on the stress factors and for predicting the optimal conditions for lutein productivity. Among the three stress variables for lutein production, stronger interactions were found between nitrate concentration and salinity than between light quality and the other two. The predicted optimal conditions for lutein production were close to the original conditions used for adaptive evolution of D. salina. This suggests that the conditions imposed during adaptive evolution may have selected for the growth optima arrived at.

Conclusions

This study shows that systematic evaluation of the relationship between abiotic environmental stresses and lutein biosynthesis can help to decipher the key parameters in obtaining high levels of lutein productivity in D. salina. This study may benefit future stress-driven adaptive laboratory evolution experiments and a strategy of applying stress in a step-wise manner can be suggested for a rational design of experiments.

The metabolism of gluconate is well characterized in prokaryotes where it is known to be degraded following phosphorylation by gluconokinase. Less is known of gluconate metabolism in humans. Human gluconokinase activity was recently identified proposing questions about the metabolic role of gluconate in humans. Here we report the recombinant expression, purification and biochemical characterization of isoform I of human gluconokinase alongside substrate specificity and kinetic assays of the enzyme catalyzed reaction. The enzyme, shown to be a dimer, had ATP dependent phosphorylation activity and strict specificity towards gluconate out of 122 substrates tested. In order to evaluate the metabolic impact of gluconate in humans we modeled gluconate metabolism using steady state metabolic network analysis. The results indicate that significant metabolic flux changes in anabolic pathways linked to the hexose monophosphate shunt (HMS) are induced through a small increase in gluconate concentration. We argue that the enzyme takes part in a context specific carbon flux route into the HMS that, in humans, remains incompletely explored. Apart from the biochemical description of human gluconokinase, the results highlight that little is known of the mechanism of gluconate metabolism in humans despite its widespread use in medicine and consumer products.

Metabolomics is a rapidly evolving analytical approach in life and health sciences. The structural elucidation of the metabolites of interest remains a major analytical challenge in the metabolomics workflow. Here, we investigate the use of ion mobility as a tool to aid metabolite identification. Ion mobility allows for the measurement of the rotationally averaged collision cross-section (CCS), which gives information about the ionic shape of a molecule in the gas phase. We measured the CCSs of 125 common metabolites using traveling-wave ion mobility-mass spectrometry (TW-IM-MS). CCS measurements were highly reproducible on instruments located in three independent laboratories (RSD < 5% for 99%). We also determined the reproducibility of CCS measurements in various biological matrixes including urine, plasma, platelets, and red blood cells using ultra performance liquid chromatography (UPLC) coupled with TW-IM-MS. The mean RSD was < 2% for 97% of the CCS values, compared to 80% of retention times. Finally, as proof of concept, we used UPLC-TW-IM-MS to compare the cellular metabolome of epithelial and mesenchymal cells, an in vitro model used to study cancer development. Experimentally determined and computationally derived CCS values were used as orthogonal analytical parameters in combination with retention time and accurate mass information to confirm the identity of key metabolites potentially involved in cancer. Thus, our results indicate that adding CCS data to searchable databases and to routine metabolomics workflows will increase the identification confidence compared to traditional analytical approaches.

Metabolic models can provide a mechanistic framework to analyze information-rich omics data sets, and are increasingly being used to investigate metabolic alternations in human diseases. An expression of the altered metabolic pathway utilization is the selection of metabolites consumed and released by cells. However, methods for the inference of intracellular metabolic states from extracellular measurements in the context of metabolic models remain underdeveloped compared to methods for other omics data. Herein, we describe a workflow for such an integrative analysis emphasizing on extracellular metabolomics data. We demonstrate, using the lymphoblastic leukemia cell lines Molt-4 and CCRF-CEM, how our methods can reveal differences in cell metabolism. Our models explain metabolite uptake and secretion by predicting a more glycolytic phenotype for the CCRF-CEM model and a more oxidative phenotype for the Molt-4 model, which was supported by our experimental data. Gene expression analysis revealed altered expression of gene products at key regulatory steps in those central metabolic pathways, and literature query emphasized the role of these genes in cancer metabolism. Moreover, in silico gene knock-outs identified unique control points for each cell line model, e.g., phosphoglycerate dehydrogenase for the Molt-4 model. Thus, our workflow is well-suited to the characterization of cellular metabolic traits based on extracellular metabolomic data, and it allows the integration of multiple omics data sets into a cohesive picture based on a defined model context.

Despite recent advances in analytical and computational chemistry, lipid identification remains a significant challenge in lipidomics. Ion-mobility spectrometry provides an accurate measure of the molecules' rotationally averaged collision cross-section (CCS) in the gas phase and is thus related to ionic shape. Here, we investigate the use of CCS as a highly specific molecular descriptor for identifying lipids in biological samples. Using traveling wave ion mobility mass spectrometry (MS), we measured the CCS values of over 200 lipids within multiple chemical classes. CCS values derived from ion mobility were not affected by instrument settings or chromatographic conditions, and they were highly reproducible on instruments located in independent laboratories (interlaboratory RSD < 3% for 98% of molecules). CCS values were used as additional molecular descriptors to identify brain lipids using a variety of traditional lipidomic approaches. The addition of CCS improved the reproducibility of analysis in a liquid chromatography-MS workflow and maximized the separation of isobaric species and the signal-to-noise ratio in direct-MS analyses (e.g., "shotgun" lipidomics and MS imaging). These results indicate that adding CCS to databases and lipidomics workflows increases the specificity and selectivity of analysis, thus improving the confidence in lipid identification compared to traditional analytical approaches. The CCS/accurate-mass database described here is made publicly available.

The temperature dependence of biological processes has been studied at the levels of individual biochemical reactions and organism physiology (e.g. basal metabolic rates) but has not been examined at the metabolic network level. Here, we used a systems biology approach to characterize the temperature dependence of the human red blood cell (RBC) metabolic network between 4 and 37 °C through absolutely quantified exo- and endometabolomics data. We used an Arrhenius-type model (Q₁₀) to describe how the rate of a biochemical process changes with every 10 °C change in temperature. Multivariate statistical analysis of the metabolomics data revealed that the same metabolic network-level trends previously reported for RBCs at 4 °C were conserved but accelerated with increasing temperature. We calculated a median Q₁₀ coefficient of 2.89 ± 1.03, within the expected range of 2-3 for biological processes, for 48 individual metabolite concentrations. We then integrated these metabolomics measurements into a cell-scale metabolic model to study pathway usage, calculating a median Q₁₀ coefficient of 2.73 ± 0.75 for 35 reaction fluxes. The relative fluxes through glycolysis and nucleotide metabolism pathways were consistent across the studied temperature range despite the non-uniform distributions of Q₁₀ coefficients of individual metabolites and reaction fluxes. Together, these results indicate that the rate of change of network-level responses to temperature differences in RBC metabolism is consistent between 4 and 37 °C. More broadly, we provide a baseline characterization of a biochemical network given no transcriptional or translational regulation that can be used to explore the temperature dependence of metabolism.

Chinese hamster ovary (CHO) cells dominate biotherapeutic protein production and are widely used in mammalian cell line engineering research. To elucidate metabolic bottlenecks in protein production and to guide cell engineering and bioprocess optimization, we reconstructed the metabolic pathways in CHO and associated them with >1,700 genes in the Cricetulus griseus genome. The genome-scale metabolic model based on this reconstruction, iCHO1766, and cell-line-specific models for CHO-K1, CHO-S, and CHO-DG44 cells provide the biochemical basis of growth and recombinant protein production. The models accurately predict growth phenotypes and known auxotrophies in CHO cells. With the models, we quantify the protein synthesis capacity of CHO cells and demonstrate that common bioprocess treatments, such as histone deacetylase inhibitors, inefficiently increase product yield. However, our simulations show that the metabolic resources in CHO are more than three times more efficiently utilized for growth or recombinant protein synthesis following targeted efforts to engineer the CHO secretory pathway. This model will further accelerate CHO cell engineering and help optimize bioprocesses.