Saccharomyces cerevisiae is a model organism in eukaryote cell research and the workhorse for the biotechnology industry. In nature and the industrial setup, environmental perturbations act as stressing factors which challenge regulation of metabolic flux and can also lead to reduced performance in industrial applications.
1. Introduction
Kinetic metabolic models are mathematical representations of a biological system that consider kinetic expressions such as rate constants. They describe the network structure, kinetic rate expressions and contain values for the parameters in these expressions[1] [6]. Thus, these descriptions are well-suited to model time-dependent dynamics. A detailed explanation of the main components in a kinetic metabolic model can be seen in Box 1. Despite the progress attained with them, a consensus version with a full coverage of CCM has not yet been achieved.
2. Glycolytic Response to Glucose Perturbations in Yeast Fermentations
Saccharomyces cerevisiae is one of the most used microorganisms in biotechnology. S. cerevisiae is a prominent cell factory involved in food, beverages, and biofuels industries
[2][3] [38,39]. On top of its favorable physiology and robustness, genetic engineering has allowed to introduce new pathways and improve existing ones, generating new strains that have widened its range of applications
[4][5][1,40]. Nonetheless, scaling up to commercial production is a challenging stage in which developed strains may emerge as inefficient
[6][41]. Long circulation times and nonideal mixing result in substrate gradients in the industrial fermenter, affecting most cell factories, including S. cerevisiae
[7][8][9][10] [3,4,42,43]. The yeast cell sees these gradients as stressing factors to which it continuously adapts, often deteriorating process yields and giving relevance to the development of stress tolerant strains
[11][12][13][14][2].
Extracellular substrate gradients alter intracellular fluxes in CCM. Carbon flux shifts between the different pathways composing CCM during these temporal transitions
[12][5]. This can become a challenge for the cell, which struggles to keep the different pathways composing CCM balanced
[13] [44], as was shown in
[14][18] for a yeast strain with a defective trehalose cycle, where sudden exposure to a high glucose concentration resulted in growth arrest. Glycolysis is found at the core of this network. This pathway digests intracellular glucose into pyruvate and produces energy in the form of ATP and glycolytic intermediates that support anabolic reactions
[15][45].
How glycolysis contributes to the metabolic processes inside the cell depends on multiple factors. The presence or absence of oxygen determines if pyruvate is used for respiration or fermentation
[16][17] [46,47]. Still, this conspicuously simple explanation is challenged at high-substrate concentrations, where the maximum respiratory capacity is reached and fermentation takes place even if oxygen is present
[18][19] [48,49], in what is known as ‘overflow metabolism’ or Crabtree effect
[14] [50]. In addition, the substrate that is used as carbon source (such as glucose or fructose) and the ability of a strain to metabolize it also affects glycolytic kinetics and process yields
[20][21][22][23][24][25][26] [
51,
52,
53,
54,
55,
56,
57]. Furthermore, the cellular state determines how glycolytic intermediates are used as biomass precursors
[15][27] [45,58]. For instance, at changing growth rates, different usage of these precursors can be observed
[28][29] [34,59]. Finally, availability of cofactors cannot always be taken for granted. A higher substrate uptake rate might be an evolutionary advantage, but it results in a demand for NADH recycling that respiration cannot achieve and thus fermentation becomes active
[30][31] [60,61].
The response of glycolysis to dynamic glucose perturbations is controlled by different regulatory layers. The first mechanism is the storage of glycogen and trehalose when glucose uptake exceeds the glycolytic processing capacity
[32][62]. On top of this, allosteric and post-translation regulation take place
[33][63]. Hexokinase (HXK) is allosterically inhibited by trehalose-6-phosphate (T6P), pyruvate kinase (PYK) is activated by fructose-1,6-bis-phosphate (FBP) and multiple metabolites act on phosphofructokinase (PFK)
[34][35][36] [27,64,65]. Simultaneously, the cAMP-protein kinase A (PKA) pathway is activated upon glucose perturbation and starts a regulation cascade in CCM
[37][66] and possible targets for Post-Translational Modifications (PTMs) have been found in multiple enzymes along the CCM
[38] [67]. Finally, to adapt to different growth conditions, yeast cells use different enzyme isoforms. For instance, hexokinases and glucokinases are balanced to adapt to different glucose concentrations
[39][68] and the regulation of intracellular pH is compartment-specific, carried out by different ATPases
[40][69].
3. The Development of Metabolic Models Has Resulted in Understanding of Key Glycolytic Properties
Many breakthroughs in metabolic modeling used genome scale models. Nonetheless, stoichiometry alone does not define function and the response to glucose perturbations is a dynamic process where stoichiometry cannot explain mechanisms that act at different time scales or the appearance of bistability, among others
[1][41][6,70]. As a result, kinetic models enable a deeper understanding of glycolytic properties. Due to the abundant data available for S. cerevisiae fermentations, models of the glycolytic networks have reached a high level of maturity for this organism.
The first kinetic models developed focused on understanding glycolytic oscillations in nongrowing yeast cells
[42][43][44][45][46][47][48][49][50][23,71,72,73,74,75,76,77,78]. Most enzymatic reactions were lumped into a few (except
[42][49][23,77]) but they acknowledged the important role of enzyme PFK and showed sensitivity to different glucose, oxygen, and acetaldehyde concentrations. Later works focused on understanding control properties and glycolytic response upon a single glucose perturbation experiment
[51][52][53][54] [19,79,80,81,82] and thanks to a progressive increase in experimental data available, more detailed models were developed
[55][56][57][58][16,17,18,22]. Much of the focus was on understanding how mutant strains lacking a functional trehalose cycle would undergo growth arrest upon the glucose perturbation
[32][35][62,64]. This was found to be due to a glycolytic imbalance between upper and lower glycolysis and attributed first to an ATP turbo metabolism
[54][
82]. Later, ref.
[57][18] explained the role that the trehalose cycle plays in the glycolytic response and highlighted how the intracellular concentrations of metabolites at a given time point modulate the outcome.
In this process, models have become more interconnected with other pathways, allowing for a more complete understanding of the glycolytic response. Ref.
[58][22] introduced glycolytic byproduct branch reactions that were necessary to reproduce the steady state. Other works modeled pathways that are directly linked to yeast glycolysis. For instance, detailed descriptions of the glycerol synthesis, trehalose cycle and PPP were developed in
[59][60][61][21,24,83], respectively. Later, a PPP model was connected to glycolysis in
[62] [33], and another model of glycolysis together with TCA was developed in
[63][20]. These networks were used to understand the control properties of glycolysis, pointing to glucose transporter (GLT) and PFK for being the enzymes with the highest controlling coefficients
[51][52][64][65][66][79,80,84,85,86] and to study the effect of genome duplications
[67][87]. For a complete overview of metabolic models developed to understand dynamic perturbations, see
Table 1.
Furthermore, the regulation exerted by cofactors has gradually become more evident, resulting in a more complex understanding of glycolysis. The depletion of inorganic phosphate concentration that was shown to be crucial in
[57][
18] had been overlooked in previous works where it was assumed to be constant over time. Simultaneously, the sum of adenosine nucleotides has been assumed to be a conserved moiety
[55][16] but under some experimental conditions this is not the case
[68][69][31,88], which can be relevant considering that controlling enzyme PFK is allosterically regulated by ATP and AMP.
Table 1. Properties of S. cerevisiae models developed to understand dynamic glucose perturbation response: glycolysis (GLYCO), tricarboxylic acid cycle (TCA), pentose phosphate pathway (PPP), trehalose cycle (TRE). Number of ‘+’ sign according to how advantageous the property is. Cofactor conservation moieties are sumAXP and sumNADX. N/A when reactions were not modeled, or data were not shown in article. Refs.
[56][63][17,20] fitted different parameter sets to multiple data sets. Other models used a unique parameter set. From the literature pool of articles obtained in the systematic reviewing process, only the works which include glycolysis are displayed.
| |
Rizzi et al.[70] | Rizzi et al. [19] |
Teusink et al. [54] | Teusink et al. [82] |
Teusink et al. [58] | Teusink et al. [22] |
van Eunen et al. [56] | van Eunen et al. [17] |
| Contribution to glycolytic understanding |
Dynamic models can accurately describe glucose perturbation. |
ATP surplus can cause the observed overactivation of initial glycolytic steps in DTps1 mutant strains. |
In vivo behavior cannot be predicted with in vitro kinetics. |
Implementation of allosteric regulation and in vivo measured parameter values is necessary to reproduce GP data. |
| GLYCO |
Individual + branch reactions (++) |
Lumped reactions (+) |
Individual + branch reactions (++) |
Individual + branch reactions (++) |
| TRE |
N/A |
N/A |
N/A |
T6P regulation (+) |
| TCA |
Individual reactions (++) |
N/A |
N/A |
N/A |
| PPP |
N/A |
N/A |
N/A |
N/A |
| Cofactors |
Conservation moiety (+) |
Conservation moiety (+) |
Conservation moiety (+) |
Conservation moiety (+) |
] but neither has been implemented.
On top of this, other variables affect individual enzyme kinetics, and have neither been considered. First, cytosolic pH decays upon extracellular glucose perturbation, affecting multiple intracellular processes, including enzyme kinetics
[35][69,94]. Second, PTMs are a fast response mechanism and multiple target sites have been found throughout CCM
[78][95]. Third, different enzyme isoforms are expressed under different growth regimes. Examples of this are the differential expression of GLK/HXK and Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) genes (
[79][96] and
[80][97], respectively).
Finally, a key challenge is the representation of variables that are not part of the carbon flux, such as cofactors. Most models have kept them constant or adopted moiety conservation cycles
[81][98], such as the sum of intracellular adenine nucleotides ([ATP] + [ADP] + [AMP] = [AXP]) or inorganic phosphate
[55][17]. Nonetheless, under intense glucose perturbations, both variables behave in a dynamic manner
[56][82][68][83] [19,26,31,99] and alter glycolytic response. An example of this is the ATP paradox, which occurs when ATP and the sum of adenine nucleotides transiently decay
[84][100]. Understanding cytosolic Pi as a dynamic variable and implementation of import from the vacuole turned out to be central in understanding the glycolytic imbalance
[57][18]. Although the availability of Pi was essential for lower glycolysis progression via GAPDH
[57][18], adenine nucleotides exert allosteric regulation on the important controlling enzyme PFK
[85][101].
5. New Intracellular Metabolomic and Fluxomic Data Boost Understanding of Glycolytic Response
Scale-down approaches have been developed to understand long-standing problems in industrial bioreactors. Although this has granted valuable knowledge, essential intracellular properties such as in vivo fluxes and kinetics have been captured with only limited resolution, constraining model development. In fact, this has become one of the main challenges in the development of high quality predictive kinetic models, since often multiple variables, such as transcriptomics, metabolomics and fluxomic data, interact to result in the final response
[1] [6].
Early works aimed to understand glycolytic oscillations did so with small datasets, reducing their range of implementation. On most occasions only extracellular data such as growth and nutrient exchange rates was available
[48][76] or a few metabolites at most
[47][75], until in vivo quantification of metabolite concentrations and fluxes became a common practice, where most cofactors, glycolytic intermediates and rates were simultaneously observable
[58] [22]. Later, a standardized dynamic glucose perturbation experimental setup with CEN-PK yeast strains was adopted (see
Table 2). This consisted of chemostat growth at dilution rate of 0.1 h
−1, followed by an external glucose perturbation, where extracellular concentration increased to 1 g L
−1. These stimulus response experiments were used to infer more physiological patterns
[82][26] and the use of Nuclear Magnetic Resonance (NMR) and Mass Spectroscopy (MS) techniques made a wide range of intracellular metabolites measurable. From only a few glycolytic concentrations, datasets gradually grew to include most metabolites in glycolysis, the trehalose cycle, the TCA cycle, and the PPP. Adenine nucleotides and NAD:NADH ratio have also been made a standard and other nucleotides and amino acids which are affected by carbon uptake dynamics are quantified in the most recent publications.
