# Dynamic Metabolic Flux Analysis Demonstrated on Cultures Where the Limiting Substrate Is Changed from Carbon to Nitrogen andVice Versa

Authors: Gaspard Lequeux; Joeri Beauprez; Jo Maertens; Ellen Van Horen; Wim Soetaert; Erick Vandamme; Peter A. Vanrolleghem
Journal: Journal of Biomedicine and Biotechnology
 (2010)

---

## Abstract

The main requirement for metabolic flux analysis (MFA) is that the
cells are in a pseudo-steady state, that there is no accumulation or depletion of intracellular metabolites. In the
past, the applications of MFA were limited to the analysis of continuous cultures. This contribution introduces the concept of
dynamic MFA and extends MFA so that it is applicable to transient cultures. Time series of concentration measurements are
transformed into flux values. This transformation involves differentiation, which typically increases the noisiness of the
data. Therefore, a noise-reducing step is needed. In this work, polynomial smoothing was used. As a test case, dynamic MFA is applied onEscherichia coli cultivations shifting from carbon limitation to nitrogen limitation and vice versa. After switching the limiting substrate from N to C, a lag phase was observed accompanied with an increase in maintenance
energy requirement. This lag phase did not occur in the C- to N-limitation case.


---

## Body

## 1. Introduction

Material balances of reactors describe how certain compounds go in the reactor and are transformed into other compounds. Metabolic models consider the cell as the reactor and, additionally to the exchange rates, also include information on the different steps (reactions) on how the input compounds are transformed to the output compounds. This should allow to gain insight in the internal cellular fluxes. Typically the number of equations (material balances of the different compounds) obtained is less than the number of unknowns (exchange rates and reaction rates) one wants to calculate.Techniques for solving metabolic models can be divided into two main groups: data-driven techniques and knowledge/assumption-driven techniques.The prime example of a data-driven technique is metabolic flux analysis (MFA) where intracellular fluxes are calculated based on measured exchange rates. The mathematical technique used was initially developed for black box modeling [1, 2], but can easily be applied to stoichiometric models [3, 4]. Exchange rates are not always sufficient to solve a stoichiometric model. In such cases, the option exists to measure intracellular fluxes via 13C enrichment analysis [5].The counterpart of MFA is flux balance analysis (FBA) which uses linear optimisation to  solve the stoichiometric model [6] and does not directly use measurements (although measurements can be included [7]). Whereas FBA typically yields one optimal solution, the calculation of elementary modes or extreme pathways gives a full view of the metabolic capabilities of an organism [8, 9]. An overview of the different extensions applied to MFA and FBA is given in [10].Both MFA and FBA have the fundamental assumption that the cell is in a pseudo-steady state, thus that there is no accumulation or depletion of intracellular metabolites. This implies that metabolic models can only be used for cells in pseudo-steady state. Models that fully describe the cellular dynamics exist, but are complex and cumbersome to implement and/or very limited in scope [11, 12].However, some examples of FBA applied in nonstationary cases can be found. Varma and Palsson [13] used flux balance analysis to describe fed batches by iteratively solving the model for maximal biomass production. At each time point, the available substrate is calculated from the results of the FBA in the previous step. This way the time profiles of cell density, glucose, and by-products could be quantitatively predicted. This concept of dynamic FBA (DFBA) is further developed by Mahadevan et al. [14]. They formalise the methodology of Varma and Palsson [13] and name it static optimisation-based DFBA. They also introduce a new methodology, called dynamic optimisation-based DFBA, in which the optimisation is done over the entire time period of interest to obtain time profiles of fluxes. The methodology of dynamic optimisation-based DFBA is combined with the concept of minimisation of metabolic adjustment (MOMA) to model the myocardial energy metabolism [15]. Minimisation of metabolic adjustment (MOMA) is an alternative goal function to better predict the behaviour of mutated strains [16]. Using this technique, it was shown that the cellular metabolism does not always remain optimal during transient perturbations [16]. Lee et al. [17] introduced integrated dynamics FBA (idFBA) that dynamically simulates cellular phenotypes arising from integrated networks. These networks combine metabolic models with genetic, regulatory, and intracellular signalling information.In this contribution, not FBA but MFA is extended so that it can be used for modelling cultures during transient conditions. This is reasonable as intracellular pseudo-steady state can also be assumed under certain dynamic conditions, because of the relatively small time constants of cellular processes, for example, mass action and metabolic adaption to novel conditions, in comparison with processes affecting the observed environmental conditions. Indeed, perturbation experiments performed inE. coli [11, 18–20] showed that intracellular metabolite pools reach a new pseudo-steady state after 20 seconds. Hence, an intracellular pseudo-steady state can be assumed during the transient.Pseudo-steady state conditions are typically encountered in chemostat operated cultures [21, 22]. Chemostats are mostly run under the concept that only one substrate is limiting. However, conditions do exist where multiple substrates can be limiting at the same time. Practical applications of such dual limitations have been described by Zinn et al. [23]. Egli [24] has described concentration zones in which both glucose and nitrogen are limiting. In this contribution, however, the limiting compound in the feed medium of steady state cultures is switched from glucose to ammonia and vice versa. No abrupt changes are thus sensed by the cells (the reactor broth is only gradually replaced), an ideal test case for applying dynamic FBA. The exchange fluxes of metabolites are determined based on their concentration in the reactor broth and are subsequently used to solve an overdetermined metabolic model, resulting in the determination of the intracellular flux values during the transient period. As exchange fluxes are based on the derivative of the concentration data, the latter ones are first smoothed to avoid amplification of the noise in the exchange fluxes.


## 2. Materials and Methods

### 2.1. Bacterial Strain

Escherichia coli MG1655 [λ-,F-,  rph–1, (del  fnr)] was obtained from the Netherlands Culture Collection of Bacteria (NCCB, Utrecht, The Netherlands).


### 2.2. Culture Conditions

#### 2.2.1. Media

Composition of the Luria Bertani Broth (LB) medium and shake flask medium (for the preculture) can be found in [25].The minimal medium used in the reactor consisted of 2.5 g/L (N-limited medium) or 5 g/L (C-limited medium)(NH4)2SO4, 2 g/L KH2PO4, 0.5 g/L NaCl, 0.5 g/L MgSO4·7H2O, 33 g/L (N-limited medium) or 16.5 g/L (C-limited medium) Glucose·H2O, 1 mL/L vitamin stock solution and 100 μL/L molybdate stock solution. Vitamin stock solution consisted of 3.6 g/L FeCl2·4H2O, 5 g/L CaCl2·2H2O, 1.3 g/L MnCl2·2H2O, 0.38 g/L CuCl2·2H2O, 0.5 g/L CoCl2·6H2O, 0.94 g/L ZnCl2, 0.0311 g/L H3BO4, 0.4 g/L Na2EDTA·2H2O, and 1.01 g/L thiamine·HCl. The molybdate stock solution contained 0.967 g/L Na2MoO4·2H2O. All components were dissolved and filter-sterilised (pore size 0.22 μm, Sartobran, Sartorius, Belgium). The pH was left at approximately 5.4.


#### 2.2.2. Cultivations

Culture conditions and inoculation procedure are described in [25].The experiment in which the limiting substrate was changed from glucose to ammonia was conducted at a dilution rate of 0.155 h-1. The experiment in which the limitation was changed from nitrogen to carbon was conducted at a dilution rate of 0.142 h-1.For each experiment, steady state samples were taken: optical density (OD), cell dry weight (CDW), and HPLC analyses. These samples were taken after waiting at least five residence times after the batch phase. To ensure that the perturbations caused by sampling did not influence steady state, switching the limiting nutrient was only done after another five residence times. The transition between the two limitations was frequently sampled for OD and metabolite data (HPLC). Sampling of the second steady state was performed again after five residence times. Furthermore, during the culture, dissolved oxygen, pH, stirrer rate, temperature, airflow and oxygen and carbon dioxide content of the off-gas were continuously logged. Two balances, one under the influent and one under the effluent barrel, were coupled to the system to allow precise measurement of the dilution rate.


### 2.3. Measurements

The sampling and measurement methodologies used are described in [25].


### 2.4. Data Analysis

Due to the amount of data that had to be combined for each experiment, a custom-made program was written in Python using the SciPy scientific library [26, 27] to process the data and apply the algorithm as explained in Section 3 (polynomial fitting, extracting the derivatives, and calculating the fluxes).Each point in the time series of the transient experiments is only measured once. This is deemed sufficient, as nonsystematical variations are smoothed out when approximating the data with polynomials (as described in Section3.2). In this polynomial approximation, error propagation calculations are not performed. For the MFA calculations, however, a covariance matrix is needed. This matrix is used to rescale the measurements so that fluxes that are known to be inaccurate do not weight as much during the data reconciliation as other, more trustworthy fluxes [1]. Therefore, a “mean” covariance matrix was calculated from data presented in [25] and used for the MFA part of this analysis. The experiments described therein are similar to those performed for this publication (also using the same equipment) and it is assumed that the utilised covariance matrix is a good approximation of the real one. The standard deviation flags in the figures are based on this covariance matrix.The MFA modelling was performed as described in [4], although error analysis on the reconciliated fluxes was not done. The model used is based on the one in [4] and is the same metabolic model as described in [25]. It contains 136 reactions (Table 1) and 150 metabolites (Table 2), of which 12 are exchangeable. There are 142 independent equations and 136+12=148 unknowns. Thus at least 6 measurements have to be performed in order to fully solve the model. Ten exchange metabolites were measured and used for solving the stoichiometric models (CO2,O2,NH3, PiOH, acetate, lactate, pyruvate, succinate, glucose, and biomass) giving 4 redundant measurements that can be balanced. The methods for balancing measurements were initially developed for black box models [1, 2, 28], but can easily be applied on stoichiometric models [3, 4].Table 1
List of reactions used in the metabolic model.
Name of reactionThe reactionPGI:G6P↔F6PPFK:ATP+F6P→ADP+FBPALD:FBP↔G3P+DHAPTPI:DHAP↔G3PG3PDH:PiOH+NAD+G3P↔NADH+H+BPGPGK:ADP+BPG↔ATP+3PGPGM:3PG↔2PGENO:2PG↔H2O+PEPPyrK:ADP+PEP→ATP+PyrPyrD:NAD+Pyr+CoA→NADH+H+AcCoA+CO2CitSY:H2O+AcCoA+OAA→CoA+CitACO:Cit↔iCitCitDH:NAD+iCit↔NADH+H+CO2+aKGAAKGDH:NAD+CoA+aKGA→NADH+H+CO2+SucCoASucCoASY:ADP+PiOH+SucCoA↔ATP+CoA+SucSucDH:FAD+Suc→FADH2+FumFumHY:H2O+Fum↔MalMalDH:NAD+Mal↔NADH+H+OAAPEPCB:H2O+PEP+CO2→PiOH+OAALacDH:NADH+H+Pyr↔NAD+LacPFLY:Pyr+CoA→AcCoA+FAAcKNLR:ADP+PiOH+AcCoA↔ATP+CoA+AcResp:1.33ADP+1.33PiOH+NADH+H+0.5O2→1.33ATP+NAD+2.33H2OH2CO3SY:H2O+CO2↔H2CO3G6PDH:NADP+G6P→NADPH+H+6PGLLAS:H2O+6PGL→6PGPGDH:NADP+6PG→NADPH+H+CO2+Rl5PPPI:Rl5P↔R5PPPE:Rl5P↔Xu5PTK1:R5P+Xu5P↔G3P+S7PTA:G3P+S7P↔F6P+E4PTK2:Xu5P+E4P↔F6P+G3PPTS:GLC+PEP→G6P+PyrPPiOHHY:PPiOH+H2O→2PiOHGluDH:NADPH+H+aKGA+NH3↔NADP+H2O+GluGluLI:ATP+NH3+Glu→ADP+PiOH+GlnAspSY:ATP+H2O+Asp+Gln→AMP+PPiOH+Asn+GluAspTA:OAA+Glu↔aKGA+AspAlaTA:Pyr+Glu↔aKGA+AlaValAT:aKIV+Glu↔aKGA+ValLeuSYLR:NAD+H2O+AcCoA+aKIV+Glu→NADH+H+CoA+CO2+aKGA+LeuaKIVSYLR:NADPH+H+2Pyr→NADP+H2O+CO2+aKIVIleSYLR:NADPH+H+Pyr+Glu+Thr→NADP+H2O+CO2+aKGA+NH3+IleProSYLR:ATP+2NADPH+2H+Glu→ADP+PiOH+2NADP+H2O+ProSerLR:NAD+H2O+3PG+Glu→PiOH+NADH+H+aKGA+SerSerTHM:Ser+THF→H2O+Gly+MeTHFH2SSYLR:2ATP+3NADPH+ThioredH2+3H+H2SO4→ADP+PPiOH+3NADP+Thiored+3H2O+H2S+PAPPAPNAS:H2O+PAP→AMP+PiOHCysSYLR:H2S+AcCoA+Ser→CoA+Cys+AcPrppSY:ATP+R5P→AMP+PRPPHisSYLR:ATP+2NAD+3H2O+Gln+PRPP→2PPiOH+PiOH+2NADH+2H+aKGA+His+AICARPheSYLR:Glu+Chor→H2O+CO2+aKGA+PheTyrSYLR:NAD+Glu+Chor→NADH+H+CO2+aKGA+TyrTrpSYLR:Gln+Ser+Chor+PRPP→PPiOH+2H2O+G3P+Pyr+CO2+Glu+TrpDhDoPHepAD:H2O+PEP+E4P→PiOH+DahpDhqSY:Dahp→PiOH+DhqDhsSYLR:Dhq→H2O+DhsShiSY:NADPH+H+Dhs↔NADP+ShiShiKN:ATP+Shi→ADP+Shi3PChorSYLR:PEP+Shi3P→2PiOH+ChorThrSYLR:ATP+H2O+HSer→ADP+PiOH+ThrMDAPSYLR:NADPH+H+Pyr+SucCoA+Glu+AspSA→NADP+CoA+aKGA+Suc+MDAPLysSY:MDAP→CO2+LysMetSYLR:H2O+SucCoA+Cys+MTHF+HSer→Pyr+CoA+Suc+NH3+Met+THFAspSASY:ATP+NADPH+H+Asp→ADP+PiOH+NADP+AspSAHSerDH:NADPH+H+AspSA↔NADP+HSerCarPSY:2ATP+H2O+H2CO3+Gln→2ADP+PiOH+Glu+CarPOrnSYLR:ATP+NADPH+H+H2O+AcCoA+2Glu→ADP+PiOH+NADP+CoA+aKGA+Orn+AcArgSYLR:ATP+Asp+Orn+CarP→AMP+PPiOH+PiOH+Fum+ArgThioredRD:NADPH+Thiored+H↔NADP+ThioredH2H2O2ox:2H2O2→2H2O+O2FAD2NAD:NAD+FADH2↔NADH+FAD+HAICARSYLR:6ATP+3H2O+CO2+Asp+2Gln+Gly+FA+PRPP→6ADP+PPiOH+6PiOH+Fum+2Glu+AICARIMPSYLR:FTHF+AICAR→H2O+THF+IMPAMPSYLR:Asp+GTP+IMP→AMP+PiOH+Fum+GDPAdKN:ATP+AMP↔2ADPADPRD:ADP+ThioredH2→Thiored+H2O+dADPdADPKN:ATP+dADP→ADP+dATPIMPDH:NAD+H2O+IMP→NADH+H+XMPGMPSY:ATP+H2O+Gln+XMP→AMP+PPiOH+Glu+GMPGuKN:ATP+GMP→ADP+GDPGDPKN:ATP+GDP→ADP+GTPGDPRD:ThioredH2+GDP→Thiored+H2O+dGDPdGDPKN:ATP+dGDP→ADP+dGTPUMPSYLR:O2+Asp+PRPP+CarP→PPiOH+PiOH+H2O+CO2+UMP+H2O2UrKN:ATP+UMP→ADP+UDPUDPKN:ATP+UDP→ADP+UTPCTPSY:ATP+H2O+Gln+UTP→ADP+PiOH+Glu+CTPCDPKN:ATP+CDP↔ADP+CTPCMPKN:ATP+CMP→ADP+CDPCDPRD:ThioredH2+CDP→Thiored+H2O+dCDPdCDPKN:ATP+dCDP→ADP+dCTPUDPRD:ThioredH2+UDP→Thiored+H2O+dUDPdUDPKN:ATP+dUDP→ADP+dUTPdUTPPPAS:H2O+dUTP→PPiOH+dUMPdTMPSY:MeTHF+dUMP→DHF+dTMPdTMPKN:ATP+dTMP→ADP+dTDPdTDPKN:ATP+dTDP→ADP+dTTPDHFRD:NADPH+H+DHF→NADP+THFFTHFSYLR:NADP+H2O+MeTHF→NADPH+H+FTHFGlyCA:NAD+Gly+THF↔NADH+H+CO2+NH3+MeTHFMeTHFRD:NADH+H+MeTHF→NAD+MTHFAcCoACB:ATP+H2O+AcCoA+CO2↔ADP+PiOH+MalCoAMalCoATA:MalCoA+ACP↔CoA+MalACPAcACPSY:MalACP→CO2+AcACPC120SY:10NADPH+10H+AcACP+5MalACP→10NADP+5H2O+5CO2+C120ACP+5ACPC140SY:12NADPH+12H+AcACP+6MalACP→12NADP+6H2O+6CO2+C140ACP+6ACPC160SY:14NADPH+14H+AcACP+7MalACP→14NADP+7H2O+7CO2+C160ACP+7ACPC181SY:15NADPH+15H+AcACP+8MalACP→15NADP+8H2O+8CO2+C181ACP+8ACPAcylTF:C160ACP+C181ACP+Go3P→2ACP+PAGo3PDH:NADPH+H+DHAP↔NADP+Go3PDGoKN:ATP+DGo→ADP+PACDPDGoSY:CTP+PA↔PPiOH+CDPDGoPSerSY:Ser+CDPDGo→CMP+PSerPSerDC:PSer→CO2+PEthAnGlnF6PTA:F6P+Gln→Glu+GA6PGlcAnMU:GA6P↔GA1PNAGUrTF:AcCoA+UTP+GA1P→PPiOH+CoA+UDPNAGLipaSYLR:ATP+2CMPKDO+2UDPNAG+C120ACP+5C140ACP→ADP+2CMP+UMP+UDP+6ACP+Lipa+2AcA5PIR:Rl5P↔Ar5PPGLCMT:G6P↔G1PCMPKDOSYLR:2H2O+PEP+Ar5P+CTP→PPiOH+2PiOH+CMPKDOADPHEPSY:ATP+S7P→PPiOH+ADPHEPUDPGlcSY:G1P+UTP→PPiOH+UDPGlcEthANPT:CMP+PEthAn↔CDPEthAn+DGoLpsSYLR:3ADPHEP+3CMPKDO+2UDPGlc+Lipa+2CDPEthAn→3ADP+3CMP+2CDP+2UDP+LpsPGSYLR:Go3P+CDPDGo→PiOH+CMP+PGCLSY:PG+CDPDGo→CMP+CLPeptidoSYLR:5ATP+NADPH+H+PEP+3Ala+MDAP+2UDPNAG→5ADP+7PiOH+NADP+UMP+UDP+PeptidoGlcgSY:ATP+G1P→ADP+PPiOH+GlcgATPHY:ATP+H2O→ADP→PiOHDNASYLR:2H2O+0.246dATP+0.254dGTP+0.254dCTP+0.246dTTP→2PiOH+DNARNASYLR:0.262ATP+2H2O+0.322GTP+0.2CTP+0.216UTP→2PiOH+RNAProtSYLR:2ATP+3H2O+0.0961Ala+0.05506Arg+0.04505Asn+0.04505Asp+0.01702Cys+0.04905Gln+0.04905Glu+0.1151Gly+0.01802His+0.05405Ile+0.08408Leu+0.06406Lys+0.02903Met+0.03504Phe+0.04104Pro+0.04004Ser+0.04705Thr+0.01101Trp+0.02603Tyr+0.07908Val+2GTP→2ADP+4PiOH+2GDP+ProtLipidSYLR:0.0266CL+0.202PG+0.7714PEthAn→LipidBiomSYLR:0.004561Glcg+0.0002663Lps+0.0008933Peptido+0.002291DNA+0.01446RNA+0.1227Prot+0.003642Lipid→BiomTable 2
List of metabolites used in the metabolic model.
MetaboliteMolecular formulaFull name2PGC3H7O7P2-phophoglycerate3PGC3H7O7P3-phophoglycerate6PGC6H13O10P6-phosphogluconate6PGLC6H11O9P6-phosphogluconolactoneAcC2H4O2AcetateAcACPC2H3OPeptAcetyl ACPAcCoAC23H34O17N7P3SAcetyl CoAACPHPeptAcyl carrier proteinADPC10H15O10N5P2Adenosine diphosphateADPHEPC17H27O16N5P2ADP-MannoheptoseAICARC9H15O8N4PAmino imidazole carboxamide ribonucleotideaKGAC5H6O5Alpha ketoglutaric acidaKIVC5H8O3Alpha-keto-isovalerateAlaC3H7O2NAlanineAMPC10H14O7N5PAdenosine monophosphateAr5PC5H11O8PArabinose-5-phosphateArgC6H14O2N4ArginineAsnC4H8O3N2AspartateAspC4H7O4NAsparagineAspSAC4H7O3NAspartate semialdehydeATPC10H16O13N5P3Adenosine triphosphateBiomCH1.63O0.392N0.244P0.021S0.00565BiomassBPGC3H8O10P21-3-biphosphoglycerateC120ACPC12H23OPeptC140ACPC14H27OPeptC160ACPC16H31OPeptC181ACPC18H33OPeptCarPCH4O5NPCarbamoyl phosphateCDPC9H15O11N3P2Citidine diphosphateCDPDGoC46H83O15N3P2CDP-diacylglycerolCDPEthAnC11H20O11N4P2CDP-ethanolamineChorC10H10O6ChorismateCitC6H8O7cisaconitateCLC77H144O16P2CardiolipinCMPC9H14O8N3PCitidine monophosphateCMPKDOC17H26O15N3PCMP-2-keto-3-deoxyoctanoateCO2CO2Carbon dioxideCoAC21H32O16N7P3SCoenzyme ACTPC9H16O14N3P3Citidine triphosphateCysC3H7O2NSCysteinedADPC10H15O9N5P2deoxy ADPDahpC7H13O10PDeoxy arabino heptulosonatedATPC10H16O12N5P3deoxy ATPdCDPC9H15O10N3P2deoxy CDPdCTPC9H16O13N3P3deoxy CTPdGDPC10H15O10N5P2deoxy GDPDGoC37H70O5Diacyl glyceroldGTPC10H16O13N5P3deoxy GTPDHAPC3H7O6PDihydroxyaceton phosphateDHFC19H21O6N7DihydrofolateDhqC7H10O6DehydroquinateDhsC7H8O5DehydroshikimateDNAC9.75H14.2O7N3.75PDNA compositiondTDPC10H16O11N2P2deoxy TDPdTMPC10H15O8N2Pdeoxy TMPdTTPC10H17O14N2P3deoxy TTPdUDPC9H14O11N2P2deoxy UDPdUMPC9H13O8N2Pdeoxy UMPdUTPC9H15O14N2P3deoxy UTPE4PC4H9O7PErythrose-4-phosphateF6PC6H13O9PFructose-6-phosphateFACH2O2Formic AcidFADC27H33O15N9P2Flavine adeninen dinucleotideFADH2C27H35O15N9P2FBPC6H14O12P2Fructose-1-6-biphosphateFTHFC20H23O7N7Formyl tetrahydrofolateFumC4H4O4FumarateG1PC6H13O9PGlucose-1-phosphateG3PC3H7O6PGlyceraldehyde-3-phosphateG6PC6H13O9PGlucose-6-phosphateGA1PC6H14O8NPD-glucosamine-6-phosphateGA6PC6H14O8NPGlucose-6-phosphateGDPC10H15O11N5P2Guanosine diphosphateGLCC6H12O6GlucoseGlcgC6H10O5GlycogenGlnC5H10O3N2GlutamineGluC5H9O4NGlutamateGlyC2H5O2NGlycineGMPC10H14O8N5PGuanosine monophosphateGo3PC3H9O6PGlycerol-3-phosphateGTPC10H16O14N5P3Guanosine triphosphateHH+HydrogeneH2CO3CH2O3BicarbonateH2OH2OWaterH2O2H2O2H2SH2SHydrogene suldeH2SO4H2O4SSulfuric acidHisC6H9O2N3HistidineHSerC4H9O3NHomoserineiCitC6H8O7IsocitratIleC6H13O2NIsoleucineIMPC10H13O8N4PInosine monophosphateLacC3H6O3LactateLeuC6H13O2NLeucineLipaC110H196O32N2P2Lipid ALipidC40.2H77.6O8.41N0.771P1.03Lipid compositionLpsC171H298O81N4P2Lipo Poly saccharideLysC6H14O2N2LysineMalC4H6O5MalateMalACPC3H3O3PeptMalonyl ACPMalCoAC24H34O19N7P3SMalonyl CoAMDAPC7H14O4N2Meso-diaminopimelateMetC5H11O2NSMethionineMeTHFC20H23O6N7Methylene tetrahydro folateMTHFC20H25O6N7Methyl tetrahydrofolateNADC21H28O14N7P2+Nicotinamide adenine dinucleotideNADHC21H29O14N7P2NADPC21H28O17N7P3+Nicotinamide adenine dinucleotide phosphateNADPHC21H29O17N7P3NH3H3NAmmoniaO2O2OxygenOAAC4H4O5OxaloacetateOrnC5H12O2N2OrnithinePAC37H71O8PPhosphatidyl acidPAPC10H15O10N5P2Phospho adenosine phosphatePEPC3H5O6PPhosphoenolpyruvatePeptidoC35H53O16N7PeptidoglycanPEthAnC39H76O8NPPhosphatidyl ethanolaminePGC40H75O9PPhosphatidyl glycerolPheC9H11O2NPhenylalaninePiOHH3O4PPhosphatePPiOHH4O7P2PyrophosphateProC5H9O2NProlineProtC4.8H7.67O1.4N1.37S0.046Protein compositionPRPPC5H13O14P35-phospho-alpha-D-ribosyl-1-pyrophosphatePSerC40H76O10NPPhosphatidyl SerinePyrC3H4O3PyruvateR5PC5H11O8PRibose-5-phosphateRl5PC5H11O8PRibulose-5-phosphateRNAC9.58H13.8O7.95N3.95PRNA compositionS7PC7H15O10PSedoheptulose-7-phosphateSerC3H7O3NSerineShiC7H10O5ShikimateShi3PC7H11O8PShikimate-3-phosphateSucC4H6O4SuccinateSucCoAC25H36O19N7P3SSuccinyl CoATHFC19H23O6N7TetrahydrofolateThioredPeptThioredoxinThioredH2H2PeptReduced thioredoxinThrC4H9O3NThreonineTrpC11H12O2N2TryptophanTyrC9H11O3NTyrosineUDPC9H14O12N2P2Uridine diphosphateUDPGlcC15H24O17N2P2UDP glucoseUDPNAGC17H27O17N3P2UDP N-acetyl glucosamineUMPC9H13O9N2PUridine monophosphateUTPC9H15O15N2P3Uridine triphosphateValC5H11O2NValineXMPC10H13O9N4PXanthosine-5-phosphateXu5PC5H11O8PXylulose-5-phosphateAlthough it is known that the biomass composition is different under C- and N-limitation [29, 30], the same model was used in both cases thus implying the same biomass composition for cells grown under N- and C-limitation. Limited testing showed that the influence of a different protein content of the biomass on the flux fluctuations was minimal.


## 2.1. Bacterial Strain

Escherichia coli MG1655 [λ-,F-,  rph–1, (del  fnr)] was obtained from the Netherlands Culture Collection of Bacteria (NCCB, Utrecht, The Netherlands).


## 2.2. Culture Conditions

### 2.2.1. Media

Composition of the Luria Bertani Broth (LB) medium and shake flask medium (for the preculture) can be found in [25].The minimal medium used in the reactor consisted of 2.5 g/L (N-limited medium) or 5 g/L (C-limited medium)(NH4)2SO4, 2 g/L KH2PO4, 0.5 g/L NaCl, 0.5 g/L MgSO4·7H2O, 33 g/L (N-limited medium) or 16.5 g/L (C-limited medium) Glucose·H2O, 1 mL/L vitamin stock solution and 100 μL/L molybdate stock solution. Vitamin stock solution consisted of 3.6 g/L FeCl2·4H2O, 5 g/L CaCl2·2H2O, 1.3 g/L MnCl2·2H2O, 0.38 g/L CuCl2·2H2O, 0.5 g/L CoCl2·6H2O, 0.94 g/L ZnCl2, 0.0311 g/L H3BO4, 0.4 g/L Na2EDTA·2H2O, and 1.01 g/L thiamine·HCl. The molybdate stock solution contained 0.967 g/L Na2MoO4·2H2O. All components were dissolved and filter-sterilised (pore size 0.22 μm, Sartobran, Sartorius, Belgium). The pH was left at approximately 5.4.


### 2.2.2. Cultivations

Culture conditions and inoculation procedure are described in [25].The experiment in which the limiting substrate was changed from glucose to ammonia was conducted at a dilution rate of 0.155 h-1. The experiment in which the limitation was changed from nitrogen to carbon was conducted at a dilution rate of 0.142 h-1.For each experiment, steady state samples were taken: optical density (OD), cell dry weight (CDW), and HPLC analyses. These samples were taken after waiting at least five residence times after the batch phase. To ensure that the perturbations caused by sampling did not influence steady state, switching the limiting nutrient was only done after another five residence times. The transition between the two limitations was frequently sampled for OD and metabolite data (HPLC). Sampling of the second steady state was performed again after five residence times. Furthermore, during the culture, dissolved oxygen, pH, stirrer rate, temperature, airflow and oxygen and carbon dioxide content of the off-gas were continuously logged. Two balances, one under the influent and one under the effluent barrel, were coupled to the system to allow precise measurement of the dilution rate.


## 2.2.1. Media

Composition of the Luria Bertani Broth (LB) medium and shake flask medium (for the preculture) can be found in [25].The minimal medium used in the reactor consisted of 2.5 g/L (N-limited medium) or 5 g/L (C-limited medium)(NH4)2SO4, 2 g/L KH2PO4, 0.5 g/L NaCl, 0.5 g/L MgSO4·7H2O, 33 g/L (N-limited medium) or 16.5 g/L (C-limited medium) Glucose·H2O, 1 mL/L vitamin stock solution and 100 μL/L molybdate stock solution. Vitamin stock solution consisted of 3.6 g/L FeCl2·4H2O, 5 g/L CaCl2·2H2O, 1.3 g/L MnCl2·2H2O, 0.38 g/L CuCl2·2H2O, 0.5 g/L CoCl2·6H2O, 0.94 g/L ZnCl2, 0.0311 g/L H3BO4, 0.4 g/L Na2EDTA·2H2O, and 1.01 g/L thiamine·HCl. The molybdate stock solution contained 0.967 g/L Na2MoO4·2H2O. All components were dissolved and filter-sterilised (pore size 0.22 μm, Sartobran, Sartorius, Belgium). The pH was left at approximately 5.4.


## 2.2.2. Cultivations

Culture conditions and inoculation procedure are described in [25].The experiment in which the limiting substrate was changed from glucose to ammonia was conducted at a dilution rate of 0.155 h-1. The experiment in which the limitation was changed from nitrogen to carbon was conducted at a dilution rate of 0.142 h-1.For each experiment, steady state samples were taken: optical density (OD), cell dry weight (CDW), and HPLC analyses. These samples were taken after waiting at least five residence times after the batch phase. To ensure that the perturbations caused by sampling did not influence steady state, switching the limiting nutrient was only done after another five residence times. The transition between the two limitations was frequently sampled for OD and metabolite data (HPLC). Sampling of the second steady state was performed again after five residence times. Furthermore, during the culture, dissolved oxygen, pH, stirrer rate, temperature, airflow and oxygen and carbon dioxide content of the off-gas were continuously logged. Two balances, one under the influent and one under the effluent barrel, were coupled to the system to allow precise measurement of the dilution rate.


## 2.3. Measurements

The sampling and measurement methodologies used are described in [25].


## 2.4. Data Analysis

Due to the amount of data that had to be combined for each experiment, a custom-made program was written in Python using the SciPy scientific library [26, 27] to process the data and apply the algorithm as explained in Section 3 (polynomial fitting, extracting the derivatives, and calculating the fluxes).Each point in the time series of the transient experiments is only measured once. This is deemed sufficient, as nonsystematical variations are smoothed out when approximating the data with polynomials (as described in Section3.2). In this polynomial approximation, error propagation calculations are not performed. For the MFA calculations, however, a covariance matrix is needed. This matrix is used to rescale the measurements so that fluxes that are known to be inaccurate do not weight as much during the data reconciliation as other, more trustworthy fluxes [1]. Therefore, a “mean” covariance matrix was calculated from data presented in [25] and used for the MFA part of this analysis. The experiments described therein are similar to those performed for this publication (also using the same equipment) and it is assumed that the utilised covariance matrix is a good approximation of the real one. The standard deviation flags in the figures are based on this covariance matrix.The MFA modelling was performed as described in [4], although error analysis on the reconciliated fluxes was not done. The model used is based on the one in [4] and is the same metabolic model as described in [25]. It contains 136 reactions (Table 1) and 150 metabolites (Table 2), of which 12 are exchangeable. There are 142 independent equations and 136+12=148 unknowns. Thus at least 6 measurements have to be performed in order to fully solve the model. Ten exchange metabolites were measured and used for solving the stoichiometric models (CO2,O2,NH3, PiOH, acetate, lactate, pyruvate, succinate, glucose, and biomass) giving 4 redundant measurements that can be balanced. The methods for balancing measurements were initially developed for black box models [1, 2, 28], but can easily be applied on stoichiometric models [3, 4].Table 1
List of reactions used in the metabolic model.
Name of reactionThe reactionPGI:G6P↔F6PPFK:ATP+F6P→ADP+FBPALD:FBP↔G3P+DHAPTPI:DHAP↔G3PG3PDH:PiOH+NAD+G3P↔NADH+H+BPGPGK:ADP+BPG↔ATP+3PGPGM:3PG↔2PGENO:2PG↔H2O+PEPPyrK:ADP+PEP→ATP+PyrPyrD:NAD+Pyr+CoA→NADH+H+AcCoA+CO2CitSY:H2O+AcCoA+OAA→CoA+CitACO:Cit↔iCitCitDH:NAD+iCit↔NADH+H+CO2+aKGAAKGDH:NAD+CoA+aKGA→NADH+H+CO2+SucCoASucCoASY:ADP+PiOH+SucCoA↔ATP+CoA+SucSucDH:FAD+Suc→FADH2+FumFumHY:H2O+Fum↔MalMalDH:NAD+Mal↔NADH+H+OAAPEPCB:H2O+PEP+CO2→PiOH+OAALacDH:NADH+H+Pyr↔NAD+LacPFLY:Pyr+CoA→AcCoA+FAAcKNLR:ADP+PiOH+AcCoA↔ATP+CoA+AcResp:1.33ADP+1.33PiOH+NADH+H+0.5O2→1.33ATP+NAD+2.33H2OH2CO3SY:H2O+CO2↔H2CO3G6PDH:NADP+G6P→NADPH+H+6PGLLAS:H2O+6PGL→6PGPGDH:NADP+6PG→NADPH+H+CO2+Rl5PPPI:Rl5P↔R5PPPE:Rl5P↔Xu5PTK1:R5P+Xu5P↔G3P+S7PTA:G3P+S7P↔F6P+E4PTK2:Xu5P+E4P↔F6P+G3PPTS:GLC+PEP→G6P+PyrPPiOHHY:PPiOH+H2O→2PiOHGluDH:NADPH+H+aKGA+NH3↔NADP+H2O+GluGluLI:ATP+NH3+Glu→ADP+PiOH+GlnAspSY:ATP+H2O+Asp+Gln→AMP+PPiOH+Asn+GluAspTA:OAA+Glu↔aKGA+AspAlaTA:Pyr+Glu↔aKGA+AlaValAT:aKIV+Glu↔aKGA+ValLeuSYLR:NAD+H2O+AcCoA+aKIV+Glu→NADH+H+CoA+CO2+aKGA+LeuaKIVSYLR:NADPH+H+2Pyr→NADP+H2O+CO2+aKIVIleSYLR:NADPH+H+Pyr+Glu+Thr→NADP+H2O+CO2+aKGA+NH3+IleProSYLR:ATP+2NADPH+2H+Glu→ADP+PiOH+2NADP+H2O+ProSerLR:NAD+H2O+3PG+Glu→PiOH+NADH+H+aKGA+SerSerTHM:Ser+THF→H2O+Gly+MeTHFH2SSYLR:2ATP+3NADPH+ThioredH2+3H+H2SO4→ADP+PPiOH+3NADP+Thiored+3H2O+H2S+PAPPAPNAS:H2O+PAP→AMP+PiOHCysSYLR:H2S+AcCoA+Ser→CoA+Cys+AcPrppSY:ATP+R5P→AMP+PRPPHisSYLR:ATP+2NAD+3H2O+Gln+PRPP→2PPiOH+PiOH+2NADH+2H+aKGA+His+AICARPheSYLR:Glu+Chor→H2O+CO2+aKGA+PheTyrSYLR:NAD+Glu+Chor→NADH+H+CO2+aKGA+TyrTrpSYLR:Gln+Ser+Chor+PRPP→PPiOH+2H2O+G3P+Pyr+CO2+Glu+TrpDhDoPHepAD:H2O+PEP+E4P→PiOH+DahpDhqSY:Dahp→PiOH+DhqDhsSYLR:Dhq→H2O+DhsShiSY:NADPH+H+Dhs↔NADP+ShiShiKN:ATP+Shi→ADP+Shi3PChorSYLR:PEP+Shi3P→2PiOH+ChorThrSYLR:ATP+H2O+HSer→ADP+PiOH+ThrMDAPSYLR:NADPH+H+Pyr+SucCoA+Glu+AspSA→NADP+CoA+aKGA+Suc+MDAPLysSY:MDAP→CO2+LysMetSYLR:H2O+SucCoA+Cys+MTHF+HSer→Pyr+CoA+Suc+NH3+Met+THFAspSASY:ATP+NADPH+H+Asp→ADP+PiOH+NADP+AspSAHSerDH:NADPH+H+AspSA↔NADP+HSerCarPSY:2ATP+H2O+H2CO3+Gln→2ADP+PiOH+Glu+CarPOrnSYLR:ATP+NADPH+H+H2O+AcCoA+2Glu→ADP+PiOH+NADP+CoA+aKGA+Orn+AcArgSYLR:ATP+Asp+Orn+CarP→AMP+PPiOH+PiOH+Fum+ArgThioredRD:NADPH+Thiored+H↔NADP+ThioredH2H2O2ox:2H2O2→2H2O+O2FAD2NAD:NAD+FADH2↔NADH+FAD+HAICARSYLR:6ATP+3H2O+CO2+Asp+2Gln+Gly+FA+PRPP→6ADP+PPiOH+6PiOH+Fum+2Glu+AICARIMPSYLR:FTHF+AICAR→H2O+THF+IMPAMPSYLR:Asp+GTP+IMP→AMP+PiOH+Fum+GDPAdKN:ATP+AMP↔2ADPADPRD:ADP+ThioredH2→Thiored+H2O+dADPdADPKN:ATP+dADP→ADP+dATPIMPDH:NAD+H2O+IMP→NADH+H+XMPGMPSY:ATP+H2O+Gln+XMP→AMP+PPiOH+Glu+GMPGuKN:ATP+GMP→ADP+GDPGDPKN:ATP+GDP→ADP+GTPGDPRD:ThioredH2+GDP→Thiored+H2O+dGDPdGDPKN:ATP+dGDP→ADP+dGTPUMPSYLR:O2+Asp+PRPP+CarP→PPiOH+PiOH+H2O+CO2+UMP+H2O2UrKN:ATP+UMP→ADP+UDPUDPKN:ATP+UDP→ADP+UTPCTPSY:ATP+H2O+Gln+UTP→ADP+PiOH+Glu+CTPCDPKN:ATP+CDP↔ADP+CTPCMPKN:ATP+CMP→ADP+CDPCDPRD:ThioredH2+CDP→Thiored+H2O+dCDPdCDPKN:ATP+dCDP→ADP+dCTPUDPRD:ThioredH2+UDP→Thiored+H2O+dUDPdUDPKN:ATP+dUDP→ADP+dUTPdUTPPPAS:H2O+dUTP→PPiOH+dUMPdTMPSY:MeTHF+dUMP→DHF+dTMPdTMPKN:ATP+dTMP→ADP+dTDPdTDPKN:ATP+dTDP→ADP+dTTPDHFRD:NADPH+H+DHF→NADP+THFFTHFSYLR:NADP+H2O+MeTHF→NADPH+H+FTHFGlyCA:NAD+Gly+THF↔NADH+H+CO2+NH3+MeTHFMeTHFRD:NADH+H+MeTHF→NAD+MTHFAcCoACB:ATP+H2O+AcCoA+CO2↔ADP+PiOH+MalCoAMalCoATA:MalCoA+ACP↔CoA+MalACPAcACPSY:MalACP→CO2+AcACPC120SY:10NADPH+10H+AcACP+5MalACP→10NADP+5H2O+5CO2+C120ACP+5ACPC140SY:12NADPH+12H+AcACP+6MalACP→12NADP+6H2O+6CO2+C140ACP+6ACPC160SY:14NADPH+14H+AcACP+7MalACP→14NADP+7H2O+7CO2+C160ACP+7ACPC181SY:15NADPH+15H+AcACP+8MalACP→15NADP+8H2O+8CO2+C181ACP+8ACPAcylTF:C160ACP+C181ACP+Go3P→2ACP+PAGo3PDH:NADPH+H+DHAP↔NADP+Go3PDGoKN:ATP+DGo→ADP+PACDPDGoSY:CTP+PA↔PPiOH+CDPDGoPSerSY:Ser+CDPDGo→CMP+PSerPSerDC:PSer→CO2+PEthAnGlnF6PTA:F6P+Gln→Glu+GA6PGlcAnMU:GA6P↔GA1PNAGUrTF:AcCoA+UTP+GA1P→PPiOH+CoA+UDPNAGLipaSYLR:ATP+2CMPKDO+2UDPNAG+C120ACP+5C140ACP→ADP+2CMP+UMP+UDP+6ACP+Lipa+2AcA5PIR:Rl5P↔Ar5PPGLCMT:G6P↔G1PCMPKDOSYLR:2H2O+PEP+Ar5P+CTP→PPiOH+2PiOH+CMPKDOADPHEPSY:ATP+S7P→PPiOH+ADPHEPUDPGlcSY:G1P+UTP→PPiOH+UDPGlcEthANPT:CMP+PEthAn↔CDPEthAn+DGoLpsSYLR:3ADPHEP+3CMPKDO+2UDPGlc+Lipa+2CDPEthAn→3ADP+3CMP+2CDP+2UDP+LpsPGSYLR:Go3P+CDPDGo→PiOH+CMP+PGCLSY:PG+CDPDGo→CMP+CLPeptidoSYLR:5ATP+NADPH+H+PEP+3Ala+MDAP+2UDPNAG→5ADP+7PiOH+NADP+UMP+UDP+PeptidoGlcgSY:ATP+G1P→ADP+PPiOH+GlcgATPHY:ATP+H2O→ADP→PiOHDNASYLR:2H2O+0.246dATP+0.254dGTP+0.254dCTP+0.246dTTP→2PiOH+DNARNASYLR:0.262ATP+2H2O+0.322GTP+0.2CTP+0.216UTP→2PiOH+RNAProtSYLR:2ATP+3H2O+0.0961Ala+0.05506Arg+0.04505Asn+0.04505Asp+0.01702Cys+0.04905Gln+0.04905Glu+0.1151Gly+0.01802His+0.05405Ile+0.08408Leu+0.06406Lys+0.02903Met+0.03504Phe+0.04104Pro+0.04004Ser+0.04705Thr+0.01101Trp+0.02603Tyr+0.07908Val+2GTP→2ADP+4PiOH+2GDP+ProtLipidSYLR:0.0266CL+0.202PG+0.7714PEthAn→LipidBiomSYLR:0.004561Glcg+0.0002663Lps+0.0008933Peptido+0.002291DNA+0.01446RNA+0.1227Prot+0.003642Lipid→BiomTable 2
List of metabolites used in the metabolic model.
MetaboliteMolecular formulaFull name2PGC3H7O7P2-phophoglycerate3PGC3H7O7P3-phophoglycerate6PGC6H13O10P6-phosphogluconate6PGLC6H11O9P6-phosphogluconolactoneAcC2H4O2AcetateAcACPC2H3OPeptAcetyl ACPAcCoAC23H34O17N7P3SAcetyl CoAACPHPeptAcyl carrier proteinADPC10H15O10N5P2Adenosine diphosphateADPHEPC17H27O16N5P2ADP-MannoheptoseAICARC9H15O8N4PAmino imidazole carboxamide ribonucleotideaKGAC5H6O5Alpha ketoglutaric acidaKIVC5H8O3Alpha-keto-isovalerateAlaC3H7O2NAlanineAMPC10H14O7N5PAdenosine monophosphateAr5PC5H11O8PArabinose-5-phosphateArgC6H14O2N4ArginineAsnC4H8O3N2AspartateAspC4H7O4NAsparagineAspSAC4H7O3NAspartate semialdehydeATPC10H16O13N5P3Adenosine triphosphateBiomCH1.63O0.392N0.244P0.021S0.00565BiomassBPGC3H8O10P21-3-biphosphoglycerateC120ACPC12H23OPeptC140ACPC14H27OPeptC160ACPC16H31OPeptC181ACPC18H33OPeptCarPCH4O5NPCarbamoyl phosphateCDPC9H15O11N3P2Citidine diphosphateCDPDGoC46H83O15N3P2CDP-diacylglycerolCDPEthAnC11H20O11N4P2CDP-ethanolamineChorC10H10O6ChorismateCitC6H8O7cisaconitateCLC77H144O16P2CardiolipinCMPC9H14O8N3PCitidine monophosphateCMPKDOC17H26O15N3PCMP-2-keto-3-deoxyoctanoateCO2CO2Carbon dioxideCoAC21H32O16N7P3SCoenzyme ACTPC9H16O14N3P3Citidine triphosphateCysC3H7O2NSCysteinedADPC10H15O9N5P2deoxy ADPDahpC7H13O10PDeoxy arabino heptulosonatedATPC10H16O12N5P3deoxy ATPdCDPC9H15O10N3P2deoxy CDPdCTPC9H16O13N3P3deoxy CTPdGDPC10H15O10N5P2deoxy GDPDGoC37H70O5Diacyl glyceroldGTPC10H16O13N5P3deoxy GTPDHAPC3H7O6PDihydroxyaceton phosphateDHFC19H21O6N7DihydrofolateDhqC7H10O6DehydroquinateDhsC7H8O5DehydroshikimateDNAC9.75H14.2O7N3.75PDNA compositiondTDPC10H16O11N2P2deoxy TDPdTMPC10H15O8N2Pdeoxy TMPdTTPC10H17O14N2P3deoxy TTPdUDPC9H14O11N2P2deoxy UDPdUMPC9H13O8N2Pdeoxy UMPdUTPC9H15O14N2P3deoxy UTPE4PC4H9O7PErythrose-4-phosphateF6PC6H13O9PFructose-6-phosphateFACH2O2Formic AcidFADC27H33O15N9P2Flavine adeninen dinucleotideFADH2C27H35O15N9P2FBPC6H14O12P2Fructose-1-6-biphosphateFTHFC20H23O7N7Formyl tetrahydrofolateFumC4H4O4FumarateG1PC6H13O9PGlucose-1-phosphateG3PC3H7O6PGlyceraldehyde-3-phosphateG6PC6H13O9PGlucose-6-phosphateGA1PC6H14O8NPD-glucosamine-6-phosphateGA6PC6H14O8NPGlucose-6-phosphateGDPC10H15O11N5P2Guanosine diphosphateGLCC6H12O6GlucoseGlcgC6H10O5GlycogenGlnC5H10O3N2GlutamineGluC5H9O4NGlutamateGlyC2H5O2NGlycineGMPC10H14O8N5PGuanosine monophosphateGo3PC3H9O6PGlycerol-3-phosphateGTPC10H16O14N5P3Guanosine triphosphateHH+HydrogeneH2CO3CH2O3BicarbonateH2OH2OWaterH2O2H2O2H2SH2SHydrogene suldeH2SO4H2O4SSulfuric acidHisC6H9O2N3HistidineHSerC4H9O3NHomoserineiCitC6H8O7IsocitratIleC6H13O2NIsoleucineIMPC10H13O8N4PInosine monophosphateLacC3H6O3LactateLeuC6H13O2NLeucineLipaC110H196O32N2P2Lipid ALipidC40.2H77.6O8.41N0.771P1.03Lipid compositionLpsC171H298O81N4P2Lipo Poly saccharideLysC6H14O2N2LysineMalC4H6O5MalateMalACPC3H3O3PeptMalonyl ACPMalCoAC24H34O19N7P3SMalonyl CoAMDAPC7H14O4N2Meso-diaminopimelateMetC5H11O2NSMethionineMeTHFC20H23O6N7Methylene tetrahydro folateMTHFC20H25O6N7Methyl tetrahydrofolateNADC21H28O14N7P2+Nicotinamide adenine dinucleotideNADHC21H29O14N7P2NADPC21H28O17N7P3+Nicotinamide adenine dinucleotide phosphateNADPHC21H29O17N7P3NH3H3NAmmoniaO2O2OxygenOAAC4H4O5OxaloacetateOrnC5H12O2N2OrnithinePAC37H71O8PPhosphatidyl acidPAPC10H15O10N5P2Phospho adenosine phosphatePEPC3H5O6PPhosphoenolpyruvatePeptidoC35H53O16N7PeptidoglycanPEthAnC39H76O8NPPhosphatidyl ethanolaminePGC40H75O9PPhosphatidyl glycerolPheC9H11O2NPhenylalaninePiOHH3O4PPhosphatePPiOHH4O7P2PyrophosphateProC5H9O2NProlineProtC4.8H7.67O1.4N1.37S0.046Protein compositionPRPPC5H13O14P35-phospho-alpha-D-ribosyl-1-pyrophosphatePSerC40H76O10NPPhosphatidyl SerinePyrC3H4O3PyruvateR5PC5H11O8PRibose-5-phosphateRl5PC5H11O8PRibulose-5-phosphateRNAC9.58H13.8O7.95N3.95PRNA compositionS7PC7H15O10PSedoheptulose-7-phosphateSerC3H7O3NSerineShiC7H10O5ShikimateShi3PC7H11O8PShikimate-3-phosphateSucC4H6O4SuccinateSucCoAC25H36O19N7P3SSuccinyl CoATHFC19H23O6N7TetrahydrofolateThioredPeptThioredoxinThioredH2H2PeptReduced thioredoxinThrC4H9O3NThreonineTrpC11H12O2N2TryptophanTyrC9H11O3NTyrosineUDPC9H14O12N2P2Uridine diphosphateUDPGlcC15H24O17N2P2UDP glucoseUDPNAGC17H27O17N3P2UDP N-acetyl glucosamineUMPC9H13O9N2PUridine monophosphateUTPC9H15O15N2P3Uridine triphosphateValC5H11O2NValineXMPC10H13O9N4PXanthosine-5-phosphateXu5PC5H11O8PXylulose-5-phosphateAlthough it is known that the biomass composition is different under C- and N-limitation [29, 30], the same model was used in both cases thus implying the same biomass composition for cells grown under N- and C-limitation. Limited testing showed that the influence of a different protein content of the biomass on the flux fluctuations was minimal.


## 3. Dynamic Metabolic Flux Analysis

### 3.1. DMFA: The Concept

In classical metabolic flux analysis, the measured fluxes are obtained by quantification of the metabolite concentrations in the reactor broth of a chemostat culture. But what should be done under transient conditions? Measurements of external metabolites cannot directly be used for fluxes. Provost and Bastin [31] applied MFA to the batch phase of the culture of Chinese Hamster Ovary (CHO) cells. The time series of the measured metabolite concentrations are approximated with a straight line, the slope of which is a direct measure of the fluxes entering and leaving the cells and is used to calculate internal reaction rates via MFA.This approach of Provost and Bastin [31] can be formalised and extended. Instead of limiting oneself to cases where a straight line is sufficient to adequately capture the dynamics of the measured metabolites during the transients, one can take the derivative in each point, without approximating the whole interval with a straight line. This derivative is, by definition, a flux.The differential equation that governs the change of a component in the reactor broth can be written as(1)dCdt=D(Cin-C)+rp,
where C is the concentration of the component, Cin the concentration of the component in the influent, rp the production rate, and D the dilution rate.From (1) the net reaction rate rp can be calculated using the following approximation for the derivative in a point n of the time series: (2)dCdt|n=Cn-Cn-1tn-tn-1.
Eventually, more data around the working point might be used to perform the linear regression. However, because of the noisiness of the data, this approach does not give satisfying results. Therefore, approximation techniques have to be used. In this work, a polynomial approximation was chosen.


### 3.2. Polynomial Fitting

The polynomial fitting algorithms found in SciPy [26] were used to approximate the data. The whole time series were covered by different polynomials. To this end, a moving window approach was used (Figure 1). The points inside the outer window, W1, are used for fitting the polynomial. The derivatives of the polynomial in the points of interest within in the inner window W2 are used for dynamic MFA. Making the inner window smaller than the outer window ensures a smooth transition between the different polynomials. However, if the inner window is too small compared to the outer window, the noise on the data is not filtered sufficiently, yielding nonrealistic fluctuations in the derivatives.Figure 1
A moving window is run through a time series of data.W1 is the polynomial fitting window, W2 is the interpolation window.Some of the available data series were impossible to smooth with polynomials. They behaved like a typical logistic curve. Hence, the logistic curve was added to the set of functions that could be used to approximate the data:(3)L(t)=p01+p1ep3t1+p2ep3t.
The parameters W1 and W2 and the maximal degree of the polynomials (or logistic curve) were manually selected by visually inspection of the generated curves. The actual degree of the polynomials could vary for each polynomial fitted (thus the curve for a time series of a certain metabolite consists of polynomials of different degrees chained together): for each W1 window the polynomial giving the least amount of outliers was selected. Data points were identified as outliers when the difference between the measured value and the estimated value was larger than 5% of the range of the response variable. Table 3 summarized the fitting parameters used for each metabolite.Table 3
Parameters for the polynomial fitting of the different time series.W1 is expressed in fraction of the time between the first measurement and the last measurement. W2 is expressed as a fraction of W1. PD is the highest degree a polynomial may have; an L is displayed if the data were fitted to a logistic curve.
C→N limitationN→C limitationW1W2PDW1W2PDGLC0.400.8750.350.453NH30.500.80211LPiOH0.700.7550.450.804O20.180.9070.110.907CO20.120.9040.120.907Ac0.700.8850.150.703Lac1170.410.807Pyr0.400.80611LSuc0.400.8060.820.995Biomass0.500.7080.370.804The polynomial or logistic functions yield a smoothed time profile of the concentrations and enable to calculate the derivative of the concentrations to time. However, when calculating this derivative directly by symbolically deriving the function and filling in the time for the point of interest, one sometimes obtains unrealistically high values (D1 in Figure 2). Therefore, the derivative was approximated by taking the slope of the point of interest and the previous point on the polynomial (D2 in Figure 2).Figure 2
Two ways of estimating the derivative in pointa: D1 is calculated by symbolically deriving the function in point a; D2 is the slope between point a and point b (where a and b are calculated from the curve). Dots represent sample points through which a curve is fitted.The derivative can subsequently be used in (1) to calculate the cellular production (or consumption) fluxes:
(4)rp=dCdt-D(Cin-C).
Once rp has been calculated, classical MFA techniques (as described in [4]) can be applied at each time point. The results of those models, evaluated at the different time points, can then be combined to obtain a time profile of internal fluxes.


### 3.3. Growth Rate Calculation

During the transient period, the assumption that the growth rate is equal to the dilution rate does not hold. However, (4) can explicitly be rewritten for the following biomass:
(5)rp,Biom=dCBiomdt+DCBiom
from which the growth rate can be calculated as (6)μ=rp,BiomCBiom=dCBiomdt1CBiom+D.


## 3.1. DMFA: The Concept

In classical metabolic flux analysis, the measured fluxes are obtained by quantification of the metabolite concentrations in the reactor broth of a chemostat culture. But what should be done under transient conditions? Measurements of external metabolites cannot directly be used for fluxes. Provost and Bastin [31] applied MFA to the batch phase of the culture of Chinese Hamster Ovary (CHO) cells. The time series of the measured metabolite concentrations are approximated with a straight line, the slope of which is a direct measure of the fluxes entering and leaving the cells and is used to calculate internal reaction rates via MFA.This approach of Provost and Bastin [31] can be formalised and extended. Instead of limiting oneself to cases where a straight line is sufficient to adequately capture the dynamics of the measured metabolites during the transients, one can take the derivative in each point, without approximating the whole interval with a straight line. This derivative is, by definition, a flux.The differential equation that governs the change of a component in the reactor broth can be written as(1)dCdt=D(Cin-C)+rp,
where C is the concentration of the component, Cin the concentration of the component in the influent, rp the production rate, and D the dilution rate.From (1) the net reaction rate rp can be calculated using the following approximation for the derivative in a point n of the time series: (2)dCdt|n=Cn-Cn-1tn-tn-1.
Eventually, more data around the working point might be used to perform the linear regression. However, because of the noisiness of the data, this approach does not give satisfying results. Therefore, approximation techniques have to be used. In this work, a polynomial approximation was chosen.


## 3.2. Polynomial Fitting

The polynomial fitting algorithms found in SciPy [26] were used to approximate the data. The whole time series were covered by different polynomials. To this end, a moving window approach was used (Figure 1). The points inside the outer window, W1, are used for fitting the polynomial. The derivatives of the polynomial in the points of interest within in the inner window W2 are used for dynamic MFA. Making the inner window smaller than the outer window ensures a smooth transition between the different polynomials. However, if the inner window is too small compared to the outer window, the noise on the data is not filtered sufficiently, yielding nonrealistic fluctuations in the derivatives.Figure 1
A moving window is run through a time series of data.W1 is the polynomial fitting window, W2 is the interpolation window.Some of the available data series were impossible to smooth with polynomials. They behaved like a typical logistic curve. Hence, the logistic curve was added to the set of functions that could be used to approximate the data:(3)L(t)=p01+p1ep3t1+p2ep3t.
The parameters W1 and W2 and the maximal degree of the polynomials (or logistic curve) were manually selected by visually inspection of the generated curves. The actual degree of the polynomials could vary for each polynomial fitted (thus the curve for a time series of a certain metabolite consists of polynomials of different degrees chained together): for each W1 window the polynomial giving the least amount of outliers was selected. Data points were identified as outliers when the difference between the measured value and the estimated value was larger than 5% of the range of the response variable. Table 3 summarized the fitting parameters used for each metabolite.Table 3
Parameters for the polynomial fitting of the different time series.W1 is expressed in fraction of the time between the first measurement and the last measurement. W2 is expressed as a fraction of W1. PD is the highest degree a polynomial may have; an L is displayed if the data were fitted to a logistic curve.
C→N limitationN→C limitationW1W2PDW1W2PDGLC0.400.8750.350.453NH30.500.80211LPiOH0.700.7550.450.804O20.180.9070.110.907CO20.120.9040.120.907Ac0.700.8850.150.703Lac1170.410.807Pyr0.400.80611LSuc0.400.8060.820.995Biomass0.500.7080.370.804The polynomial or logistic functions yield a smoothed time profile of the concentrations and enable to calculate the derivative of the concentrations to time. However, when calculating this derivative directly by symbolically deriving the function and filling in the time for the point of interest, one sometimes obtains unrealistically high values (D1 in Figure 2). Therefore, the derivative was approximated by taking the slope of the point of interest and the previous point on the polynomial (D2 in Figure 2).Figure 2
Two ways of estimating the derivative in pointa: D1 is calculated by symbolically deriving the function in point a; D2 is the slope between point a and point b (where a and b are calculated from the curve). Dots represent sample points through which a curve is fitted.The derivative can subsequently be used in (1) to calculate the cellular production (or consumption) fluxes:
(4)rp=dCdt-D(Cin-C).
Once rp has been calculated, classical MFA techniques (as described in [4]) can be applied at each time point. The results of those models, evaluated at the different time points, can then be combined to obtain a time profile of internal fluxes.


## 3.3. Growth Rate Calculation

During the transient period, the assumption that the growth rate is equal to the dilution rate does not hold. However, (4) can explicitly be rewritten for the following biomass:
(5)rp,Biom=dCBiomdt+DCBiom
from which the growth rate can be calculated as (6)μ=rp,BiomCBiom=dCBiomdt1CBiom+D.


## 4. Results

### 4.1. Polynomial Fitting

The polynomials fitted through the data (Figures3 and 4) are used to obtain the concentration, C, and the derivative of the concentration to the time, dC/dt, in each point. From those values, the conversion rate can be calculated (4). Biomass concentrations were calculated from OD values, after generating a calibration curve between cell dry weight (CDW) and OD values. This calibration curve was generated at the end of each experiment. As the curves for both experiments were very close to each other, no distinct calibration was used for the C-limited and N-limited conditions.Polynomial fit of some metabolites (expressed in g/L) for the experiment where carbon-limiting medium was replaced with nitrogen-limiting medium. The switch is made at time zero, when the cells are in carbon-limited steady state. Steady state values of five residence times before the switch and five residence times after the switch are not shown.
(a)(b)(c)(d)(e)(f)(g)(h)Polynomial fit of some metabolites (expressed in g/L) for the experiment where nitrogen-limiting medium was replaced with carbon-limiting medium. The switch is made at time zero, when the cells are in nitrogen-limited steady state. Steady state values of five residence times before the switch and five residence times after the switch are not shown.
(a)(b)(c)(d)(e)(f)(g)(h)The steady state concentration ofNH3 during N-limitation in the C- to N-limitation experiment is around 0.2 g/L (data are not shown on Figure 3, as this value falls out the time frame shown), while it is 0.1 g/L in the N to C-limitation case. These values are high compared to what is found in the literature for the residual nitrogen concentration in similar strains [29] and the nitrogen measurement kit used is probably not suitable for accurately measuring very low NH3 concentrations.


### 4.2. Growth Rate and Biomass Production

At time zero, the influent medium was simply switched. No cells were extracted and resuspended in the other medium, as done in [32]. Thus at the switch of the limiting compound, there is an abundance of both, as the medium limited in carbon is saturated with nitrogen while the medium limited in nitrogen is saturated with glucose. In the C to N-limitation experiment, it seems as if N-limitation is reached not much later than when carbon abundance starts (Figure 3). In the N to C-limitation case, nitrogen abundance is reached before all the glucose is depleted: after 5 hours, both glucose and nitrogen are not limiting. This can be seen in Figure 4, in the glucose profile. The glucose concentration becomes only zero after 12 hours while the nitrogen concentration already reaches steady state after 5 hours. At the medium switch, as there is excess nitrogen, the growth rate could increase, but apparently  E. coli needs some time to adapt to the new conditions. The biomass flux even decreases as nitrogen is added and becomes almost zero (top of Figure 5). Only seven hours after the medium switch, the growth rate increases sharply to attain almost 0.6 h-1 (right part of Figure 5). This increased growth rate coincides with the uptake of all glucose left in the broth (Figure 4).Fluxes in mol/L/h of the biomass production. Open symbols are the values as derived from the polynomials with formula 4; closed symbols are values obtained after flux balancing (top). The growth rate in 1/h of the cells during the transients (bottom). Left:C-limitation to N-limitation; right: N-limitation to C-limitation. The first point left of each figure is the steady state value before the medium switch. The last point right on each figure is the steady state value after at least 50 hours. Error bars represent the standard deviations.
(a)(b)(c)(d)Figure5 depicts the different growth rates during the transient conditions occurring when switching the limiting substrate. The same trends as in the biomass flux (top of Figure 5) can be observed. For the C- to N-limitation case, a sharp increase in growth rate is observed in the beginning, while for the N- to C-limitation the growth rate decreases to almost zero before increasing to 0.6 h-1, not far from the maximal growth rate of 0.7 h-1. In the carbon to nitrogen-limitation case, the maximal possible growth rate is never achieved and after a sharp increase in the beginning to 0.35 h-1 the growth rate decreases.The sharp increase in growth rate of the cells at the beginning of the medium switch in the C- to N-limitation experiment (left part of Figure5) is not accompanied by an increased flux through the Krebs cycle (Figure 8, the first hour after the steady state point, labeled with a triangle). An increase in biomass flux does not imply that intracellular fluxes relative to the biomass flux should increase as well. On the contrary, in the beginning the cells are more energy-efficient than during the carbon-limited steady state (left subfigure of Figure 7 shows that the ATP hydrolysis flux relative to the biomass flux is lower in the beginning implying that the maintenance energy requirement is lower). However, very quickly the flux through the glycolysis increases.During the switch of limiting substrate from carbon to nitrogen, less oxygen is consumed and lessCO2 is produced than during steady state (Figure 6). In the beginning of the medium switch, a sharp decline in O2 consumption and CO2 production is observed. This correlates with the increased growth rate and biomass production rate (Figure 5). For the experiment where nitrogen limitation is switched to carbon limitation, two parts can be observed: in the first period (until five hours), the O2 consumption rate and CO2 production rate increase. Then they suddenly decrease (around 7 hours), to steadily increase again until they reach the steady state level. In both experiments, the steady state value of oxygen consumption and carbon dioxide production is around 16 g/mol Biomass/h. This is in accordance with values obtained with a similar strain [4].Amount of oxygen consumed (top) and carbon dioxide produced (bottom). Left:C-limitation to N-limitation; right: N-limitation to C-limitation. The first point left of each figure is the steady state value before the medium switch. The last point right on each figure is the steady state value after at least 50 hours. Error bars represent the standard deviations.
(a)(b)(c)(d)Fluxes through the ATP hydrolysis reaction (ATPHY reaction): it can be considered as a measure for the maintenance requirement of the cells, as it combines all the futile cycles and ATP hydrolysed in nonspecific reactions. Left: C-limitation to N-limitation; right: N-limitation to C-limitation. The first point left of each figure is the steady state value before the medium switch. The last point right on each figure is the steady state value after at least 50 hours. Error bars represent the standard deviations.
(a)(b)Figure 8
Fluxmap of the glycolysis, penthose phosphate pathway, and citric acid cycle for the experiment where carbon-limiting medium is changed to nitrogen-limiting medium at time zero. The ordinate on each graph represents the flux expressed in mol/mol Biomass/h while the abscissa represents time. The first point left of each figure is the steady state value before the medium switch. The last point right on each figure is the steady state value after at least 50 hours. Error bars represent the standard deviations.The ATP hydrolysis reaction (Figure7) lumps all the maintenance requirements and futile cycles (as including all the futile cycles in the metabolic model would yield parallel pathways and the system of equations would not be solvable without measuring intracellular fluxes). It can be seen that during the lag phase of the experiment switching nitrogen to carbon-limitation (right part of Figure 7), the ATP maintenance requirement is around 0.6 mol ATP/mol Biom/h, which is high compared to the values found in a similar  E. coli [4]. It seems even more unusually high because there is minimal growth during this period. This high ATP requirement in the beginning of the lag phase correlates with a high uptake of glucose (PTS reaction in Figure 9) and an increased flux in the TCA cycle (Figure 9). The somewhat higher steady state ATP flux at the end of the C- to N-limitation (left part of Figure 7) may be due to multiple physiological steady states [33].Figure 9
Fluxmap of the glycolysis, penthose phosphate pathway, and citric acid cycle for the experiment where nitrogen-limiting medium is changed to carbon-limiting medium at time zero. The ordinate on each graph represents the flux expressed in mol/mol Biomass/h while the abscissa represents time. The first point left of each figure is the steady state value before the medium switch. The last point right on each figure is the steady state value after at least 50 hours. Error bars represent the standard deviations.


### 4.3. Central Carbon Metabolism

During the transient in the C- to N-limitation experiment, the fluxes through the glycolysis increase while those through the Krebs cycle decrease. Also an increase in acetate excretion is observed (Figure8). The opposite is observed in the N- to C-limitation experiment: fluxes through the glycolysis and acetate production decrease while the fluxes through the Krebs cycle increase (Figure 9). The decrease is not linear: at 15 hours the fluxes through the glycolysis reach a minimum to increase again until the steady state values are obtained after the subsequent 5 hours. In the Krebs cycle, the fluxes gradually increase the first five hours (coinciding with the lag phase observed in the biomass production, Figure 5) after which a sharp decline is observed, immediately followed by a gradual increase, peaking at around 15 hours, and then decrease to steady state level. The succinate flux sharply decreases in the beginning of the switch to gradually increase during the lag phase (5 hours) at which point it reaches a value that does not change anymore. The PEPCB flux, which is responsible for replenishing to the Krebs cycle, follows (or predates) this succinate excretion rate. The excretion of pyruvate decreases during the lag phase and then sharply increases again to reach a lower steady state than under N limitation.The pentose phosphate pathway is entirely driven by the demands in biomass components, and as each flux is expressed relatively to the biomass, to remove the influence of biomass on the flux values, the fluxes of the pentose phosphate pathway do not change, which results in seemingly huge error bars.


## 4.1. Polynomial Fitting

The polynomials fitted through the data (Figures3 and 4) are used to obtain the concentration, C, and the derivative of the concentration to the time, dC/dt, in each point. From those values, the conversion rate can be calculated (4). Biomass concentrations were calculated from OD values, after generating a calibration curve between cell dry weight (CDW) and OD values. This calibration curve was generated at the end of each experiment. As the curves for both experiments were very close to each other, no distinct calibration was used for the C-limited and N-limited conditions.Polynomial fit of some metabolites (expressed in g/L) for the experiment where carbon-limiting medium was replaced with nitrogen-limiting medium. The switch is made at time zero, when the cells are in carbon-limited steady state. Steady state values of five residence times before the switch and five residence times after the switch are not shown.
(a)(b)(c)(d)(e)(f)(g)(h)Polynomial fit of some metabolites (expressed in g/L) for the experiment where nitrogen-limiting medium was replaced with carbon-limiting medium. The switch is made at time zero, when the cells are in nitrogen-limited steady state. Steady state values of five residence times before the switch and five residence times after the switch are not shown.
(a)(b)(c)(d)(e)(f)(g)(h)The steady state concentration ofNH3 during N-limitation in the C- to N-limitation experiment is around 0.2 g/L (data are not shown on Figure 3, as this value falls out the time frame shown), while it is 0.1 g/L in the N to C-limitation case. These values are high compared to what is found in the literature for the residual nitrogen concentration in similar strains [29] and the nitrogen measurement kit used is probably not suitable for accurately measuring very low NH3 concentrations.


## 4.2. Growth Rate and Biomass Production

At time zero, the influent medium was simply switched. No cells were extracted and resuspended in the other medium, as done in [32]. Thus at the switch of the limiting compound, there is an abundance of both, as the medium limited in carbon is saturated with nitrogen while the medium limited in nitrogen is saturated with glucose. In the C to N-limitation experiment, it seems as if N-limitation is reached not much later than when carbon abundance starts (Figure 3). In the N to C-limitation case, nitrogen abundance is reached before all the glucose is depleted: after 5 hours, both glucose and nitrogen are not limiting. This can be seen in Figure 4, in the glucose profile. The glucose concentration becomes only zero after 12 hours while the nitrogen concentration already reaches steady state after 5 hours. At the medium switch, as there is excess nitrogen, the growth rate could increase, but apparently  E. coli needs some time to adapt to the new conditions. The biomass flux even decreases as nitrogen is added and becomes almost zero (top of Figure 5). Only seven hours after the medium switch, the growth rate increases sharply to attain almost 0.6 h-1 (right part of Figure 5). This increased growth rate coincides with the uptake of all glucose left in the broth (Figure 4).Fluxes in mol/L/h of the biomass production. Open symbols are the values as derived from the polynomials with formula 4; closed symbols are values obtained after flux balancing (top). The growth rate in 1/h of the cells during the transients (bottom). Left:C-limitation to N-limitation; right: N-limitation to C-limitation. The first point left of each figure is the steady state value before the medium switch. The last point right on each figure is the steady state value after at least 50 hours. Error bars represent the standard deviations.
(a)(b)(c)(d)Figure5 depicts the different growth rates during the transient conditions occurring when switching the limiting substrate. The same trends as in the biomass flux (top of Figure 5) can be observed. For the C- to N-limitation case, a sharp increase in growth rate is observed in the beginning, while for the N- to C-limitation the growth rate decreases to almost zero before increasing to 0.6 h-1, not far from the maximal growth rate of 0.7 h-1. In the carbon to nitrogen-limitation case, the maximal possible growth rate is never achieved and after a sharp increase in the beginning to 0.35 h-1 the growth rate decreases.The sharp increase in growth rate of the cells at the beginning of the medium switch in the C- to N-limitation experiment (left part of Figure5) is not accompanied by an increased flux through the Krebs cycle (Figure 8, the first hour after the steady state point, labeled with a triangle). An increase in biomass flux does not imply that intracellular fluxes relative to the biomass flux should increase as well. On the contrary, in the beginning the cells are more energy-efficient than during the carbon-limited steady state (left subfigure of Figure 7 shows that the ATP hydrolysis flux relative to the biomass flux is lower in the beginning implying that the maintenance energy requirement is lower). However, very quickly the flux through the glycolysis increases.During the switch of limiting substrate from carbon to nitrogen, less oxygen is consumed and lessCO2 is produced than during steady state (Figure 6). In the beginning of the medium switch, a sharp decline in O2 consumption and CO2 production is observed. This correlates with the increased growth rate and biomass production rate (Figure 5). For the experiment where nitrogen limitation is switched to carbon limitation, two parts can be observed: in the first period (until five hours), the O2 consumption rate and CO2 production rate increase. Then they suddenly decrease (around 7 hours), to steadily increase again until they reach the steady state level. In both experiments, the steady state value of oxygen consumption and carbon dioxide production is around 16 g/mol Biomass/h. This is in accordance with values obtained with a similar strain [4].Amount of oxygen consumed (top) and carbon dioxide produced (bottom). Left:C-limitation to N-limitation; right: N-limitation to C-limitation. The first point left of each figure is the steady state value before the medium switch. The last point right on each figure is the steady state value after at least 50 hours. Error bars represent the standard deviations.
(a)(b)(c)(d)Fluxes through the ATP hydrolysis reaction (ATPHY reaction): it can be considered as a measure for the maintenance requirement of the cells, as it combines all the futile cycles and ATP hydrolysed in nonspecific reactions. Left: C-limitation to N-limitation; right: N-limitation to C-limitation. The first point left of each figure is the steady state value before the medium switch. The last point right on each figure is the steady state value after at least 50 hours. Error bars represent the standard deviations.
(a)(b)Figure 8
Fluxmap of the glycolysis, penthose phosphate pathway, and citric acid cycle for the experiment where carbon-limiting medium is changed to nitrogen-limiting medium at time zero. The ordinate on each graph represents the flux expressed in mol/mol Biomass/h while the abscissa represents time. The first point left of each figure is the steady state value before the medium switch. The last point right on each figure is the steady state value after at least 50 hours. Error bars represent the standard deviations.The ATP hydrolysis reaction (Figure7) lumps all the maintenance requirements and futile cycles (as including all the futile cycles in the metabolic model would yield parallel pathways and the system of equations would not be solvable without measuring intracellular fluxes). It can be seen that during the lag phase of the experiment switching nitrogen to carbon-limitation (right part of Figure 7), the ATP maintenance requirement is around 0.6 mol ATP/mol Biom/h, which is high compared to the values found in a similar  E. coli [4]. It seems even more unusually high because there is minimal growth during this period. This high ATP requirement in the beginning of the lag phase correlates with a high uptake of glucose (PTS reaction in Figure 9) and an increased flux in the TCA cycle (Figure 9). The somewhat higher steady state ATP flux at the end of the C- to N-limitation (left part of Figure 7) may be due to multiple physiological steady states [33].Figure 9
Fluxmap of the glycolysis, penthose phosphate pathway, and citric acid cycle for the experiment where nitrogen-limiting medium is changed to carbon-limiting medium at time zero. The ordinate on each graph represents the flux expressed in mol/mol Biomass/h while the abscissa represents time. The first point left of each figure is the steady state value before the medium switch. The last point right on each figure is the steady state value after at least 50 hours. Error bars represent the standard deviations.


## 4.3. Central Carbon Metabolism

During the transient in the C- to N-limitation experiment, the fluxes through the glycolysis increase while those through the Krebs cycle decrease. Also an increase in acetate excretion is observed (Figure8). The opposite is observed in the N- to C-limitation experiment: fluxes through the glycolysis and acetate production decrease while the fluxes through the Krebs cycle increase (Figure 9). The decrease is not linear: at 15 hours the fluxes through the glycolysis reach a minimum to increase again until the steady state values are obtained after the subsequent 5 hours. In the Krebs cycle, the fluxes gradually increase the first five hours (coinciding with the lag phase observed in the biomass production, Figure 5) after which a sharp decline is observed, immediately followed by a gradual increase, peaking at around 15 hours, and then decrease to steady state level. The succinate flux sharply decreases in the beginning of the switch to gradually increase during the lag phase (5 hours) at which point it reaches a value that does not change anymore. The PEPCB flux, which is responsible for replenishing to the Krebs cycle, follows (or predates) this succinate excretion rate. The excretion of pyruvate decreases during the lag phase and then sharply increases again to reach a lower steady state than under N limitation.The pentose phosphate pathway is entirely driven by the demands in biomass components, and as each flux is expressed relatively to the biomass, to remove the influence of biomass on the flux values, the fluxes of the pentose phosphate pathway do not change, which results in seemingly huge error bars.


## 5. Discussion

### 5.1. Fluxes

When going from an nitrogen-limited medium to a carbon-limited one, a lag phase was observed (Figure5). A hypothesis could be that when more nitrogen becomes available, the cell starts to break down the accumulated glycogen which is redirected to energy metabolism for protein production. As the biomass composition is assumed to be constant in the model, this increased respiration is measured in the off-gas (O2 and CO2 fluxes, Figure 6) and propagated in the model to the ATPHY flux (Figure 7). However, the question remains why the cell would degrade glycogen while there is still plenty of carbon available in the environment. Furthermore, the carbon balance data (not shown) do not reflect an increase in carbon output compared to the carbon input.Another explanation is that ammonia abundance for cells previously grown under ammonia limiting conditions causes a futile cycle of ammonia import and export. Such cycles have been described previously in micro-organisms, although not under the experimental conditions described in this work [34]. Normally bacterial cells are highly resistant to ammonium, and the negative effects of high NH4+ concentrations are due to an enhanced osmolarity or increased ionic strength of the medium and not caused by ammonium itself [35]. But when changing the medium from nitrogen-limitation to carbon-limitation, the environment of the bacteria switches from low NH4+ concentrations to high ones, and it could be that  E. coli needs a certain adaptation period. The scavenging active during N-limitation could still be active in the beginning of the transition, resulting in a high influx of ammonia. To counter this, active efflux is needed, draining ATP and explaining the high energetic cost in the beginning of the addition of ammonium (right part of Figure 7).Such a mechanism has been observed in plants, whereNH4+ toxicity is due to the inability of root cells to limit the influx of ammonium. The high cytosolic NH4+ concentration activates the high-capacity and energy demanding ammonium efflux systems. This ammonium efflux can constitute as much as 80% of primary influx, resulting in a futile cycle of nitrogen across the plasma membrane of root cells. This futile cycle imposes a high energetic cost on the cell that is independent of N metabolism and is accompanied by a decline in growth. Plants that are resistant to high ammonium concentrations (e.g., Oriza sativa) limit the influx of ammonium by lowering the polarisation of the cellular membrane. This lowering of membrane polarisation is not observed in  Hordeum vulgare, sensitive to ammonium [36].The lack of growth resulting from adding ammonium to cells that are adapted to ammonium limitation has previously only been observed in cells lackingglnE (giving the gene product ATase, responsible for activating/inactivating glutamine synthetase) while the wild-type  E. coli did not show an impaired growth when shifted from ammonium limitation to excess [37, 38]. Furthermore, the lack of growth in  glnE mutants could be alleviated in cells constitutively expressing only low levels of glutamine synthetase [37], suggesting that the toxicity of ammonium in this case is not due to NH4+ itself, but maybe to accumulation of glutamine and/or glutamate. Normally, when shifting from ammonium limitation to carbon-limitation, the amount of active glutamine synthetase decreases instantly [39]. Possibly due to the lack of ATase, the enzymatic glutamine synthetase regulation is lost, and only genetic regulation is present. If this would be the case for the strain used in this study, it could explain the long lag phase after the medium switch. The cells survive until enough glutamine synthetase is degraded before they are able to grow again.Thus it seems asE. coli MG1655, the strain used in this study, has some anomalies in its nitrogen metabolism. It is known that multiple strains are denominated under this name [40]. The strain used in this work originated from the Coli Genetic Stock Center (CGSC) (personal communications from the Netherlands Culture Collection of Bacteria). The NCCB acquired this strain from CGSC prior to 1991. At that time, the strain was sent out as being  fnr-. Later it was found that the strain was a mixture of  fnr+ and  fnr- strains (personal communication from Mary Berlyn). For the experiments in this work, it turned out that the fnr- strain was used. More research is needed to find out why this strain does not grow well when NH4+ is added to a nitrogen-limiting environment.The lag phase found in the N-limitation to C-limitation experiment stands in sharp contrast with the experiment in which the carbon-limiting medium is replaced with the nitrogen-limiting one and thus carbon is added in the beginning of the switch while nitrogen is not yet limiting. The growth rate increases instantly (top of Figure5). Death and Ferenci [41] demonstrated that during carbon-limitation, sugar regulons are upregulated by internally synthesised sugars. Lequeux et al. [4] observed that under glucose limitation  ptsG is upregulated. The results of the experiments described here show that this up-regulation permits the cells to almost instantly increase their growth rate when carbon is added to a carbon-limited culture. The cells grow as fast as possible, until all nitrogen is consumed and thus the glucose concentration in the reactor broth increases around the same time as the nitrogen decreases (Figure 3).Hence the peak of succinate (and the smaller peak of pyruvate) occurring around 5 hours after the medium switch (Figure8) could be due to the depletion of ammonium (Figure 3). The high flux through the glycolysis is not stopped quickly enough and the carbon is excreted as pyruvate or diverted to PEPCB to be excreted as succinate.


### 5.2. Dynamic MFA

The extension of metabolic flux analysis to allow to predict also intracellular fluxes for cells growing under nonsteady state conditions was successfully applied to transient cultures. As raw measurement data are too noisy, smoothing had to be performed prior to the calculation of the exchange rates, needed for solving the metabolic model. An improvement for this smoothing process would be to incorporate more than two points in the calculation of the derivative, to avoid steep derivatives (and thus unrealistic fluxes) caused by the quirks of the higher-order polynomials.A possible improvement to the model is to allow fluctuations of the biomass composition. However, this should be accompanied with the appropriate measurements such as protein, RNA, DNA, lipid... content of the biomass. The biomass reaction can then be removed from the model and instead the different biomass compounds can be considered as excreted by the cell.


## 5.1. Fluxes

When going from an nitrogen-limited medium to a carbon-limited one, a lag phase was observed (Figure5). A hypothesis could be that when more nitrogen becomes available, the cell starts to break down the accumulated glycogen which is redirected to energy metabolism for protein production. As the biomass composition is assumed to be constant in the model, this increased respiration is measured in the off-gas (O2 and CO2 fluxes, Figure 6) and propagated in the model to the ATPHY flux (Figure 7). However, the question remains why the cell would degrade glycogen while there is still plenty of carbon available in the environment. Furthermore, the carbon balance data (not shown) do not reflect an increase in carbon output compared to the carbon input.Another explanation is that ammonia abundance for cells previously grown under ammonia limiting conditions causes a futile cycle of ammonia import and export. Such cycles have been described previously in micro-organisms, although not under the experimental conditions described in this work [34]. Normally bacterial cells are highly resistant to ammonium, and the negative effects of high NH4+ concentrations are due to an enhanced osmolarity or increased ionic strength of the medium and not caused by ammonium itself [35]. But when changing the medium from nitrogen-limitation to carbon-limitation, the environment of the bacteria switches from low NH4+ concentrations to high ones, and it could be that  E. coli needs a certain adaptation period. The scavenging active during N-limitation could still be active in the beginning of the transition, resulting in a high influx of ammonia. To counter this, active efflux is needed, draining ATP and explaining the high energetic cost in the beginning of the addition of ammonium (right part of Figure 7).Such a mechanism has been observed in plants, whereNH4+ toxicity is due to the inability of root cells to limit the influx of ammonium. The high cytosolic NH4+ concentration activates the high-capacity and energy demanding ammonium efflux systems. This ammonium efflux can constitute as much as 80% of primary influx, resulting in a futile cycle of nitrogen across the plasma membrane of root cells. This futile cycle imposes a high energetic cost on the cell that is independent of N metabolism and is accompanied by a decline in growth. Plants that are resistant to high ammonium concentrations (e.g., Oriza sativa) limit the influx of ammonium by lowering the polarisation of the cellular membrane. This lowering of membrane polarisation is not observed in  Hordeum vulgare, sensitive to ammonium [36].The lack of growth resulting from adding ammonium to cells that are adapted to ammonium limitation has previously only been observed in cells lackingglnE (giving the gene product ATase, responsible for activating/inactivating glutamine synthetase) while the wild-type  E. coli did not show an impaired growth when shifted from ammonium limitation to excess [37, 38]. Furthermore, the lack of growth in  glnE mutants could be alleviated in cells constitutively expressing only low levels of glutamine synthetase [37], suggesting that the toxicity of ammonium in this case is not due to NH4+ itself, but maybe to accumulation of glutamine and/or glutamate. Normally, when shifting from ammonium limitation to carbon-limitation, the amount of active glutamine synthetase decreases instantly [39]. Possibly due to the lack of ATase, the enzymatic glutamine synthetase regulation is lost, and only genetic regulation is present. If this would be the case for the strain used in this study, it could explain the long lag phase after the medium switch. The cells survive until enough glutamine synthetase is degraded before they are able to grow again.Thus it seems asE. coli MG1655, the strain used in this study, has some anomalies in its nitrogen metabolism. It is known that multiple strains are denominated under this name [40]. The strain used in this work originated from the Coli Genetic Stock Center (CGSC) (personal communications from the Netherlands Culture Collection of Bacteria). The NCCB acquired this strain from CGSC prior to 1991. At that time, the strain was sent out as being  fnr-. Later it was found that the strain was a mixture of  fnr+ and  fnr- strains (personal communication from Mary Berlyn). For the experiments in this work, it turned out that the fnr- strain was used. More research is needed to find out why this strain does not grow well when NH4+ is added to a nitrogen-limiting environment.The lag phase found in the N-limitation to C-limitation experiment stands in sharp contrast with the experiment in which the carbon-limiting medium is replaced with the nitrogen-limiting one and thus carbon is added in the beginning of the switch while nitrogen is not yet limiting. The growth rate increases instantly (top of Figure5). Death and Ferenci [41] demonstrated that during carbon-limitation, sugar regulons are upregulated by internally synthesised sugars. Lequeux et al. [4] observed that under glucose limitation  ptsG is upregulated. The results of the experiments described here show that this up-regulation permits the cells to almost instantly increase their growth rate when carbon is added to a carbon-limited culture. The cells grow as fast as possible, until all nitrogen is consumed and thus the glucose concentration in the reactor broth increases around the same time as the nitrogen decreases (Figure 3).Hence the peak of succinate (and the smaller peak of pyruvate) occurring around 5 hours after the medium switch (Figure8) could be due to the depletion of ammonium (Figure 3). The high flux through the glycolysis is not stopped quickly enough and the carbon is excreted as pyruvate or diverted to PEPCB to be excreted as succinate.


## 5.2. Dynamic MFA

The extension of metabolic flux analysis to allow to predict also intracellular fluxes for cells growing under nonsteady state conditions was successfully applied to transient cultures. As raw measurement data are too noisy, smoothing had to be performed prior to the calculation of the exchange rates, needed for solving the metabolic model. An improvement for this smoothing process would be to incorporate more than two points in the calculation of the derivative, to avoid steep derivatives (and thus unrealistic fluxes) caused by the quirks of the higher-order polynomials.A possible improvement to the model is to allow fluctuations of the biomass composition. However, this should be accompanied with the appropriate measurements such as protein, RNA, DNA, lipid... content of the biomass. The biomass reaction can then be removed from the model and instead the different biomass compounds can be considered as excreted by the cell.


## 6. Conclusion

This work introduced the concept of dynamic metabolic flux analysis. We demonstrated how to transform extracellular measurement data from dynamic experiments to flux values. It is not always clear in transient experiments whether the decrease in concentration of the extracellular metabolites is due to the cells stopping production and the remaining product diluting out or whether the cells are actively taking up the metabolites. Transforming the extracellular time series of concentrations to flux values can give the answer to that question. Furthermore, those flux values can help to get insight into the intracellular reactions.The transformation of time series of concentration measurements to flux values is based on differentiation (in the mathematical sense: finding the derivative) of those time series. Differentiation typically amplifies the noise on the data; therefore a noise-reducing step is needed prior to the differentiation. In this work, this was done by polynomial filtering. Selecting the right parameters for this filtering was done manually and produced satisfying results for the case presented. However, a statistically more sound approach may be desirable.We subsequently applied dynamic metabolic flux analysis onE. coli cultures in which the limiting compound in the medium was changed from nitrogen to glucose or from glucose to nitrogen.A lag phase occurred when cells adapted to a nitrogen-limiting environment were supplied with nitrogen. During this lag phase of several hours, the ATP demand was high. No clear physiological reason could be given. The only similar case described in literature forE. coli is a  glnE knock out, in which case the toxicity was not due to ammonium itself, but to the accumulation of glutamine/glutamate because the down-regulation of glutamine synthetase is not operative (glutamine synthetase is down-regulated by the gene product of glnE). A possible explanation for the long lag phase observed in the N-to-C-limitation experiment could be that sufficient glutamine synthetase has to be degraded so that the accumulation of glutamine/glutamate does not impair growth anymore.No lag phase occurred when supplying cells adapted to carbon limitation, with carbon. Instead the growth rate increased immediately. However, the highest growth rates were attained in the N to C-limitation experiment and not in the C- to N-limitation one.



---

*Source: 621645-2010-08-23.xml*