In industrial large-scale bioreactors, microorganisms encounter heterogeneous substrate concentration conditions, which can impact growth or product formation. Here we carried out an extended (12 h) experiment of repeated glucose pulsing with a 10-min period to simulate fluctuating glucose concentrations with Aspergillus niger producing glucoamylase, and investigated its dynamic response by rapid sampling and quantitative metabolomics. The 10-min period represents worst-case conditions, as in industrial bioreactors the average cycling duration is usually in the order of 1 min. We found that cell growth and the glucoamylase productivity were not significantly affected, despite striking metabolomic dynamics. Periodical dynamic responses were found across all central carbon metabolism pathways, with different time scales, and the frequently reported ATP paradox was confirmed for this A. niger strain under the dynamic conditions. A thermodynamics analysis revealed that several reactions of the central carbon metabolism remained in equilibrium even under periodical dynamic conditions. The dynamic response profiles of the intracellular metabolites did not change during the pulse exposure, showing no significant adaptation of the strain to the more than 60 perturbation cycles applied. The apparent high tolerance of the glucoamylase producing A. niger strain for extreme variations in the glucose availability presents valuable information for the design of robust industrial microbial hosts.

In a 54 m³ large-scale penicillin fermentor, the cells experience substrate gradient cycles at the timescales of global mixing time about 20–40 s. Here, we used an intermittent feeding regime (IFR) and a two-compartment reactor (TCR) to mimic these substrate gradients at laboratory-scale continuous cultures. The IFR was applied to simulate substrate dynamics experienced by the cells at full scale at timescales of tens of seconds to minutes (30 s, 3 min and 6 min), while the TCR was designed to simulate substrate gradients at an applied mean residence time ((Formula presented.)) of 6 min. A biological systems analysis of the response of an industrial high-yielding P. chrysogenum strain has been performed in these continuous cultures. Compared to an undisturbed continuous feeding regime in a single reactor, the penicillin productivity (q_PenG) was reduced in all scale-down simulators. The dynamic metabolomics data indicated that in the IFRs, the cells accumulated high levels of the central metabolites during the feast phase to actively cope with external substrate deprivation during the famine phase. In contrast, in the TCR system, the storage pool (e.g. mannitol and arabitol) constituted a large contribution of carbon supply in the non-feed compartment. Further, transcript analysis revealed that all scale-down simulators gave different expression levels of the glucose/hexose transporter genes and the penicillin gene clusters. The results showed that q_PenG did not correlate well with exposure to the substrate regimes (excess, limitation and starvation), but there was a clear inverse relation between q_PenG and the intracellular glucose level.

In the present work, by performing chemostat experiments at 400 and 600 RPM, two typical power inputs representative of industrial penicillin fermentation (P/V, 1.00 kW/m³ in more remote zones and 3.83 kW/m³ in the vicinity of the impellers, respectively) were scaled-down to bench-scale bioreactors. It was found that at 400 RPM applied in prolonged glucose-limited chemostat cultures, the previously reported degeneration of penicillin production using an industrial Penicillium chrysogenum strain was virtually absent. To investigate this, the cellular response was studied at flux (stoichiometry), residual glucose, intracellular metabolite and transcript levels. At 600 RPM, 20% more cell lysis was observed and the increased degeneration of penicillin production was accompanied by a 22% larger ATP gap and an unexpected 20-fold decrease in the residual glucose concentration (C_s,out). At the same time, the biomass specific glucose consumption rate (q_s) did not change but the intracellular glucose concentration was about sixfold higher, which indicates a change to a higher affinity glucose transporter at 600 RPM. In addition, power input differences cause differences in the diffusion rates of glucose and the calculated Batchelor diffusion length scale suggests the presence of a glucose diffusion layer at the glucose transporting parts of the hyphae, which was further substantiated by a simple proposed glucose diffusion-uptake model. By analysis of calculated mass action ratios (MARs) and energy consumption, it indicated that at 600 RPM glucose sensing and signal transduction in response to the low C_s,out appear to trigger a gluconeogenic type of metabolic flux rearrangement, a futile cycle through the pentose phosphate pathway (PPP) and a declining redox state of the cytosol. In support of the change in glucose transport and degeneration of penicillin production at 600 RPM, the transcript levels of the putative high-affinity glucose/hexose transporter genes Pc12g02880 and Pc06g01340 increased 3.5- and 3.3-fold, respectively, and those of the pcbC gene encoding isopenicillin N-synthetase (IPNS) were more than twofold lower in the time range of 100–200 hr of the chemostat cultures. Summarizing, changes at power input have unexpected effects on degeneration and glucose transport, and result in significant metabolic rearrangements. These findings are relevant for the industrial production of penicillin, and other fermentations with filamentous microorganisms.

Oxygen limitation is regarded as a useful strategy to improve enzyme production by mycelial fungus like Aspergillus niger. However, the intracellular metabolic response of A. niger to oxygen limitation is still obscure. To address this, the metabolism of A. niger was studied using multi-omics integrated analysis based on the latest GEMs (genome-scale metabolic model), including metabolomics, fluxomics and transcriptomics. Upon sharp reduction of the oxygen supply, A. niger metabolism shifted to higher redox level status, as well as lower energy supply, down-regulation of genes for fatty acid synthesis and a rapid decrease of the specific growth rate. The gene expression of the glyoxylate bypass was activated, which was consistent with flux analysis using the A. niger GEMs iHL1210. The increasing flux of the glyoxylate bypass was assumed to reduce the NADH formation from TCA cycle and benefit maintenance of the cellular redox balance under hypoxic conditions. In addition, the relative fluxes of the EMP pathway were increased, which possibly relieved the energy demand for cell metabolism. The above multi-omics integrative analysis provided new insights on metabolic regulatory mechanisms of A. niger associated with enzyme production under oxygen-limited condition, which will benefit systematic design and optimization of the A. niger microbial cell factory.

A powerful approach for the optimization of industrial bioprocesses is to perform detailed simulations integrating large-scale computational fluid dynamics (CFD) and cellular reaction dynamics (CRD). However, complex metabolic kinetic models containing a large number of equations pose formidable challenges in CFD-CRD coupling and computation time afterward. This necessitates to formulate a relatively simple but yet representative model structure. Such a kinetic model should be able to reproduce metabolic responses for short-term (mixing time scale of tens of seconds) and long-term (fed-batch cultivation of hours/days) dynamics in industrial bioprocesses. In this paper, we used Penicillium chrysogenum as a model system and developed a metabolically structured kinetic model for growth and production. By lumping the most important intracellular metabolites in 5 pools and 4 intracellular enzyme pools, linked by 10 reactions, we succeeded in maintaining the model structure relatively simple, while providing informative insight into the state of the organism. The performance of this 9-pool model was validated with a periodic glucose feast–famine cycle experiment at the minute time scale. Comparison of this model and a reported black box model for this strain shows the necessity of employing a structured model under feast–famine conditions. This proposed model provides deeper insight into the in vivo kinetics and, most importantly, can be straightforwardly integrated into a computational fluid dynamic framework for simulating complete fermentation performance and cell population dynamics in large scale and small scale fermentors. Biotechnol. Bioeng. 2017;114: 1733–1743.

The trajectories, referred to as lifelines, of individual microorganisms in an industrial scale fermentor under substrate limiting conditions were studied using an Euler-Lagrange computational fluid dynamics approach. The metabolic response to substrate concentration variations along these lifelines provides deep insight in the dynamic environment inside a large-scale fermentor, from the point of view of the microorganisms themselves. We present a novel methodology to evaluate this metabolic response, based on transitions between metabolic “regimes” that can provide a comprehensive statistical insight in the environmental fluctuations experienced by microorganisms inside an industrial bioreactor. These statistics provide the groundwork for the design of representative scale-down simulators, mimicking substrate variations experimentally. To focus on the methodology we use an industrial fermentation of Penicillium chrysogenum in a simplified representation, dealing with only glucose gradients, single-phase hydrodynamics, and assuming no limitation in oxygen supply, but reasonably capturing the relevant timescales. Nevertheless, the methodology provides useful insight in the relation between flow and component fluctuation timescales that are expected to hold in physically more thorough simulations. Microorganisms experience substrate fluctuations at timescales of seconds, in the order of magnitude of the global circulation time. Such rapid fluctuations should be replicated in truly industrially representative scale-down simulators.

Background: Aspergillus niger is widely used for enzyme production and achievement of high enzyme production depends on the comprehensive understanding of cell's metabolic regulation mechanisms. Results: In this paper, we investigate the metabolic differences and regulation mechanisms between a high glucoamylase-producing strain A. niger DS03043 and its wild-type parent strain A. niger CBS513.88 via an integrated isotope-assisted metabolomics and ¹³C metabolic flux analysis approach. We found that A. niger DS03043 had higher cell growth, glucose uptake, and glucoamylase production rates but lower oxalic acid and citric acid secretion rates. In response to above phenotype changes, A. niger DS03043 was characterized by an increased carbon flux directed to the oxidative pentose phosphate pathway in contrast to reduced flux through TCA cycle, which were confirmed by consistent changes in pool sizes of metabolites. A higher ratio of ATP over AMP in the high producing strain might contribute to the increase in the PP pathway flux as glucosephosphate isomerase was inhibited at higher ATP concentrations. A. niger CBS513.88, however, was in a higher redox state due to the imbalance of NADH regeneration and consumption, resulting in the secretion of oxalic acid and citric acid, as well as the accumulation of intracellular OAA and PEP, which may in turn result in the decrease in the glucose uptake rate. Conclusions: The application of integrated metabolomics and ¹³C metabolic flux analysis highlights the regulation mechanisms of energy and redox metabolism on flux redistribution in A. niger. Graphical abstract An integrated isotope-assisted metabolomics and ¹³C metabolic flux analysis was was firstly systematically performed in A. niger. In response to enzyme production, the metabolic flux in A. niger DS03043 (high-producing) was redistributed, characterized by an increased carbon flux directed to the oxidative pentose phosphate pathway as well as an increased pool size of pentose. The consistency in ¹³C metabolic flux analysis and metabolites quantification indicated that an imbalance of NADH formation and consumption led to the accumulation and secretion of organic acids in A. niger CBS513.88 (wild-type)