Team:ETH Zurich/Modeling/Chemostat Selection

From 2008.igem.org

(Difference between revisions)
(Competition in the chemostat)
 
(35 intermediate revisions not shown)
Line 12: Line 12:
==Chemostat Selection==
==Chemostat Selection==
-
As mentioned in the [[Team:ETH Zurich/Project/Overview|Project Overview]] page, one of the main concepts in our project is the possibility to select for organisms with reduced genome size. To achieve this we have thought of letting our system develop in a continuous environment with nucleotide limitation and base the selection method on the difference in growthrate. The aim of this model is to show the applicability of this concept and to determine the parameters for the operation of the chemostat in the lab.
+
As mentioned in the [[Team:ETH Zurich/Project/Overview|Project Overview]] page, one of the main concepts in our project is the possibility to select for organisms with reduced genome size. To achieve this we have thought of letting our system develop in a continuous environment with nucleotide limitation and base the selection method on the difference in growth rate. The aim of this model is to show the applicability of this concept and to determine the parameters for the operation of the chemostat in the lab.
-
The questions we are trying to adress with this model are:
+
The questions we are trying to address with this model are:
* What are the ideal settings of nutrient concentration in the feed and dilution rate in order to select for reduced strains?
* What are the ideal settings of nutrient concentration in the feed and dilution rate in order to select for reduced strains?
-
* Which is the sensitivity regarding growth rate selection?
+
* What is the sensitivity regarding growth rate selection?
-
* What are the optimal timing parameters in order to induce two subsequent rounds of restriction enzyme expression?
+
* What are the optimal timing parameters in order to induce two subsequent rounds of restriction enzyme gene expression?
-
We obtain the results for the calculated growth rates from the flux balance analysis performed in section [[Team:ETH Zurich/Modeling/Genome-Scale Model|Genome-Scale Model]] for the random deletions. These simulations performed on the different strains in our chemostat are used to simulate the development of the population concentrations in continuous culture. The goal is to build a model in which we can vary parameters like substrate feed concentration and dilution rate to optimize the operation of the continuous culture in the lab as described in [[Team:ETH Zurich/Wetlab/Chemostat Selection|Chemostat Selection]].
+
We obtain the results for the calculated growth rates from the flux balance analysis performed in section [[Team:ETH Zurich/Modeling/Genome-Scale Model|Genome-Scale Model]] for the random deletions. These simulations performed in chemostat on the different strains are used to simulate the development of the population concentrations in continuous culture. The goal is to build a model in which we can vary parameters like substrate feed concentration and dilution rate to optimize the operation of the continuous culture in the lab as described in [[Team:ETH Zurich/Wetlab/Chemostat Selection|Chemostat Selection]].
Line 27: Line 27:
== Cemostat basics ==
== Cemostat basics ==
-
A chemostat allows for continuous fermentation in a constant volume. It consists of a culture vessel with an input aperture for sterile medium and an overflow aperture to mantain the volume constant. This allows to mantain a constant population in the chemostat but still have continuous cell division.  
+
A chemostat allows for continuous fermentation in a constant volume. It consists of a culture vessel with an input aperture for sterile medium and an overflow aperture to mantain the volume constant. This allows to maintain a constant population in the chemostat but still have continuous cell division.  
Line 33: Line 33:
-
Here '''X''' is the concentration of biomass in kg/m3, '''S0''' the feed concentration in kg/m3 and '''F''' Feed in l/h.
+
Here '''X''' is the concentration of biomass in kg/m3, '''S0''' the feed concentration in kg/m3 and '''F''' Feed in m3/h.
By varying parameters like the dilution rate (D) and the concentration of nutrients in the medium, different steady states can be obtained, thus allowing to fine tune the selection process. In our case, the feed rate can be set externally and cell growth is limited by a selected nutrient (in our case nucleotides).
By varying parameters like the dilution rate (D) and the concentration of nutrients in the medium, different steady states can be obtained, thus allowing to fine tune the selection process. In our case, the feed rate can be set externally and cell growth is limited by a selected nutrient (in our case nucleotides).
Line 41: Line 41:
(Rate of accumulation)=(Input rate) - (Output rate) + (Production rate)
(Rate of accumulation)=(Input rate) - (Output rate) + (Production rate)
-
The equation for biomass and substrate are given by[[Image:Gleichung_nCompetitors.jpg|center|400px]].
+
The equation for biomass and substrate are given by
-
[[Image:Nomenclature.jpg|center|700px]]
+
-
[[Image:substrate1.png]]
+
-
Where '''X''' is the biomass in kg/m3, '''D''' the dilution rate in 1/h given by the Monod equation, '''mu''' is the growth rate in 1/h, '''Sf''' is the feed concentration, '''S''' is the substrate concentration and '''rs''' substrate rate.
+
 
 +
[[Image:Gleichung_nCompetitors.jpg|center|300px]].
 +
 
 +
[[Image:Nomenclature.jpg|center|500px]]
== Competition in the chemostat ==
== Competition in the chemostat ==
-
We have extended the basic chemostat model to allow for selection of strains with reduced genome sizes by including mass balances for the biomass of more than one strain. We couple the differential equations in order to let the strains compete for the available limiting substrate. By changing parameters like the different growth rates of the different populations or the dilution rate we can variate the parameters and optimize the selection mechanism. We are interested in the effects of differences in growth rate and dilution rates on the sensitivity of our method to validate our assumptions.
+
We have extended the basic chemostat model to allow for selection of strains with reduced genome sizes by including mass balances for the biomass of more than one strain. We couple the differential equations in order to let the strains compete for the available limiting substrate. By changing parameters like dilution rate we can variate optimize the selection mechanism. We are interested in the effects of differences in growth rate and dilution rates on the sensitivity of our method to validate our assumptions.
Line 60: Line 61:
* Strain 2 has a higher biomass concentration at the beginning, but due to the slightly smaller growth rate is eventually outgrown by strain 1.
* Strain 2 has a higher biomass concentration at the beginning, but due to the slightly smaller growth rate is eventually outgrown by strain 1.
* Strain 3 has a growth rate that is smaller than the dilution rate and thus is washed out irrespective of the initial concentration.
* Strain 3 has a growth rate that is smaller than the dilution rate and thus is washed out irrespective of the initial concentration.
 +
 +
The system of differential equations is described as follows:<br>
 +
[[Image:GleichungenChemostat.jpg|400px]]
The simulation results confirm our assumptions that a selection for strains showing different growth rates is possible under adequate conditions.
The simulation results confirm our assumptions that a selection for strains showing different growth rates is possible under adequate conditions.
 +
 +
To estimate the the parameters needed to further develop our model we used the model for flux balance analysis, as presented in section [[Team:ETH Zurich/Modeling/Genome-Scale Model|Genome-Scale Model]]. In our model the uptake rate of thymidine is 0.016 mmol/h (this is equivalent to 0.0038 g/h). With this parameter we calculate the yield for our initial population:
-
The simulation files can be found here.
+
[[Image:yield.jpg|center]]
 +
 
 +
 
 +
The yield can be directly computed with this model, as the growth rate and the thymidine uptake rate can be directly used to calculate these parameters.
 +
 
 +
 
 +
[[Image:FBA.jpg|center|300px]]
 +
 
 +
The simulation files can be found [[Team:ETH Zurich/Modeling/Download|here]].
== Pulse ==
== Pulse ==
-
The next step in our model is to simulate the effect of pulsing the expression of restriction enzyme in our populations. We start the culture with a wildtype strain and let it grow to steady state. By applying a pulse as desribed in the [[Team:ETH Zurich/Wetlab/Genome Reduction|Genome Reduction]] section we generate mutants. For this model we have assumed, that the expression of restriction enzimes will be very agressive, killing allmost everything in the culture vessel and generating a small amount of mutants with reduced genomes (this assumption can be easily changed in our model). As this step is used to estimate the sensitivity of the selection mechanism, we have varied the growth rate of the mutants by a random change in growthrate of 0.03.  
+
The next step in our model is to simulate the effect of pulsing the expression of restriction enzyme in our populations. We start the culture with a wild type strain and let it grow to steady state. By applying a pulse as described in the [[Team:ETH Zurich/Wetlab/Genome Reduction|Genome Reduction]] section we generate mutants. For this model we have assumed that the expression of restriction enzymes will be very aggressive, killing almost everything in the culture vessel and generating a small amount of mutants with reduced genomes (this assumption can be easily changed in our model). As this step is used to estimate the sensitivity of the selection mechanism, we have varied the growth rate of the mutants by a random change in growth rate of 0.03.  
-
This allows to estimate the settings for dilution rate and substrate feed concentration as the diference in growthrate, as this difference in growth rate corresponds with the sensitivity described in the literature [1].
+
This allows to estimate the settings for dilution rate and substrate feed concentration as the difference in growth rate, as this difference in growth rate corresponds with the sensitivity described in the literature [https://2008.igem.org/Team:ETH_Zurich/Modeling/Chemostat_Selection#References [1]].
In the next section we will present the results with data obtained from the [[Team:ETH Zurich/Modeling/Genome-Scale Model|Genome-Scale Model]].
In the next section we will present the results with data obtained from the [[Team:ETH Zurich/Modeling/Genome-Scale Model|Genome-Scale Model]].
-
[[Image:chemo3.png|center|frame|none|Figure 3: Selection of strains after multiple pulses]]
+
[[Image:chemo3.png|center|frame|none|Figure 3: Growth simulation of different strains after multiple pulses(schematic)]]
-
The differential equation for these three strains are described as follows
+
 
-
[[Image:GleichungenChemostat.jpg|thumb|left|500px]]
+
The preceding figure shows the possibility to select for strains after a variable number of reduction rounds. We observe, that with increasing number of deletion rounds the sensitivity of the selection method is the limiting factor. Small differences in growth lead to an increase in necessary time to separate the different strains. In this first estimate, the model shows that the strategy to be followed is an increase of selection time between pulses. This is accounted for in the framework where we include the computed growth rates for the reduced mutants.
The preceding figure shows the possibility to select for strains after a variable number of reduction rounds. We observe, that with increasing number of deletion rounds the sensitivity of the selection method is the limiting factor. Small differences in growth lead to an increase in necessary time to separate the different strains. In this first estimate, the model shows that the strategy to be followed is an increase of selection time between pulses. This is accounted for in the framework where we include the computed growth rates for the reduced mutants.
-
The simulation files can be found here.
+
The differential equation for these three strains are described as follows
== Results ==
== Results ==
 +
In Figure 4, you see the chemostat simulation results. The parameters were set as follows:<br>
 +
Ks = 0.1 kg/m^3 <br>
 +
Sf = 10 kg/m^3 <br>
 +
S0 = 8kg/m^3 <br>
 +
D = 0.55 1/h <br>
 +
Y = 157 <br>
 +
For this simulation restriction enzyme ''SmiI'' was used. After five times pulsing with restriction enzymes(where the time between two pulses was set to 100h), the growth rates of mutants vary between 0.6114 and 0.6191.
-
[[Image:5steps_Thy0_016_D5_5.jpg|center|frame|none|Figure 4: Growth Simulations of reduced-genome mutants.
+
[[Image:ForWiki 5 100h Y157 GR0 6.jpg|center|frame|Figure 4: Growth Simulations of reduced-genome mutants.
 +
The black points denote a pulse of restriction enzymes which was given every 100 h. The uptake rate of thymidine was set to 0.016 1/h. The other curves describe the growth of reduce-genome
 +
mutants. In our model, each pulse of restriction enzymes reduces the current population dramatically and leads
 +
to emergence of new mutants. The simulation assumes that from each strain maximally two mutants can be generated and their concentration is proportional to the amount of this strain before pulsing with restriction enzymes. During the whole simulation(4 pulses with duration of 100 h) the whole population is assumed to be present in the chemostat.  ]]
 +
 
 +
 
 +
<!--[[Image:5steps_Thy0_016_D5_5.jpg|center|frame|Figure 4: Growth Simulations of reduced-genome mutants.
The black points denote a pulse of restriction enzyme which was given every 100 h.  Dashed red line denotes
The black points denote a pulse of restriction enzyme which was given every 100 h.  Dashed red line denotes
-
the concentration of substrat, namely amount of Thymidin. The uptake rate of thymidin was set to 0.016 1/h. The other curvea describe the growth of reduce-genome
+
the concentration of substrate, namely amount of thymidine. The uptake rate of thymidine was set to 0.016 1/h. The other curves describe the growth of reduce-genome
mutants. In our model, each pulse of restriction enzyme reduces the current population dramatically and leads
mutants. In our model, each pulse of restriction enzyme reduces the current population dramatically and leads
-
to emergence of mutants. The simulation assumes that from each strain maximally two mutants can be generated and their conentration is proportinal to the amount of this strain before pulsing with restriction enzymes. During the whole
+
to emergence of mutants. The simulation assumes that from each strain maximally two mutants can be generated and their concentration is proportional to the amount of this strain before pulsing with restriction enzymes. During the whole
-
simulation(4 pulses with duration of 100 h) the whole population is assumed to be present in the chemostat.  ]]
+
simulation(4 pulses with duration of 100 h) the whole population is assumed to be present in the chemostat.  ]]-->
-
 
+
The simulation files can be found [[Team:ETH Zurich/Modeling/Download|here]].<br>
-
In Figure 4, you see the Simulation of Chemostat which was modeled as follows:
+
-
where the parameters are ....
+
-
(K_s = 0.08, Sf = 10, D = 0.55, Y = 2, S0 =8)
+
-
This simulation results were obtained using restriction enzyme ''SmiI''. Every 100h
+
-
a pulse of this restriction enzyme was given to the current population. The growth rates of mutants
+
-
vary between 0.6114 and 0.6174.
+
-
We have shown that a selection mechanism based on growthrate, as predicted by the Genome-Scale Model, is theoretically prossible.
+
We have shown that a selection mechanism based on growthrate, as predicted by the Genome-Scale Model, is theoretically possible.
-
We have further estimated a first feasable range for the principal parameters, including dilution rate and feed concentration. In our case and for a difference in growth rate of 0.03 we find values of D= and Sf=
+
We showed that a continuous culture is potentially a suitable environment to select for the desired organisms. Parameter ranges for the timing of pulses, dilution rate and feed concentration were identified. We developed a tool that can be refined by experimental data to predict the characteristics of the selection method.
-
We further have a first estimation of the regeneration time between pulses. Here we have also a first set of parameters for the time between pulses of 150 hours, with regeneration times twice as high every 5 to 8 pulses.
+
<!--We further have a first estimation of the regeneration time between pulses. Here we have also a first set of parameters for the time between pulses of 150 hours, with regeneration times twice as high every 5 to 8 pulses.-->
== References ==
== References ==

Latest revision as of 05:30, 30 October 2008

Contents


Chemostat Selection

As mentioned in the Project Overview page, one of the main concepts in our project is the possibility to select for organisms with reduced genome size. To achieve this we have thought of letting our system develop in a continuous environment with nucleotide limitation and base the selection method on the difference in growth rate. The aim of this model is to show the applicability of this concept and to determine the parameters for the operation of the chemostat in the lab.

The questions we are trying to address with this model are:

  • What are the ideal settings of nutrient concentration in the feed and dilution rate in order to select for reduced strains?
  • What is the sensitivity regarding growth rate selection?
  • What are the optimal timing parameters in order to induce two subsequent rounds of restriction enzyme gene expression?

We obtain the results for the calculated growth rates from the flux balance analysis performed in section Genome-Scale Model for the random deletions. These simulations performed in chemostat on the different strains are used to simulate the development of the population concentrations in continuous culture. The goal is to build a model in which we can vary parameters like substrate feed concentration and dilution rate to optimize the operation of the continuous culture in the lab as described in Chemostat Selection.


In the following sections we give a short introduction to continuous cultures and proceed to develop the used model.

Cemostat basics

A chemostat allows for continuous fermentation in a constant volume. It consists of a culture vessel with an input aperture for sterile medium and an overflow aperture to mantain the volume constant. This allows to maintain a constant population in the chemostat but still have continuous cell division.


Figure 1: Schematic of a chemostat


Here X is the concentration of biomass in kg/m3, S0 the feed concentration in kg/m3 and F Feed in m3/h.

By varying parameters like the dilution rate (D) and the concentration of nutrients in the medium, different steady states can be obtained, thus allowing to fine tune the selection process. In our case, the feed rate can be set externally and cell growth is limited by a selected nutrient (in our case nucleotides).

The specific growth rate in the chemostat can be determined by the material mass balances of biomass and substrate.

(Rate of accumulation)=(Input rate) - (Output rate) + (Production rate)

The equation for biomass and substrate are given by


Gleichung nCompetitors.jpg
.
Nomenclature.jpg

Competition in the chemostat

We have extended the basic chemostat model to allow for selection of strains with reduced genome sizes by including mass balances for the biomass of more than one strain. We couple the differential equations in order to let the strains compete for the available limiting substrate. By changing parameters like dilution rate we can variate optimize the selection mechanism. We are interested in the effects of differences in growth rate and dilution rates on the sensitivity of our method to validate our assumptions.


Figure 2: Three different strains competing in a chemostat


In this figure we observe three different strains competing in the chemostat.

  • In the first case, Strain 1 has the highest growth rate and as expected after some time reaches the highest biomass in the reactor, although the initial concentration is the smallest.
  • Strain 2 has a higher biomass concentration at the beginning, but due to the slightly smaller growth rate is eventually outgrown by strain 1.
  • Strain 3 has a growth rate that is smaller than the dilution rate and thus is washed out irrespective of the initial concentration.

The system of differential equations is described as follows:
GleichungenChemostat.jpg

The simulation results confirm our assumptions that a selection for strains showing different growth rates is possible under adequate conditions.

To estimate the the parameters needed to further develop our model we used the model for flux balance analysis, as presented in section Genome-Scale Model. In our model the uptake rate of thymidine is 0.016 mmol/h (this is equivalent to 0.0038 g/h). With this parameter we calculate the yield for our initial population:

Yield.jpg


The yield can be directly computed with this model, as the growth rate and the thymidine uptake rate can be directly used to calculate these parameters.


FBA.jpg

The simulation files can be found here.

Pulse

The next step in our model is to simulate the effect of pulsing the expression of restriction enzyme in our populations. We start the culture with a wild type strain and let it grow to steady state. By applying a pulse as described in the Genome Reduction section we generate mutants. For this model we have assumed that the expression of restriction enzymes will be very aggressive, killing almost everything in the culture vessel and generating a small amount of mutants with reduced genomes (this assumption can be easily changed in our model). As this step is used to estimate the sensitivity of the selection mechanism, we have varied the growth rate of the mutants by a random change in growth rate of 0.03.

This allows to estimate the settings for dilution rate and substrate feed concentration as the difference in growth rate, as this difference in growth rate corresponds with the sensitivity described in the literature [1].

In the next section we will present the results with data obtained from the Genome-Scale Model.

Figure 3: Growth simulation of different strains after multiple pulses(schematic)

The preceding figure shows the possibility to select for strains after a variable number of reduction rounds. We observe, that with increasing number of deletion rounds the sensitivity of the selection method is the limiting factor. Small differences in growth lead to an increase in necessary time to separate the different strains. In this first estimate, the model shows that the strategy to be followed is an increase of selection time between pulses. This is accounted for in the framework where we include the computed growth rates for the reduced mutants.

The differential equation for these three strains are described as follows

Results

In Figure 4, you see the chemostat simulation results. The parameters were set as follows:
Ks = 0.1 kg/m^3
Sf = 10 kg/m^3
S0 = 8kg/m^3
D = 0.55 1/h
Y = 157
For this simulation restriction enzyme SmiI was used. After five times pulsing with restriction enzymes(where the time between two pulses was set to 100h), the growth rates of mutants vary between 0.6114 and 0.6191.

Figure 4: Growth Simulations of reduced-genome mutants. The black points denote a pulse of restriction enzymes which was given every 100 h. The uptake rate of thymidine was set to 0.016 1/h. The other curves describe the growth of reduce-genome mutants. In our model, each pulse of restriction enzymes reduces the current population dramatically and leads to emergence of new mutants. The simulation assumes that from each strain maximally two mutants can be generated and their concentration is proportional to the amount of this strain before pulsing with restriction enzymes. During the whole simulation(4 pulses with duration of 100 h) the whole population is assumed to be present in the chemostat.


The simulation files can be found here.

We have shown that a selection mechanism based on growthrate, as predicted by the Genome-Scale Model, is theoretically possible.

We showed that a continuous culture is potentially a suitable environment to select for the desired organisms. Parameter ranges for the timing of pulses, dilution rate and feed concentration were identified. We developed a tool that can be refined by experimental data to predict the characteristics of the selection method.


References

(1) "Selection in chemostats." ,D E Dykhuizen and D L Hartl, Journal of Bacteriology, Microbiol Rev. 1983 June; 47(2): 150–168.
(2) "Biological Reaction Engineering Dynamic Modelling Fundamentals with Simulation Examples", Dunn, Irving J. / Heinzle, Elmar / Ingham, John / Prenosil, Jiri E., Wiley-VCH Verlag GmbH,. Weinheim, 2003