Bacterial Growth in Pyruvate-Based Media
The optical density (OD) measurements at 600 nm were taken at various time points to assess the bacterial growth in three separate flasks (Fig. 1). At the start of the experiment (0 hours), the initial OD values were low, with Flask 1 at 0.108, Flask 2 at 0.090, and Flask 3 at 0.091, indicating a low bacterial density. After 3 hours, the OD values increased slightly, with Flask 1 measuring 0.23, Flask 2 at 0.223, and Flask 3 at 0.229, demonstrating early bacterial growth. At 9 hours, just prior to the addition of IPTG for protein induction, the OD values had further increased to approximately 0.55 in all flasks. Following IPTG induction at 9 hours, the OD values at 10 hours (1 hour post-induction) showed a small but consistent rise, with values around 0.55 in all flasks. By 24 hours (14 hours post-induction), the cultures reached their peak OD readings, with Flask 1 at 0.742, Flask 2 at 0.782, and Flask 3 at 0.772. At 26 hours (16 hours post-induction), the OD values remained stable around 0.78 for all flasks, indicating continued growth. However, at the 36-hour time point (26 hours post-induction), a slight decrease in OD was observed across all flasks, with values dropping to 0.679 in Flask 1, 0.667 in Flask 2, and 0.700 in Flask 3. This decline suggests that the bacterial cultures may have entered the stationary phase or experienced nutrient depletion in the medium.
Metabolite Profiling in the Culture Medium
The bar plots in Fig. 1 illustrate the temporal concentration changes of key metabolites involved in central metabolic pathways during bacterial growth and recombinant protein production. These metabolites were measured at various time points: 0h, 3h, 9h (IPTG induction), 10h, 24h, 26h, and 36h, capturing metabolic dynamics from the lag phase through the exponential and stationary phases of bacterial growth. The analysis of metabolite concentrations across different time points revealed significant changes in key intermediates involved in central carbon metabolism during bacterial growth and recombinant protein production. These metabolic shifts were observed both before and after IPTG induction, reflecting the cells' metabolic response to the high availability of pyruvate as the primary carbon source.
Pyruvate (Fig. 1A), the main carbon source in this experiment, exhibited a rapid decrease starting well before IPTG induction, with its levels already substantially reduced by the 9-hour mark. This early decline in pyruvate suggests that bacterial cells are actively consuming pyruvate through multiple metabolic pathways to support energy production and biosynthesis. The sustained reduction of pyruvate post-9h further indicates that the cells continue to use pyruvate intensively during recombinant protein production. The early depletion of pyruvate is likely a result of its diversion into multiple pathways, including both the TCA cycle for energy production and overflow metabolic pathways. Concurrent with the depletion of pyruvate, both acetate (Fig. 1B) and lactate (Fig. 1C) showed substantial increases in the culture media prior to IPTG induction. This observation suggests that bacterial cells, faced with an abundance of pyruvate, are converting excess pyruvate into acetate and lactate as a means of mitigating the buildup of intracellular pyruvate. These metabolites are likely excreted into the medium as waste products to avoid toxic intracellular accumulation. The high concentrations of acetate and lactate in the culture media reflect overflow metabolism, a process whereby excess carbon from glycolysis (in the form of pyruvate) is diverted into pathways that produce acetate and lactate, which are then secreted from the cells. This process is typical of bacterial cultures experiencing rapid growth and high glycolytic flux, where the TCA cycle becomes saturated and alternative metabolic routes are used to handle the surplus carbon.
After IPTG induction at 9 hours, the concentrations of acetate and lactate in the medium remained elevated. The sustained high levels of acetate and lactate post-IPTG likely reflect the fact that once these metabolites are expelled into the culture medium, their reuptake and reintegration into intracellular metabolism is limited. Acetate and lactate, once secreted, may no longer serve as useful substrates for energy production or biosynthesis. Instead, their accumulation in the medium suggests that the bacterial cells are unable, or inefficient, at reclaiming these byproducts. Citrate (Fig. 1D), a key intermediate of the TCA cycle, showed a gradual increase in concentration, particularly after IPTG induction. This suggests that while pyruvate is being converted into acetate and lactate, a portion of it continues to fuel the TCA cycle to meet the energy demands and biosynthetic needs of the cells. Citrate's rise is indicative of active TCA cycle flux, particularly after recombinant protein production begins, as the cells require sustained energy and carbon skeletons for the synthesis of amino acids and other macromolecules. Together, the patterns observed in pyruvate, acetate, lactate, and citrate concentrations highlight the metabolic flexibility of bacterial cells in response to pyruvate availability. Before IPTG induction, excess pyruvate is channeled into overflow pathways, resulting in the secretion of acetate and lactate into the medium. Following IPTG induction, this pattern continues, but the rise in citrate suggests that a portion of the pyruvate is also directed into the TCA cycle to support the increased metabolic demands of recombinant protein production. The secretion of acetate and lactate into the culture media throughout the experiment underscores the cells' need to expel excess metabolites and maintain intracellular metabolic balance under conditions of high pyruvate flux.
Glycine (Fig. 1E) shows stable concentrations before IPTG induction, followed by a modest increase during the later time points. This suggests that glycine metabolism is linked to protein biosynthesis during the extended growth phases. Branched-chain amino acids (BCAAs), including valine (Fig. 1F), leucine (Fig. 1G), and isoleucine (Fig. 1H), show moderate increases post-IPTG induction, indicative of their essential roles in protein synthesis. The upregulation of BCAA biosynthesis following IPTG induction reflects the heightened demand for amino acids during recombinant protein production. The concentration of maleate (Fig. 1I) increases gradually, peaking at 24 hours, suggesting its involvement in energy production pathways as bacterial cells transition into the stationary phase. ATP (Fig. 1J) levels rise significantly after IPTG induction, peaking at 10 hours, indicating an increased demand for energy during the early stages of recombinant protein synthesis. This is followed by a gradual decline in ATP levels, reflecting a stabilization of energy metabolism as the bacterial culture approaches stationary phase. The ADP (Fig. 1L) profile mirrors that of ATP, peaking at 10 hours and gradually decreasing, indicating high energy turnover during protein production.
Aspartate (Fig. 1K) and arginine (Fig. 1L) show continuous increases in concentration after IPTG induction, with aspartate peaking at 24 hours. This suggests their critical roles in supporting protein synthesis and nitrogen metabolism, particularly as amino acid demand rises during recombinant protein production. Similarly, betaine (Fig. 1K) concentrations increase over time, peaking at 24 hours. Betaine is known to function as an osmoprotectant, and its accumulation may reflect bacterial stress responses during extended periods of growth and protein synthesis. Adenine (Fig. 1M), a purine nucleotide, increases steadily following IPTG induction, peaking at 36 hours, indicating heightened nucleotide turnover associated with increased DNA/RNA synthesis during recombinant protein production. Finally, 2-hydroxyisovalerate (Fig. 1N), a metabolite linked to branched-chain amino acid metabolism, shows a gradual increase over time, peaking at 36 hours. This suggests its involvement in secondary metabolic pathways, particularly as bacterial cells shift toward the stationary phase and metabolic processes become focused on maintenance rather than growth.
Taken together, these results reveal a clear metabolic reprogramming in response to IPTG induction and recombinant protein production. Pyruvate consumption, acetate overflow, and shifts in amino acid and nucleotide metabolism reflect the bacterial cells’ need to balance energy production with the synthesis of macromolecules for protein production. The most significant metabolic changes occur between 9 and 10 hours, shortly after IPTG induction, indicating the critical metabolic adjustments required to meet the demands of protein synthesis. As the bacterial culture transitions to the stationary phase (24–36 hours), metabolic activity stabilizes, as evidenced by the gradual leveling off of key metabolite concentrations.
Multivariate Analysis of Metabolite Profiles
The PCA scores plot (Figure 3A) highlights distinct separation of metabolite profiles across time points, with the first two principal components (PC1 and PC2) explaining 49.9% and 22.4% of the total variance, respectively. The samples collected at 0h and 3h cluster tightly, indicating minimal metabolic variation during the early growth (lag) phase when bacterial metabolism is primarily focused on adjusting to the medium. However, a significant shift occurs between 3h and 9h, corresponding to the addition of IPTG at the 9-hour mark to induce recombinant protein production.
Post-induction, the 10h samples (1 hour after IPTG induction) show a distinct shift along PC1, indicating rapid metabolic changes following protein induction. This is likely due to the metabolic burden of recombinant protein synthesis, which alters the demand for energy and metabolic intermediates derived from pyruvate. Notably, the 24h and 26h time points (14 and 16 hours after IPTG induction, respectively) form closely clustered groups, suggesting that the bacteria have reached a more stable metabolic state during this extended period of protein production. This clustering is consistent with the bacteria entering the stationary phase, during which growth slows and metabolism stabilizes, with a potential shift toward metabolite recycling and waste product excretion (e.g., acetate, lactate).
By 36h (26 hours post-induction), the metabolite profile has shifted slightly compared to earlier post-IPTG time points, likely due to nutrient depletion or accumulation of metabolic byproducts in the culture medium. The confidence ellipses show tight clustering of replicate samples at later time points (24h, 26h, 36h), indicating reproducibility of metabolic states during and after recombinant protein production. The PCA results clearly indicate that IPTG induction at 9h causes a dramatic metabolic shift, with the most pronounced changes occurring between 9h and 10h. These shifts reflect the metabolic adaptations required for the bacteria to support increased protein synthesis, including changes in energy metabolism, amino acid synthesis, and metabolite excretion. The pairwise scatter plots of principal components (Figure 3B) provide a more detailed view of the variance captured by the top five components. The first two principal components (PC1 and PC2) capture most of the variance, as indicated by the clear clustering of samples according to time points, especially around IPTG induction (9h-10h). Time points immediately following IPTG induction, such as 10h, show clear separation from pre-induction time points, illustrating the significant metabolic impact of initiating recombinant protein synthesis. Minor components, such as PC3 (4.9%) and PC4 (3.9%), capture additional subtle metabolic changes that may reflect variations in the levels of specific metabolites involved in later phases of growth and protein production. For instance, these components might be influenced by the accumulation of waste products like acetate or lactate as bacterial growth slows. In summary, the PCA analysis reveals significant metabolic shifts corresponding to the addition of IPTG at 9h, with marked differences in metabolite profiles between pre- and post-induction phases. The clustering of time points from 24h to 36h suggests that bacterial metabolism stabilizes during the later stages of protein production, possibly due to nutrient limitation or the plateauing of recombinant protein synthesis.
To further investigate the relationships between the metabolite profiles at different time points during recombinant protein production, a hierarchical clustering analysis was performed (Figure 3C). The clustering dendrogram groups samples based on their metabolic similarity, revealing clear temporal patterns in the data. The samples collected at 0h, 3h, and 9h (prior to or at the time of IPTG induction) form a distinct cluster, indicating that the metabolic profiles during the early phases of growth (lag and early exponential phases) are highly similar. The proximity of these time points suggests minimal metabolic reprogramming in the absence of IPTG and during the initial phase of bacterial growth in the pyruvate-supplemented medium. The 9h time point represents the moment IPTG is added to induce recombinant protein production, yet it remains metabolically similar to earlier time points, reflecting a delay in the metabolic response to IPTG. In contrast, the samples collected 1 hour after IPTG induction (10h) form a separate, highly distinct cluster from the earlier time points. This separation underscores the significant metabolic changes that occur rapidly following IPTG induction, likely driven by the metabolic demands associated with recombinant protein synthesis. These changes may include shifts in energy metabolism, amino acid biosynthesis, and central carbon metabolism, all necessary to support the increased biosynthetic activity. The 24h, 26h, and 36h time points (14, 16, and 26 hours post-IPTG induction, respectively) form a tight, well-defined cluster. This clustering indicates that once the bacteria have adapted to the protein production process, their metabolism stabilizes. The minimal differences between these time points suggest that bacterial cells have reached a stationary phase, with metabolic activities largely focused on maintaining cellular homeostasis and recycling metabolites. This stabilization is likely due to nutrient depletion and the plateauing of protein synthesis as bacterial growth slows down. The close grouping of replicates within these time points further demonstrates the reproducibility of the metabolic changes observed.
Interestingly, hierarchical clustering also reveals that some of the 9h replicates show partial overlap with 10h replicates, suggesting a transitional phase in metabolic reprogramming as cells switch from regular growth to recombinant protein production. This gradual transition is expected, as the induction of recombinant protein synthesis imposes a substantial metabolic burden, requiring time for bacterial cells to adjust their metabolic pathways to meet the demands of protein synthesis. Overall, the hierarchical clustering analysis provides a clear temporal overview of the metabolic shifts that occur during bacterial growth and recombinant protein production. The clustering of time points prior to IPTG induction, followed by distinct metabolic states post-induction, highlights the significant metabolic reprogramming that accompanies protein production. This analysis confirms that bacterial metabolism stabilizes during the later stages of growth, likely in response to nutrient limitation and the completion of protein synthesis. To complement the unsupervised analysis provided by PCA, Partial Least Squares Discriminant Analysis (PLS-DA) was performed to further explore the metabolic changes in response to IPTG induction and identify metabolites that contribute most to the observed variations across time points. The PLS-DA scores plot (Figure 3D) highlights the discriminative power of the PLS-DA model, which maximizes the separation between metabolite profiles at different time points. In this plot, Component 1 explains 48.4% of the variance, and Component 2 explains 22.3%, together capturing a significant portion of the total variability. The clear separation of the time points along these two components reflects the dynamic changes in bacterial metabolism over time, particularly in response to IPTG induction at the 9-hour mark. Pre-induction time points, 0h, 3h, and 9h, form distinct clusters, with 9h showing more separation from earlier time points, as expected, given that IPTG induction occurs at this time. Interestingly, the 9h and 10h samples (1 hour post-IPTG induction) show a sharp metabolic shift, as these two clusters are well separated, indicating a rapid metabolic reprogramming following IPTG addition. This shift is likely due to the metabolic demands imposed by recombinant protein production, which alters the cellular utilization of pyruvate and other metabolites involved in central carbon metabolism and amino acid biosynthesis. The later time points, 24h, 26h, and 36h, form tightly clustered groups, with some overlap, indicating that the bacterial metabolism reaches a more stable state after several hours of protein production. This stabilization likely reflects the adaptation of bacterial cells to the metabolic burden of protein production and the onset of the stationary phase, where growth has slowed, and nutrient consumption is reduced. The tighter clustering at these time points suggests that while metabolic activity persists, the variation between these samples is minimal, consistent with a plateau in protein synthesis and bacterial growth. The confidence ellipses around each time point cluster indicate the consistency of replicate samples, particularly for the post-induction time points, showing that the metabolic state of the cultures is highly reproducible after the initial response to IPTG. This distinct separation of time points in the PLS-DA plot, particularly around the 9h and 10h time points, confirms that IPTG induction is a major driver of metabolic variation in the dataset. The pairwise scatter plots of the first five components of the PLS-DA model (Figure 3E) provide a more detailed breakdown of how each component contributes to the separation of time points. As shown in the plot, Component 1(48.4%) and Component 2 (22.3%) account for most of the variance, with clear separation of early (0h, 3h) and post-induction time points. This further supports the observation that IPTG induction introduces significant metabolic shifts, which are well captured by the first two components. Other components, such as Component 3 (4.4%) and Component 4 (2.8%), capture more subtle variations that may be related to secondary metabolic processes or changes in less abundant metabolites that contribute to the overall metabolic profile. These components likely reflect more nuanced variations in the cellular metabolism, such as changes in specific amino acid levels or byproduct excretion (e.g., acetate, lactate) as the bacteria adapt to prolonged recombinant protein production. The pairwise component plots confirm that the metabolic changes induced by IPTG are most strongly captured by Components 1 and 2, with minimal variance explained by subsequent components. This suggests that the major metabolic reprogramming in response to IPTG induction can be largely described by the shifts along these two axes, with smaller components reflecting minor metabolic adjustments as the cultures reach the stationary phase.The PLS-DA analysis provides a supervised view of the temporal metabolic shifts in bacterial cultures during recombinant protein production. The analysis confirms that IPTG induction at 9h is a major driver of metabolic variation, leading to rapid and distinct metabolic changes as bacteria adapt to the demands of protein synthesis. The stabilization of metabolite profiles at later time points, as seen in both the scores plot and pairwise component plots, reflects the transition into the stationary phase, where nutrient consumption slows, and the metabolic state becomes more consistent.
To further understand the metabolites driving the separation of time points in the PLS-DA analysis, we calculated the Variable Importance in Projection (VIP) scores (Figure 3F). VIP scores rank metabolites based on their contribution to the PLS-DA model, identifying those most influential in distinguishing between time points, particularly in response to IPTG induction and the subsequent phases of bacterial growth and protein production.
The metabolites with the highest VIP scores are displayed in Figure 3F, indicating their importance in explaining the metabolic shifts across the time points. The top-ranking metabolites include pyruvate, acetate, glycine, citrate, isoleucine, and valine, which are directly or indirectly involved in energy metabolism, amino acid biosynthesis, and central carbon metabolism. These metabolites are particularly important in the context of bacterial growth under minimal medium conditions, where pyruvate serves as the sole carbon source.
Pyruvate, the primary carbon source used in this study, ranks as the most important metabolite, as expected. Pyruvate is not only a key energy source but also a central intermediate in multiple metabolic pathways, including the citric acid cycle (TCA cycle) and gluconeogenesis. Its importance reflects the fact that the metabolic shifts occurring after IPTG induction and during recombinant protein production are tightly linked to how pyruvate is utilized and converted into essential biomolecules. Acetate, a common byproduct of bacterial metabolism, also ranks highly in the VIP score. Its accumulation is typically a sign of overflow metabolism, where excess carbon from pyruvate is converted into acetate, especially under conditions where the TCA cycle is saturated or energy demands exceed the capacity of oxidative metabolism. This aligns with the later time points, where bacterial growth slows, and acetate may be secreted as a waste product during the stationary phase. Glycine, an amino acid involved in protein biosynthesis and the folate cycle, also contributes significantly to the separation of time points. Its importance in the VIP score suggests that shifts in glycine metabolism are closely linked to the demands of protein production, especially after IPTG induction, when the need for amino acids increases. Citrate, a key metabolite in the TCA cycle, plays a critical role in energy production and biosynthesis. The high VIP score for citrate indicates that variations in its levels are essential in distinguishing between time points, reflecting the dynamic shifts in energy metabolism as bacteria transition through different growth phases, particularly in response to the metabolic burden of recombinant protein synthesis. Isoleucine and valine, both branched-chain amino acids (BCAAs), are also among the top-ranked metabolites. These amino acids are not only essential for protein synthesis but are also linked to energy metabolism. Bacteria may upregulate BCAA biosynthesis or alter their catabolism in response to the increased demand for protein production following IPTG induction. The significant contribution of these amino acids to the PLS-DA model suggests that their levels are highly responsive to the metabolic state of the culture and the progression through growth phases. To complement the VIP score rankings, the heatmap on the right side of Figure 3F shows the relative concentrations of these key metabolites across the different time points. This visualization highlights the dynamic changes in metabolite levels over time, with warmer colors (red) indicating higher concentrations and cooler colors (blue) indicating lower concentrations. For example, pyruvate levels decrease after IPTG induction, suggesting its rapid utilization in metabolic pathways during protein production. In contrast, acetate levels increase over time, consistent with its role as a byproduct of overflow metabolism, especially as cells transition into the stationary phase. Glycine, citrate, and BCAAs (isoleucine, valine) also exhibit time-dependent fluctuations, reflecting their roles in protein synthesis and energy metabolism. The heatmap provides a clear visualization of how the concentration of these metabolites changes in response to bacterial growth and protein production. The most substantial shifts are seen post-IPTG induction, particularly at the 10h, 24h, and 36h time points, corresponding to the exponential and stationary phases of bacterial growth. The VIP score analysis identifies key metabolites that are critical in explaining the metabolic shifts observed during the bacterial growth phases and in response to IPTG induction. Pyruvate, acetate, glycine, citrate, and branched-chain amino acids play pivotal roles in supporting energy metabolism and protein biosynthesis under minimal medium conditions with pyruvate as the carbon source. The heatmap further illustrates the dynamic changes in these metabolites, providing insights into how bacterial metabolism adapts to the metabolic demands of recombinant protein production.
To understand the relationships between metabolites during recombinant protein production, we performed a Pearson correlation analysis, followed by hierarchical clustering to group metabolites based on their co-regulation patterns. These analyses provide insights into how metabolites interact with one another during different phases of bacterial growth. The heatmap in Figure 4A shows the Pearson correlation coefficients between metabolites across all time points, providing a comprehensive view of metabolite interactions. The color scale ranges from dark red (strong positive correlation, close to +1) to dark blue (strong negative correlation, close to -1), with white indicating no correlation. Several distinct clusters of highly correlated metabolites emerge, suggesting that these groups are co-regulated during bacterial growth and recombinant protein production. Pyruvate, a central metabolite in glycolysis and the citric acid cycle, shows strong positive correlations with several key energy-related metabolites such as acetate, succinate, and citrate, reflecting its pivotal role in energy metabolism and carbon flow through these pathways. The tight clustering of these metabolites suggests that they are part of a coordinated metabolic response to the use of pyruvate as the sole carbon source. In contrast, metabolites such as NADPH and NADH, which are involved in redox reactions and cellular energy balance, display negative correlations with certain amino acids like isoleucine and valine, possibly reflecting a metabolic trade-off between energy production and biosynthesis during recombinant protein production. This is consistent with the metabolic burden that bacteria face when diverting resources toward protein synthesis. The metabolites involved in amino acid metabolism, such as leucine, valine, and isoleucine (branched-chain amino acids, BCAAs), form a highly correlated cluster. These BCAAs play critical roles in protein synthesis and are tightly regulated during recombinant protein production, particularly after IPTG induction. Their strong correlations with other amino acids, such as serine and glycine, suggest that the demand for amino acids during the growth and production phases drives coordinated metabolic shifts.
Additionally, several metabolites related to amino acid biosynthesis and protein production, such as isoleucine, valine, and citrate, form another cluster. The clustering of these metabolites suggests that they are closely linked to the metabolic demands of recombinant protein synthesis, especially following IPTG induction when the demand for amino acids sharply increases. This cluster also includes succinyl-CoA, an intermediate in the TCA cycle, further indicating the importance of energy production and carbon flow through central metabolism to support protein biosynthesis.
The hierarchical clustering dendrogram (Figure 4B) further elucidates the relationships between metabolites based on their correlation profiles. The dendrogram groups metabolites into clusters according to their similarity in correlation patterns, with the vertical axis representing the distance (or dissimilarity) between clusters. Metabolites with similar roles or functions are grouped together, indicating shared regulatory mechanisms or metabolic pathways. At the top of the dendrogram, pyruvate and acetate form a distinct cluster, reflecting their central roles in energy metabolism and overflow metabolism, respectively. The strong correlation between pyruvate and acetate likely reflects the bacterial metabolism's response to the high flux of carbon through glycolysis and the citric acid cycle, leading to acetate secretion as a byproduct during later stages of growth. Interestingly, metabolites involved in nucleotide metabolism, such as adenine and uracil, cluster together, reflecting their role in DNA and RNA synthesis, which may be essential for supporting bacterial growth and protein production during the exponential phase. This cluster also includes metabolites like thymine, emphasizing the coordinated regulation of nucleotide metabolism in response to the energy and biosynthetic demands imposed by recombinant protein production. Overall, the hierarchical clustering reveals that metabolite profiles are organized into functional groups based on their roles in central carbon metabolism, amino acid biosynthesis, and energy production. The tight clustering of key metabolites involved in pyruvate metabolism and BCAA biosynthesis highlights the importance of these pathways in supporting bacterial growth and protein production under minimal medium conditions with pyruvate as the sole carbon source. The Pearson correlation heatmap and hierarchical clustering dendrogram together reveal distinct metabolic modules that are co-regulated during bacterial growth and recombinant protein production. Metabolites involved in energy metabolism, amino acid biosynthesis, and nucleotide metabolism show strong correlations, reflecting the bacterial cells' need to balance energy production with biosynthesis to support protein production. These findings provide a comprehensive overview of how central metabolism is reprogrammed in response to the metabolic demands imposed by recombinant protein synthesis.