Original Research ARTICLE
Physiological Response of Corynebacterium glutamicum to Increasingly Nutrient-Rich Growth Conditions
- Institute of Biochemical Engineering, University of Stuttgart, Stuttgart, Germany
To ensure economic competitiveness, bioprocesses should achieve maximum productivities enabled by high growth rates (μ) and equally high substrate consumption rates (qS) as a prerequisite of sufficient carbon-to-product conversion. Both traits were investigated and improved via bioprocess engineering approaches studying the industrial work horse Corynebacterium glutamicum. Standard minimal medium CGXII with glucose as sole carbon source was supplemented with complex brain-heart-infusion (BHI) or amino acid (AA) cocktails. Maximum μ of 0.67 h-1 was exclusively observed in 37 g BHI L-1 whereas only minor growth stimulation was found after AA supplementation (μ = 0.468 h-1). Increasing glucose consumption rates (qGlc) were solely observed in certain dosages of BHI (1–10 g L-1), while 37 g BHI L-1 and AA addition revealed qGlc below the reference experiments. Moreover, BHI supplementation revealed Monod-type saturation kinetics of μ (KBHI = 2.73 g BHI L-1) referring to the preference of non-AAs as key boosting nutrients. ATP-demands under reference, 1 g BHI L-1, and AA conditions were nearly constant but halved in BHI concentrations above 5 g L-1 reflecting the energetic advantage of consuming complex nutrient components in addition to “simple” building blocks such as AAs. Furthermore, C. glutamicum revealed maximum biomass per carbon yields of about 18 gCDW C-mol-1 irrespective of the medium. In AA supplementation experiments, simultaneous uptake of 17 AAs was observed, maximum individual consumption rates determined, and L-asparagine and L-glutamine were distinguished as compounds with the highest consumption rates. Employment of the expanded stoichiometric model iMG481 successfully reproduced experimental results and revealed the importance of C. glutamicum’s transaminase network to compensate needs of limiting AA supply. Model-based sensitivity studies attributed the highest impact on μ to AAs with high ATP and NADPH demands such as L-tryptophan or L-phenylalanine.
Introduction and Motivation
Corynebacterium glutamicum is a Gram-positive, non-sporulating, facultative anaerobic prokaryote (Liebl, 2006; Nishimura et al., 2007), generally recognized as safe (GRAS), and best known for its production of L-glutamate (Kinoshita et al., 1957). Revelation of the C. glutamicum ATCC 13032 wild type (WT) genome by Ikeda and Nakagawa (2003) and Kalinowski et al. (2003) enabled application of metabolic engineering techniques and the construction of genome scale metabolic models (Kjeldsen and Nielsen, 2009; Shinfuku et al., 2009). Thereby, targeted modification of the strain’s genome, often predicted by said in silico models, enhanced C. glutamicum’s abilities to produce lactate, succinate, and pyruvate (Dominguez et al., 1993; Okino et al., 2005; Wendisch et al., 2006; Wieschalka et al., 2012), to form non-natural compounds such as diamins (Kind et al., 2010) or heterologous proteins (Scheele et al., 2013), and to emerge in the field of bio-based chemicals, materials, and fuels (Becker and Wittmann, 2012). Besides, the advantageous physiological properties, summarized by Vertès et al. (2012), further empower the microorganism to endure harsh production conditions as simulated and investigated in several scale-down studies mimicking production-like environments (Buchholz et al., 2014b; Käß et al., 2014a,b). Accordingly, the strain is well accepted as industrial platform and predominately employed for the production of the amino acids (AAs) L-glutamate, L-lysine, and L-valine (Hermann, 2003) serving demands of food, feed, cosmetic, and pharmaceutical markets (Takors et al., 2007). More than 2 million tons per year of L-glutamate (Becker and Wittmann, 2012) are likely to be produced by the strain in large scale bioreactors reaching volumes of 500 m3 and more.
Regarding industrial AA production, C. glutamicum competes with Escherichia coli. Whereas the first offers advantages in robustness also including minimized phage contamination vulnerability, E. coli naturally possesses transhydrogenases which allow quick equilibration of NADH/NAD+ and NADPH/NADP+ pools (Sauer et al., 2004). The latter is of particular importance to engineer top producers of NADPH demanding compounds such as L-lysine requiring 4 mol NADPH in the synthesis of 1 mol L-lysine (Marx et al., 1997). Besides, C. glutamicum reveals only second-best properties in terms of growth kinetics: maximum growth rates of about 0.58 h-1 have been measured for C. glutamicum strain ATCC 17965 in batch and continuous cultivations with minimal medium containing glucose (Cocaign-Bousquet et al., 1996; Gourdon et al., 2003) and in turbidostatic processes using C. glutamicum ATCC 13032 (Bäumchen et al., 2007). Grünberger et al. (2013) achieved even slightly faster growth rates of 0.60–0.64 h-1 in highly diluted microliter-scale conditions using CGXII minimal medium with glucose as sole carbon source. However, the same medium only shows approximately 0.4 h-1 in conventional bioreactor setups (Blombach et al., 2013; Grünberger et al., 2013). On the contrary, E. coli achieves maximum growth rates of up to 1.04 h-1 in glucose supplemented minimal medium and 1.73 h-1 in complex broth (Cox, 2004). In general, low growth rates inherently limit the productivity of the production process: first by achieving the desired biomass concentration, second by limiting product formation in growth coupled kinetics. By analogy, low biomass specific glucose uptake rates represent an intrinsic disadvantage of the producing cell as they may cause low biomass specific productivities, too.
Consequently, to increase the attraction of C. glutamicum for industrial application, improvement of fundamental kinetics would be highly appreciated. Strain and process engineering measures could be applied. The first, however, is beyond the scope of this contribution. Instead, properly applied medium supplements may improve productivities of existing production processes without necessary strain engineering. Accordingly, our contribution investigates process engineering means, in particular medium supplementation, for increasing growth and glucose consumption rates of C. glutamicum WT as a prerequisite of improving volumetric productivities. Industrial bioprocess development aims at reducing complex media such as yeast extract or corn steep liquor to reduce medium costs and variabilities. Of course, key performance criteria of the fermentation should still be achieved. Accordingly, the identification of so-called “goodies” of medium constituents is welcomed. However, the general physiological response of C. glutamicum to complex medium (compounds) was not studied yet. To date, supplementation investigations with the strain mostly aimed at media optimization (Coello et al., 2002; Jeon et al., 2013) or usage of alternative carbon and/or nitrogen sources (Lee et al., 2014; Wendisch et al., 2016; Lange et al., 2017). For E. coli, in contrast, profound knowledge exists based on physiological investigations in Luria–Bertani (LB) broth (Sezonov et al., 2007) and due to comparable transcriptional expression analyses of late-exponential E. coli cultures grown in minimal and LB medium (Tao et al., 1999). This approach was further expanded in a broad gene expression study analyzing the different growth phases of the enterobacterium grown in the same complex medium (Baev et al., 2006a,b,c).
Thus, we performed a two-staged analysis to study how supplementation of standard growth minimal medium CGXII impacts the physiology of C. glutamicum WT in bioreactor batch cultivations. In the first part, complex brain-heart-infusion (BHI) was added to CGXII to investigate and identify the maximum metabolic capacity of the cell regarding growth. As a step toward the identification of optimum synthetic growth medium, AA supplementation was studied in the second part which also allows the analysis of theoretical impacts via stoichiometric modeling. Until now, AA consumption was only investigated in transporter characterization studies (e.g., Burkovski et al., 1996; Trötschel et al., 2008; Shang et al., 2013) or under anaerobic conditions (Michel et al., 2015). In summary, the manuscript focuses on process engineering approaches to accelerate growth of C. glutamicum.
Materials and Methods
Corynebacterium glutamicum ATCC 13032 WT cells from working cell bank glycerol stocks were spread on 2 x tryptone-yeast extract (2 x TY) medium (Sambrook, 2001) agar plates and incubated for 48–60 h at 30°C. Colonies were used to inoculate 5 mL of 2 x TY medium in glass reaction tubes which were incubated for 8 h at 30°C and 120 rpm on a bench-top rotary shaker (Infors HT, Bottmingen, Switzerland). The content of one glass reaction tube was used to inoculate 50 mL modified CGXII minimal medium (Buchholz et al., 2014a) supplemented with 4% (w/v) glucose. The shaking flask precultures were incubated overnight at 30°C and 120 rpm.
Cultivation in Bioreactors
Bioreactor cultivations were carried out in two different bioreactor systems. Experiments with BHI-supplementation were conducted in biological triplicates in a parallel triplex glass bioreactor system (HWS Labortechnik, Mainz, Germany) with a working volume of 250 mL. CGXII minimal medium (Buchholz et al., 2014a) was supplemented with 2% (w/v) glucose (reference) and optionally 1, 5, 10, or 37 g BHI (1, 5, 10, 37 BHI) medium per liter. Culture broth from preculture shaking flasks was harvested, centrifuged at 4°C for 10 min at 4000 × g (5430 R, Eppendorf, Hamburg, Germany) and resuspendend in sterile 0.9% (w/v) NaCl-solution. The differently supplemented CGXII media were inoculated with the washed preculture installing initial optical densities (OD600) of 1 in 200 mL reaction volume. The parallel reactor system was equipped with dissolved oxygen (pO2)- and pH-probes (Mettler Toledo GmbH, Albstadt, Germany). The temperature was kept constant at 30°C, pO2 was kept >30% by increasing the impeller speed by 50 rpm (initial speed: 200 rpm) and the aeration rate by 20 mL min-1 (initial rate: 20 mL min-1) in a cascadic fashion. pH was kept constant by adding 25% (v/v) NH4OH. Antifoam agent (Struktol®J 647, Schill + Seilacher, Hamburg, Germany) was added manually when needed.
For AA supplementation experiments, biological triplicates were conducted in a 3 L steel-tank bioreactor (Bioengineering, Wald, Switzerland) ensuring high cell densities and enabling off-gas analysis. CGXII minimal cultivation media supplemented with 4% (w/v) glucose and optionally 5 mmol of the following AAs per liter were used: L-alanine (Ala), L-asparagine (Asn), L-aspartate (Asp), L-arginine (Arg), L-glutamine (Gln), L-glutamate (Glu), glycine (Gly), L-histidine (His), L-isoleucine (Ile), L-leucine (Leu), L-lysine (Lys), L-methionine (Met), L-phenylalanine (Phe), L-proline (Pro), L-serine (Ser), L-threonine (Thr), L-tryptophan (Trp), and L-valine (Val); 100 mL washed shaking flask preculture was used to inoculate supplemented CGXII media to give a starting OD600 of 3 and a total of 1.5 L reaction volume. The bioreactor was equipped with a six-blade Rushton impeller, a pO2-, and pH-probe and operated at a total pressure of 1.5 bar at 30°C. pO2 was kept > 30% by increasing the impeller speed by 50 rpm (initial speed: 250 rpm) and the aeration rate by 0.3 L min-1 (initial rate: 0.3 L min-1) in a cascadic fashion. pH was kept constant by adding 25% (v/v) NH4OH. Exhaust gas analysis of molar O2- and CO2-fractions was performed with non-dispersive (photometric) infrared gas analyzers (BCPO2 and BCPCO2, BlueSens, Herten, Germany). Antifoam agent was added manually when needed.
During the BHI-supplementation experiments, biosuspension samples were withdrawn hourly to measure the biomass concentration and to obtain cell-free samples. Determination of biomass was performed in technical replicates by following the optical density (OD600) of the cell suspension with a spectrophotometer (DR 3900, Dr. Lange, Berlin, Germany) at 600 nm. During the AA supplementation experiments, biomass was additionally measured gravimetrically: 1 mL biosuspension was centrifuged at 20,173 × g and 4°C for 2 min (5430 R, Eppendorf, Hamburg, Germany), washed twice with deionized water, and dried at 105°C for at least 24 h in a convection oven (Heraeus, Hanau, Germany). Correlating the OD600 and cell-dry-weight (CDW) measurements resulted in a correlation factor of 0.27 g biomass L-1. Cell-free samples were obtained by applying syringe filters with 0.2 μm pore size (Rotilabo®, Carl Roth, Karlsruhe, Germany) and storing the filtrates at -20°C until glucose-, total organic carbon- (TOC), and AA-quantification was performed. Glucose concentration was determined with an enzymatic assay following the instructions of the manufacturer (R-biopharma, Darmstadt, Germany). Following the protocols given by Buchholz et al. (2014a), AA concentrations were determined by means of high-pressure liquid chromatography, and total inorganic carbon (TIC)-, and TOC-amounts in biosuspension samples, cell-free filtrates, and biomass samples were measured with total carbon (TC) analysis. Therefore, a TC analyzer (Multi N/C 2100s, Analytik Jena, Jena, Germany) was employed to analyze samples regarding (i) the TIC content by acidification with 10% ortho-phosphoric acid to free carbonate-derived CO2 and (ii) the TC content by total combustion of organic and inorganic carbon species to CO2 at 750°C in a furnace system. TOC amounts were obtained with the differential detection mode of the device (TOC = TC – TIC). In the AA supplementation experiments, harvesting frequency was doubled after 4 h, and additional sampling enabled measuring of TOC and TIC values to check CO2 exhaust-gas analysis as proposed by Buchholz et al. (2014a).
Determination of Kinetic Parameters
Exponential Growth Rate
Exponential growth rates (μ) were determined using linear regression of logarithmic biomass concentration curves over the process time as described by Grünberger et al. (2013) for batch cultivations.
Consumption Rates and Yields
Biomass-specific glucose, AA, and TOC consumption rates (qGlc, qAA, qC) were determined by dividing the exponential growth rate (μ) by the biomass substrate yield (YXS) previously calculated by linear regression of substrate (S, representing glucose or carbon) concentrations curves over the biomass (X) concentration.
Respiratory Rates in Amino Acid Supplementation Experiments
Biomass-specific respiratory rates (qO2, qCO2) from AA supplementation experiments were determined by dividing the volumetric oxygen consumption or carbon dioxide emission rate by the biomass concentration. qCO2 deduced from the off-gas analysis was corrected by accounting for dissolved carbon species in the liquid volume (TIC) in the early growth phase (Buchholz et al., 2014a).
Estimation of Respiratory Rates and ATP-Biomass Yield in BHI Supplementation Experiments
To estimate the produced CO2 during the BHI-supplementation experiments, molar (n) carbon (C) balancing was performed assuming that the consumed molar amount of TOC (nC(TOC)) in the time frame of the exponential growth phase was either converted to biomass or to CO2 neglecting TIC species accumulating in the fermentation broth in the early growth phase
Determination of TOC in biomass and cell-free samples by TC-analysis allowed calculation of the theoretically produced amount of CO2. Thereof, the CO2-biomass yield was obtained and converted to qCO2,est by multiplication with μ. The corresponding oxygen consumption rate qO2,est was deduced from qCO2,est by assuming a respiratory quotient (RQ) of 1, which was supported by the experimental rates of the reference and AA-supplementation experiments (c.f. Tables 1, 2) in the steel-tank bioreactor
TABLE 1. Summary of kinetic parameters determined and estimated (index est) for C. glutamicum ATCC 13032 grown in modified CGXII minimal medium supplemented with 2% (w/v) glucose (Glc, reference) and additionally 1, 5, 10, or 37 g BHI L-1 (abbr. to 1, 5, 10, 37 BHI).
TABLE 2. Summary of kinetic parameters determined for C. glutamicum grown in modified CGXII minimal medium supplemented with 4% (w/v) glucose (Glc) or additionally 5 mmol of 18 amino acids (AAs) per L (abbr. to 5 AA).
ATP production and consumption were assumed to be equilibrated. The first was deduced from qO2,est assuming a P/O-ratio of 1.2 as a conservative estimate. Accordingly, the biomass-specific ATP production rate (qATP) could be calculated as
In turn, division of qATP by μ resulted in the ATP-biomass yield YATPX. Standard deviations (SDs) for C-balancing are based on experimental results (biological triplicates) of carbon measurement of biomass and TOC, while SDs in all further calculation steps (cf. Eqs 3 and 4) were calculated by Gaussian error propagation.
Fitting of Monod-Type Parameters
The exponential growth rate μ of a microorganism can be calculated based on Monod-type kinetics (Monod, 1949)
where cS is the substrate concentration of the only growth limiting substrate, KS is the half-saturation constant (which equals the substrate concentration leading to the half maximum growth rate), and μmax is the maximum growth rate observed for the used substrate. In this study, Monod-kinetics were applied to examine the effect of different BHI concentrations (1–37 g L-1) on the corresponding growth rate (μBHI) by calculating the difference (Δμ) between growth rates obtained from experiments without BHI [μRef, 2% (w/v) glucose as only carbon source] and μBHI
Inserting Eq. (6) into Eq. (5) depicted the influence of BHI-concentration (cBHI) on the resulting growth difference
To determine the unknown parameters Δμmax, that is the difference between the reference growth rate and the maximum growth rate observed under 37 g BHI L-1, and the saturation constant of BHI (KBHI), which is the BHI-concentration leading to half-maximum growth, the Curve Fitting Toolbox of Matlab (MATLAB Release 2014a, The MathWorks, Inc., Natick, MA, United States) was used. Thereby, non-linear least squares fitting to the experimental data based on the Monod-equation (7) were performed making use of the Levenberg–Marquadt algorithm and setting the confidence level to 90%.
Simulations With Stoichiometric Network Model iMG481
Curation of Model iEZ475 and Expansion to iMG481
The stoichiometric model curated and expanded in this work is based on model iEZ475 of Zelle et al. (2015). Recently, another C. glutamicum reconstruction was crafted by Zhang et al. (2017), but was not considered in this investigation. After examination of AA biosynthesis and transporter reactions in iEZ475, the following alterations and expansions were performed according to published knowledge on the topics. iMG481 is attached in the Supplementary Material as sbml- and Excel-file.
Regarding AA biosynthesis, the main route for Ala synthesis via alanine-transferase AlaT (Marienhagen and Eggeling, 2008) using pyruvate (Pyr) and Glu as precursors was added as reaction “alaT” yielding Ala and oxo-glutarate as products. Recently, Buerger et al. (2016) annotated gene glsK (cg2728) encoding the only known glutaminase in C. glutamicum. The corresponding hydrolyzation of Gln to Glu and NH3 was added as irreversible reaction “glsK.” Concerning the transport of AAs into the in silico cell, the direction of already implemented reactions “SER_t_NA” and “THR_t_NA,” formerly used for export of the AAs, was reversed to enable Ser and Thr import with Na+-symport (Palmieri et al., 1996) since the export of both AAs is independent of sodium ions (Palmieri et al., 1996). The Na+-flux, formerly in antiport with Gln, was reversed in reaction “GLN_t_NA” according to Siewe et al. (1995). Lys import can be facilitated by antiport with Ala, Ile, and Val (Bröer and Krämer, 1990) and reaction “LYS_t_LEU” which falsely uses Leu was corrected to “LYS_t_ILE” transporting Ile.
Incorporating C. glutamicum-specific information on AA transporters, several reactions were added: (i) Asn transport, hitherto only possible via diffusion, was expanded by a Na+-symporter (“ASN_t_NA”) making use of (supposedly) the same transporter as for Gln (Siewe et al., 1995); (ii) Asp uptake was expanded by ATP-dependency (“ASP_t_ATP”) mirroring the same transporter type as for Glu uptake (Krämer et al., 1990; Marin and Krämer, 2007); and (iii) Met transport facilitated by Na+-symport (Trötschel et al., 2008) was added (“MET_t_NA”) to the existing ATP-dependent system. No Arg-transporter was implemented in model iEZ475 because Arg transport or a specific Arg transporter is not yet characterized or validated for C. glutamicum. However, Arg consumption was observed during the experiments (cf. section “Amino Acid Consumption and Amino Acid-Biomass Yield”), and therefore an ATP-dependent importer (“ARG_t_ATP”) was added based on experimental evidence for E. coli (Gong et al., 2003).
Flux Balance Analysis
Growth rates, AA consumption rates, and metabolic flux distributions of the in silico cell were calculated using the optimizeCbModel function of “The COBRA Toolbox” (Heirendt et al., 2017) with solver set “glpk” and Matlab software. Biomass production was set as objective function, while experimentally determined glucose consumption rates (qGlc), AA consumption rates (qAA), and respiratory rates (qO2, qCO2) were used as constraints. To simulate abundant AA supply, AA uptake was unconstrained by setting the upper flux bound to 1000.
To study the impact of AA supplementation on the growth rate of the in silico cell, sensitivity analyses were performed. Thus, FBAs were conducted where the glucose uptake rate was constrained to 3.5 mmol g-1 h-1, representing the operation mode of exponential growth (cf. section “Amino Acid Consumption and Amino Acid-Biomass Yield”), and that of one single AA (qAA) to 3.5 nmol g-1 h-1 with the objective to maximize growth (μSim). FBA with only glucose uptake and without AA uptake yielded the reference growth rate μRef. The sensitivity of each AA was thus calculated as:
BHI Supplementation Experiments
The impact of BHI medium on growth behavior and kinetic parameters of C. glutamicum ATCC 13032 was studied in triplicate bioreactor cultivations. Therefore, 1, 5, 10, or 37 g BHI L-1 (1, 5, 10, 37 BHI) was added to modified CGXII minimal medium [2% (w/v) glucose] and the determined kinetic properties were compared to reference cultivations with glucose as sole carbon source. Calculated exponential growth rates μ obtained from each condition showed an increase of μ with increasing BHI amounts (μ = 0.44–0.67 h-1; Table 1). Under reference conditions, μ = 0.32 ± 0.02 h-1 was observed and a maximum rate of 0.67 ± 0.03 h-1 in 37 BHI. Pairwise t-tests between conditions with the most similar μ confirmed significant differences in growth mirroring the different BHI supplements. Besides, BHI stimulated the glucose consumption under 1–10 BHI conditions, revealing qGlc of 0.0277 ± 0.0016 C-mol g-1 h-1 in 1 BHI, which was 29% higher than the reference. However, a decline of qGlc with increasing growth rate promoted by increasing BHI amounts was observed (Figure 1, open circles), finally reaching a value of 0.0191 ± 0.0019 C-mol g-1 h-1 in 37 BHI that was 11% smaller than in reference conditions. In contrast to glucose consumption, the TC consumption rate qC increased proportionally with μ as illustrated in Figure 1 (filled black circles). Thereby, a constant biomass-carbon yield (YXC) of 18.47 ± 0.55 g C-mol-1 could be calculated by linear regression of biomass production versus carbon consumption.
FIGURE 1. Illustration of biomass-specific glucose (qGlc) and total organic carbon (qC) consumption rates obtained in experiments with C. glutamicum WT grown in modified CGXII minimal medium with 2% (w/v) glucose (Glc) and brain-heart-infusion (BHI)-supplemented media (1, 5, 10, or 37 g BHI L-1) over the corresponding exponential growth rate (μ). A proportional increase of qC in BHI-supplemented conditions with μ indicated a constant biomass-organic carbon yield (YXC, cf. Table 1). Values are deduced from three biological replicates ± standard deviation.
The supplementation of CGXII minimal medium with different amounts of BHI induced a diversification of the growth rate of C. glutamicum which motivated us to investigate the influence of BHI concentration on μ. Therefore, differences in growth (Δμ) of BHI-supplemented and non-supplemented conditions (reference) were calculated (Eq. 6) and plotted against related BHI concentrations (1–37 g BHI L-1), illustrated in Figure 2 (filled black circles). A Monod-like saturation curve was observed with a steep increase of Δμ for low BHI additions (1–5 g BHI L-1) leveling out at BHI-concentrations above 10 g L-1 and approximating the maximal growth rate difference (Δμmax). Therefore, non-linear least squares fitting of the experimental data was performed, using the Monod-type equation (7) with 90% confidence level of the fit. Δμmax and the saturation constant of BHI (KBHI), that is the BHI concentration leading to half-Δμmax, were thus estimated with 0.359 h-1 and 2.73 g BHI L-1, respectively. Plugging both values into Eq. 7 results in the estimation curve (blue solid line) pictured in Figure 2 with the dotted curves above and below as a measure for the confidence level.
FIGURE 2. Exponential growth rate differences (Δμ) were calculated for C. glutamicum WT grown in modified CGXII minimal medium supplemented with 2% (w/v) glucose and with additional 1, 5, 10, or 37 g brain-heart-infusion (BHI) L-1 medium (three biological replicates ± standard deviation). Illustration of Δμ over the corresponding amount of BHI revealed a Monod-like saturation curve. The half-velocity constant for BHI (KBHI) and the maximum growth rate difference (Δμmax) given in the figure were deduced from non-linear least-squares fitting (blue solid line) with a confidence level of 90% (blue dotted lines).
Estimated ATP-Biomass Yield
The ATP consumption of C. glutamicum cultivated in supplemented CGXII media was estimated to evaluate the energy demand of the cell under minimal and increasingly rich conditions. Since no intracellular ATP pools were measured, the theoretical amount of ATP produced and consumed by the strain was estimated as given in Eqs 2–4. Therefore, carbon balancing was performed to estimate the produced molar amount of CO2 which was converted to the (estimated) carbon emission rate qCO2,est (Table 1). Assuming a RQ of 1 led to the corresponding oxygen consumption rate qO2,est, which was then used to calculated the ATP consumption rate by assuming a P/O ratio of 1.2. The resulting ATP-biomass yield YATPX is plotted against growth rates in Figure 3. Under reference conditions (glucose as sole carbon source), the highest ATP-yield of 26.03 ± 2.08 mmol g-1 was observed, whereas the yield with 1 g BHI L-1 was in the same range (23.9 ± 3.04 mmol g-1). However, elevation of BHI-amounts above 1 g L-1 lead to a new ATP-level of around 14 mmol g-1 considering the SDs of the respective conditions (5, 10, and 37 BHI). This equals a drop of the ATP-biomass yield by ca. 46% in comparison to reference conditions.
FIGURE 3. Illustration of estimated ATP-biomass yields (YATPX) over exponential growth rates (μ) from BHI-supplementation experiments with C. glutamicum WT grown in modified CGXII minimal medium with 2% (w/v) glucose (Glc) and brain-heart-infusion (BHI)-supplemented media (1, 5, 10, or 37 g BHI L-1). YATPX-estimation was based on carbon balancing to calculate the theoretical respiratory rates and assuming a P/O-ratio of 1.2 (cf. Eqs 2–4 and Table 1). Values are deduced from three biological replicates and error propagation was used to determine corresponding standard deviations pictured as error bars in the figure.
Amino Acid Supplementation Experiments
To unravel whether the stimulative effect of BHI on growth rate and TC consumption of C. glutamicum WT was due to AAs, which are components of BHI in free or peptide-bound form, AA supplementation fermentations were performed. Thus, 5 mM of all essential AAs (except for L-cysteine and L-tyrosine) were added to modified CGXII minimal medium [4% (w/v) glucose] to comparably study the kinetic properties of the strain cultivated in the “simulated” complex medium. Impacts of elevated biomass concentrations were studied in the pressurized bioreactor to prevent oxygen limitation. To ensure comparability, reference cultivations with glucose as only carbon source were repeated in this bioreactor. Deduced kinetic parameters (Table 2) showed good conformity with corresponding values obtained in the previous system (Table 1). Accordingly, all results can be discussed irrespective of the reactor system. Addition of AAs to CGXII minimal medium increased μ by about 40% from 0.337 ± 0.033 (reference) to 0.468 ± 0.003 h-1 (Table 2) which was within the range of 1–5 BHI supplementation (μ = 0.44–0.54 h-1; Table 1). Interestingly, the maximum observed growth rate (0.67 ± 0.03 h-1) in 37 BHI was not reached by AA supplementation. As previously observed in 37 BHI, qGlc determined under AA-supplemented conditions was within the range of the reference considering the deviance (AA: 0.0205 ± 0.0006 C-mol g-1 h-1; reference: 0.0213 ± 0.0020 C-mol g-1 h-1). In contrast, TC consumption was about 21% higher (0.0257 ± 0.0006 C-mol g-1 h-1) making it comparable to the values observed in 1 BHI (0.0255 ± 0.0021 C-mol g-1 h-1). Determination of the biomass-carbon yield YXC in AA supplementation revealed a value of 18.00 ± 0.52 g C-mol-1 which fits to the linear regression of YXC from the BHI-experiments (18.47 ± 0.55 g C-mol-1). Calculation of respiratory rates deduced from exhaust gas analysis displayed an increase of about 45% compared to the reference (qO2 = 6.06 ± 0.48 mmol g-1 h-1; qCO2 = 6.17 ± 0.65 mmol g-1 h-1) which equals the increase of growth (about 40%).
Amino Acid Consumption and Amino Acid-Biomass Yield
Regarding the consumption of the 18 added AAs in the AA supplementation experiments, all AAs except for Asp were consumed (cf. Supplementary Figure S1) by C. glutamicum WT. Simultaneous decrease of AA concentrations in cell-free filtrates was observed but for Asn which was consumed only after Gln was nearly depleted (3–4 h, Supplementary Figure S1). Related AA consumption rates were calculated for each AA and are summarized in Table 3. The highest rate of 0.694 ± 0.005 mmol g-1 h-1 was obtained for Gln, whereas Trp was consumed with the smallest rate (0.022 ± 0.001 mmol g-1 h-1). Next, AA-biomass yields (YAA) were determined and compared to three biomass AA-compositions of different C. glutamicum species (Table 4): Micrococcus glutamicus No. 541 (C. glutamicum ATCC 13058, Kimura, 1963), C. glutamicum ATCC 17965 (Cocaign-Bousquet et al., 1996), and C. glutamicum MH20-22B (lysine producer based on C. glutamicum ATCC 13032, Marx et al., 1996). Mean biomass compositions were derived to check whether experimental YAA mirror the apparent needs for biomass build-up (last column in Table 4). Thereby, three groups could be identified: (i) AAs fitting well with the reference data (Arg, His, Ile, Leu, Met, Phe, Ser, Thr, Trp), (ii) AAs taken up in higher (Gln, Asn, Pro), or (iii) in lower amounts (Ala, Glu, Gly, Lys, Val, Asp) than expected from composition. Figure 4A illustrates the observation by color-coding AAs of group (i) in white, AAs of group (ii) in yellow, and AAs of group (iii) in red. Accordingly, consumed amounts of AAs of group (ii) exceed the biomass need approximately by factor three. AAs from group (iii) poorly satisfied the composition-demand between 0 (Asp) and 52% (Glu).
TABLE 3. Amino acid (AA) consumption rates (qAA) of C. glutamicum cultivated in CGXII minimal medium supplemented with 4% (w/v) glucose and 5 mM of all AAs except L-cysteine and L-tyrosine.
TABLE 4. Comparison of published amino acid (AA)-compositions of different C. glutamicum strains determined by (A) Kimura, 1963, (B) Cocaign-Bousquet et al., 1996, (C) Marx et al., 1996, and experimentally determined AA-yields (YAA).
FIGURE 4. (A) Illustration of percentage amino acid (AA)-biomass yields (mmol g-1) calculated from AA supplementation experiments in relation to the mean AA-composition of biomass from different C. glutamicum strains (cf. Table 4). AA-amounts consumed during the experiments which were higher (above 100%) or lower (below 100%) than the published values are highlighted in yellow and red, respectively. (B) Results of the sensitivity analyses indicating the impact of AA-uptake on the growth rate of the in silico cell (cf. Eq. 8). The experimental biomass-glucose yield (YXGlc) observed for the strain grown in CGXII [4% (w/v) glucose, Table 2] is given as reference. Color-coding is adapted from panel (A). (C) Results of unbound AA-uptake (qmax,AA) by the in silico cell. The experimental glucose consumption rate (qGlc) observed for growth in CGXII [4% (w/v) glucose, Table 2] is given as reference. Color-coding is adapted from panel (A).
In silico AA-Uptake and -Consumption Studies With iMG481
Validation of stoichiometric network model iMG481
Besides studying the in vivo AA consumption of C. glutamicum WT during bioprocesses, in silico AA uptake and its influence on growth was further analyzed with the expanded stoichiometric network model iMG481. To check the suitability of the model, FBA was performed for reference conditions using qGlc = 3.55 mmol g-1 h-1 as input, leaving qO2 and qCO2 non-constrained. μSim = 0.337 h-1 was the output which agrees well with the experimental value (Table 2) and with μSim obtained from the predecessor model iEZ475 (Zelle et al., 2015). Because iMG481 was derived from iEZ475 by the implementation of additional AA importers, measured AA consumption rates, in particular qArg, qSer, and qThr, could be used as additional inputs yielding μSim = 0.481 h-1, which was slightly higher than the experimental μ (0.468 ± 0.003 h-1). Constraining respiration of the in silico cell to the experimental values reduced the predicted growth rate by ca. 20% in all tested cases. Therefore, following investigations were performed with unconstrained oxygen consumption or carbon dioxide emission.
Comparing predicted intracellular fluxes of iEZ475 with iMG481 under reference conditions revealed that the new reaction “alaT” (using Glu as amino donor for Ala-synthesis) was used by iMG481 instead of the alternative route “NCgl0388” using Val as amino donor. This result was in accordance with Marienhagen and Eggeling (2008) describing the former reaction as the main synthesis route for Ala in C. glutamicum. Furthermore, under AA supplementation conditions, the new Gln hydrolysis reaction “glsK” was used by the in silico cell instead of reactions “glmS” and “pyrG_2gln” which, likewise to “glsK,” produce Glu from exceeding amounts of Gln.
Apart from the two new synthesis pathways, iMG481 has a broader transport variability than iEZ475 by incorporation of additional transporters for Ala, Asp, and Met. Accordingly, different uptake options for Ala, Asn, Asp, Glu, Ile, Leu, Met, Phe, and Val were available which were tested with regard to the preferred uptake system of the in silico cell. When AA uptake was unconstrained and at least one of the said AA was assumed to be active, the highest μSim was always achieved by the least-energy demanding uptake system. The ranking found that diffusion was preferred to ion exchange assigning ATP-requiring uptake the least attraction. For further investigations, all transporters which were supposed to be active under non-limited substrate conditions were considered (Supplementary Table S1).
Sensitivity analysis: impact of AA uptake on μ
The impact of AA supplementation on the simulated growth rate μSim was studied by means of sensitivity analysis. Substrate consumption rates under exponential growth conditions were taken as the reference. The impact on growth increase was investigated by adding infinitesimal amounts of singular AAs with 3.5 nmol g-1 h-1 which represented 1/1000 of the glucose uptake rate constraint. Thereby, the AA-related sensitivity on growth induced by each of the 18 AAs was obtained employing Eq. (8). Figure 4B ranks the growth sensitivities of these AAs giving the experimental biomass-glucose yield (YXGlc, Table 2) as reference (dashed line). As a result, Trp, Phe, Ile, Met, Leu, Arg, and His showed sensitivities above the reference and therefore had the highest impact on μSim. Interestingly, this group of AAs was already identified in the experimental approach as AAs being consumed according to the AA-biomass composition (cf. Table 4), as highlighted by the color scheme adopted from Figure 4A. Vice versa, AAs exerting non-growth coupled impact on μSim (e.g., Pro, Glu, Gln, etc.) were either consumed in lower or higher than expected amounts taking cellular AA-biomass composition as threshold.
Complementary to the sensitivity analysis, the consequence of extreme AA uptake has also been studied. The simulations considered constant glucose uptake (3.5 mmol g-1 h-1) and eliminated any upper limit of AA uptake. The results pictured in Figure 4C show that Asn, Asp, Glu, and Gln surpassed the molar uptake of glucose (dashed line) by up to 85 times. Calculating the AA-biomass yield with the help of the simulated AA uptake and corresponding growth rate according to Eq. (1) (data not shown) further revealed that apart from the four mentioned AAs, Ala, Leu, Ser, Thr, and Val were taken up in higher amounts than expected from the model’s AA-biomass composition (cf. reaction “PROTEIN_a”). Gly, however, was the only AA with a taken-up amount below the composition value. Comparison of these AAs to the experimental AA-grouping and color scheme of Figure 4A identifies these AA (apart from Ser) as AAs exhibiting non-proportional uptake in the experiments.
The study was motivated by the fundamental question to elucidate the physiological response of C. glutamicum ATCC 13032 on increasingly nutrient-rich growth conditions. Therefore, a two-staged supplementation analysis was performed. In the first part, standard minimal medium CGXII with glucose as sugar source was supplemented with varying amounts of BHI to identify the cellular maximum of metabolic capacity. Since AAs are crucial components of BHI (in free or peptide-bound form), their influence on growth was analyzed in the second part of the study by adding an AA-cocktail comprising 18 AAs with individual concentrations of 5 mM. Using AA supplementation in synthetic medium enabled further in silico FBA studies. Both BHI and AA supplementation are discussed in the following.
Impact of BHI and AAs on Growth and Carbon Consumption of C. glutamicum
Previous growth studies by Grünberger et al. (2013) showed that the strain achieved maximum growth rates of 0.42 ± 0.03 h-1 when cultivated in lab-scale bioreactors using CGXII. Thereof, Unthan et al. (2014) identified the supplement protocatechuic acid (PCA) in CGXII as growth-accelerating compound. Strictly speaking, it is not yet fully unraveled whether PCA serves as additional carbon source, affects iron transport, enhances growth by a combinatorial effect, or may even support metabolism by other means. Accordingly, CGXII minimal medium without PCA was used in our investigations. As expected, growth rates of 0.32 ± 0.02 h-1 were observed in non-supplemented conditions which could be increased to 0.45 ± 0.01 h-1 when PCA was used (data not shown). In the following, CGXII always implies PCA-free minimal medium with glucose as standard carbon source.
Adding varying amounts of BHI (1–37 g BHI L-1) to CGXII [2% (w/v) glucose] significantly increased μ from 0.32 ± 0.02 h-1 to a maximum of 0.67 ± 0.03 h-1 revealing Monod-type dependence of the growth rate increase (Δμmax) on BHI-concentration. Interesting enough, KBHI of 2.73 g L-1 indicates the existence of a growth boosting component in BHI that is crucial for growth acceleration in low amounts but achieves saturating levels in BHI > 5 g L-1.
To elucidate whether the growth enhancement was due to AAs present in BHI, 18 AAs were added to CGXII [4% (w/v) glucose]. To ensure maximum AA uptake rates, abundant AA concentrations of 5 mM were installed which were sufficiently higher than published KM-values of the strain’s AA transporters (cf. Supplementary Table S1). Hence, an equally high growth rate of about 0.6 h-1 as induced by 37 g BHI L-1 was expected, if AAs were the only growth enhancing factors in BHI. Surprisingly, a remarkably smaller exponential μmax of 0.468 ± 0.003 h-1 was observed which equaled the addition of 2 g BHI L-1 (Figure 2). However, quantification of free and peptide-bound AAs in BHI (data not shown) revealed that only about 8 mM AAs were available in 2 g BHI L-1 which was significantly lower than the total amount used in the AA experiments (90 mM). Since these findings indicated that AAs were not the only growth enhancing factors in BHI, we evaluated the findings by further energetic analysis.
Therefore, we qualified the substrates by estimating the ATP-biomass yields of the different nutrient scenarios. Assuming that total consumed carbon (supplied by glucose and BHI) was used to produce biomass, CO2, and no by-products, CO2 emission rates and O2 consumption rates were calculated (Table 1) and compared to experimentally determined rates (Table 2) of the glucose reference conditions. ATP formation was calculated assuming a growth medium-independent P/O-ratio of 1.2 which represents a realistic estimation based on own observations. Moreover, we presumed that ATP formation equaled consumption. Accordingly, the ATP-yield of about 26 mmol g-1 for reference conditions was observed (Table 1). The estimation fits well with findings of Cocaign-Bousquet et al. (1996) (29.1 mmol g-1), but differs from the 60 mmol g-1 determined by Pfelzer (2016). For comparison, E. coli shows demands of 71 mmol g-1 (Varma et al., 1993). However, estimations strictly depend on the choice of P/O-ratio. For example, choosing the P/O ratio of two yields 43 mmol g-1. Accordingly, estimated YATPX values only show tendencies allowing comparison of the different media.
Almost identical YATPX of 23.9 ± 3.04 and 24.3 ± 2.2 mmol g-1 were observed for the addition of 1 g BHI L-1 and AA supplementation, respectively. The observation somewhat mirrors the Monod-type growth increase after BHI addition (Figure 2) compared to the growth rise after AA supplementation. Both revealed similar growth benefits, apparently resembling similar energetic status of the cells. However, adding BHI ≥ 5 g L-1, steadily decreased YATPX to about 14 mmol g-1 which corresponds to a 46% drop with respect to the reference. Several hypothesis can be formulated: (i) under reference, 1 g BHI L-1-, and AA-supplemented conditions, ATP demanding futile cycles existed which disappeared under BHI-rich conditions; (ii) P/O ratios are not constant but increase under BHI-rich conditions which would compensate the YATPX drop; and (iii) BHI-rich media provide crucial precursors, thereby preventing own ATP-demanding cellular synthesis. Hypothesis (i) appears unlikely qualifying a 46% share of ATP production for futile cycling production as too high. Hypothesis (ii) requires for detailed (13C based) metabolic flux analysis which is beyond the scope of this study. Ruling out options (i) and (ii), hypothesis (iii) was taken as the guiding idea for subsequent analysis.
Accordingly, net ATP costs for AA biosynthesis starting from central metabolism precursors were estimated modifying the approach of Kaleta et al. (2013) who calculated the net ATP costs for AA synthesis routes in E. coli by considering ATP costs in the synthesis pathways but also accounted for the amount of produced ATP by precursor synthesis (cf. Table 4 and Supplementary Tables S2, S3). In contrast to Kaleta et al. (2013), NADPH was not considered as energy (ATP) yielding precursor in respiration of C. glutamicum because the strain does not possess nicotinamide nucleotide-dependent transhydrogenases like E. coli (Kabus et al., 2007). As a result, ATP costs of about 4 mmol for 1 g biomass were calculated that represents less than 15% of the net ATP formation from glucose (about 26 mmol g-1, Table 1) considering the P/O ratio of 1.2. Accordingly, enriching the medium by AAs did create only minor energetic benefit as shown for the AA supplemented medium compared to the reference (difference of approximately 2 mmol g-1). Studies of Ingraham et al. (1983) analyzing ATP demands of E. coli grown in minimal and rich medium further support the observation: ATP demands of protein formation outnumber those of AA formation by far. Similar to our findings for C. glutamicum, E. coli revealed 40% reduced YATPX values when cultivated in rich medium compared to minimal medium (Ingraham et al., 1983). Apparently, the reduction of 46% observed for C. glutamicum agrees well. Most likely, other ATP-demanding building blocks, such as fatty acids or nucleotides, are consumed in highly concentrated BHI environments which relieves the cellular need to create ATP thereby reducing the glucose uptake accordingly. Therefore, the growth phenotype resulting from different nutrient supplies can be explained as follows: C. glutamicum benefits from the abundant addition of AAs by omitting the related AA biosynthesis which yields a moderate μ increase. Growth advantages of AA supplementation may be rooted in the elimination of kinetics limits. For instance, external supply of AAs may fuel intracellular transaminases thereby preventing own “slow” biosynthesis. Alternatively, external supply may meet short-term needs for charging tRNAs as a prerequisite of fast translation for protein formation (Nieß et al., 2017). Adding high amounts of BHI further released the cell of ATP-costly fatty acid and nucleotide synthesis enhancing μ to its apparent limit of 0.67 h-1. Necessary ATP production was reduced to about 14 mmol ATP per g biomass resembling the remaining cellular needs for maintenance and biosynthesis. The strong growth boosting effect of 1 g BHI L-1 (Figure 2) anticipates that other growth stimulating components beside AAs (about 4 mM, cf. above) play a role, too. Likely candidates are trace elements or vitamins. The latter have been tested (supplementation of 0.2 mg thiamine-HCl L-1 in CGXII) and an elevated growth rate of approximately 0.4 h-1 was determined (data not shown).
Biomass specific glucose uptake rates qGlc also showed particularities for different nutrient-rich media. Whereas 1 g BHI L-1 caused qGlc increase, further BHI supplementation yielded decreasing rates (Table 1) reaching reference-like levels in 37 g BHI L-1 and in the AA supplementation experiments. Notably, TC consumption increased proportionally with μ showing a biomass-carbon yield YXC of 18.47 g C-mol-1. The observations agree with findings of Dominguez et al. (1997) who found YXC of 17.6 g C-mol-1 for the consumption of different glucose-fructose mixtures. Wendisch et al. (2000) observed 12 g C-mol-1 for the consumption of glucose-acetate mixtures and Frunzke et al. (2008) measured 13 g C-mol-1 for glucose, gluconate, and glucose–gluconate mixtures in shaking flask cultivations. Our study somewhat differs from the others by considering not only carbon uptake alone but also the parallel consumption of nitrogen containing AAs. Consequently, the maximum value of 18.47 g C-mol-1 may reflect the total cellular capacity limit of processing carbon and nitrogen (via AAs) through metabolism and translation, whereas YXC of 17.6 g C-mol-1 mirrors limits of sugar uptake and metabolism. Interestingly, both values are similar which hints to optimized carbon-to-biomass conversion irrespective of the medium composition and to equally adjusted capacities for AA processing, i.e., protein formation.
Amino Acid Consumption of the in vivo and in silico Cell
The addition of an abundant supply of 18 AAs in minimal medium revealed three groups of AAs with distinctive consumption characteristics. AAs of the first group, consisting of Arg, His, Ile, Leu, Met, Phe, Ser, Thr, and Trp, were consumed according to AA-biomass composition. Biosynthesis pathways of most of these AAs, e.g., Met, Trp, or Phe, are energy-intensive, i.e., they have a positive net ATP-cost (Supplementary Table S3). Consequently, the cell safes energy when the AAs are consumed instead of being synthesized. The experimental observation agrees well with in silico analyses: growth rate-sensitivity studies revealed that the same AAs (apart from Thr and Ser) had the strongest influence on μ. Except for serving as AA building block like the other AAs, Ser may serve as carbon source in the absence of other available C-sources (Netzer et al., 2004). It is a precursor for synthesis of activated C1-compounds and for L-cysteine synthesis (Wada and Takagi, 2006). The latter was not supplemented in the experiments.
The second group consisted of Asn, Gln, and Pro showing consumption rates more than threefold higher than expected from the AA-biomass composition. In contrast, the third group of AAs (Ala, Asp, Glu, Gly, Lys, Val) was consumed remarkably below the expected threshold. The consumption rate of Gln was the highest of all AAs (0.694 mmol g-1 h-1) and approximated the maximum uptake rate of the Gln transporter (Siewe et al., 1995; Supplementary Table S1). Notably, Gln may serve carbon and nitrogen demands at the same time being decomposed by hydrolysis or deamination via glutaminase Glsk or comparable enzymes (Buerger et al., 2016). Glu and ammonia are the products which correlates increased Gln uptake with the reduced need of Glu and ammonia consumption. In case of Glu, the inadequate consumption of the AA might have been caused by the presence of glucose in the medium causing catabolite repression (Krämer et al., 1990; Kronemeyer et al., 1995). Likewise, the aforementioned compensation via Gln uptake may explain the observed low Glu consumption rates. Taking Gln instead of Glu also saves ATP which would have been needed when Gln had to be synthesized from Glu via Gln synthase. Irrespective whether Glu is produced from Gln, taken up, or synthesized from alpha ketoglutarate via glutamate dehydrogenase, Glu may serve as amino donor in transaminase reactions such as aspartate transaminase (AspB) or alanine transferase (AlaT, Marienhagen and Eggeling, 2008) thereby forming Asp or Ala, respectively. Accordingly, Asp and Ala uptake rates are likely to be reduced as observed in the experiments. Ala itself may serve as amino donor for Val production via AvtA-reaction (Marienhagen and Eggeling, 2008) which finally links low Val, Ala, Asp, and Glu uptakes with tremendously high Gln consumption. For completion, Ala uptake via MetP (Trötschel et al., 2008; Supplementary Table S1) may have been hampered further by competition with Met, while Lange et al. (2003) showed that C. glutamicum WT cannot use Val as sole carbon or nitrogen source if supplemented in CGXII shaking flask experiments. Similar to Gln, the immense Asn uptake likewise reduced the uptake of other AAs of the Asp-family. Aparaginase AnsA which is suspected to act as an overflow enzyme (Mesas et al., 1990) might have been active to hydrolyze Asn to Asp and ammonia thereby minimizing the need for Asp uptake. Besides, Ziert (2014) showed that Asp is a poor nitrogen source, is even less attractive as carbon source, and was therefore not consumed in batch Asp supplementation experiments conducted in microflower plates. Interesting enough, Asp is a precursor for Lys biosynthesis which may explain why immense Asn uptake and its conversion to Asp may finally reduce Lys uptake to 26% of the biomass needs. The corresponding consumption rate (0.025 mmol g-1 h-1) was in the same range as the Lys uptake determined in isotopic studies by Seep-Feldhaus et al. (1991) (0.015 mmol g-1 h-1). Both rates are significantly lower than the rate needed to satisfy the composition demand at the measured growth rate of 0.468 h-1 (0.099 mmol g-1 h-1). Summarizing, the immense uptake of Gln and Asn was likely to minimize the uptakes of Ala, Asp, Glu, Lys, and Val, which is also reflected by FBA assuming unlimited uptake of one single AA. In particular, Asn, Asp, Gln, and Glu, the protagonists of the transaminase network, revealed maximum (theoretical) uptake rates in addition to glucose (Figure 4C).
So far, no clear experimental observations give evidence on crucial Pro consumption (Bott and Niebisch, 2003; Jensen and Wendisch, 2013) although putA (cg0129) has been annotated in the genome of C. glutamicum encoding proline dehydrogenase which would ultimately provide Glu and NADH (Ling et al., 1994). Besides, the low consumption of Gly (14% with respect to the biomass composition) might be explained by its intracellular production via serine hydroxymethyltransferase being a by-product of 5,10-methylenetetrahydrofolate (MTHF) synthesis (Simic et al., 2002). FBA simulations supported the hypothesis: unlimited Gly uptake as a co-substrate of limited glucose uptake revealed that only 60% of the consumed Gly replenished the in silico cell’s demand for Gly, while the remaining amount was available as by-product of the MTHF pathway.
In this two-staged study, physiological investigations were conducted to improve growth and nutrient consumption of the industrial production host C. glutamicum. Supplementation of minimal medium with complex compounds (BHI) relieved the strain of energy intensive biosynthesis, most likely of proteins and fatty acids, subsequently reducing the ATP demand of the cell. In turn, growth was accelerated to the WT strain’s apparent maximum of 0.67 h-1 in 37 g BHI L-1, while the TC consumption was proportionally increased. However, glucose consumption was only stimulated slightly under certain circumstances (1–10 g BHI L-1) and revealed falling tendencies in rich environments. Besides, biomass carbon yields deduced from BHI and AA supplementation experiments revealed nutrient-irrespective limits of about 18 g C-mol-1 that may be a target of further strain engineering studies. Experimental consumption rates of 18 AAs under process-relevant conditions were reported, and Asn and Gln were identified as main protagonists due to their role in the transaminase network of C. glutamicum. FBA simulations with iMG481 confirmed the experimental findings. Summarizing, the revealed AA uptake kinetics, the ranking of individual AA needs, observed AA exchange abilities, the properties of BHI supplementation, and the unraveled apparent limit of carbon-to-biomass conversion offer a highly fruitful basis for improving process productivities and for metabolic engineering. Such studies may expand the application of C. glutamicum not only as an AAs producer but also for the production of rather growth coupled products such as technical proteins.
MG designed the study, carried out the bioreactor experiments, analyzed the datasets, and drafted the manuscript. MG and JZ expanded the stoichiometric network model and performed the flux balance analyses. TH and AN analyzed the datasets and corrected the manuscript. BB and RT conceived the whole study and corrected the manuscript. All authors read and approved the final manuscript.
The authors gratefully acknowledge the funding of this work by the German Federal Ministry of Education and Research (BMBF; Grant 031A302A).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank Anna-Lena Mayer, Salaheddine Laghrami, and Andreas Freund for excellent support with bioreactor fermentations, as well as Mira Lenfers-Lücker for assistance with the HPLC analyses. Moreover, we thank all members of the 0.6 plus project group for valuable comments on the topic and for their great cooperation.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.02058/full#supplementary-material
Baev, M. V., Baev, D., Radek, A. J., and Campbell, J. W. (2006a). Growth of Escherichia coli MG1655 on LB medium: monitoring utilization of sugars, alcohols, and organic acids with transcriptional microarrays. Appl. Microbiol. Biotechnol. 71, 310–316. doi: 10.1007/s00253-006-0317-6
Baev, M. V., Baev, D., Radek, A. J., and Campbell, J. W. (2006b). Growth of Escherichia coli MG1655 on LB medium: monitoring utilization of amino acids, peptides, and nucleotides with transcriptional microarrays. Appl. Microbiol. Biotechnol. 71, 317–322. doi: 10.1007/s00253-005-0310-5
Baev, M. V., Baev, D., Radek, A. J., and Campbell, J. W. (2006c). Growth of Escherichia coli MG1655 on LB medium: determining metabolic strategy with transcriptional microarrays. Appl. Microbiol. Biotechnol. 71, 323–328. doi: 10.1007/s00253-006-0392-8
Bäumchen, C., Knoll, A., Husemann, B., Seletzky, J., Maier, B., Dietrich, C., et al. (2007). Effect of elevated dissolved carbon dioxide concentrations on growth of Corynebacterium glutamicum on D-glucose and L-lactate. J. Biotechnol. 128, 868–874. doi: 10.1016/j.jbiotec.2007.01.001
Becker, J., and Wittmann, C. (2012). Bio-based production of chemicals, materials and fuels–Corynebacterium glutamicum as versatile cell factory. Curr. Opin. Biotechnol. 23, 631–640. doi: 10.1016/j.copbio.2011.11.012
Blombach, B., Buchholz, J., Busche, T., Kalinowski, J., and Takors, R. (2013). Impact of different CO2/HCO3- levels on metabolism and regulation in Corynebacterium glutamicum. J. Biotechnol. 168, 331–340. doi: 10.1016/j.jbiotec.2013.10.005
Buchholz, J., Graf, M., Freund, A., Busche, T., Kalinowski, J., Blombach, B., et al. (2014b). CO 2/HCO 3- perturbations of simulated large scale gradients in a scale-down device cause fast transcriptional responses in Corynebacterium glutamicum. Appl. Microbiol. Biotechnol. 98, 8563–8572. doi: 10.1007/s00253-014-6014-y
Buerger, J., Rehm, N., Grebenstein, L., and Burkovski, A. (2016). Glutamine metabolism of Corynebacterium glutamicum: role of the glutaminase GlsK. FEMS Microbiol. Lett. 363:fnw230. doi: 10.1093/femsle/fnw230
Burkovski, A., Weil, B., and Krämer, R. (1996). Characterization of a secondary uptake system for L-glutamate in Corynebacterium glutamicum. FEMS Microbiol. Lett. 136, 169–173. doi: 10.1111/j.1574-6968.1996.tb08044.x
Cocaign-Bousquet, M., Guyonvarch, A., and Lindley, N. D. (1996). Growth rate-dependent modulation of carbon Flux through central metabolism and the kinetic consequences for glucose-limited chemostat cultures of Corynebacterium glutamicum. Appl. Environ. Microbiol. 62, 429–436.
Coello, N., Montiel, E., Concepcion, M., and Christen, P. (2002). Optimisation of a culture medium containing fish silage for L-lysine production by Corynebacterium glutamicum. Bioresour. Technol. 85, 207–211. doi: 10.1016/S0960-8524(02)00084-6
Cox, R. A. (2004). Quantitative relationships for specific growth rates and macromolecular compositions of Mycobacterium tuberculosis, Streptomyces coelicolor A3(2) and Escherichia coli B/r: an integrative theoretical approach. Microbiology 150, 1413–1426. doi: 10.1099/mic.0.26560-0
Dominguez, H., Cocaign-Bousquet, M., and Lindley, N. D. (1997). Simultaneous consumption of glucose and fructose from sugar mixtures during batch growth of Corynebacterium glutamicum. Appl. Microbiol. Biotechnol. 47, 600–603. doi: 10.1007/s002530050980
Dominguez, H., Nezondet, C., Lindley, N. D., and Cocaign, M. (1993). Modified carbon flux during oxygen limited growth of Corynebacterium glutamicum and the consequences for amino acid overproduction. Biotechnol. Lett. 15, 449–454. doi: 10.1007/BF00129316
Frunzke, J., Engels, V., Hasenbein, S., Gätgens, C., and Bott, M. (2008). Co-ordinated regulation of gluconate catabolism and glucose uptake in Corynebacterium glutamicum by two functionally equivalent transcriptional regulators, GntR1 and GntR2. Mol. Microbiol. 67, 305–322. doi: 10.1111/j.1365-2958.2007.06020.x
Gong, S., Richard, H., and Foster, J. W. (2003). YjdE (AdiC) is the arginine: agmatine antiporter essential for arginine-dependent acid resistance in Escherichia coli. J. Bacteriol. 185, 4402–4409. doi: 10.1128/JB.185.15.4402-4409.2003
Gourdon, P., Raherimandimby, M., Dominguez, H., Cocaign-Bousquet, M., and Lindley, N. D. (2003). Osmotic stress, glucose transport capacity and consequences for glutamate overproduction in Corynebacterium glutamicum. J. Biotechnol. 104, 77–85. doi: 10.1016/S0168-1656(03)00165-2
Grünberger, A., van Ooyen, J., Paczia, N., Rohe, P., Schiendzielorz, G., Eggeling, L., et al. (2013). Beyond growth rate 0.6: Corynebacterium glutamicum cultivated in highly diluted environments. Biotechnol. Bioeng. 110, 220–228. doi: 10.1002/bit.24616
Heirendt, L., Arreckx, S., Pfau, T., Mendoza, S. N., Richelle, A., Heinken, A., et al. (2017). Creation and analysis of biochemical constraint-based models: the COBRA Toolbox v3.0. [Preprint]. Available at: https://arxiv.org/abs/1710.04038
Jeon, J. M., Song, E. J., and Lee, H. W. (2013). Media optimization of Corynebacterium glutamicum for succinate production under oxygen-deprived condition. J. Microbiol. Biotechnol. 23, 211–217. doi: 10.4014/jmb.1206.06057
Kabus, A., Georgi, T., Wendisch, V. F., and Bott, M. (2007). Expression of the Escherichia coli pntAB genes encoding a membrane-bound transhydrogenase in Corynebacterium glutamicum improves L-lysine formation. Appl. Microbiol. Biotechnol. 75, 47–53. doi: 10.1007/s00253-006-0804-9
Kalinowski, J., Bathe, B., Bartels, D., Bischoff, N., Bott, M., Burkovski, A., et al. (2003). The complete Corynebacterium glutamicum ATCC 13032 genome sequence and its impact on the production of L-aspartate-derived amino acids and vitamins. J. Biotechnol. 104, 5–25. doi: 10.1016/S0168-1656(03)00154-8
Käß, F., Hariskos, I., Michel, A., Brandt, H. J., Spann, R., Junne, S., et al. (2014a). Assessment of robustness against dissolved oxygen/substrate oscillations for C. glutamicum DM1933 in two-compartment bioreactor. Bioprocess Biosyst. Eng. 37, 1151–1162. doi: 10.1007/s00449-013-1086-0
Käß, F., Junne, S., Neubauer, P., Wiechert, W., and Oldiges, M. (2014b). Process inhomogeneity leads to rapid side product turnover in cultivation of Corynebacterium glutamicum. Microb. Cell Fact. 13:6. doi: 10.1186/1475-2859-13-6
Kind, S., Jeong, W. K., Schröder, H., Zelder, O., and Wittmann, C. (2010). Identification and elimination of the competing N-acetyldiaminopentane pathway for improved production of diaminopentane by Corynebacterium glutamicum. Appl. Environ. Microbiol. 76, 5175–5180. doi: 10.1128/AEM.00834-10
Kjeldsen, K. R., and Nielsen, J. (2009). In silico genome-scale reconstruction and validation of the Corynebacterium glutamicum metabolic network. Biotechnol. Bioeng. 102, 583–597. doi: 10.1002/bit.22067
Kronemeyer, W., Peekhaus, N., Krämer, R., Sahm, H., and Eggeling, L. (1995). Structure of the gluABCD cluster encoding the glutamate uptake system of Corynebacterium glutamicum. J. Bacteriol. 177, 1152–1158. doi: 10.1128/jb.177.5.1152-1158.1995
Lange, C., Rittmann, D., Wendisch, V. F., Bott, M., and Sahm, H. (2003). Global expression profiling and physiological characterization of Corynebacterium glutamicum grown in the presence of L-valine. Appl. Environ. Microbiol. 69, 2521–2532. doi: 10.1128/AEM.69.5.2521-2532.2003
Lange, J., Müller, F., Bernecker, K., Dahmen, N., Takors, R., and Blombach, B. (2017). Valorization of pyrolysis water: a biorefinery side stream, for 1, 2-propanediol production with engineered Corynebacterium glutamicum. Biotechnol. Biofuels 10:277. doi: 10.1186/s13068-017-0969-8
Lee, J., Sim, S. J., Bott, M., Um, Y., Oh, M. K., and Woo, H. M. (2014). Succinate production from CO2-grown microalgal biomass as carbon source using engineered Corynebacterium glutamicum through consolidated bioprocessing. Sci. Rep. 4:5819. doi: 10.1038/srep05819
Ling, M., Allen, S. W., and Wood, J. M. (1994). Sequence analysis identifies the proline dehydrogenase and Δ1-pyrroline-5-carboxylate dehydrogenase domains of the multifunctional Escherichia coli PutA protein. J. Mol. Biol. 243, 950–956. doi: 10.1006/jmbi.1994.1696
Marienhagen, J., and Eggeling, L. (2008). Metabolic function of Corynebacterium glutamicum aminotransferases AlaT and AvtA and impact on L-valine production. Appl. Environ. Microbiol. 74, 7457–7462. doi: 10.1128/AEM.01025-08
Marin, K., and Krämer, R. (2007). “Amino acid transport systems in biotechnologically relevant bacteria,” in Amino Acid Biosynthesis - Pathways, Regulation and Metabolic Engineering, ed. V. F. Wendisch (Berlin: Springer), 289–325.
Marx, A., de Graaf, A. A., Wiechert, W., Eggeling, L., and Sahm, H. (1996). Determination of the fluxes in the central metabolism of Corynebacterium glutamicum by nuclear magnetic resonance spectroscopy combined with metabolite balancing. Biotechnol. Bioeng. 49, 111–129. doi: 10.1002/(SICI)1097-0290(19960120)49:2<111::AID-BIT1>3.0.CO;2-T
Marx, A., Striegel, K., de Graaf, A. A., and Eggeling, L. (1997). Response of central metabolism of Corynebacterium glutamicum to different flux burdens. Biotechnol. Bioeng. 56, 168–180. doi: 10.1002/(SICI)1097-0290(19971020)56:2
Mesas, J. M., Gil, J. A., and Martín, J. F. (1990). Characterization and partial purification of L-asparaginase from Corynebacterium glutamicum. Microbiology 136, 515–519. doi: 10.1099/00221287-136-3-515
Michel, A., Koch-Koerfges, A., Krumbach, K., Brocker, M., and Bott, M. (2015). Anaerobic growth of Corynebacterium glutamicum via mixed-acid fermentation. Appl. Environ. Microbiol. 81, 7496–7508. doi: 10.1128/AEM.02413-15
Netzer, R., Peters-Wendisch, P., Eggeling, L., and Sahm, H. (2004). Cometabolism of a nongrowth substrate: L-serine utilization by Corynebacterium glutamicum. Appl. Environ. Microbiol. 70, 7148–7155. doi: 10.1128/AEM.70.12.7148-7155.2004
Nieß, A., Failmezger, J., Kuschel, M., Siemann-Herzberg, M., and Takors, R. (2017). Experimentally validated model enables debottlenecking of in vitro protein synthesis and identifies a control shift under in vivo conditions. ACS Synth. Biol. 6, 1913–1921. doi: 10.1021/acssynbio.7b00117
Nishimura, T., Vertès, A. A., Shinoda, Y., Inui, M., and Yukawa, H. (2007). Anaerobic growth of Corynebacterium glutamicum using nitrate as a terminal electron acceptor. Appl. Microbiol. Biotechnol. 75, 889–897. doi: 10.1007/s00253-007-0879-y
Okino, S., Inui, M., and Yukawa, H. (2005). Production of organic acids by Corynebacterium glutamicum under oxygen deprivation. Appl. Microbiol. Biotechnol. 68, 475–480. doi: 10.1007/s00253-005-1900-y
Palmieri, L., Berns, D., Krämer, R., and Eikmanns, M. (1996). Threonine diffusion and threonine transport in Corynebacterium glutamicum and their role in threonine production. Arch. Microbiol. 165, 48–54. doi: 10.1007/s002030050295
Sauer, U., Canonaco, F., Heri, S., Perrenoud, A., and Fischer, E. (2004). The soluble and membrane-bound transhydrogenases UdhA and PntAB have divergent functions in NADPH metabolism of Escherichia coli. J. Biol. Chem. 279, 6613–6619. doi: 10.1074/jbc.M311657200
Scheele, S., Oertel, D., Bongaerts, J., Evers, S., Hellmuth, H., Maurer, K. H., et al. (2013). Secretory production of an FAD cofactor-containing cytosolic enzyme (sorbitol–xylitol oxidase from Streptomyces coelicolor) using the twin-arginine translocation (Tat) pathway of Corynebacterium glutamicum. Microb. Biotechnol. 6, 202–206. doi: 10.1111/1751-7915.12005
Seep-Feldhaus, A. H., Kalinowski, J., and Pühler, A. (1991). Molecular analysis of the Corynebacterium glutamicum gene involved in lysine uptake. Mol. Microbiol. 5, 2995–3005. doi: 10.1111/j.1365-2958.1991.tb01859.x
Shang, X., Zhang, Y., Zhang, G., Chai, X., Deng, A., Liang, Y., et al. (2013). Characterization and molecular mechanism of AroP as an aromatic amino acid and histidine transporter in Corynebacterium glutamicum. J. Bacteriol. 195, 5334–5342. doi: 10.1128/JB.00971-13
Shinfuku, Y., Sorpitiporn, N., Sono, M., Furusawa, C., Hirasawa, T., and Shimizu, H. (2009). Development and experimental verification of a genome-scale metabolic model for Corynebacterium glutamicum. Microb. Cell Fact. 8:43. doi: 10.1186/1475-2859-8-43
Simic, P., Willuhn, J., Sahm, H., and Eggeling, L. (2002). Identification of glyA (encoding serine hydroxymethyltransferase) and its use together with the exporter ThrE to increase L-threonine accumulation by Corynebacterium glutamicum. Appl. Environ. Microbiol. 68, 3321–3327. doi: 10.1128/AEM.68.7.3321-3327.2002
Takors, R., Bathe, B., Rieping, M., Hans, S., Kelle, R., and Huthmacher, K. (2007). Systems biology for industrial strains and fermentation processes—example: amino acids. J. Biotechnol. 129, 181–190. doi: 10.1016/j.jbiotec.2007.01.031
Tao, H., Bausch, C., Richmond, C., Blattner, F. R., and Conway, T. (1999). Functional genomics: expression analysis of Escherichia coli growing on minimal and rich media. J. Bacteriol. 181, 6425–6440.
Trötschel, C., Follmann, M., Nettekoven, J. A., Mohrbach, T., Forrest, L. R., Burkovski, A., et al. (2008). Methionine uptake in Corynebacterium glutamicum by MetQNI and by MetPS, a novel methionine and alanine importer of the NSS neurotransmitter transporter family. Biochemistry 47, 12698–12709. doi: 10.1021/bi801206t
Unthan, S., Grünberger, A., van Ooyen, J., Gätgens, J., Heinrich, J., Paczia, N., et al. (2014). Beyond growth rate 0.6: what drives Corynebacterium glutamicum to higher growth rates in defined medium. Biotechnol. Bioeng. 111, 359–371. doi: 10.1002/bit.25103
Wendisch, V. F., Bott, M., and Eikmanns, B. J. (2006). Metabolic engineering of Escherichia coli and Corynebacterium glutamicum for biotechnological production of organic acids and amino acids. Curr. Opin. Microbiol. 9, 268–274. doi: 10.1016/j.mib.2006.03.001
Wendisch, V. F., Brito, L. F., Lopez, M. G., Hennig, G., Pfeifenschneider, J., Sgobba, E., et al. (2016). The flexible feedstock concept in industrial biotechnology: metabolic engineering of Escherichia coli, Corynebacterium glutamicum, Pseudomonas, Bacillus and yeast strains for access to alternative carbon sources. J. Biotechnol. 234, 139–157. doi: 10.1016/j.jbiotec.2016.07.022
Wendisch, V. F., de Graaf, A. A., Sahm, H., and Eikmanns, B. J. (2000). Quantitative determination of metabolic fluxes during coutilization of two carbon sources: comparative analyses with Corynebacterium glutamicum during growth on acetate and/or glucose. J. Bacteriol. 182, 3088–3096. doi: 10.1128/JB.182.11.3088-3096.2000
Wieschalka, S., Blombach, B., and Eikmanns, B. J. (2012). Engineering Corynebacterium glutamicum for the production of pyruvate. Appl. Microbiol. Biotechnol. 94, 449–459. doi: 10.1007/s00253-011-3843-9
Zelle, E., Nöh, K., and Wiechert, W. (2015). “Growth and production capabilities of Corynebacterium glutamicum: interrogating a genome-scale metabolic network model”, in Corynebacterium glutamicum: from Systems Biology to Biotechnological Applications, ed. A. Burkovski (Norfolk, VA: Caister Academic Press), 39–54. doi: 10.21775/9781910190050.04
Zhang, Y., Cai, J., Shang, X., Wang, B., Liu, S., Chai, X., et al. (2017). A new genome-scale metabolic model of Corynebacterium glutamicum and its application. Biotechnol. Biofuels 10:169. doi: 10.1186/s13068-017-0856-3
Ziert, C. (2014). Metabolic Engineering of Corynebacterium glutamicum for the Production of L-aspartate and its Derivatives β-Alanine and Ectoine. Doctoral dissertation, Bielefeld, Universität Bielefeld
Keywords: Corynebacterium glutamicum, growth rate, substrate consumption, carbon and nitrogen consumption, amino acid uptake, stoichiometric network model iMG481, sensitivity analyses
Citation: Graf M, Zieringer J, Haas T, Nieß A, Blombach B and Takors R (2018) Physiological Response of Corynebacterium glutamicum to Increasingly Nutrient-Rich Growth Conditions. Front. Microbiol. 9:2058. doi: 10.3389/fmicb.2018.02058
Received: 29 May 2018; Accepted: 13 August 2018;
Published: 29 August 2018.
Edited by:Dimitris G. Hatzinikolaou, National and Kapodistrian University of Athens, Greece
Reviewed by:Stefan Junne, Technische Universität Berlin, Germany
Armen Trchounian, Yerevan State University, Armenia
Copyright © 2018 Graf, Zieringer, Haas, Nieß, Blombach and Takors. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ralf Takors, email@example.com