Team:Imperial College/Modelling PLAdeg
From 2013.igem.org
Contents |
Polylactic acid (PLA) Degradation Module
Introduction
The efficiency of Polylactic acid is important for the performance of MAPLE. Therefore, we chose a strong enzyme Proteinase K which is suggested by many literature as an efficient PLA degrading enzyme. We performed several assays in order to determine the kinetic properties of Proteinase K. The PLA degradation model was then built based on the experimental results.
Objective and Design
1. With the defined kinetic properties of Proteinase K. The model can predict the time needed to degrade a certain amount of PLA. The information is important as it estimates the efficiency of the MAPLE system when degrading the large amount of plastic.Therefore, further improvements can be made if the system is not efficient.
2. The model also involves the gene expression model of the Proteinase K with an inducible promoter. Therefore, gene expression can be regulated by adjusting the inducer concentration. The overall plastic degradation can be regulated by the gene expression.
3. The secretion model is also contained in the model which assumes the efficiency of enzyme secretion to the culture.
Specifications of the model
Polylactic acid (PLA) degradation involves proteinase K:
enzyme | source organism | biobrick | reference |
---|---|---|---|
Proteinase K | Engyodontium album | [http://parts.igem.org/Part:BBa_K1149002 BBa_K1149002] | [http://www.ncbi.nlm.nih.gov/nuccore/X14689.1] |
The overall PLA degradation model is shown as below:
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 Proteinase K
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 (Proteinase K) | 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 | - | - |
[Proteinase K] | Concentration of Proteinase K | - | mM | - | - |
1.Derivation of the maximal expression rate,β
- Average molecular weight (Mw) of a base pair = 660g/mol[http://www.geneinfinity.org/sp/sp_dnaprop.html][http://www.lifetechnologies.com/uk/en/home/references/ambion-tech-support/rna-tools-and-calculators/dna-and-rna-molecular-weights-and-conversions.html]
- Average mass of a base pair = 660g/mol x 1.66x10-24 = 1.1x10-21g
- Volume of an E.coli cell = 1µm3[http://kirschner.med.harvard.edu/files/bionumbers/fundamentalBioNumbersHandout.pdf] = 1x10-15L
- BioBrick assembly plasmid pSB1C3 is a high copy number plasmid (100-300 copies per cell)[http://parts.igem.org/Part:pSB1C3?title=Part:pSB1C3]
- assume 200 copies per cell
- ∴ concentration of the gene per cell = N x 200 x 1.66x10-6mM, where N = number of base pairs
- ∴ concentration of the gene BDH2 (N = 768) in the volume of an E.coli cell is = 0.25mM
- Transcription rate in E.coli = 80bp/s[http://kirschner.med.harvard.edu/files/bionumbers/fundamentalBioNumbersHandout.pdf] = 80 x 1.66x10-6mM/s = 80 x 1.66x10-6 x 60mM/min = 7.97x10-3mM/min
- ∴ Rate of mRNA_Proteinase K production under the control of pBAD = 7.97x10-3 ÷ 0.25 = 0.032/min
- Average molecular weight(Mw) of an amino acid(aa)= 110g/mol[http://www.genscript.com/conversion.html][http://www.promega.com/~/media/Files/Resources/Technical%20References/Amino%20Acid%20Abbreviations%20and%20Molecular%20Weights.pdf]
- Average mass of an amino acid = 110g/mol x 1.66x10-24=1.83x10-22g/L
- Translation rate = 20aa/s = (20 x 1.66x10-5 x 60)mM/min = 0.020mM/min
- BDH2 comprises of 256aa[http://www.uniprot.org/uniprot/Q2PEN2&format=html]
- ∴concentration of BDH2's aa in the volume of an E.coli = 1.66x10-5mM x 256 = 4.25x10-3mM
- ∴ Rate of protein production = 0.020 ÷ 4.25x10-3 = 4.7/min
Degradation
The reaction equation of the PLA degradation is:
[Polylactic acid]+[Proteinase K]= 660 [lactic acid] + [Proteinase K]
Assumptions:
We assumed 1 mole of polylactic acid can produce 660 moles of lactic acid.
The molecular weight of a single polylactic acid monomer is 90 g/mol [http://www.chemspider.com/Chemical-Structure.592.html]whereas the molecular weight of the solid polylactic acid is around 59500 g/mol.[http://www.sciencedirect.com/science/article/pii/S0014305711003582]Therefore, the short-chain polylactic acid consists approximately 660 monomers. 660 molecules of lactic acid will be produced by degrading one chain of polymer.
We also assumed a simple Michaelis-Menten mechanism for proteinase K
The parameters for the kinetic equations are:
Parameter | Description | Value | Units | Sources |
---|---|---|---|---|
Km | Michaelis constant | 0.032 | mM | - |
Vmax | Maximum velocity | 2.472 | mM/min | - |
Proteinase K | Concentration of the enzyme in assays | 30 | mg/L | - |
Kcat | turnover number | 2348.4 | 1/min | See derivation below |
Mw | Molecular weight of proteinase K | 28.5 | KDa | [http://www.sigmaaldrich.com/etc/medialib/docs/Sigma/Datasheet/2/p4850dat.Par.0001.File.tmp/p4850dat.pdf] |
Derivation:Turnover number (Kcat) = Vmax*Mw/[proteinase K] = 2348.4
The efficiency of Secretion is assumed to be 90% secretion over 2 hours.[http://www.ncbi.nlm.nih.gov/pmc/articles/PMC1251600/] The rate of secretion in the model is therefore:
rate of secretion = 0.9[concentration of Proteinase K]/120 (mM/mins)
Simulations and Results
[File:PLAD.png]