Team:Imperial College/Modelling PURdeg
From 2013.igem.org
Line 74: | Line 74: | ||
<td>Please see derivation 1 below.</td> | <td>Please see derivation 1 below.</td> | ||
</tr> | </tr> | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
<tr> | <tr> | ||
<th>K</th> | <th>K</th> | ||
Line 108: | Line 101: | ||
<tr> | <tr> | ||
<th>k<sub>2</sub></th> | <th>k<sub>2</sub></th> | ||
- | <td>Protein production rate ( | + | <td>Protein production rate (PudA)</td> |
<td>4.7</td> | <td>4.7</td> | ||
<td>1/min</td> | <td>1/min</td> | ||
Line 114: | Line 107: | ||
<td>Please see derivation 2 below.</td> | <td>Please see derivation 2 below.</td> | ||
</tr> | </tr> | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
<tr> | <tr> | ||
<th>[Arabinose]</th> | <th>[Arabinose]</th> | ||
Line 139: | Line 125: | ||
</tr> | </tr> | ||
<tr> | <tr> | ||
- | <th>[ | + | <th>[PudA]</th> |
- | <td>Concentration of | + | <td>Concentration of PudA</td> |
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
<td>-</td> | <td>-</td> | ||
<td>mM</td> | <td>mM</td> | ||
Line 162: | Line 132: | ||
<td>-</td> | <td>-</td> | ||
</tr> | </tr> | ||
+ | |||
+ | |||
</table> | </table> | ||
<h4>PUR degradation module</h4> | <h4>PUR degradation module</h4> |
Revision as of 21:14, 2 October 2013
Contents |
Polyurethane (PUR) Degradation Module
Introduction
Both efficiencies of Polyurethane (PUR) degradation and ethylene glycol production are important for the performance of MAPLE system. We built a mathematical and deterministic model based on MATLAB extension Simbiology for Polyurethane degradation. The model contains the kinetic property of degradation enzymes that is helpful for the design of assays. As we scaled up the initial concentrations of all substrates to meet the conditions for a bio-reactor, the model can provide preliminary simulations and predictions for the MAPLE system.
Design
Objective
Here are some specific objectives for the model to achieve:
1. The model should contain the gene expression model of the degradation enzymes because the enzyme concentration determines the rate of plastic degradation. In our case for PUR degradation, we used [http://parts.igem.org/Part:BBa_K206000 pBAD strong promoter K206000] for most enzymes. We built the gene expression model based on inducible pBAD promoter, which gene expression rate can be regulated by inducer concentration.
2. The model should show the efficiency of the enzyme secretion to the culture from the cells. It's also important because the enzyme concentration in the culture depends on it. Here we used pelB secretion tag for most enzymes in order to achieve a high efficiency.
3. The model basically predict how long will take to degrade a known concentration of soluble polyurethane. It is assumed that the enzyme in our assays has the same kinetic properties as the enzyme used in the literature. The model can suggest a suitable concentration of the plastic to use in order to get good results from the assays.
4. It is known that ethylene glycol is toxic to E.coli. However, it has no clear effect on the growth of our MG1655 strain when the concentration of ethylene glycol is below 200mM. Therefore, the model should suggest a safe range of PUR concentration to avoid a high concentration (>200mM) of ethylene glycol produced.
The Model
Polyurethane (PUR) degradation involves 5 different degradation enzymes:
enzyme | source organism | biobrick | reference |
---|---|---|---|
EstCS2 | uncultured unknown bacterium (GU256649.1) | [http://parts.igem.org/Part:BBa_K1149002 BBa_K1149002] | [http://www.microbialcellfactories.com/content/10/1/41 Kang et.al 2011] |
pueA | Pseudomonas chlororaphis | [http://parts.igem.org/Part:BBa_K1149003 BBa_K1149003] | [http://onlinelibrary.wiley.com/doi/10.1111/j.1574-6968.2000.tb09056.x/pdf Stern et al., 2000] |
pueB | Pseudomonas chlororaphis | [http://parts.igem.org/Part:BBa_K1149004 BBa_K1149004] | [http://www.sciencedirect.com/science/article/pii/S0964830501000427 Howard et al., 2001] |
pudA | Comamonas acidovorans | [http://parts.igem.org/Part:BBa_K1149005 BBa_K1149005] | [http://www.sciencedirect.com/science/article/pii/S0964830598000663 Allen et al. 1999] |
pulA | Pseudomonas fluorescens | [http://parts.igem.org/Part:BBa_K1149006 BBa_K1149006] | [http://www.sciencedirect.com/science/article/pii/S0964830598000687 Vega et al., 1999] |
However, 4 of the 5 enzymes are not well characterised before, so we could't find enough kinetic data from the literature. The only well-characterised PUR degradation enzyme PudA is used in the model as an illustration of all PUR degradation enzymes. The model will be more complete when the kinetic data of the other enzymes are defined. The finished PUR degradation model is shown as below:
Key
There are two compartments which represents cells and the culture from left to right. The "cell" compartment contains the gene expression module whereas the "culture" compartment contains the degradation module. The "secretion" block that connects two compartments is the secretion module.
Parameters and assumptions
Gene expression module of PudA
Parameter | Description | Value | Units | Sources | Assumptions |
---|---|---|---|---|---|
β | maximum rate of transcription | 0.032 | mM/min | Please see derivation 1 below. | Please see derivation 1 below. |
K | Activation coefficient | 0.0031 | mM | [http://parts.igem.org/Part:BBa_K206000:Characterization] | Taking the "switch point" as the activation coefficient |
dmRNA | mRNA degradation rate | 0.10 | 1/min | [http://www.ncbi.nlm.nih.gov/pmc/articles/PMC420366/pdf/07X.pdf] | Taking the value of mRNA half-life in E.coli strain MG1655 as 6.8min. rate = ln2/half-life = ln2/6.8 = 0.10 |
dprotein | Protein degradation rate | 0.050 | 1/min | [http://jb.asm.org/content/189/23/8746.full] | There is no active degradation pathway and that dilution is the dominant way by which the protein level decreases in a cell. Rate = 1/doubling time, where doubling time = 20min. Assuming steady-state growth in LB broth as presented in paper. |
k2 | Protein production rate (PudA) | 4.7 | 1/min | Please see derivation 2 below. | Please see derivation 2 below. |
[Arabinose] | Concentration of arabinose | Initial: 0.008 | mM | ||
[mRNA] | Concentration of mRNA | - | mM | - | - |
[PudA] | Concentration of PudA | - | mM | - | - |
PUR degradation module
Simulations and Results