Team:DTU-Denmark/Reactor Model

From 2013.igem.org

(Difference between revisions)
(The Reactor)
(The Reactor)
Line 26: Line 26:
-
For a CSTR perfect mixing is assumed, which means that the concentrations of all compounds are the same everywhere in the reactor and also at the exit point. If the system is run in steady state the so called design equation can be used to calculate the necessary volume of the reactor:
+
For a CSTR perfect mixing is assumed, which means that the concentrations of all compounds are the same everywhere in the reactor and also at the exit point. If the system is run in steady state the so called <b>design equation</b> can be used to calculate the necessary volume of the reactor:
[[File:DTU_mod_2e.png|190px|center]]
[[File:DTU_mod_2e.png|190px|center]]

Revision as of 12:37, 4 October 2013

Contents


The Reactor

Summary

We want to get a first estimate of how large a reactor filled with our transformed E. coli cells must be to treat ammonia polluted wastewater. Because Mutant 2 requires anaerobic conditions, we plan to build two reactors in series and find with the help of modeling that the second should be around 10 times the size of the first. This is under the assumption that all transformed proteins will be expressed at the same level, so expressing proteins at different levels is something we should look into if we want to avoid a large difference in reactor sizes.

Methods

With our genetically engineered E. coli we want to treat ammonia polluted waste water. The fastest way to do this would be in a continuous process like a continuous stirred tank reactor (CSTR). Ideally this process would generate an amount of nitrous oxide that is so large that the energy needed for running the system can be delivered by catalytically decomposing the nitrous oxide. We are planing to test a simple CSTR system but hold the engineered E. coli back in the reactor, so that the outgoing water is free of bacteria. Also we will need to use two CSTRs in sequence because the second mutant requires anaerobic conditions. A picture of two CSTR reactors in series can be found in Figure 1. Another interesting solution would be to grow a bacterial biofilm on solid carriers as shown in Figure 2, a technology known as moving bed biofilm reactor.

DTU modeling 2CSTR.png

Figure 1: Two CSTRs in series.


DTU mod 2c.png

Figure 2: In moving bed biofilm reactors the biofilm is grown on carriers as shown above. The picture is taken from http://www.hydrotech-group.com/en/products/industrial-wwtp/mbbr/.

A CSTR can be balanced using the following equation [4]:

DTU mod 2a.pngDTU mod 2d.png


For a CSTR perfect mixing is assumed, which means that the concentrations of all compounds are the same everywhere in the reactor and also at the exit point. If the system is run in steady state the so called design equation can be used to calculate the necessary volume of the reactor:

DTU mod 2e.png

We can split the molar flow rate of species j into concentration of species j times flow. The reaction rate is dependent on the substrate concentration and can be described with Michaelis-Menten kinetics in our case: DTU modeling Rj.png. So we can calculate the volume of the first reactor for series of different flow rates and substrate concentrations. To include the right enzyme concentration we need to estimate how many cells per litre we will have in the reactor and how many enzymes each cell will contain. It was estimated that the cells take up 40% of the volume in the reactor and that they contain 1000 copies of each protein of the pathway (see kinetic modeling) resulting in a enzyme concentration of 8.456 nmol/L. However this equation does not allow a complete removal of the substrate ammonia from the system, therefore we look at substrate outlet values around 10% of the substrate inlet values (this corresponds to 90% of NH3 conversion to NO2 ). The matlab function describing the design equation and the script to evaluate it for different flows and substrate concentrations can be found under scripts below. The result of evaluating the design equation for different flows and substrate concentrations is shown in Figure 6.

The second reactor takes the effluent water of the first reactor as feed water. Assuming 90% of the NH3is converted to NO2 this means that the feed water of the second reactor will contain 5.2848 mM NO2 (if the NH3 concentration in the feed water of the first reactor was 5.872 mM as described above). Toreach an initial estimation of the necessary reactor size we assume that the conversion of NO to N2O happens instantaneously. This is reasonable because NO is toxic to the cells and should therefore not accumulate, further the reaction of NO2 to NO is slower than that of NO to N2O. Whith this assumption we can calculate the necessary volume of the second reactor in the same way as for the fist. The result is shown in Figure

Results and Discussion

References

[1] WK Keener and DJ Arp. Kinetic studies of ammonia monooxygenase inhibition in Nitrosomonas europaea by hydrocarbons and halogenated hydrocarbons in an optimized whole-cell assay. Applied and Environmental Microbiology, 59(8): 2501–2510, 1993.

[2] Serena Rinaldo. Biology of the Nitrogen Cycle. Francesca Cutruzzola, 2007 (37-55).

[3] SW Snyder and TC Hollocher. Purification and some characteristics of nitrous oxide reductase from paracoccus denitrificans. Journal of Biological Chemistry, 262: 6515–6525, 1987.

[4] H. Scott Fogler. Essentials of Chemical Reaction Engineering. Prentice Hall, 2010.