During batch fermentation, a variety of compounds are synthesized, as microorganisms undergo distinct growth phases: lag, exponential, growth-no-growth transition, stationary, and decay. A detailed understanding of the metabolic pathways involved in these phases is crucial for optimizing the production of target compounds. Dynamic flux balance analysis (dFBA) offers insight into the dynamics of metabolic pathways. However, explaining secondary metabolism remains a challenge. A multiphase and multi-objective dFBA scheme (MPMO model) has been proposed for this purpose. However, its formulation is discontinuous, changing from phase to phase; its accuracy in predicting intracellular fluxes is hampered by the lack of a mechanistic link between phases; and its simulation requires considerable computational effort. To address these limitations, we combine a novel model with a genome-scale model to predict the distribution of intracellular fluxes throughout batch fermentation. This integrated multiphase continuous model (IMC) has a unique formulation over time, and it incorporates empirical regulatory descriptions to automatically identify phase transitions and incorporates the hypotheses that yeasts might vary their cellular objective over time to adapt to the changing environment. We validated the predictive capacity of the IMC model by comparing its predictions with intracellular metabolomics data for Saccharomyces uvarum during batch fermentation. The model aligns well with the data, confirming its predictive capabilities. Notably, the IMC model accurately predicts trehalose accumulation, which was enforced in the MPMO model. We further demonstrate the generalizability of the IMC model, explaining the dynamics of primary and secondary metabolism of three Saccharomyces species. The model provides biological insights consistent with the literature and metabolomics data, establishing it as a valuable tool for exploring the dynamics of novel fermentation processes. IMPORTANCE: This work presents an integrated multiphase continuous dynamic genome-scale model (IMC model) for batch fermentation, a crucial process widely used in industry to produce biofuels, enzymes, pharmaceuticals, and food products or ingredients. The IMC model integrates a continuous kinetic model with a genome-scale model to address the critical limitations of existing dynamic flux balance analysis schemes, such as the difficulty of explaining secondary metabolism, the lack of mechanistic links between growth phases, or the high computational demands. The model also introduces the hypothesis that cells adapt the FBA objective over time. The IMC improves the accuracy of intracellular flux predictions and simplifies the implementation process with a unique dFBA formulation over time. Its ability to predict both primary and secondary metabolism dynamics in different Saccharomyces species underscores its versatility and robustness. Furthermore, its alignment with empirical metabolomics data validates its predictive power, offering valuable insights into metabolic processes during batch fermentation. These advances pave the way for optimizing fermentation processes, potentially leading to more efficient production of target compounds and novel biotechnological applications.

Saccharomyces non‐cerevisiae yeasts are gaining momentum in wine fermentation due to their potential to reduce ethanol content and achieve attractive aroma profiles. However, the design of the fermentation process for new species requires intensive experimentation. The use of mechanistic models could automate process design, yet to date, most fermentation models have focused on primary metabolism. Therefore, these models do not provide insight into the production of secondary metabolites essential for wine quality, such as aromas. In this work, we formulate a continuous model that accounts for the physiological status of yeast, that is, exponential growth, growth under nitrogen starvation and transition to stationary or decay phases. To do so, we assumed that nitrogen starvation is associated with carbohydrate accumulation and the induction of a set of transcriptional changes associated with the stationary phase. The model accurately described the dynamics of time series data for biomass and primary and secondary metabolites obtained for various yeast species in single culture fermentations. We also used the proposed model to explore different process designs, showing how the addition of nitrogen could affect the aromatic profile of wine. This study underlines the potential of incorporating yeast physiology into batch fermentation modelling and provides a new means of automating process design.

Yeast batch fermentation is widely used in industrial biotechnology, yet its performance is strongly influenced by temperature and nitrogen availability, which affect growth kinetics and metabolite production. The development of predictive models that accurately describe these effects is essential for automating and optimizing fermentation design, reducing trial-and-error experimentation, and improving process efficiency and product quality. However, most mathematical models focus on primary metabolism and lack a systematic approach to integrate the effects of temperature. Existing models often rely on empirical corrections with limited predictive power beyond specific experimental conditions. Furthermore, there is no unified framework for optimizing fermentation processes while accounting for the temperature-dependent metabolic responses. We addressed these gaps by developing a temperature-dependent kinetic model for nitrogen-limited batch fermentation by Saccharomyces cerevisiae. The modeling approach is based on advanced systems identification, integrating identifiability analyses (structural and practical), multi-experiment parameter estimation, and automated model selection to determine the most appropriate temperature dependencies for key metabolic processes. Validated across five industrial S. cerevisiae strains in an illustrative example related to wine fermentation, the model exhibited strong predictive performance (NRMSE $<10.5\%$, median R²$>0.95$) and enabled simulation-based process optimization, including nitrogen-supplementation strategies and strain selection for improved fermentation outcomes. By providing a systematic modeling framework that accounts for temperature effects, this work bridges a critical gap in predictive modeling and advances the rational design and control of industrial fermentation processes.

Batch fermentation is a biotechnological dynamic process that produces various products by employing microorganisms that undergo different growth phases: lag, exponential, growth-non-growth, stationary, and decay. Genome-scale constrained-based models are commonly used to explore the phenotypic potential of these microorganisms. Previous studies have primarily used dynamic Flux Balance Analysis (dFBA) to elucidate the metabolism during the exponential phase. However, this approach falls short in addressing the multi-phase nature of the process and secondary metabolism, posing significant challenges to our understanding of batch fermentation. A recent attempt at a solution was a discontinuous, multi-phase, multi-objective dFBA implementation. However, this approximation lacks the mechanistic connection between phases, limiting its applicability in predicting intracellular fluxes during batch fermentation. To overcome these limitations, we combined a novel continuous model with a genome-scale model to predict the distribution of intracellular fluxes throughout the batch fermentation process. The proposed model includes empirical descriptions of regulation that automatically identify the transition between phases. Its application to explain primary and secondary metabolism of Saccharomyces species in batch fermentation results in biological insights that are in good agreement with the previous literature. The ability to account for all process phases and explain secondary metabolism makes this model a valuable and easy-to-use tool for exploring novel fermentation processes. IMPORTANCE: This research proposes a novel dynamic genome-scale modelling approach for batch fermentation, a crucial process widely used to produce a diverse range of products such as biofuels, enzymes, pharmaceuticals, and food products or ingredients. The proposed approach automatically accounts for the transitions between different phases of the fermentation process (lag, exponential, growth-no-growth, and stationary). This is a significant advancement over previous methods that required different model formulations for different phases. We have successfully applied this modelling approach to explore the primary and secondary metabolism of three yeast species under batch fermentation conditions. The model accurately explained experimental data and provided biological insights consistent with previous research findings, instilling confidence in its reliability and accuracy. The ability of this modelling approach to explain primary and secondary metabolism makes it a valuable tool for designing novel, more efficient, and effective fermentation processes, which could have far-reaching implications in industrial biotechnology.