Introduction

This model was designed to investigate the effect of induction on overall bacterial cellulose (BC) yield and it was used to assist the wetlab by proposing means for improving BC yield within an experiment time frame. The ATCC53582 Gluconacetobacter can be induced by isopropyl β-D-1-thiogalactopyranoside (IPTG) to alter its enzymatic activities, increasing the proportion of nutrients diverted to BC production. After being induced, the ratio of nutrients directed to BC production to that to cell growth will increase and a variation in induction time will inevitably result in a change in overall BC yield. With a set of ordinary differential equations (ODEs), the model predicts the maximum BC yield and cell population profiles for different induction time points. As a result, the model provides an effective approach to approximate the optimal induction time under specific experimental conditions, which significantly enhances the efficiency of BC production (please refer to E. coli, induction experiments for further information).

The model consists of two parts: the first part approximately predicts the profiles for BC production, substrate (glucose) consumption, BC-producing cell population and mutant (dormant) cell population at each induction time. The model simulation is experiment-specific as optimal induction time varies with both initial substrate concentration and the time frame allowed for the experiment. As an output, part one of the model will generate 1) a plot of maximum BC yield achievable within the experiment time frame against induction time, and 2) a plot for the time necessary to achieve maximum BC yield against induction time.

Based on the estimation provided by part one of the model, the second part takes advantage of the integrated form of the ODEs and generates a series of more precisely calibrated profiles at a given induction time (i.e. the optimal induction time approximated by model part one). The profiles of substrate consumption, cell population and BC production can be easily interpreted and inform the wetlab experiments with a refined estimation so that the efficiency of BC production process can be better understood.

Mathematical Background

Model part one (ODE simulation):

• Dependent Variables:
  • P_plus(t): population of BC producing cells at time t
  • P_minus(t): population of mutant cells (dormant cells producing no bacterial cellulose) at time t
  • BC(t): BC yield at time t
  • S(t): substrate (glucose) concentration at time t
  • O(t): oxygen level at time t
• Parameters:
  • Kd: mutation rate from BC producing cells to mutant (dormant) cells
  • Vp: maximum rate of cell population change achieved by the system, parameter in Michaelis-Menten kinetics
  • Kp: Michaelis constant in rate equation for cell population change, parameter in Michaelis-Menten kinetics describing the enzymatic activities involved in substrate consumption
  • K_O2: Michaelis constant in rate equation for cell population change, parameter in Michaelis-Menten kinetics describing the enzymatic activities involved in oxygen consumption
  • Kd_plus: death rate of BC producing cells
  • Vc: maximum rate of BC yield achieved by the system, parameter in Michaelis-Menten kinetics
  • Kc: Michaelis constant in rate equation for BC production, parameter in Michaelis-Menten kinetics describing the enzymatic activities involved in substrate consumption
  • K_up: together with O(t), the term (K_up+O(t))/K_up captures the upregulating effect of oxygen on bacterial cellulose yield
  • Kd_minus: death rate of mutant (dormant) cells
  • K_growth: substrate consumption rate of cell population
  • K_BC: substrate consumption rate of BC production
• ODEs:
  • Rate of change of BC producing cells population:
  • Rate of change of mutant (dormant) cells production:
  • Rate of BC production:
  • Rate of substrate (glucose) consumption:

Model part two (integrated form):

The simulation of part one of the induction model provides an approximation of the optimal induction time for which overall BC yield is maximised. This induction time will then serve as an input to the second part of the model, where an integrated form of the rate equations are used to predict more detailed profiles for BC production, substrate consumption and cell population.

• Solutions:
  • 
  • 
  • 
  • 
• Dependent Variables:
  • P(t): total cell population, all cells produced since starting time
  • P_plus(t) : BC producing cell population.
• Parameters
  • K_maintain: maintenance coefficient, a measurement of the substrate consumption rate to maintain the BC producing cell population (viable cells). This is an added feature in model part two.
  • Experiments were set up to determine P₀, P_max, µ_BC, µ_P, P_plus,0, P_plus,max, µ_plus, K_growth, K_BC, and K_maintain where these parameters obtained from experiments characterise the strain ATCC53582 specifically. The measurements taken to characterise the parameters include: the dry mass of bacterial cellulose being produced, glucose (substrate) level just below the pellicle and cell concentration.

Results and Conclusions

Model part one (ODE simulation):

*Key: +cell population = BC producing cell population; -cell population = dormant (mutant) cell population (not producing BC)

Figure 1 demonstrates the change in four dependent variable profiles as per induction time, simulated for 50 time points starting after 60 minutes with one hour as the time interval. The time frame simulated in the model is 50 hours long, with the y axis labelled with arbitrary units since the model is intended to provide relative values, rather than absolute values, which will be calculated in second part of the model. This part of the model (part one) provides an approximation of the optimal induction time by comparing maximum BC yield (within experiment time frame) at each induction time, and the corresponding time needed to achieve maximum yield.

As shown in figure 2, a maximum BC yield is exhibited when induced after around four hours of growth, during the growth (log) phase of cell cycle. If induced earlier than this time point, as illustrated in the movie above, cell growth will be delayed during lag phase and therefore a competent cell population will not be established to yield satisfactory amount of BC within the experiment time frame, which in this simulation is 50 hours. The maximum BC yield stays at a constant level when induction occurred after 10 hours and it is probably due to the fact that nutrient (substrate) is being depleted after around 10 hours.

Figure 3 illustrates the variation in time needed to achieve maximum BC production within the experiment time frame as per induction time. Similarly, the fact that time needed to achieve maximum BC yield stabilises after 10 hours of growth is also resulted from the depletion of nutrient (substrate).

Figure 4 illustrates the efficiency of BC production as a trade-off between maximum BC yield and corresponding time needed to achieve such yield at each induction time. It is observed that an induction occurred after seven hours will not only result in a relatively high yield but also save time in achieving this high yield. Again, due to depletion of nutrient, any induction occurred after 10 hours converges at the same point as the one induced at 10 hour time point.

Model part two (integrated form):

The above presented simulations are also adaptable to larger time scales (i.e. days). With appropriate boundary conditions, the provided model can be used to simulate the dependent variable profiles at a given induction time (e.g. the approximated optimal induction time point). This is explained further below:

For the cell concentration measurements, no difference has been made between BC producing cells and dormant (mutant) cells, as simulated in part one of the model. Total cell population thus considers an accumulation of all cells produced since starting point, and the viable cell population is an approximation of cell population in the oxygenated layer that are producing BC.

As described in the mathematical background section, the nutrient profile presented here also takes cell maintenance consumption into consideration, which is a further refined feature specific for the second part of the model.

References

1. M. Hornung; M. Ludwig; A.M. Gerrard; H.-P. Schmauder Optimizing the Production of Bacterial Cellulose in Surface Culture: Evaluation of Substrate Mass Transfer Influences on the Bioreaction (Part 1) Eng. Life Sci. 2006, 6, No. 6, 537–545

Appendix: Code

The code for this model can be found here: Induction_model_imperial_igem_2014.m