Metabolic flux analysis and the NAD(P)H/NAD(P)+ ratios in chemostat cultures of Azotobacter vinelandii
© The Author(s) 2018
Received: 10 August 2017
Accepted: 15 January 2018
Published: 22 January 2018
Azotobacter vinelandii is a bacterium that produces alginate and polyhydroxybutyrate (P3HB); however, the role of NAD(P)H/NAD(P)+ ratios on the metabolic fluxes through biosynthesis pathways of these biopolymers remains unknown. The aim of this study was to evaluate the NAD(P)H/NAD(P) + ratios and the metabolic fluxes involved in alginate and P3HB biosynthesis, under oxygen-limiting and non-limiting oxygen conditions.
The results reveal that changes in the oxygen availability have an important effect on the metabolic fluxes and intracellular NADPH/NADP+ ratio, showing that at the lowest OTR (2.4 mmol L−1 h−1), the flux through the tricarboxylic acid (TCA) cycle decreased 27.6-fold, but the flux through the P3HB biosynthesis increased 6.6-fold in contrast to the cultures without oxygen limitation (OTR = 14.6 mmol L−1 h−1). This was consistent with the increase in the level of transcription of phbB and the P3HB biosynthesis. In addition, under conditions without oxygen limitation, there was an increase in the carbon uptake rate (twofold), as well as in the flux through the pentose phosphate (PP) pathway (4.8-fold), compared to the condition of 2.4 mmol L−1 h−1. At the highest OTR condition, a decrease in the NADPH/NADP+ ratio of threefold was observed, probably as a response to the high respiration rate induced by the respiratory protection of the nitrogenase under diazotrophic conditions, correlating with a high expression of the uncoupled respiratory chain genes (ndhII and cydA) and induction of the expression of the genes encoding the nitrogenase complex (nifH).
We have demonstrated that changes in oxygen availability affect the internal redox state of the cell and carbon metabolic fluxes. This also has a strong impact on the TCA cycle and PP pathway as well as on alginate and P3HB biosynthetic fluxes.
Azotobacter vinelandii is a Gram negative bacterium that produces two polymers of biotechnological importance, poly-3-hydroxybutyrate (P3HB), an intracellular polyester of the polyhydroxyalkanoates (PHAs) family, and alginate, an extracellular polysaccharide [1, 2]. This bacterium fixes nitrogen under aerobic conditions, being able to protect its oxygen sensitive nitrogenase complex from damage by respiratory protection . Due to this characteristic, A. vinelandii exhibits a high respiratory activity, especially when exposed to high oxygen concentrations . Therefore, when it is grown without control of the dissolved oxygen tension (DOT), the cultures operate at a DOT close to zero (microaerophilic conditions). Under this condition, the maximum oxygen transfer rate (OTRmax) and/or the maximum oxygen uptake rate (OURmax) have been used to evaluate the respiratory metabolism of A. vinelandii [5, 6]. Several studies been focused on how the agitation rate and the OTRmax affect alginate and P3HB production [7–10], showing that the OTRmax is positively correlated with alginate biosynthesis, whereas P3HB biosynthesis follows an inverse relation.
Several studies have reported the relationship between the NAD(P)H+/NAD(P) cofactor levels and the metabolic flux distributions in aerobic cultures. For example, genetic manipulations aiming at perturbing NADH cofactor levels and/or regeneration rates were employed as a tool for the metabolic engineering of Escherichia coli [11–13], Lactococcus lactis , Bacillus subtilis [15, 16] and Saccharomyces cerevisiae . These studies showed that changes in the ratio of NADH/NAD+ determined the metabolic products. Other studies have shown that strategies increasing the NADPH cofactor levels improved, for example, penicillin formation , methylenomycin biosynthesis  and P3HB biosynthesis [20–23].
To our knowledge, there have been no studies about the relationship between NADH and NADPH levels and the metabolic behavior of A. vinelandii. However, it is well known that the NADPH is consumed during PHB biosynthesis and the alginate biosynthesis produces NADH in this bacterium. For example, the limiting step and the control point for alginate biosynthesis is the activity of the GDP-mannose dehydrogenase enzyme, which is involved in the irreversible oxidation of GDP-mannose to GDP-mannuronic acid , and is NAD+-dependent . This enzyme performs a double oxidation, in which two molecules of NADH are released for each GDP-mannuronic acid synthesized; therefore, the biosynthesis of alginate could be affected by changes in the intracellular redox state. So far, there is no evidence pointing at a possible effect of the reduction power on the production of this polymer. On the other hand, the accumulation of reducing power (NADH and NADPH) in cells growing under oxygen-limiting conditions is probably involved in the allosteric down-regulation of some TCA cycle enzymes such as citrate synthase and isocitrate dehydrogenase, decreasing the flux of acetyl-CoA into this cycle, thereby increasing its availability for P3HB biosynthesis . This is because a molecule of NADHP being required for each monomer of 3-hydroxybutyryl formed during P3HB biosynthesis.
Recently, in A. vinelandii growing in shake flasks it was demonstrated through metabolic flux analysis  that changes in oxygen availability have an important impact on the metabolic fluxes, which was also reflected on the alginate and P3HB yields. However, the NADH/NAD + and NADPH/NADP + ratios were not quantified, therefore preventing dissection of the interplay between the rate of redox cofactor regeneration, the intracellular carbon flux distribution and the rate of carbon substrate uptake. In the present study, the relationship of the intracellular NAD(P)H/NAD(P) + ratios and metabolic flux analysis through the central carbon metabolism were analyzed in A. vinelandii cultures grown under oxygen-limiting and non-limiting conditions.
DOT and OTR under oxygen-limiting and non-limiting conditions
Kinetic parameters during steady-state of A. vinelandii ATCC 9046, grown in continuous cultures (D = 0.08 h−1) under oxygen-limiting (300, 500 rpm) and non-limiting oxygen (700 rpm) conditions
Agitation rate (rpm)
OTRmax (mmol L−1 h−1)
Yprotein/glucose (g g−1)
qO2 (mmol g prot −1 h−1)
2.4 ± 0.32
0.15 ± 0.02
3.7 ± 0.16
0.7 ± 0.08
0.75 ± 0.04
6.2 ± 0.57
0.11 ± 0.01
8.26 ± 0.24
1.18 ± 0.18
0.34 ± 0.05
14.3 ± 1.2
0.08 ± 0.02
20 ± 0.34
1.07 ± 0.15
0.05 ± 0.002
Glucose consumption, biomass, alginate and P3HB production
The residual glucose concentrations in the steady-state were 3.5 ± 0.3, 5.5 ± 0.35 and 1.3 ± 0.2 g L−1, for the OTR of 2.4, 6.2 and 14.3 mmol L−1 h−1, respectively. Under similar conditions, Díaz-Barrera et al.  reported an affinity constant (Ks) for sucrose of 0.1 g L−1. Therefore, the cultures were not limited by the carbon source. In addition, the ammonium concentration in the steady-state under oxygen-limiting conditions (2.4 and 6.2 mmol L−1 h−1) was higher (Fig. 2e) than the K S described for ammonium in A. vinelandii (0.11 g L−1) . Nevertheless, in cultures developed at 14.3 mmol L−1 h−1, the ammonium concentration was close to zero (Fig. 2e), suggesting that these cultures were probably fixing atmospheric nitrogen. As shown in Fig. 2b, when the OTR was higher (14.3 mmol L−1 h−1), the q G increased (twofold) to 6.12 ± 0.35 mmol g−1 h−1, while the q G at the lowest OTR condition was 2.9 ± 0.2 mmol g−1 h−1, showing a close correlation with the values of qO 2 (Table 1).
The alginate and P3HB productions under the steady-state at the different OTR tested are shown in Table 1. Under oxygen-limiting conditions, the alginate production increased when the OTR was higher, from 0.7 ± 0.08 g L−1 at an OTR of 2.4 mmol L−1 h−1, to 1.18 ± 0.18 g L−1 when the cultures were conducted at 6.2 mmol L−1 h−1. A further increase in the OTR under non-limiting oxygen conditions (14.3 mmol L−1 h−1) did not increase alginate production (Table 1). Interestingly, the specific alginate production rate (q Alg ) was similar (0.61 ± 0.04 mmol g protein −1 h−1) at the high OTR values (6.2 and 14.3); and this value was 35% higher than that achieved in cultures at 2.4 mmol L−1 h−1 (0.44 ± 0.03 mmol g protein −1 h−1).
The P3HB production in A. vinelandii has been reported to be accumulated mainly under oxygen-limiting conditions [2, 30, 31]. Consistently, in this study the P3HB production was negatively correlated by the increase in OTR, as shown in Table 1. The highest P3HB concentration (0.75 ± 0.04 g L−1) was obtained in cultures developed at 2.4 mmol L−1 h−1, whereas a value of 0.05 ± 0.002 g L−1 was obtained at the highest OTR (14.3 mmol L−1 h−1).
Internal redox state, metabolic flux distribution and its corresponding gene expression
% C-mole usage in A. vinelandii ATCC 9046, grown in continuous cultures (D = 0.08 h−1) under oxygen-limiting (2.4 and 6.2 mmol L−1 h−1) and non-limiting oxygen (14.3 mmol L−1 h−1) conditions
OTRmax (mmol L−1 h−1)
% C-mol biomass
% C-mol alginate
% C-mol P3HB
% C-mol CO2
2.4 ± 0.32
6.2 ± 0.57
14.3 ± 1.2
The results of the present study show, for the first time, that the OTR affects metabolic fluxes through the central carbon metabolism network and the NAD(P)H/NAD(P) + ratios, thereby pointing at a relationship of the redox state with the alginate and P3HB production in chemostat cultures of A. vinelandii. Previous studies with this bacterium conducted in shaken flasks, under high and low aeration conditions, showed that changes in oxygen availability had a considerable impact on the growth profiles, the alginate and P3HB yields and metabolic fluxes . However, shaken flasks are not the most suitable growth system for the evaluation of metabolic fluxes, because several non-controlled culture parameters (e.g. pH, DOT) and most of the relevant variables, substrate, product, or biomass concentrations continuously vary. Therefore, in this study, the continuous culture mode was chosen, because it allows the specific growth rate to be fixed, allowing a strict comparison among different conditions at the steady-state and also to evaluate the concentrations of NADH, NADPH and metabolic fluxes.
The cellular protein yield was twofold higher in those cultures conducted under oxygen limitation (2.4 and 6.2 mmol L−1 h−1) compared to the cultivations conducted at the highest (non-limiting) OTR (14.3 mmol L−1 h−1), although the highest q G and qO 2 were obtained under this last condition (Table 1). These differences might be related to the protein synthesis (under high OTR) of the alternative uncoupled respiratory chain of A. vinelandii, which has been reported to be active under high oxygen concentrations  and in diazotrophic conditions, allowing the oxygen sensitive nitrogenase complex to be protected. The results of the present study are in agreement with the functioning of the uncoupled respiratory chain of A. vinelandii , because under the highest OTR condition, some of its protein components, such as NADH oxidoreductase II and cytochrome bd oxidases were clearly induced (Fig. 5), probably as a response to the high oxygen concentration. This alternative respiratory chain has a high affinity for oxygen, but a low net ATP production . These results, together with the analysis of carbon usage (% C-mol), suggest that the carbon source oxidized at higher oxygen consumption would be released as CO2 (Table 2), consume less energy for cell growth, in agreement with our results for qO 2 , YProtein/Glucose (Table 1) and q G (Fig. 2b). On the other hand, this respiratory path is required under conditions of nitrogen fixation. This is in agreement with the requirement imposed by the NH3 limitation observed in the cultures under no oxygen-limited conditions at 14.3 mmol L−1 h−1 OTR (Fig. 2e) where the ammonium concentration was close to zero (Fig. 2e), suggesting that these cultures were probably fixing atmospheric nitrogen. Moreover, under this OTR condition, the transcription level of nifH, which is part of the nifHDK operon that encodes structural components of the nitrogenase system, was induced up to 60-fold with respect to its expression level at the lowest OTR condition. These results support the fact that A. vinelandii responds rapidly to nitrogen depletion under this condition by inducing the expression of nitrogenase; this is also related to the observed increase in the expression of the alternative uncoupled respiratory electron transport system that allows the nitrogenase to be protected under high oxygen transfer rates. As shown in Fig. 5, the relative transcription levels of cydA and ndhII (genes that encode the proteins cytochrome b and NADH dehydrogenase II of the uncoupled respiratory chain) increased up to threefold at 14.3 mmol L−1 h−1 compared to the cultivations developed at low OTR. It is also important to highlight that the enzymatic reduction of molecular nitrogen to ammonia requires high amounts of energy and reducing power (ATP and NADH) . This explanation is supported by the fact that a higher relative metabolic flux through the TCA cycle and CO2 production was determined in the cultivations conducted at 14.3 mmol L−1 h−1 (Fig. 4).
In addition, as shown in Fig. 4, other metabolic pathways affected by the availability of oxygen were: the pyruvate and acetyl-CoA generation points, and the pentose phosphate (PP) pathway. Moreover, these differences were further reflected on the fluxes through the alginate and P3HB biosynthesis pathways. However, we found that the expression of genes encoding central metabolism (idh, akdh and zwf) were similar under all of the conditions evaluated (Fig. 5). These results were inconsistent with the differences in carbon flux towards the central metabolism (PP pathway and TCA cycle). In the cultures at higher oxygen consumption (14.3 mmol L−1 h−1), the relative flux through the PP pathway was nearly 1.8- and 4.8-fold higher than the flux achieved in cultures at 6.2 and 2.4 mmol L−1 h−1, respectively. A similar increase was observed in the metabolic fluxes towards alginate production when this bacterium was grown at high OTR (6.2 and 14.3 mmol L−1 h−1) reaching the highest fluxes at 14.3 mmol L−1 h−1. Previous studies have shown that when A. vinelandii is grown in continuous cultures, without DOT control (oxygen-limited), the alginate concentration can also be increased by increasing qO 2 [29, 35]. When the OTR, and therefore, the qO 2 increased, the carbon fluxes through alginate biosynthesis pathway also increased, from 1.1 at 2.4 mmol L−1 h−1 to 2.6 when the OTR was 14.3 mmol L−1 h−1. These changes at the metabolic flux level are associated with the increase in q Alg values, which increased from 0.45 mmol g−1 h−1 in cultures developed at the lowest OTR, to 0.6 mmol g−1 h−1 when the OTR was 14.3 mmol L−1 h−1. However, this value is the same in the cultures developed at 6.2 mmol L−1 h−1 (Fig. 2c). These results suggest that changes in the cellular respiration (in the range from 6.2 to 14.3 mmol L−1 h−1) did not affect alginate biosynthesis. In contrast, Sabra et al.  reported that the alginate could serve as a barrier to protect the nitrogenase system at high OTR and under diazotrophic conditions. One possible explanation for this behavior could be that, under non-limiting oxygen conditions (14.3 mmol L−1 h−1), the energy demand (in the form of ATP) increased, for nitrogen fixation purposes, as well for the mechanism known as respiratory protection instead of being used to increase alginate biosynthesis . The induction of algD gene transcription is a key point in regulation of the alginate synthesis pathway and is mediated by alginate switching and regulatory genes. In this study, we found that the transcription of the algD gene was higher at 6.2 mmol L−1 h−1 (Fig. 5) and these results were consistent with the higher alginate production under this condition. However, the alginate production was similar in cultures at higher oxygen consumption (14.3 mmol L−1 h−1) when the algD transcription level decreased.
On the other hand, the carbon fluxes towards P3HB production were considerably higher (sixfold) at the lowest OTR condition, in contrast to the flux through this pathway observed for those cultivations conducted at 14.3 mmol L−1 h−1. These findings correlate with the q P3HB , which was 49-fold higher in those cultures conducted with an OTR of 2.4 mmol L−1 h−1 compared to those developed at the highest OTR. Our results are consistent with previous studies showing that P3HB biosynthesis in A. vinelandii is a response to oxygen limitation and carbon excess . As shown in Fig. 5, the phbB transcription levels increased at 2.4 mmol L−1 h−1 in correspondence with the higher production of P3HB observed under that condition. Several authors have proposed that NADPH concentration increases under oxygen-limiting conditions, thereby inhibiting the activities of the tricarboxylic acid (TCA) cycle enzymes, citrate synthase and isocitrate dehydrogenase. Subsequently, the inhibition of these TCA enzymes would cause an increase in the acetyl-CoA available to flux through the P3HB synthesis pathway [2, 37]. On the other hand, as shown schematically in Fig. 6, we observed a slight decrease in gene expression levels of the TCA enzymes ICTDH and aKGDH, in those cultivations conducted at the lowest OTR, as compared to these genes transcription levels at high OTR in agreement with the fluxes discussed.
In the present study, we observed that the internal redox state (reflected by the NADH/NAD+ and NADPH/NADP+ ratios) was affected by changes in the OTR (Fig. 3). Additionally, it was observed that independently of the OTR tested, the NADH/NAD+ ratio was always significantly lower than the NADPH/NADP+ ratio. Under oxygen-limiting conditions (OTR = of 2.4 mmol L−1 h−1), the NADPH/NADP ratio was half of that obtained in the cultures performed at 6.2 mmol L−1 h−1. There is probably more NADPH regenerated at the intermediate oxygen-limiting condition than at extreme oxygen limitation, as we can see higher flux through the PP pathway (an important NADPH generation pathway) was higher (2.5-fold) when the OTR was 6.2 mmol L−1 h−1.
It is important to point out that the NADH cofactor is mainly generated in the central metabolic pathways [Entner–Doudoroff (ED) pathway, TCA cycle and alginate biosynthesis], being consumed mainly in the respiratory chains to produce energy in the form of ATP. Although the values of the NADH/NAD+ ratios achieved under the different conditions changed slightly with OTR (Fig. 3), they were always low. This could be due to fact that the NADH is rapidly consumed by the respiratory chains, and a high proportion of the oxidized form NAD+ has to be preserved in order to favor oxidation of the carbon source; therefore, this intracellular ratio is low. Nevertheless, there were significant differences in the NADH/NAD ratios among the different cultivation conditions evaluated. It is known, that when the oxygen consumption rate increases, the total dinucleotide pool generated during the oxidation of glucose is consumed by the first protein complex of the respiratory chain . Therefore, the intracellular ratio of this cofactor is expected to decrease as the OTR increases. It is important to point out that at an OTR of 6.2 mmol L−1 h−1, the highest NADH/NAD ratio (0.23) was achieved, which is in agreement with an increase in the carbon fluxes through the alginate biosynthesis and TCA cycle, both of which are NADH regeneration pathways (Fig. 4). Therefore, it is possible that under this condition, alginate biosynthesis contributes to increase the NADH availability. On the other hand, when the OTR increased up to 14.3 mmol L−1 h−1, the A. vinelandii cells were not oxygen-limiting, but were exposed to nitrogen limitation conditions. Under these condition, although the carbon fluxes through the TCA cycle, and therefore for NADH production, were higher than under the other conditions evaluated, the electrons of cofactors (NADH and NADPH) could be consumed faster in the respiratory chains of A. vinelandii (the coupled and uncoupled) as part of the protection of the nitrogenase system under high oxygen concentrations , and it was reflected as a decrease in the NADH/NAD+ ratio (Fig. 6).
In A. vinelandii, the main NADPH producing pathways are the ED pathway, the isocitrate dehydrogenase step of the TCA cycle and PPP; these pathways allow a high NADPH/NADP+ ratio to be achieved . The reduced form NADPH has a significant contribution during reductive biosynthetic reactions (fatty acids, amino acids and P3HB), and when A. vinelandii is exposed to high oxygen concentrations and nitrogen fixation, this cofactor could be directed through the uncoupled respiratory chain by NADH dehydrogenase II . The above fact suggests that, the lower NADPH/NADP+ ratio observed in those cultivations conducted at the highest OTR could be due to an increase in NADH dehydrogenase II activity.
The NADH/NAD+ and NADPH/NADP+ cofactor pairs are involved in many biochemical reactions. In the case of A. vinelandii, it has been shown that changes in the oxygen availability have an important effect on the metabolic fluxes. Our results show that at the lowest OTR (2.4 mmol L−1 h−1), the NAD(P)H levels were lower than those achieved at 6.2 mmol L−1 h−1, probably due to a decrease in the carbon fluxes towards the central metabolism (PP pathway an TCA cycle) and to an increase in P3HB biosynthesis. However, the decrease of these reduced cofactor pairs under non-limiting oxygen conditions (14.3 mmol L−1 h−1) compared to the condition at 6.2 mmol L−1 h−1, might be due to the process of nitrogen fixation and the respiratory protection of A. vinelandii that this work shows being induced at high oxygen availability, which was associated with the nitrogen limitation that occurs under this condition.
Bacterial strain and culture medium
Azotobacter vinelandii wild-type strain ATCC 9046 was used. The strain was conserved at 29 °C in Burk’s nitrogen-free salts supplemented with 20 g L−1 of glucose . The composition of the culture medium was as follows (in g L−1): K2HPO4 0.66, KH2PO4 0.16, MOPS 1.42, CaSO4·2H2O 0.05, NaCl 0.2, MgSO4 0.2, NaMoO4·2H2O 0.0029, FeSO4·7H2O 0.027 and 0.8 (NH4)2SO4. The glucose concentration used for chemostat stage cultivation, was of 10 g L−1. The initial pH was adjusted to 7.2 using 2 N NaOH.
Continuous cultures were carried out in a 3 L Applikon bioreactor (Schiedam, Netherlands) with 2 L working volume. The bioreactor, was equipped with two Rushton turbines and aerated at 2 L min−1 (1.0 vvm). The cultivations were conducted at 300, 500 and 700 rpm and the pH was kept constant at 7.2 and it was controlled by the addition of 2 N NaOH. Temperature was maintained at 29 °C. The DOT was measured using a polarographic oxygen probe (Applikon Schiedam, Netherlands). The continuous cultures were operated a dilution rate of 0.08 h−1, this value corresponded to 80% of μmax at 300 rpm. The values of μmax were determined in batch cultivations at different agitation rates using the same culture conditions (pH 7.2, 29 °C, 1.0 vvm). The steady-state condition was achieved after three residence times, when the protein and the glucose concentration remained constant (≤ 5% variation). Samples of cultures (50 mL) were withdrawn from the bioreactor for analytical measurements. All the experiments were conducted by triplicate and the mean value of these results was calculated and reported.
Cellular protein, glucose and ammonium concentration assessments
Microbial growth was evaluated through protein measurements using the Lowry method , with bovine serum albumin as standard, the residual biomass was estimated by subtracting the P3HB content from the total biomass. Glucose concentration was determined by the dinitrosalicylic acid (DNS) reagent . Ammonium concentration was measured by the phenol-hypochlorite method, as described by Kaplan .
Poly(3-hydroxybutyrate) (P3HB) and alginate quantification
P3HB content was quantified by HPLC, after its conversion into crotonic acid. Firstly, the biomass was dried under vacuum at 60 °C. For each sample, 3 mg of biomass was weighed in a 1.5 mL Eppendorf tube; 1 mL of H2SO4 was added, and the sample was heated at 90 °C for 1 h. Subsequently, the sample was cooled at room temperature and diluted with Milli-Q water to concentrations within the range of the calibration curve. The P3HB quantification was performed using an HPLC system with an Aminex HPX-87H column (300 × 7.8 mm) (Bio-Rad, Hercules, CA, USA) at 50 °C and using H2SO4 (7 mM) as the eluent with a flow rate of 0.65 ml min−1. Crotonic acid was quantified using UV absorption at 220 nm. Commercial P3HB (Sigma-Aldrich), treated identically as the samples was used as standard. Alginate concentration was quantified gravimetrically by precipitation with 3 volumes of ice propanol. The precipitate was filtered (0.22 μm Millipore filters) and dried at 85 °C to constant weight .
Determination of the OTR/OUR and the specific oxygen uptake rate
Measurements of the intracellular cofactors concentrations
Amounts of NAD+, NADH, NADP+, and NADPH were quantified by enzymatic methods . These cofactors, were extracted and assayed using the EnzyChrom™ assay kit following the supplier’s instructions (BioAssay Systems, Hayward, CA, USA). 12.6 mg of wet cells from cultures were immediately received in methanol (70% v/v) at − 50 °C for a rapid inactivation of the cellular metabolism . The cell pellet was washed with cold PBS and resuspended with base or acid buffer (BioAssay Systems, Hayward, CA, USA) to extract the reduced or oxidized pyridine nucleotides . The assays utilized glucose-6-phosphate dehydrogenase and lactate dehydrogenase for NAD(H) and NADP(H) quantification respectively, at 565 nm (BioAssay Systems, Hayward, CA, USA).
Flux balance analysis
The components of the metabolic network used for flux calculation were as follows
GAP + DHA
ACCOA + CO2
PYR + CO2
PEP + CO2
P5P + CO2
GAP + S7P
GAP + S7P
GAP + S7P
cE4P + F6P
cE4P + F6P
GAP + S7P
cE4P + P5P
GAP + F6P
GAP + F6P
cE4P + P5P
OAA + ACCOA
CO2 + aKG
SUC + CO2
SUC + GOX
ACCOA + GOX
PYR + GAP
0.0228 * c3PG_bio + 0.0024 * PEP_bio + 0.0012 * cE4P_bio + 7.0E−4 * G6P_bio + 0.0194 * OAA_bio + 0.2768 * ACCOA_bio + 0.0062 * P5P_bio + 0.0512 * PYR_bio + 0.0267 * aKG_bio
Primers used for the quantitative RT-PCR
Forward primers (5′–3′)
Reverse primers (5′–3′)
The specific glucose uptake rate (q G ) and specific alginate and P3HB production rate (q P ) in steady-state conditions were calculated as follows:
Specific glucose uptake rate
Specific alginate or P3HB production rate
Carbon usage determination
The % C-mol was calculated as the ratio of the C-mol of biomass and products and the C-mol of glucose assuming that the carbon source consumed was 100%. The C-mol of biomass results from the biomass elemental composition described as CH2O0.52N0.16 (24.6 g C-mol−1). For alginate and PHB, the monomeric composition was considered: hydroxybutyric acid for PHB and mannuronic acid and acetyl groups for alginate. The following values were used: 24.6 g CDW = 1C-mol for % C-mol biomass; 26.0 g hydroxybutyric acid = 1C-mol for % C-molPHB; 32.3 g mannuronic = 1C-mol for % C-molAlg; 21.5 g acetyl = 1C-mol for % C-molAcetyl. The % C-mol of CO2 was estimated as the difference between the 100% C-mol and the % C-mol used for the biomass and products.
AG performed the experimental work, data analysis and manuscript writing. PF and JA were involved in the analysis of metabolism fluxes and the critical revision of the manuscript. TM and DS were involved in the critical revision of the manuscript. CP, supervised experimental work, data analysis and he was involved in the critical revision of the manuscript. All authors read and approved the final manuscript.
The authors acknowledge the financial support obtained from CONACyT (28535) and DGAPA-UNAM (IT100513). We thank Dra. Celia Flores for technical assistance.
The authors declare that they have no competing interests.
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
- Galindo E, Peña C, Nuñez C, Segura D, Espín G. Molecular and bioengineering strategies to improve alginate and polydydroxyalkanoate production by Azotobacter vinelandii. Microb Cell Fact. 2007;6:7.View ArticleGoogle Scholar
- Peña C, Galindo E, Büchs J. The viscosifying power, degree of acetylation and molecular mass of the alginate produced by Azotobacter vinelandii in shake flasks are determined by the oxygen transfer rate. Process Biochem. 2011;46:230–57.View ArticleGoogle Scholar
- Oelze J. Respiratory protection of nitrogenase in Azotobacter species: is a widely held hypothesis unequivocally supported by experimental evidence. FEMS Microbiol Rev. 2000;24:321–33.View ArticleGoogle Scholar
- Post E, Kleiner D, Oelze J. Whole cell respiration and nitrogenase activities in Azotobacter vinelandii growing in oxygen-controlled continuous culture. Arch Microbiol. 1983;134:68–72.View ArticleGoogle Scholar
- Richard A, Margaritis A. Rheology, oxygen transfer and molecular weight characteristics of poly(glutamic acid) fermentation by Bacillus subtilis. Biotechnol Bioeng. 2003;82:299–305.View ArticleGoogle Scholar
- Anderlei T, Zang W, Papaspyrou M, Büchs J. Online respiration activity measurement (OTR, CTR, RQ) in shake flasks. J Biochem Eng. 2004;17:187–94.View ArticleGoogle Scholar
- Díaz-Barrera A, Peña C, Galindo E. The oxygen transfer rate influences the molecular mass of the alginate produced by Azotobacter vinelandii. Appl Microb Biotechnol. 2007;76:903–10.View ArticleGoogle Scholar
- Díaz-Barrera A, Silva P, Avalos R, Acevedo F. Alginate molecular mass produced by Azotobacter vinelandii in response to changes of the O2 transfer rate in chemostat cultures. Biotechnol Lett. 2009;31(6):903–10.View ArticleGoogle Scholar
- Lozano E, Galindo E, Peña C. The quantity and molecular mass of the alginate produced by Azotobacter vinelandii under oxygen-limited and non oxygen-limited conditions are determined by the maximal oxygen transfer rate (OTRmax). Microb Cell Fact. 2011;10:1–13.View ArticleGoogle Scholar
- Gómez-Pazarín K, Flores C, Castillo T, Buchs J, Galindo E, Peña C. Molecular weight and viscosifying power of alginates produced in Azotobacter vinelandii cultures in shake flasks under low power input. J Chem Technol Biotechnol. 2016;91:1485–92.View ArticleGoogle Scholar
- Berrios-Rivera SJ, Bennett GN, San KY. Metabolic engineering of Escherichia coli: increase of NADH availability by overexpressing an NAD+-dependent formate dehydrogenase. Metab Eng. 2002;4:217–29.View ArticleGoogle Scholar
- Berrios-Rivera SJ, Bennett GN, San KY. The effect of manipulating NADH availability on the redistribution of metabolic fluxes in Escherichia coli chemostat cultures. Metab Eng. 2002;4:230–7.View ArticleGoogle Scholar
- Sanchez AM, Bennett GN, San KY. Effect of different levels of NADH availability on metabolic fluxes of Escherichia coli chemostat cultures in defined medium. J Biotechnol. 2005;117:395–405.View ArticleGoogle Scholar
- De Felipe FL, Kleerebezem M, de Vos WM, Hugenholtz J. Cofactor engineering: a novel approach to metabolic engineering in Lactococcus lactis by controlled expression of NADH oxidase. J Bacteriol. 1998;180:3804–8.Google Scholar
- Fu J, Wang Z, Chen T, Liu W, Shi T, Wang G, Tang YJ, Zhao X. NADH plays the vital role for chiral pure d-(−)-2,3-butanediol production of Bacillus subtilis under limited oxygen conditions. Biotechnol Bioeng. 2014;111:2126–31.View ArticleGoogle Scholar
- Bao T, Zhang X, Zhao X, Rao Z, Yang T, Yang S. Regulation of the NADH pool and NADH/NADPH ratio redistributes acetoin and 2,3-butanediol proportion in Bacillus subtilis. Biotechnol J. 2015;10:1298–306.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
- Henriksen CM, Christensen LH, Nielsen J, Villadsen J. Growth energetics and metabolic fluxes in continuous cultures of Penicillium chrysogenum. J Biotechnol. 1996;45:149–64.View ArticleGoogle Scholar
- Obanye AIC, Hobbs G, Gardner DCJ, Oliver SG. Correlation between carbon flux through the pentose phosphate pathway and production of the antibiotic methylenomycin in Streptomyces coelicolor A3 (2). Microbiology. 1996;142:133–7.View ArticleGoogle Scholar
- Lim SJ, Jung YM, Shin HD, Lee YH. Amplification of the NADPH-related genes zwf and gnd for the oddball biosynthesis of PHB in an E. coli transformant harboring a cloned phbCAB operon. J Biosci Bioeng. 2002;93:543–9.View ArticleGoogle Scholar
- Choi JC, Shin HD, Lee YH. Modulation of 3-hydroxyvalerate molar fraction in poly(3-hydroxybutyrate-3-hydroxyvalerate) using Ralstonia eutropha transformant co-amplifying phbC and NADPH generation-related zwf genes. Enzyme Microb Technol. 2003;32:178–85.View ArticleGoogle Scholar
- Lee JN, Shin HD, Lee YH. Metabolic engineering of pentose phosphate pathway in Ralstonia eutropha for enhanced biosynthesis of poly-β-hydroxybutyrate. Biotechnol Prog. 2003;19:1444–9.View ArticleGoogle Scholar
- Centeno-Leija S, Utrilla J, Flores N. Metabolic and transcriptional response of Escherichia coli with a NADP+-dependent glyceraldehyde 3-phosphate dehydrogenase from Streptococcus mutans. Antonie Van Leeuwenhoek. 2013;104(6):913–24.View ArticleGoogle Scholar
- Tatnell P, Russell N, Gacesa P. GDP-mannose dehydrogenase in the key regulatory enzyme in alginate biosynthesis in Pseudomonas aeruginosa: evidence form metabolite studies. Microbiology. 1994;140:1745–54.View ArticleGoogle Scholar
- Snook CF, Tipton PA, Beamer LJ. Crystal structure of GDP-mannose dehydrogenase: a key enzyme of alginate biosynthesis in Pseudomonas aeruginosa. Biochemistry. 2003;42:4658–68.View ArticleGoogle Scholar
- Senior PJ, Beech GA, Ritchie GAF, Dawes EA. The role of oxygen limitation in the formation of poly-3-hydroxybutyrate during batch and continuous culture of Azotobacter beijerinckii. Biochem J. 1972;128:1193–201.View ArticleGoogle Scholar
- Castillo MT, Heinzle E, Peifer S, Schneider K, Peña C. Oxygen supply strongly influences metabolic fluxes, the production of poly(3-hydroxybutyrate) and alginate, and the degree of acetylation of alginate in Azotobacter vinelandii. Process Biochem. 2013;48:995–1003.View ArticleGoogle Scholar
- Maier U, Losen M, Büchs J. Advances in understanding and modeling the gas–liquid mass transfer in shake flasks. Biochem Eng J. 2004;17:155–67.View ArticleGoogle Scholar
- Díaz-Barrera A, Soto E, Altamirano C. Alginate production and alg8 gene expression by Azotobacter vinelandii in continuous cultures. J Ind Microbiol Biotechnol. 2012;39:613–21.View ArticleGoogle Scholar
- Noguez R, Segura D, Moreno S, Hernández A, Juárez K, Espín G. Enzyme INtr, NPr and IIANtr are involved in regulation of the poly-β-hydroxybutyrate biosynthetic genes in Azotobacter vinelandii. J Mol Microbiol Biotechnol. 2008;15:244–54.View ArticleGoogle Scholar
- Castillo MT, Galindo E, Peña C. The acetylation degree of alginates in Azotobacter vinelandii ATCC9046 is determined by the dissolved oxygen and the specific growth rate: studies in glucose-limited chemostat cultivations. J Ind Microbiol Biotechnol. 2013;40:715–23.View ArticleGoogle Scholar
- Bertsova YV, Bogachev AV, Skulachev VP. Noncoupled NADH: ubiquinone oxidoreductase of Azotobacter vinelandii is required for diazotrophic growth at high oxygen concentrations. J Bacteriol. 2001;183:6869–74.View ArticleGoogle Scholar
- Linkerhägner K, Oelze J. Nitrogenase activity and regeneration of the cellular ATP pool in Azotobacter vinelandii adapted to different oxygen concentrations. J Bacteriol. 1996;179:1362–7.View ArticleGoogle Scholar
- Sabra W, Zeng AP, Deckwer WD. Bacterial alginate: physiology, product quality and process aspects. Appl Microbiol Biotechnol. 2001;56:315–25.View ArticleGoogle Scholar
- Díaz-Barrera A, Maturana N, Pacheco-Leyva I, Martínez I, Altamirano C. Different responses in the expression of alginases, alginate polymerase and acetylation genes during alginate production by Azotobacter vinelandii under oxygen-controlled conditions. J Ind Microbiol Biotechnol. 2017;1:1–11.Google Scholar
- Inomura K, Bragg J, Follows MJ. A quantitative analysis of the direct and indirect costs of nitrogen fixation: a model based on Azotobacter vinelandii. ISME J. 2017;11:166–75.View ArticleGoogle Scholar
- Sabra W, Zeng AP, Sabry S, Omar S, Deckwer WD. Effect of phosphate and oxygen concentrations on alginate production and stoichiometry of metabolism of Azotobacter vinelandii under microaerobic conditions. Appl Microbiol Biotechnol. 1999;52:773–80.View ArticleGoogle Scholar
- Bertsova YV, Bogachev AV, Skulachev VP. Two NADH: ubiquinone oxidoreductases of Azotobacter vinelandii and their role in respiratory protection. Biochim Biophys Acta. 1998;1363:125–33.View ArticleGoogle Scholar
- Spaans SK, Weusthuis RA, van der Oost J, Kengen SW. NADPH-generating systems in bacteria and archaea. Front Microbiol. 2015;6:742.View ArticleGoogle Scholar
- Lowry OH, Rosebrough NJ, Farr AL, Randall RJ. Protein measurement with the Folin phenol reagent. J Biol Chem. 1951;193:265–75.Google Scholar
- Miller G. Use of dinitrosalicylic acid reagent for determination of reducing sugar. Anal Chem. 1959;31:426–8.View ArticleGoogle Scholar
- Kaplan A. The determination of urea, ammonia and urease. Methods Biochem Anal. 1969;17:311–24.Google Scholar
- Peña C, Trujillo-Roldán M, Galindo E. Influence of dissolved oxygen tension and agitation speed on alginate production and its molecular weight in cultures of Azotobacter vinelandii. Enzyme Microb Technol. 2000;27:390–8.View ArticleGoogle Scholar
- Zeng AP, Byun T, Posten C, Deckwer WD. Use of the respiratory quotient as a control parameter for optimum oxygen supply and scale-up of 2,3-butanediol production under microaerobic conditions. Biotechnol Bioeng. 1994;19:1107–14.View ArticleGoogle Scholar
- Stanbury PF, Whitaker A, Hall S. Principles of fermentation technology. Oxford, United Kingdom: Butterworth Heinemann; 2003. p. 360.Google Scholar
- Bergmeyer H. Methods of enzymatic analysis, vol. 2. 3rd ed. Weinheim: VCH Verlagsgesellschaft mbH; 1985. p. 1047–64.Google Scholar
- Rocha I, Maia P, Evangelista P, Vilaça P, Soares S. OptFlux: an open-source software platform for in silico metabolic engineering. BMC Syst Biol. 2010;4:45–56.View ArticleGoogle Scholar
- Flores C, Moreno S, Espín G, Peña C, Galindo E. Expression of alginases and alginate polymerases in response to oxygen, and their relationship with the alginate molecular weight in Azotobacter vinelandii. Enzyme Microbiol Technol. 2013;53:85–91.View ArticleGoogle Scholar
- Barry T, Geary S, Hannify S, MacGearailt C, Shalloo M, Heery D. Rapid minipreparations of total RNA from bacteria. Nucleic Acids Res. 1992;20:4940.View ArticleGoogle Scholar
- Livak KJ, Shmittgen TD. Analysis of relative gene expression data using realtime quantitative PCR and the 2-CT method. Methods. 2001;25:402–8.View ArticleGoogle Scholar