Flux Balance Analysis Research Articles

‘Social’ versus ‘asocial’ cells—dynamic competition flux balance analysis

In multicellular organisms cells compete for resources or growth factors. If any one cell type wins, the co-existence of diverse cell types disappears. Existing dynamic Flux Balance Analysis (dFBA) does not accommodate changes in cell density caused by competition. Therefore we here develop ‘dynamic competition Flux Balance Analysis’ (dcFBA). With total biomass synthesis as objective, lower-growth-yield cells were outcompeted even when cells synthesized mutually required nutrients. Signal transduction between cells established co-existence, which suggests that such ‘socialness’ is required for multicellularity. Whilst mutants with increased specific growth rate did not outgrow the other cell types, loss of social characteristics did enable a mutant to outgrow the other cells. We discuss that ‘asocialness’ rather than enhanced growth rates, i.e., a reduced sensitivity to regulatory factors rather than enhanced growth rates, may characterize cancer cells and organisms causing ecological blooms. Therapies reinforcing cross-regulation may therefore be more effective than those targeting replication rates.

A role for fermentation in aerobic conditions as revealed by computational analysis of maize root metabolism during growth by cell elongation.

The root is a well-studied example of cell specialisation, yet little is known about the metabolism that supports the transport functions and growth of different root cell types. To address this, we used computational modelling to study metabolism in the elongation zone of a maize lateral root. A functional-structural model captured the cell-anatomical features of the root and modelled how they changed as the root elongated. From these data, we derived constraints for a flux balance analysis model that predicted metabolic fluxes of the 11 concentric rings of cells in the root. We discovered a distinct metabolic flux pattern in the cortical cell rings, endodermis and pericycle (but absent in the epidermis) that involved a high rate of glycolysis and production of the fermentation end-products lactate and ethanol. This aerobic fermentation was confirmed experimentally by metabolite analysis. The use of fermentation in the model was not obligatory but was the most efficient way to meet the specific demands for energy, reducing power and carbon skeletons of expanding cells. Cytosolic acidification was avoided in the fermentative mode due to the substantial consumption of protons by lipid synthesis. These results expand our understanding of fermentative metabolism beyond that of hypoxic niches and suggest that fermentation could play an important role in the metabolism of aerobic tissues.

Model predictive control of a fermenter using dynamic flux balance analysis coupled with convolutional neural networks

Flux balance analysis-based models are increasingly used in bioprocess control and optimization. Unlike unstructured models, flux balance analysis-based models investigate the genome-scale network reconstructions of the microorganisms under study. Although these models are accurate, they pose computational cost challenges in rigorous optimization tasks or online optimal control strategies. In this work, we develop low-computational cost hybrid models by using deep convolutional neural networks to surrogate the problem/model. Furthermore, to address the computational challenges associated with optimal control of flux balance analysis-based models, we propose a successive linearization scheme that incorporates a Laguerre function-based model predictive control strategy coupled with a Luenberger-like observer. To investigate the effectiveness of the proposed method, optimal control of a fed-batch process is considered. Results show the acceptable accuracy of the proposed hybrid model and control scheme while reducing the computational cost significantly.

Balancing biomass reaction stoichiometry and measured fluxes in flux balance analysis.

Flux balance analysis (FBA) is widely recognized as an important method for studying metabolic networks. When incorporating flux measurements of certain reactions into an FBA problem, it is possible that the underlying linear program may become infeasible, e.g. due to measurement or modeling inaccuracies. Furthermore, while the biomass reaction is of central importance in FBA models, its stoichiometry is often a rough estimate and a source of high uncertainty. In this work, we present a method that allows modifications to the biomass reaction stoichiometry as a means to (i) render the FBA problem feasible and (ii) improve the accuracy of the model by corrections in the biomass composition. Optionally, the adjustment of the biomass composition can be used in conjunction with a previously introduced approach for balancing inconsistent fluxes to obtain a feasible FBA system. We demonstrate the value of our approach by analyzing realistic flux measurements of E.coli. In particular, we find that the growth-associated maintenance (GAM) demand of ATP, which is typically integrated with the biomass reaction, is likely overestimated in recent genome-scale models, at least for certain growth conditions. In light of these findings, we discuss issues related to the determination and inclusion of GAM values in constraint-based models. Overall, our method can uncover potential errors and suggest adjustments in the assumed biomass composition in FBA models based on inconsistencies between the model and measured fluxes. The developed method has been implemented in our software tool CNApy available from https://github.com/cnapy-org/CNApy.

SL-scan identifies synthetic lethal interactions in cancer using metabolic networks

Exploiting synthetic lethality is a promising strategy for developing targeted cancer therapies. However, identifying clinically significant synthetic lethal (SL) interactions among a large number of gene combinations is a challenging computational task. In this study, we developed the SL-scan pipeline based on metabolic network modeling to discover SL interaction. The SL-scan pipeline identifies the association between simulated Flux Balance Analysis knockout scores and mutation data across cancer cell lines and predicts putative SL interactions. We assessed the concordance of the SL pairs predicted by SL-scan with those of obtained from analysis of the CRISPR, shRNA, and PRISM datasets. Our results demonstrate that the SL-scan pipeline outperformed existing SL prediction approaches based on metabolic networks in identifying SL pairs in various cancers. This study emphasizes the importance of integrating multiple data sources, particularly mutation data, when identifying SL pairs for targeted cancer therapies. The findings of this study may lead to the development of novel targeted cancer therapies.

The rise of the sourdough: Genome-scale metabolic modeling-based approach to design sourdough starter communities with tailored-made properties

Sourdough starters harbor microbial consortia that benefit the final product's aroma and volume. The complex nature of these spontaneously developed communities raises challenges in predicting the fermentation phenotypes. Herein, we demonstrated for the first time in this field the potential of genome-scale metabolic modeling (GEMs) in the study of sourdough microbial communities. Broad in-silico modeling of microbial growth was applied on communities composed of yeast (Saccharomyces cerevisiae) and different Lactic Acid Bacteria (LAB) species, which mainly predominate in sourdough starters. Simulations of model-represented communities associated specific bacterial compositions with sourdough phenotypes. Based on ranking the phenotypic performances of different combinations, Pediococcus spp. – Lb. sakei group members were predicted to have an optimal effect considering the increase in S. cerevisiae growth abilities and overall CO2 secretion rates. Flux Balance Analysis (FBA) revealed mutual relationships between the Pediococcus spp. – Lb. sakei group members and S. cerevisiae through bidirectional nutrient dependencies, and further underlined that these bacteria compete with the yeast over nutrients to a lesser extent than the rest LAB species. Volatile compounds (VOCs) production was further modeled, identifying species-specific and community-related VOCs production profiles. The in-silico models' predictions were validated by experimentally building synthetic sourdough communities and assessing the fermentation phenotypes. The Pediococcus spp. – Lb. sakei group was indeed associated with increased yeast cell counts and fermentation rates, demonstrating a 25 % increase in the average leavening rates during the first 10 fermentation hours compared to communities with a lower representation of these group members. Overall, these results provide a possible novel strategy towards the de-novo design of sourdough starter communities with tailored-made characterizations, including a shortened leavening period.

Maximizing multi-reaction dependencies provides more accurate and precise predictions of intracellular fluxes than the principle of parsimony.

Intracellular fluxes represent a joint outcome of cellular transcription and translation and reflect the availability and usage of nutrients from the environment. While approaches from the constraint-based metabolic framework can accurately predict cellular phenotypes, such as growth and exchange rates with the environment, accurate prediction of intracellular fluxes remains a pressing problem. Parsimonious flux balance analysis (pFBA) has become an approach of choice to predict intracellular fluxes by employing the principle of efficient usage of protein resources. Nevertheless, comparative analyses of intracellular flux predictions from pFBA against fluxes estimated from labeling experiments remain scarce. Here, we posited that steady-state flux distributions derived from the principle of maximizing multi-reaction dependencies are of improved accuracy and precision than those resulting from pFBA. To this end, we designed a constraint-based approach, termed complex-balanced FBA (cbFBA), to predict steady-state flux distributions that support the given specific growth rate and exchange fluxes. We showed that the steady-state flux distributions resulting from cbFBA in comparison to pFBA show better agreement with experimentally measured fluxes from 17 Escherichia coli strains and are more precise, due to the smaller space of alternative solutions. We also showed that the same principle holds in eukaryotes by comparing the predictions of pFBA and cbFBA against experimentally derived steady-state flux distributions from 26 knock-out mutants of Saccharomyces cerevisiae. Furthermore, our results showed that intracellular fluxes predicted by cbFBA provide better support for the principle of minimizing metabolic adjustment between mutants and wild types. Together, our findings point that other principles that consider the dynamics and coordination of steady states may govern the distribution of intracellular fluxes.

Unveiling abundance-dependent metabolic phenotypes of microbial communities.

Constraint-based modeling has risen as an alternative for characterizing metabolism of communities. Adaptations of flux balance analysis have been proposed to model metabolic interactions, which in most cases consider the maximization of biomass production as their objective. In nature, novel essential functions are not directly related to cell growth force communities to display suboptimal growth rates. These suboptimal states allow a degree of plasticity in their metabolism, thus allowing quick shifts between alternative flux distributions as an initial response to environmental changes. In this work, we introduce the abundance-growth space as a representation of metabolic phenotypes of a community. This space is defined by the composition of a community, represented by its members' relative abundances, and their growth rate. The analysis of this space allows us to pinpoint how critical reactions respond to shifts of the environment, showing where changes in community plasticity occur. Interestingly, it highlights the relevance of the relative abundance of its members in the lost or gain of plasticity. This method is applied to two simple communities that exchange metabolites. A synthetic community of two mutant Escherichia coli strains and an environmental bioleaching community composed by Acidithiobacillus ferrooxidans Wenelen and Sulfobacillus thermosulfidooxidans Cutipay, where only Cutipay consumes organic matter disposed of by the community.IMPORTANCEIn nature, organisms live in communities and not as isolated species, and their interactions provide a source of resilience to environmental disturbances. Despite their importance in ecology, human health, and industry, understanding how organisms interact in different environments remains an open question.In this work, we provide a novel approach that, only using genomic information, studies the metabolic phenotype exhibited by communities, where the exploration of suboptimal growth flux distributions and the composition of a community allows to unveil its capacity to respond to environmental changes, shedding light of the degrees of metabolic plasticity inherent to the community.

Cardinality optimization in constraint-based modelling: application to human metabolism.

Several applications in constraint-based modelling can be mathematically formulated as cardinality optimization problems involving the minimization or maximization of the number of nonzeros in a vector. These problems include testing for stoichiometric consistency, testing for flux consistency, testing for thermodynamic flux consistency, computing sparse solutions to flux balance analysis problems and computing the minimum number of constraints to relax to render an infeasible flux balance analysis problem feasible. Such cardinality optimization problems are computationally complex, with no known polynomial time algorithms capable of returning an exact and globally optimal solution. By approximating the zero-norm with nonconvex continuous functions, we reformulate a set of cardinality optimization problems in constraint-based modelling into a difference of convex functions. We implemented and numerically tested novel algorithms that approximately solve the reformulated problems using a sequence of convex programs. We applied these algorithms to various biochemical networks and demonstrate that our algorithms match or outperform existing related approaches. In particular, we illustrate the efficiency and practical utility of our algorithms for cardinality optimization problems that arise when extracting a model ready for thermodynamic flux balance analysis given a human metabolic reconstruction. Open source scripts to reproduce the results are here https://github.com/opencobra/COBRA.papers/2023_cardOpt with general purpose functions integrated within the COnstraint-Based Reconstruction and Analysis toolbox: https://github.com/opencobra/cobratoolbox.

Can Genome Sequencing Coupled to Flux Balance Analyses Offer Precision Guidance for Industrial Strain Development? The Lessons from Carbon Trafficking in Corynebacterium glutamicum ATCC 21573.

Systems biology tools offer new prospects for industrial strain selection. For bacteria that are significant for industrial applications, whole-genome sequencing coupled to flux balance analysis (FBA) can help unpack the complex relationships between genome mutations and carbon trafficking. This work investigates the l-tyrosine (l-Tyr) overproducing model system Corynebacterium glutamicum ATCC 21573 with an eye to more rational and precision strain development. Using genome-wide mutational analysis of C. glutamicum, we identified 27,611 single nucleotide polymorphisms and 479 insertion/deletion mutations. Mutations in the carbon uptake machinery have led to phosphotransferase system-independent routes as corroborated with FBA. Mutations within the central carbon metabolism of C. glutamicum impaired the carbon flux, as evidenced by the lower growth rate. The entry to and flow through the tricarboxylic acid cycle was affected by mutations in pyruvate and α-ketoglutarate dehydrogenase complexes, citrate synthase, and isocitrate dehydrogenase. FBA indicated that the estimated flux through the shikimate pathway became larger as the l-Tyr production rate increased. In addition, protocatechuate export was probabilistically impossible, which could have contributed to the l-Tyr accumulation. Interestingly, aroG and cg0975, which have received previous attention for aromatic amino acid overproduction, were not mutated. From the branch point molecule, prephenate, the change in the promoter region of pheA could be an influential contributor. In summary, we suggest that genome sequencing coupled with FBA is well poised to offer rational guidance for industrial strain development, as evidenced by these findings on carbon trafficking in C. glutamicum ATCC 21573.

Elevated nocturnal respiratory rates in the mitochondria of CAM plants: current knowledge and unanswered questions.

Crassulacean acid metabolism (CAM) is a metabolic adaptation that has evolved convergently in 38 plant families to aid survival in water-limited niches. Whilst primarily considered a photosynthetic adaptation, CAM also has substantial consequences for nocturnal respiratory metabolism. Here, we outline the history, current state and future of nocturnal respiration research in CAM plants, with a particular focus on the energetics of nocturnal respiratory oxygen consumption. Throughout the 20th century, research interest in nocturnal respiration occurred alongside initial discoveries of CAM, although the energetic and mechanistic implications of nocturnal oxygen consumption and links to the operation of the CAM cycle were not fully understood. Recent flux balance analysis (FBA) models have provided new insights into the role that mitochondria play in the CAM cycle. Several FBA models have predicted that CAM requires elevated nocturnal respiratory rates, compared to C3 species, to power vacuolar malic acid accumulation. We provide physiological data, from the genus Clusia, to corroborate these modelling predictions, thereby reinforcing the importance of elevated nocturnal respiratory rates for CAM. Finally, we outline five unanswered questions pertaining to nocturnal respiration which must be addressed if we are to fully understand and utilize CAM plants in a hotter, drier world.

Flux balance analysis-based metabolic modeling of microbial secondary metabolism: Current status and outlook.

In microorganisms, different from primary metabolism for cellular growth, secondary metabolism is for ecological interactions and stress responses and an important source of natural products widely used in various areas such as pharmaceutics and food additives. With advancements of sequencing technologies and bioinformatics tools, a large number of biosynthetic gene clusters of secondary metabolites have been discovered from microbial genomes. However, due to challenges from the difficulty of genome-scale pathway reconstruction and the limitation of conventional flux balance analysis (FBA) on secondary metabolism, the quantitative modeling of secondary metabolism is poorly established, in contrast to that of primary metabolism. This review first discusses current efforts on the reconstruction of secondary metabolic pathways in genome-scale metabolic models (GSMMs), as well as related FBA-based modeling techniques. Additionally, potential extensions of FBA are suggested to improve the prediction accuracy of secondary metabolite production. As this review posits, biosynthetic pathway reconstruction for various secondary metabolites will become automated and a modeling framework capturing secondary metabolism onset will enhance the predictive power. Expectedly, an improved FBA-based modeling workflow will facilitate quantitative study of secondary metabolism and in silico design of engineering strategies for natural product production.

The Quantitative Measurement of Peptidoglycan Components Obtained from Acidic Hydrolysis in Gram-Positive and Gram-Negative Bacteria via Hydrophilic Interaction Liquid Chromatography Coupled with Mass Spectrometry.

The high throughput in genome sequencing and metabolic model (MM) reconstruction has democratised bioinformatics approaches such as flux balance analysis. Fluxes' prediction accuracy greatly relates to the deepness of the MM curation for a specific organism starting from the cell composition. One component is the cell wall, which is a functional barrier (cell shape, exchanges) with the environment. The bacterial cell wall (BCW), including its thickness, structure, and composition, has been extensively studied in Escherichia coli but poorly described for other organisms. The peptidoglycan (PG) layer composing the BCW is usually thinner in Gram- bacteria than in Gram+ bacteria. In both bacteria groups, PG is a polymeric mesh-like structure of amino acids and sugars, including N-acetylglucosamine, N-acetylmuramic acid, and amino acids. In this study, we propose a high-throughput method to characterise and quantify PG in Gram-positive and Gram-negative bacteria using acidic hydrolysis and hydrophilic interaction liquid chromatography coupled with mass spectrometry (HILIC-MS). The method showed a relatively short time frame (11 min analytical run), low inter- and intraday variability (3.2% and 4%, respectively), and high sensitivity and selectivity (limits of quantification in the sub mg/L range). The method was successfully applied on two Gram-negative bacteria (Escherichia coli K12 MG1655, Bacteroides thetaiotaomicron DSM 2079) and one Gram-positive bacterium (Streptococcus salivarius ssp. thermophilus DSM20259). The PG concentration ranged from 1.6% w/w to 14% w/w of the dry cell weight. The results were in good correlation with previously published results. With further development, the PG concentration provided by this newly developed method could reinforce the curation of MM.

Microbiome distribution modeling using gradient descent strategies for mock, in vitro and clinical community distributions.

The human gut is home to a complex array of microorganisms interacting with the host and each other, forming a community known as the microbiome. This community has been linked to human health and disease, but understanding the underlying interactions is still challenging for researchers. Standard studies typically use high-throughput sequencing to analyze microbiome distribution in patient samples. Recent advancements in meta-omic data analysis have enabled computational modeling strategies to integrate this information into an in silico model. However, there is a need for improved parameter fitting and data integration features in microbial community modeling. This study proposes a novel alternative strategy utilizing state-of-the-art dynamic flux balance analysis (dFBA) to provide a simple protocol enabling accurate replication of abundance data composition through dynamic parameter estimation and integration of metagenomic data. We used a recurrent optimization algorithm to replicate community distributions from three different sources: mock, in vitro, and clinical microbiome. Our results show an accuracy of 98% and 96% when using in vitro and clinical bacterial abundance distributions, respectively. The proposed modeling scheme allowed us to observe the evolution of metabolites. It could provide a deeper understanding of metabolic interactions while taking advantage of the high contextualization features of GEM schemes to fit the study case. The proposed modeling scheme could improve the approach in cases where external factors determine specific bacterial distributions, such as drug intake.

Integrative teaching of metabolic modeling and flux analysis with interactive python modules.

The modeling of rates of biochemical reactions-fluxes-in metabolic networks is widely used for both basic biological research and biotechnological applications. A number of different modeling methods have been developed to estimate and predict fluxes, including kinetic and constraint-based (Metabolic Flux Analysis and flux balance analysis) approaches. Although different resources exist for teaching these methods individually, to-date no resources have been developed to teach these approaches in an integrative way that equips learners with an understanding of each modeling paradigm, how they relate to one another, and the information that can be gleaned from each. We have developed a series of modeling simulations in Python to teach kinetic modeling, metabolic control analysis, 13C-metabolic flux analysis, and flux balance analysis. These simulations are presented in a series of interactive notebooks with guided lesson plans and associated lecture notes. Learners assimilate key principles using models of simple metabolic networks by running simulations, generating and using data, and making and validating predictions about the effects of modifying model parameters. We used these simulations as the hands-on computer laboratory component of a four-day metabolic modeling workshop and participant survey results showed improvements in learners' self-assessed competence and confidence in understanding and applying metabolic modeling techniques after having attended the workshop. The resources provided can be incorporated in their entirety or individually into courses and workshops on bioengineering and metabolic modeling at the undergraduate, graduate, or postgraduate level.

On paradoxes between optimal growth, metabolic control analysis, and flux balance analysis

In Microbiology it is often assumed that growth rate is maximal. This may be taken to suggest that the dependence of the growth rate on every enzyme activity is at the top of an inverse-parabolic function, i.e. that all flux control coefficients should equal zero. This might seem to imply that the sum of these flux control coefficients equals zero. According to the summation law of Metabolic Control Analysis (MCA) the sum of flux control coefficients should equal 1 however. And in Flux Balance Analysis (FBA) catabolism is often limited by a hard bound, causing catabolism to fully control the fluxes, again in apparent contrast with a flux control coefficient of zero. Here we resolve these paradoxes (apparent contradictions) in an analysis that uses the ‘Edinburgh pathway’, the ‘Amsterdam pathway’, as well as a generic metabolic network providing the building blocks or Gibbs energy for microbial growth. We review and show that (i) optimization depends on so-called enzyme control coefficients rather than the ‘catalytic control coefficients’ of MCA's summation law, (ii) when optimization occurs at fixed total protein, the former differ from the latter to the extent that they may all become equal to zero in the optimum state, (iii) in more realistic scenarios of optimization where catalytically inert biomass is compensating or maintenance metabolism is taken into consideration, the optimum enzyme concentrations should not be expected to equal those that maximize the specific growth rate, (iv) optimization may be in terms of yield rather than specific growth rate, which resolves the paradox because the sum of catalytic control coefficients on yield equals 0, (v) FBA effectively maximizes growth yield, and for yield the summation law states 0 rather than 1, thereby removing the paradox, (vi) furthermore, FBA then comes more often to a ‘hard optimum’ defined by a maximum catabolic flux and a catabolic-enzyme control coefficient of 1. The trade-off between maintenance metabolism and growth is highlighted as worthy of further analysis.

Characterizing cancer metabolism from bulk and single-cell RNA-seq data using METAFlux

Cells often alter metabolic strategies under nutrient-deprived conditions to support their survival and growth. Characterizing metabolic reprogramming in the tumor microenvironment (TME) is of emerging importance in cancer research and patient care. However, recent technologies only measure a subset of metabolites and cannot provide in situ measurements. Computational methods such as flux balance analysis (FBA) have been developed to estimate metabolic flux from bulk RNA-seq data and can potentially be extended to single-cell RNA-seq (scRNA-seq) data. However, it is unclear how reliable current methods are, particularly in TME characterization. Here, we present a computational framework METAFlux (METAbolic Flux balance analysis) to infer metabolic fluxes from bulk or single-cell transcriptomic data. Large-scale experiments using cell-lines, the cancer genome atlas (TCGA), and scRNA-seq data obtained from diverse cancer and immunotherapeutic contexts, including CAR-NK cell therapy, have validated METAFlux’s capability to characterize metabolic heterogeneity and metabolic interaction amongst cell types.

Heuristic-enabled active machine learning: A case study of predicting essential developmental stage and immune response genes in Drosophila melanogaster.

Computational prediction of absolute essential genes using machine learning has gained wide attention in recent years. However, essential genes are mostly conditional and not absolute. Experimental techniques provide a reliable approach of identifying conditionally essential genes; however, experimental methods are laborious, time and resource consuming, hence computational techniques have been used to complement the experimental methods. Computational techniques such as supervised machine learning, or flux balance analysis are grossly limited due to the unavailability of required data for training the model or simulating the conditions for gene essentiality. This study developed a heuristic-enabled active machine learning method based on a light gradient boosting model to predict essential immune response and embryonic developmental genes in Drosophila melanogaster. We proposed a new sampling selection technique and introduced a heuristic function which replaces the human component in traditional active learning models. The heuristic function dynamically selects the unlabelled samples to improve the performance of the classifier in the next iteration. Testing the proposed model with four benchmark datasets, the proposed model showed superior performance when compared to traditional active learning models (random sampling and uncertainty sampling). Applying the model to identify conditionally essential genes, four novel essential immune response genes and a list of 48 novel genes that are essential in embryonic developmental condition were identified. We performed functional enrichment analysis of the predicted genes to elucidate their biological processes and the result evidence our predictions. Immune response and embryonic development related processes were significantly enriched in the essential immune response and embryonic developmental genes, respectively. Finally, we propose the predicted essential genes for future experimental studies and use of the developed tool accessible at http://heal.covenantuniversity.edu.ng for conditional essentiality predictions.

Cytosolic malic enzyme and glucose-6-phosphate dehydrogenase modulate redox balance in NSCLC with acquired drug resistance.

Lung cancer cells often show elevated levels of reactive oxygen species (ROS) and nicotinamide adenine dinucleotide phosphate (NADPH). However, the connections between deregulated redox homeostasis in different subtypes of lung cancer and acquired drug resistance in lung cancer have not yet been fully established. Herein, we analyzed different subtypes of lung cancer data reported in the Cancer Cell Line Encyclopedia (CCLE) database, the Cancer Genome Atlas program (TCGA), and the sequencing data obtained from a gefitinib-resistant non-small-cell lung cancer (NSCLC) cell line (H1975GR). Using flux balance analysis (FBA) model integrated with multiomics data and gene expression profiles, we identified cytosolic malic enzyme 1 (ME1) and glucose-6-phosphate dehydrogenase as the major contributors to the significantly upregulated NADPH flux in NSCLC tissues as compared with normal lung tissues, and gefitinib-resistant NSCLC cell line as compared with the parental cell line. Silencing the gene expression of either of these two enzymes in two osimertinib-resistant NSCLC cell lines (H1975OR and HCC827OR) exhibited strong antiproliferative effects. Our findings not only underscored the pivotal roles of cytosolic ME1 and glucose-6-phosphate dehydrogenase in regulating redox states in NSCLC cells but also provided novel insights into their potential roles in drug-resistant NSCLC cells with disturbed redox states.

A genome-scale metabolic model of the effect of dissolved oxygen on 1,3-propanediol fermentation by Klebsiella pneumoniae.

Although 1,3-propanediol (1,3-PD) is usually considered an anaerobic fermentation product from glycerol by Klebsiella pneumoniae, microaerobic conditions proved to be more conducive to 1,3-PD production. In this study, a genome-scale metabolic model (GSMM) specific to K. pneumoniae KG2, a high 1.3-PD producer, was constructed. The iZY1242 model contains 2090 reactions, 1242 genes and 1433 metabolites. The model was not only able to accurately characterise cell growth, but also accurately simulate the fed-batch 1,3-PD fermentation process. Flux balance analyses by iZY1242 was performed to dissect the mechanism of stimulated 1,3-PD production under microaerobic conditions, and the maximum yield of 1,3-PD on glycerol was 0.83mol/mol under optimal microaerobic conditions. Combined with experimental data, the iZY1242 model is a useful tool for establishing the best conditions for microaeration fermentation to produce 1,3-PD from glycerol in K. pneumoniae.