Integrated isotope-assisted metabolomics and 13C metabolic flux analysis reveals metabolic flux redistribution for high glucoamylase production by Aspergillus niger
© Lu et al. 2015
Received: 27 May 2015
Accepted: 31 August 2015
Published: 17 September 2015
Aspergillus niger is widely used for enzyme production and achievement of high enzyme production depends on the comprehensive understanding of cell’s metabolic regulation mechanisms.
In this paper, we investigate the metabolic differences and regulation mechanisms between a high glucoamylase-producing strain A. niger DS03043 and its wild-type parent strain A. niger CBS513.88 via an integrated isotope-assisted metabolomics and 13C metabolic flux analysis approach. We found that A. niger DS03043 had higher cell growth, glucose uptake, and glucoamylase production rates but lower oxalic acid and citric acid secretion rates. In response to above phenotype changes, A. niger DS03043 was characterized by an increased carbon flux directed to the oxidative pentose phosphate pathway in contrast to reduced flux through TCA cycle, which were confirmed by consistent changes in pool sizes of metabolites. A higher ratio of ATP over AMP in the high producing strain might contribute to the increase in the PP pathway flux as glucosephosphate isomerase was inhibited at higher ATP concentrations. A. niger CBS513.88, however, was in a higher redox state due to the imbalance of NADH regeneration and consumption, resulting in the secretion of oxalic acid and citric acid, as well as the accumulation of intracellular OAA and PEP, which may in turn result in the decrease in the glucose uptake rate.
Aspergillus niger, one of the most important filamentous fungi, has been widely used for the production of organic acids [1–3], such as citric and gluconic acid, and for the production of enzymes [4, 5] including glucoamylase, proteases, cellulase and glucose oxidase due to its excellent capacity in expression and secretion of proteins.
Despite the wide usage of A. niger for enzyme production, metabolic regulation mechanisms for high enzyme production by A. niger are still largely unclear. The complexity in metabolic regulation for high enzyme production calls for the introduction of multi-omics analysis. Metabolomics and fluxomics analysis, as an important complement to the commonly used genomics, transcriptomics and proteomics analysis, have undergone a quick development during the past years [6, 7]. With the aid of metabolomics analysis, the bottleneck for the target product synthesis can be identified , thus the addition of the pivotal limiting precursors can largely improve the strain performance. In the perturbed experiments, the profiles in dynamic changes of metabolites pool sizes can be used to assess the kinetic information of intracellular enzymes and finally estimate the limiting steps for the cell metabolism [9, 10]. Furthermore, the metabolomics play an important role in revealing the effects of energy metabolism and accumulation of polysaccharides on the strain physiology .
13C metabolic flux analysis has become a powerful tool to study the strain metabolic properties due to its capability of providing precise in vivo flux data [12–14]. For fructofuranosidase production by A. niger, Driouch et al.  reported that the production of recombinant enzyme led to a flux redistribution and in the recombinant strain the cytosolic pentose phosphate pathway (PPP) and mitochondrial malic enzyme was activated, indicating that the supply of NADPH is essential for efficient production of this protein. Klein et al.  observed that there is a positive relation between the lipid content and protein secretion by Schizosaccharomyces pombe. The latest 13C metabolic flux analysis also revealed that the imbalance in regeneration and utilization of cofactors result in the secretion of by-products .
To comprehensively understand the metabolic regulation mechanism for high glucoamylase production by A. niger, an integrated isotope-assisted metabolomics and 13C metabolic flux analysis was performed in this work. The integrative analysis of 13C MFA and metabolomics made then possible to identify the metabolic regulation mechanisms for high glucoamylase production by A. niger from the clear distinctions between the two strains in energy and redox metabolism.
Physiology of A. niger DS03043 and CBS513.88
Kinetic parameters and carbon recovery for A. niger DS03043 and CBS513.88
A. niger DS03043
A. niger CBS513.88
0.14 ± 0.01
0.10 ± 0.01
qS (mmol/gDCW h)
1.83 ± 0.01
1.46 ± 0.01
qGA (μmol/gDCW h)a
0.40 ± 0.03
0.05 ± 0.00
qCIT (mmol/gDCW h)
0.02 ± 0.00
0.11 ± 0.01
qOA (mmol/gDCW h)
0.11 ± 0.01
0.35 ± 0.02
qO2 (mmol/gDCW h)
2.94 ± 0.03
2.67 ± 0.02
qCO2 (mmol/gDCW h)
3.21 ± 0.02
2.79 ± 0.05
0.49 ± 0.03
0.47 ± 0.02
0.12 ± 0.01
0.02 ± 0.00
0.01 ± 0.00
0.07 ± 0.01
0.02 ± 0.00
0.07 ± 0.01
0.30 ± 0.02
0.32 ± 0.03
Estimation of intracellular metabolic flux and metabolites pool sizes
Pool sizes of amino acids for A. niger DS03043 and CBS513.88
A. niger DS03043
A. niger CBS513.88
P value of T test
Pool sizes of nucleotides and coenzymes in A. niger DS03043 and CBS513.88 during batch cultivations
A. niger DS03043
A. niger CBS513.88
Central carbon metabolism
The central carbon metabolism, including EMP, PP and TCA cycle pathways, exhibited obvious differences between two strains. Firstly, consistent with a larger specific glucose uptake rate, most of the absolute fluxes in A. niger DS03043 were relatively higher than those of A. niger CBS513.88 (Fig. 4). In the following comparison, all fluxes are normalized to the specific glucose uptake rate (Fig. 3). The relative flux channeled into the PP pathway was obviously larger in A. niger DS03043 than in A. niger CBS513.88. Accordingly, the flux through EMP in A. niger DS03043 was decreased. Exceeding the demand for biomass synthesis, a fraction of pentose was recycled back into the EMP pathway in the form of F6P and GAP. As the main pathway for NADH formation, the flux through TCA cycle in A. niger DS03043 was only half of that in wild type strain (Fig. 3). Such results agree with previous studies involving Aspergilli [15, 19] and the detailed comparison could be found in Additional file 3.
The pool sizes of metabolites from the core carbon metabolic network showed a relationship with the flux distribution for both strains. Due to the increased flux from G6P to 6PG, the pool sizes of R5P and S7P in A. niger DS03043 were clearly larger (P value <0.05) than those in A. niger CBS513.88 (Fig. 4). In the upper EMP pathway, the pool sizes of G6P and F6P were nearly constant for the two strains, indicating that their pool sizes were strictly controlled. However, in the lower EMP pathway, except for PYR, the pool sizes of 3PG and PEP were largely increased in A. niger CBS513.88 (P value <0.05). It was also found that the intracellular OAA, CIT, isCIT and AKG pool sizes were obviously bigger in A. niger CBS513.88 than those in A. niger DS03043 (P value <0.05). On the contrary, in the TCA cycle of CBS513.88, the pool sizes of MAL and FUM were slightly smaller than A. niger DS03043. It could be speculated that there exists a close relationship between the accumulation of intermediates in the TCA cycle and the excretion of large amounts of oxalic acid and citric acid in A. niger CBS513.88.
Amino acids metabolism
Energy and redox metabolism
Estimated production and consumption of NADPH, NADH, ATP ( mmol/gDCW h)
A. niger DS03043
A. niger CBS513.88
A. niger DS03043
A. niger CBS513.88
A. niger DS03043
A. niger CBS513.88
Amino acid synthesis
AcCoA formation by Ac
Among the three pathways to generate ATP, including EMP, TCA cycle and the oxidative phosphorylation, the oxidative phosphorylation was the main source of ATP in both strains and the total specific generation rate of ATP were 12.38 and 11.22 mmol ATP/gDCW h respectively (Table 4) in A. niger DS03043 and CBS513.88. As shown in Table 4, most ATP was used for the strain growth and the fractions of ATP for biomass formation were about 68 and 57 % of total ATP generated in A. niger DS03043 and CBS513.88. The specific rate of ATP being used for glucoamylase synthesis was increased from 0.17 mmol ATP/gDCW h in A. niger CBS513.88 to 1.3 mmol ATP/gDCW h in A. niger DS03043. The rest ATP was potentially consumed for the non-growth maintenance requirements, which was calculated to be 3.73 and 5.68 mmol ATP/gDCW h respectively in A. niger DS03043 and CBS 513.88.
Amino acid metabolism
The pool sizes of intracellular free amino acids are influenced by multiple factors , including the energy availability, the flux distribution, as well as the synthesis and degradation of specific proteins. In this work, we found that there existed no direct relations between the pool sizes of amino acids and the synthesis of glucoamylase (data not shown), except for Ala, Asp, Glu and Pro (Table 2), which may be due to the fact that during the batch phase, the main part of carbon substrate (about 80 %) was used for strain growth and aerobic respiration. Similar to previous reports , the overall distribution of amino acid pools were slightly affected by the enzyme production levels, indicating the robustness of strain metabolism. It has been reported that , for Pichia pastoris, the single amino acid pool size was correlated with overall cell protein amino acid compositional variations. However, in this study, we found that the amino acids pool sizes were directly related to the precursor metabolites. The average pool size of amino acids from the families of aromatic and histidine amino acids were lower on the whole (Fig. 5). Therefore, we speculated that the synthesis of these specific amino acids might be the bottleneck for further improvement in enzyme production.
Mechanism for PP pathway regulation
TCA cycle pathway and energy metabolism
The TCA cycle pathway is the main source for NADH formation , as evident from the observation that the generation rate of NADH from the TCA cycle pathway was nearly twofold of that from the EMP pathway (Table 4). Via oxidative phosphorylation, NADH was transformed into ATP to meet the energy demand for maintenance, growth and product formation. For protein synthesis, the changes in the TCA cycle flux might depend on the strain. For Schizosaccharomyces pombe , it has been reported that to compensate for metabolic burden for heterologous protein synthesis, the flux through TCA cycle pathway was increased. As a result, there exists a positive correlation between the protein yield and the flux directed into the TCA cycle pathway. In contrast, for A. niger producing fructofuranosidase , the flux of the TCA cycle pathway dropped. Similarly, we found that this flux changed in the same direction in A. niger DS03043 (Fig. 3) and, more strikingly, the pool sizes of TCA cycle intermediates, like OAA, CIT, is CIT and AKG were also reduced (Fig. 4), which was accompanied with the reduction in by-products (oxalic acid and citric acid) formation (Table 1). The reported m_ATP for A. niger was 1.9 mmol ATP/gDCW h , while the calculated maintenance energy for both strains in this study were obviously higher, especially for A. niger CBS513.88. In general, the ATP cost for non-growth maintenance depends on product and by-products formation and the type of substrate used. It has been reported  that the cellular maintenance ATP consumption in engineered E. coli cells producing fatty acids was nearly twofold of that in control strain due to an increased osmotic stress resulting from product accumulation. Chung and co-workers  also reported that maintenance requirements for P. pastoris growing on glucose and glycerol–methanol mixed media were 2.3 and about 6 mmol ATP/gDCW h, respectively. Moreover, the transmembrane transport of organic acids was also accompanied with loss of energy . Therefore, it could be concluded here that A. niger CBS513.88 needs more non-growth maintenance energy due to the accumulation and secretion of oxalic and citric acids, which, to some extent, decreased the utilization efficiency of energy in A. niger CBS513.88.
How metabolites modulate metabolic flux ?
Among thousands of metabolites, cofactors participate in 300 oxidation and reduction reactions  and play an important role in regulating the cell metabolic activities. Generally, the formation of by-products are correlated with the cofactors metabolism, for example, the glycerol synthesis could contribute to the recycle of NADH . It has been reported that the substrate uptake rate was depressed under higher intracellular ratio of NADH/NAD+ . However, the mechanisms for the modulation in fluxes via cofactors are generally not clear. The physiological parameters of A. niger CBS513.88 and DS03043, as found in this study, indicated that with large amounts of by-products formation the specific growth rate of the wild type strain was decreased. Interestingly, we found that the ratios of NADH/NAD+ and NADPH/NADP+ in A. niger CBS513.88 were higher than that in A. niger DS03043 (Fig. 6b, c), meaning that the A. niger CBS513.88 was in high redox condition. Accompanied with the high ratio of NADH/NAD+, the intermediates in the TCA cycle (CIT, is CIT, AKG and OAA) were accumulated in A. niger CBS513.88 (Fig. 4). The accumulation of organic acid from the TCA cycle also happened under oxygen limited conditions when the NADH recycle was hindered . As reported by Heux et al.  and Hou et al. , decreasing the intracellular NADH level by introducing NADH oxidase or alternative oxidase lowered the by-products formation greatly. So it could be speculated that the higher ratio of NADH/NAD+ contributes to the increase in pool sizes of these intermediates from TCA cycle. As PEP was linked with the TCA cycle through anaplerotic reactions, the pool size of PEP in the cytoplasm was also obviously increased accordingly in A. niger CBS513.88 (Fig. 4).
According to Ogawa et al. , the accumulation of PEP could regulate the glucose uptake and glycolytic flux in cells by inhibiting glucokinase and glucose-6-phosphate isomerase. For mammalian cells culture, it was found that the high lactate level can exert an inhibitory effect on glycolysis enzyme phosphofructokinase, thus reducing the glycolytic flux . So for A. niger CBS513.88 in this work, the lower absolute glycolytic flux (Fig. 4) might be caused by the inhibitory effect of organic acid, including PEP and 3PG, which in turn results in a lower specific growth rate (Table 1). As the accumulation of organic acids is related to the cell redox condition, the regulation of glycolytic flux actually seems to depend on the redox condition of the strain. In a comprehensive view, the possible regulation mechanism in flux and metabolites pools by the ratio of NADH/NAD+ could be summarized by Fig. 7b. With integrative analysis of metabolomics and 13C MFA, this work provides new clues to illustrate the interactive effects of flux and metabolites, which could be applied in other strains or process studies as well.
In this study, the combination of 13C metabolic flux analysis and quantitative metabolomics was adopted for the first time to investigate the metabolic changes in A. niger correlated with high glucoamylase production, initially revealing the metabolic regulation mechanisms between intracellular flux redistribution and metabolites pool sizes. The consistency in 13C metabolic flux analysis and metabolites quantification indicated that an imbalance of NADH formation and consumption led to the accumulation and secretion of organic acids in A. niger CBS513.88. The significantly reduced by-products formation rates resulted in an increased energy utilization efficiency and specific growth rate in A. niger DS03043. To enable a more efficient enzyme production and faster cell growth, the total amino acids pool size, the ratio of ATP to AMP and the flux channeled into the PP pathway were obviously improved in A. niger DS03043, which bring new insights on metabolic engineering targets for improved enzyme production using A. niger.
Strain and cultivations
Two A. niger strains with different glucoamylase production levels were kindly donated by DSM Corporation (Delft, the Netherlands). A. niger DS03043, a high glucoamylase-producing strain, contains seven copies of glucoamylase glaA gene, while A. niger CBS513.88, as a wild-type parent strain, contains only one glaA gene copy .
To obtain spores, Petri dishes containing PDA (Potato Dextrose Agar) medium were incubated with spores from a frozen stock (stored in 50 % glycerin at −80 °C). The seed cultivation was conducted in 500 mL shaking flasks with three baffles. Each flask was inoculated with 107 spores per 100 mL broth. The batch cultivations were conducted in a 5 L fermenter with an electronic balance, inoculated with the seed broth from the shake flasks after 24 h seed culture at 34 °C, with the speed of shaking table controlled at 150 rpm. The working volume for 5 L fermenter was 3.0 L. The agitation and aeration in the 5 L fermenter were kept at 375 rpm and 1 vvm respectively, with overpressure maintained at 0.05 MPa and temperature at 34 °C during the whole process. The pH was automatically maintained at pH 4.5 with 5 % (w/w) NH4·OH. The cultivations for each strain were performed in three independent cultivation replicates.
The medium for seed culture contained (g/L) glucose·H2O 22 and spray dried corn steep liquor 20. Before sterilization, the initial pH was adjusted to 6.5 using 3 mol/L NaOH. A chemically modified defined medium was used for all cultivations in this study, which contained (g/L) glucose·H2O 20, KH2PO4 3, NaH2PO4·H2O 1.5, (NH4)2SO4 3, MgSO4·7H2O 1.0, CaCl2·2H2O 0.1, ZnCl2 0.02, CuSO4·5H2O 0.015, CoCl2·6H2O 0.015, MnSO4·H2O 0.04, FeSO4·7H2O 0.3.
13C labeling experiments
Considering the high cost, the 13C labeling experiments were conducted in a reactor system consisting of a 5 L bioreactor and a 250 mL mini bioreactor equipped with pH and oxygen sensors. Before the formal experiment, the agitation and aeration in the mini bioreactor were carefully optimized to guarantee that the strain had similar physiological parameters in two scales of fermenters. The initial glucose concentration in the 5 L bioreactor was 10 g/L. When the glucose concentration decreased to approximately 5 g/L, about 170 mL broth was rapidly transferred from the 5 L bioreactor into the mini bioreactor. Then the labeled glucose including 0.56 g [U–13C]–Glc and 1.41 g [1–13C]–Glc (isotopic enrichment 98–99 %, Cambridge Isotope Laboratories, Inc.) were added as a pulse into the mini bioreactor.
Quantification of biomass and enzyme activity
10 mL broth was filtered through pre-dried and pre-weighed suction filter paper. Before filtering, the filter should be dried to a constant weight (24 h at 80 °C). To remove solutes, the samples were rinsed three times with deionized water. Then the wet filters with the biomass were put in the 80 °C oven and dried for 24 h. The dried filter was re-weighed immediately. A previous protocol was adopted to determine the enzyme activity in all samples .
Analysis of extracellular glucose and organic acids
A glucose analyzer (Shandong Academy of Sciences, China) was used for residual glucose concentration analysis. The amounts of extracellular organic acids (acetic, citric, oxalic, fumaric, malic, pyruvic, and succinic acids) were analyzed using a high-performance liquid chromatography (HPLC) system. An ion-exclusion column (Hi-Plex H, Agilent) was eluted at 50 °C with 10 mM H2SO4 at a flow rate of 0.5 mL/min and connected with an absorbance detector spectrophotometer at 210 nm. No acetic acid was detected under these conditions of cultivation.
Sampling and quantification of intracellular metabolites concentration
The protocol for quantification of intracellular metabolites pool size was adapted from the works of Douma et al. . With a dedicated fast-sample equipment, 1–2 mL broth was rapidly withdrawn from sensor bioreactors in tubes with 10 mL precooled quench solution (−27.6 °C 40 % v/v methanol solution) at 30, 60, 90, 120, 180, 240 and 300 min since the start of 13C labelled experiments. The tubes were weighed before and after sampling procedure in order to estimate the exact amount of broth taken. Subsequently, fast filtration with a vacuum pump was conducted to remove the extracellular metabolites as soon as possible. Then, 120 mL precooled quench solution was used to wash the filter cake. To accurately determine the metabolites concentration, Isotope Dilution Mass Spectrometry (IDMS)  was adopted in this study. During the extraction process, the washed filter cake and 13C internal standard solution were together added into the pre-warmed 25 mL 75 % (v/v) ethanol solution and the extraction continued for 3 min at 95 °C. The metabolites pools were determined with the UPLC-MS/MS (Thermo Fisher Scientific Coporation).
For analysis of the intracellular pool size of NADH, NAD+, NADPH, and NADP+, the EnzyFluo™ Assay Kit (BioAssay Systems Coporation) was used. Firstly, 1 mL broth was quenched with 10 mL quench solution. The supernatant was discarded after centrifugation at −4 °C and 6000 rpm. The remaining processes were referred to the standard specifications (BioAssay Systems Coporation).
Analysis of amino acids enrichment with GC–MS
To obtain the enrichment of intracellular free amino acids during the 13C labelled experiments, the extracting solution was concentrated to 100 μL of supernatant before being transferred to a glass vial. After lyophilization, 80 μL acetonitrile and 80 μL of N-methyl-N-(tert-butyldimethylsilyl) trifluoroacetamide (MTBSTFA, Thermo Scientific) were added and the vial was incubated for 2 h at 60 °C. After that, the sample was centrifuged (10,000g, 2 min) and 160 μL of the supernatant was transferred to a GC glass vial with an insert. The sample was then analyzed by GC–MS instrument coupled to a 5975 C MSD single quadrupole mass spectrometer (Agilent, Santa Clara, CA, USA).
Calculation of metabolic fluxes based on 13C-labeled experiments
HL carried out bioreactor cultivations, 13C labelled experiments, 13C flux calculation and metabolomics analysis. XL participated in the measurement of metabolites concentration. MH and JC participated in the design and coordination of the study. JX, YZ, SZ and NH helped to draft the manuscript. All authors read and approved the final manuscript.
This work was financially supported by Royal DSM (Delft, the Netherlands) and partially supported by National Basic Research Program (973 Program 2013CB733600), NWO-MoST Joint Program (2013DFG32630), National High Technology Research and Development Program of China (863 Program 2012AA021201) and National Key Technology R&D Program (2012BAI44G01).
Compliance with ethical guidelines
Competing interests The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Andersen MR, Nielsen J. Current status of systems biology in Aspergilli. Fungal Genet Biol. 2009;46:180–90.View ArticleGoogle Scholar
- Meyer V, Wu B, Ram AFJ. Aspergillus as a multi-purpose cell factory: current status and perspectives. Biotechnol Lett. 2011;33:469–76.View ArticleGoogle Scholar
- Knuf C, Nielsen J. Aspergilli: systems biology and industrial applications. Biotechnol J. 2012;7:1147–55.View ArticleGoogle Scholar
- Ward OP. Production of recombinant proteins by filamentous fungi. Biotechnol Adv. 2012;30:1119–39.View ArticleGoogle Scholar
- Sun H, Zhao P, Ge X, Xia Y, Hao Z, Liu J, et al. Recent advances in microbial raw starch degrading enzymes. Appl Biochem Biotechnol. 2010;160:988–1003.View ArticleGoogle Scholar
- Sévin DC, Kuehne A, Zamboni N, Sauer U. Biological insights through nontargeted metabolomics. Curr Opin Biotechnol. 2015;34:1–8.View ArticleGoogle Scholar
- Ellis DI, Goodacre R. Metabolomics-assisted synthetic biology. Curr Opin Biotechnol. 2012;23:22–8.View ArticleGoogle Scholar
- Xia M, Huang D, Li S, Wen J, Jia X, Chen Y. Enhanced FK506 production in Streptomyces tsukubaensis by rational feeding strategies based on comparative metabolic profiling analysis. Biotechnol Bioeng. 2013;110:2717–30.View ArticleGoogle Scholar
- Taymaz-Nikerel H, De Mey M, Baart G, Maertens J, Heijnen JJ, van Gulik W. Changes in substrate availability in Escherichia coli lead to rapid metabolite, flux and growth rate responses. Metab Eng. 2013;16:115–29.View ArticleGoogle Scholar
- de Jonge LP, Buijs NAA, ten Pierick A, Deshmukh A, Zhao Z, Kiel JAKW, et al. Scale-down of penicillin production in Penicillium chrysogenum. Biotechnol J. 2011;6:944–58.View ArticleGoogle Scholar
- Jordà J, Rojas H, Carnicer M, Wahl A, Ferrer P, Albiol J. Quantitative metabolomics and instationary 13C-metabolic flux analysis reveals impact of recombinant protein production on trehalose and energy metabolism in Pichia pastoris. Metabolites. 2014;4:281–99.View ArticleGoogle Scholar
- Zamboni N, Fendt S-M, Ruhl M, Sauer U. 13C-based metabolic flux analysis. Nat Protoc. 2009;4:878–92.View ArticleGoogle Scholar
- Antoniewicz MR. 13C metabolic flux analysis: optimal design of isotopic labeling experiments. Curr Opin Biotechnol. 2013;24:1116–21.View ArticleGoogle Scholar
- Niedenführ S, Wiechert W, Nöh K. How to measure metabolic fluxes: a taxonomic guide for 13C fluxomics. Curr Opin Biotechnol. 2015;34:82–90.View ArticleGoogle Scholar
- Driouch H, Melzer G, Wittmann C. Integration of in vivo and in silico metabolic fluxes for improvement of recombinant protein production. Metab Eng. 2012;14:47–58.View ArticleGoogle Scholar
- Klein T, Lange S, Wilhelm N, Bureik M, Yang T-H, Heinzle E, et al. Overcoming the metabolic burden of protein secretion in Schizosaccharomyces pombe—a quantitative approach using 13C-based metabolic flux analysis. Metab Eng. 2014;21:34–45.View ArticleGoogle Scholar
- Nie Y, Huang M, Lu J, Qian J, Lin W, Chu J, et al. Impacts of high β-galactosidase expression on central metabolism of recombinant Pichia pastoris GS115 using glucose as sole carbon source via 13C metabolic flux analysis. J Biotechnol. 2014;187:124–34.View ArticleGoogle Scholar
- Young JD. INCA: a computational platform for isotopically non-stationary metabolic flux analysis. Bioinformatics. 2014;30:1333–5.View ArticleGoogle Scholar
- Pedersen H, Carlsen M, Nielsen J. Identification of enzymes and quantification of metabolic fluxes in the wild type and in a recombinant Aspergillus oryzae strain. Appl Environ Microbiol. 1999;65:11–9.Google Scholar
- Ball WJ, Atkinson DE. Adenylate energy charge in Saccharomyces cerevisiae during starvation. J Bacteriol. 1975;121:975–82.Google Scholar
- Chapman AG, Fall L, Atkinson DE. Adenylate energy charge in Escherichia coli during growth and starvation. J Bacteriol. 1971;108:1072–86.Google Scholar
- Vrabl P, Mutschlechner W, Burgstaller W. Dynamics of energy charge and adenine nucleotides during uncoupling of catabolism and anabolism in Penicillium ochrochloron. Mycol Res. 2009;113:1422–32.View ArticleGoogle Scholar
- Carnicer M, ten Pierick A, van Dam J, Heijnen JJ, Albiol J, van Gulik W, et al. Quantitative metabolomics analysis of amino acid metabolism in recombinant Pichia pastoris under different oxygen availability conditions. Microb Cell Fact. 2012;11:83–93.View ArticleGoogle Scholar
- Bordbar A, Monk JM, King ZA, Palsson BO. Constraint-based models predict metabolic and associated cellular functions. Nat Rev Genet. 2014;15:107–20.View ArticleGoogle Scholar
- Heyland J, Fu J, Blank LM, Schmid A. Carbon metabolism limits recombinant protein production in Pichia pastoris. Biotechnol Bioeng. 2011;108:1942–53.View ArticleGoogle Scholar
- Wasylenko TM, Ahn WS, Stephanopoulos G. The oxidative pentose phosphate pathway is the primary source of NADPH for lipid overproduction from glucose in Yarrowia lipolytica. Metab Eng. 2015;30:27–39.View ArticleGoogle Scholar
- Chen H, Hao G, Wang L, Wang H, Gu Z, Liu L, et al. Identification of a critical determinant that enables efficient fatty acid synthesis in oleaginous fungi. Sci Rep. 2015;5:11247.View ArticleGoogle Scholar
- Park SH, Kim HU, Kim TY, Park JS, Kim S-S, Lee SY. Metabolic engineering of Corynebacterium glutamicum for l-arginine production. Nat Commun. 2014;5:4618.Google Scholar
- Ghorbaniaghdam A, Henry O, Jolicoeur M. An in silico study of the regulation of CHO cells glycolysis. J Theor Biol. 2014;357:112–22.View ArticleGoogle Scholar
- Takama M, Nosoh Y. Effect of ATP on glucose-6-phosphate isomerase from Bacillus caldotenax. Biochim Biophys Acta. 1982;705:127–30.View ArticleGoogle Scholar
- Zhou J, Liu L, Shi Z, Du G, Chen J. ATP in current biotechnology: regulation, applications and perspectives. Biotechnol Adv. 2009;27:94–101.View ArticleGoogle Scholar
- David H, Åkesson M, Nielsen J. Reconstruction of the central carbon metabolism of Aspergillus niger. Eur J Biochem. 2003;270:4243–53.View ArticleGoogle Scholar
- He L, Xiao Y, Gebreselassie N, Zhang F, Antoniewicz MR, Tang YJ, et al. Central metabolic responses to the overproduction of fatty acids in Escherichia coli based on 13C-metabolic flux analysis. Biotechnol Bioeng. 2014;111:575–85.View ArticleGoogle Scholar
- Chung B, Selvarasu S, Camattari A, Ryu J, Lee H, Ahn J, et al. Genome-scale metabolic reconstruction and in silico analysis of methylotrophic yeast Pichia pastoris for strain improvement. Microb Cell Fact. 2010;9:50–65.View ArticleGoogle Scholar
- van Maris AJA, Winkler AA, Porro D, van Dijken JP, Pronk JT. Homofermentative lactate production cannot sustain anaerobic growth of engineered Saccharomyces cerevisiae: possible consequence of energy-dependent lactate export. Appl Environ Microbiol. 2004;70:2898–905.View ArticleGoogle Scholar
- Förster J, Famili I, Fu P, Palsson BØ, Nielsen J. Genome-scale reconstruction of the Saccharomyces cerevisiae metabolic network. Genome Res. 2003;13:244–53.View ArticleGoogle Scholar
- Rigoulet M, Aguilaniu H, Avéret N, Bunoust O, Camougrand N, Grandier-Vazeille X, et al. Organization and regulation of the cytosolic NADH metabolism in the yeast Saccharomyces cerevisiae. Mol Cell Biochem. 2004;256:73–81.View ArticleGoogle Scholar
- Heux S, Cachon R, Dequin S. Cofactor engineering in Saccharomyces cerevisiae: expression of a H2O-forming NADH oxidase and impact on redox metabolism. Metab Eng. 2006;8:303–14.View ArticleGoogle Scholar
- Diano A, Peeters J, Dynesen J, Nielsen J. Physiology of Aspergillus niger in oxygen-limited continuous cultures: influence of aeration, carbon source concentration and dilution rate. Biotechnol Bioeng. 2009;103:956–65.View ArticleGoogle Scholar
- Hou J, Lages NF, Oldiges M, Vemuri GN. Metabolic impact of redox cofactor perturbations in Saccharomyces cerevisiae. Metab Eng. 2009;11:253–61.View ArticleGoogle Scholar
- Ogawa T, Mori H, Tomita M, Yoshino M. Inhibitory effect of phosphoenolpyruvate on glycolytic enzymes in Escherichia coli. Res Microbiol. 2007;158:159–63.View ArticleGoogle Scholar
- Mulukutla BC, Gramer M, Hu W-S. On metabolic shift to lactate consumption in fed-batch culture of mammalian cells. Metab Eng. 2012;14:138–49.View ArticleGoogle Scholar
- van Dijck PWM, Selten GCM, Hempenius RA. On the safety of a new generation of DSM Aspergillus niger enzyme production strains. Regul Toxicol Pharm. 2003;38:27–35.View ArticleGoogle Scholar
- Lu H, Li C, Tang W, Wang Z, Xia J, Zhang S, et al. Dependence of fungal characteristics on seed morphology and shear stress in bioreactors. Bioprocess Biosyst Eng. 2015;38:917–28.View ArticleGoogle Scholar
- Douma RD, de Jonge LP, Jonker CTH, Seifar RM, Heijnen JJ, van Gulik WM. Intracellular metabolite determination in the presence of extracellular abundance: application to the penicillin biosynthesis pathway in Penicillium chrysogenum. Biotechnol Bioeng. 2010;107:105–15.View ArticleGoogle Scholar
- Wu L, Mashego MR, van Dam JC, Proell AM, Vinke JL, Ras C, et al. Quantitative analysis of the microbial metabolome by isotope dilution mass spectrometry using uniformly C-13-labeled cell extracts as internal standards. Anal Biochem. 2005;336:164–71.View ArticleGoogle Scholar
- Antoniewicz MR, Kelleher JK, Stephanopoulos G. Elementary metabolite units (EMU): a novel framework for modeling isotopic distributions. Metab Eng. 2007;9:68–86.View ArticleGoogle Scholar
- Antoniewicz MR, Kelleher JK, Stephanopoulos G. Determination of confidence intervals of metabolic fluxes estimated from stable isotope measurements. Metab Eng. 2006;8:324–37.View ArticleGoogle Scholar