Team:Imperial College/BioPlastic Recycling: PHB

From 2013.igem.org

Revision as of 02:35, 5 October 2013 by Jps09 (Talk | contribs)


Module 2: Plastic Fantastic

Module 2 summary.jpg





Plastic Fantastic is a complete P(3HB) bioplastic recycling platform, where P(3HB) is degraded into monomeric form and then re-polymerised back into de novo P(3HB) for future applications.

  • Overview

  • Specifications

  • Design

  • Modelling

  • Testing & results

  • Future work

Overview

testing

We have developed a system by which PHB can be recycled when products made from it come to the end of their life. In order to do this, we have engineered E.coli to break down P3HB extracellularly to release the monomers of 3-hydroxybutyrate (3-HB). We have designed the metabolic connection between P3HB degradation and synthesis so that 3-HB can be used as a carbon source. We are testing a putative 3-HB permease to increase the rate of uptake of the molecule into the cell. An intracellular 3-HB dehydrogenase will then form acetoacetate from 3-HB which the cell can utilise. We have thoroughly characterised the phaZ1 PHB depolymerase enzyme and demonstrated that it degrades PHB. Our models predict that use of the enzyme on a bioreactor scale will be effective in degrading bioplastic.

We are also establishing a second bioplastic recycling platform for polylactic acid (PLA)

This aspect of our module is in collaboration with the Yale iGEM team.

Yale are working to biologically synthesise PLA. We have designed E.coli to secrete proteinase K to form lactic acid monomers, forming half of a second bioplastic recycling system. Lactic acid is used in the chemical synthesis of PLA and so this system could also be implemented in an industrial setting.

Specification

Degradation of P(3HB)

Our bacteria should resist any toxic effects that are associated with P(3HB) or 3HB

The bacteria should secrete a functional (P3HB) degrading enzyme

Synthesis of P(3HB)

The bacteria should take up and internalise 3HB from the surrounding media

The bacteria should be able to utilise P(3HB) as a sole carbon source

Degradation of PLA

Our bacteria should resist any toxic effects that are associated with lactic acid, the breakdown product of PLA

Our bacteria should secrete a functional PLA degrading enzyme

P(3HB) Degradation: Pathway

The phaZ1 enzyme from Ralstonia eutropha can break down PHB. This bacterium uses PHB as storage molecule and expresses PhaZ1 when it needs to mobilise the nutrients from the polymer inside the cell. In our system, the enzyme will be secreted.

Module2 4-10.jpg

P(3HB) Degradation: Biobrick design

BBa_K1149010: Extracellular expression of phaZ1, PHB depolymerase enzyme. The construct design incorporates the RBS 0034 and is regulated by a strong xylose-inducible promoter BBa_I741018. For extracellular secretion, we fused pelB secretion tag BBa_J32015 to the N terminus of the protein.

PhaZ1_brick_diagram.jpg

Synthesis of P(3HB) from 3HB: Pathway

For the Plastic Fantastic recycling module, we designed a modifications to the PHB synthesis pathway of the Resourceful Waste module. The additional steps connect 3-Hydroxy-Butyrate to the downstream metabolite of the synthesis: acetoacetyl-CoA where from the next two steps are carried out by phaB and phaA enzymes, as described in Module1. This pathway contains both new enzymes that we are introducing to E.coli and also endogenous enzymes, such as AtoA/D. Because of this, the pathway is in connection with the cells normal metabolism and our models were of great use for predicting the flux through the pathway and informing the design.

Plastic synthesis from 3HB monomers. 3HB is taken up into the cell by the putative permease. The 3HB is then converted to P(3HB) by the pathways shown. The genes in green are those we have expressed within MG1655. Our system takes advantage of the endogenous ato operon.

Synthesis of P(3HB) from 3HB: Biobrick design

3HB Permease:

For 3HB import, we designed a biobrick for the expression of a Putative permease, found in the literature [1]. We optimised the sequence for E.coli and planned experiments with the construct.

Construct for bdh2 dehydrogenase expression

BBa_K1149050: Intracellular expression of bdh2, 3HB dehydrogenase enzyme. We have used a strong arabinose inducible promoter BBa_K206000 in front of the operon for controlled expression of the enzyme; this enzyme is only required by the cell if 3HB is present. We have included superfolder GFP BBa_K515005 in the operon with phaZ1 so that we can monitor gene expression from the promoter via fluorescence measurements. (see our corresponding data under results)

Bdh2_brick_diagram.jpg

Reaction_bdh.jpg

PLA degradation: pathway and biobrick design

Poly(lactic acid) (PLA) depolymerases are categorized into protease-type and lipase-type (7). We examined the literature and found sequences for one enzyme from each type.

Protease-type

  • Proteinase K sequence from Tritirachium album, hydrolyzes poly(L-lactic acid) (PLLA) only.

Lipase-type

  • Cutinase-like enzyme (CLE) sequence from Cryptococcus sp. strain S-2 preferentially hydrolyzes poly(D-lactic acid) (PDLA), but can also cleave poly(L-lactic acid) (PLLA).

We considered the safety aspects of using these two enzymes, since their sequences were originally from Risk Group 2 (RG2) organisms. We submitted our safety information to iGEM for review. Safety forms were approved on October 2nd, 2013 by the iGEM Safety Committee.

We did not use these organisms only the parts and the final sequences were codon optimised for expression in E. coli.


PLA Biobricks.jpg
Proteinase K and CLE, Poly lactic acid (PLA) degradation mechanisms. Adapted from Kawai et. al 2011

1.0 P(3HB) Degradation Model

The P(3HB) degradation model was built and simulated in MATLAB Simbiology in the similar fashion as that seen in the modelling section of Module 1. This model encompasses gene expression of PhaZ1 (P(3HB) depolymerase) controlled by a xylose inducible promoter[2]and the kinetics of PhaZ1 degrading the polymer P(3HB) into shorter oligomers or 3HB monomers.

Imperial College Phaz circ.JPG

Imperial College Phaz key.JPG

1.1 Ordinary differential equations

Genetic expression and regulation of PhaZ1
Imperial College PhaZ1 gene.JPG
Values, sources and assumptions
Parameter Description Value Units Sources Assumptions/Notes
β maximum rate of transcription 0.019 mM/min Please see derivation 1 below. Please see derivation 1 below.
n Hill coefficient 2.0 dimensionless Estimation. Not on parts registry [3] Didn't manage to find an appropriate source with such value. Assume cooperativity of degree 2.
K Activation coefficient 23.5 mM Personal conversation with Christopher Hirst, Department of Bioengineering, Imperial College London Taken from data for xylF in MOPS with 0.4% glucose with extra xylR on the plasmid, however the plasmid is low copy and also the presence of the xylR gene on the plasmid definitely alters the expression level. Data fitted to Hill equation. This value was the most accurate we managed to obtain, though a rough one as 3% glucose in LB was used in the expreiment.
dmRNA mRNA degradation rate 0.10 1/min [4] Taking the value of mRNA half-life in E.coli strain MG1655 as 6.8min. rate = ln2/half-life = ln2/6.8 = 0.10
k2 Protein production rate 2.9 1/min Please see derivation 2 below. Please see derivation 2 below.
kprotein Protein degradation rate 0.050 1/min [5] Assume the only mechanism by which protein degradation can occur is dilution due to cell division (doubling time of E. coli MG1655 is 20 mins so degradation rate 1/20mins). Assume enzyme doesn't degrade extracellularly via any means throughout the time course of the simulation. Despite the portion degraded inside of the cell, 100% leaks out from the periplasm to the outside.
[Xylose] Concentration of xylose Initial: mM - -
[mRNA] Concentration of mRNA - mM - -
[PhaZ1] Concentration of PhaZ1 - mM - -
Derivations

1.Derivation of the maximal expression rate,β

  • Average molecular weight (Mw) of a base pair = 660g/mol[6][7]
  • Average mass of a base pair = 660g/mol x 1.66x10-24 = 1.1x10-21g
  • Volume of an E.coli cell = 1µm3[8] = 1x10-15L
    • ∴Mass concentration = Imperial College Mass conc cal beta.JPG
    • ∴Molar concentration of 1 base pair in the volume of E.coli = Imperial College Molar conc cal beta.JPG = 1.66x10-6 mM
  • BioBrick assembly plasmid pSB1C3 is a high copy number plasmid (100-300 copies per cell)[9]
    • 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 phaZ1 (N = 1245)[10] in the volume of an E.coli cell is = 0.41mM
  • Transcription rate in E.coli= 80bp/s[11] = 80 x 1.66x10-6mM/s = 80 x 1.66x10-6 x 60mM/min = 7.97x10-3mM/min
  • ∴ Rate of mRNA_phaZ1 production under the control of pBAD = 7.97x10-3 ÷ 0.41 = 0.019/min

2.Protein production rate of phaZ1, k2

  • Average molecular weight(Mw) of an amino acid(aa)= 110g/mol[12][13]
  • Average mass of an amino acid = 110g/mol x 1.66x10-24=1.83x10-22g/L
    • ∴Mass concentration of one aa in the volume of an E.coli = Imperial College Mass conc cal beta aa.JPG = 1.83x10-6g/L
    • ∴Molar concentration of one aa = Imperial College Molar conc cal aa.JPG = 1.66x10-5mM
  • Translation rate = 20aa/s = (20 x 1.66x10-5 x 60)mM/min = 0.020mM/min
  • PhaZ1 comprises of 414aa[14]
    • ∴concentration of PhaZ1's aa in the volume of an E.coli= 1.66x10-5mM x 414 = 6.87x10-3mM
  • ∴ Rate of protein production = 0.020 ÷ 6.87x10-3 = 2.9/min


1.2 Enzyme kinetics

PhaZ1 reaction with P(3HB)
Imperial College PhaZ1 kinetics.JPG
Values, sources and assumptions
Parameter Description Value Units Source Assumptions/Notes
kcat,PhaZ1 Turnover number of PhaZ1 with P(3HB) as substrate 8.3x105 1/min [15] P(3HB) depolymerase A purified from Pseudomonas lemoignei . Took specific activity with substrate P(3HB) from data in paper: "19.1 U/mg protein x1000" = 19100 µmol/(min.mg)=(19100x10-6x1000)mol/(min.g)= 19.1 mol/(min.g). ∴ kcat = 19.1 mol/(min.g) x 43496 g/mol = 8.3x105, where 43496 = Mw of PhaZ1[16]
Km Michaelis constant 4.3x10-5 mM [17][18] Assume the substrate P(3HB) has a molecular weight (Mw) of 2.3x106Da. Using the Km given in mg/ml (value=0.1), Km in mM = (0.1x1000)/2.3x106
[PhaZ1] Concentration of PhaZ1 - mM - -
[P(3HB)] Concentration of P(3HB) - mM - -

1.3 Simulation results

Due to the fact that the model built was a single-cell model, the raw output from Simbiology would show the degradation of P(3HB) into 3HB by one cell. Please refer to explanations and assumptions below the graphs for how the simulation was scaled up.

(i) Single-cell simulation (units: mM)

Imperial College 1phaz.JPG

(ii) Single-cell simulation (units: g/L)

Imperial College 2phaz.JPG

(iii)Scaled-up macroscopic simulation (units: g/L)

Imperial College 3phaz.JPG

Graphs (i) to (ii) unit conversion
  • To convert from mM to g/L:
    • Molecular weight (Mw) of 3HB = 104.1g/mol[19]
    • ∴3HB concentration in mM x 10-3 x 104.1g/mol gives concentration in g/L
Graphs (iii) scaling-up assumptions
  • Assume 1.5x1012cells in a 1L broth, calculated from value obtained from BioNumbers. A range was given: 1-2x109 cells/ml, so 1.5x109 was taken. [20]
    • Conditions and considerations accompanying this value are:
      • growth medium: LB broth
      • Cultivation temperature: 37°C
      • wild-type E. coli K-12 strain MG1655 (NB: though not exactly our engineered one, it is the correct strain)
  • Assume stationary growth phase so that the number of cells remain constant throughout the time course of simulation.
  • Assume that the polymer is thin and small such that that the rate of degradation is omly minimally affected by the dimensions of the substrate.
  • Note the simulation time was increased to one day to show how the P(3HB) would degrade on a longer time scale

P(3HB) degradation model:Download

Additional script:Download

2.0 P(3HB) production from 3HB monomers

In module 2, the main objective is the bioplastic recycling. We have a detailed explanation of the P(3HB) synthesis model and the optimization in P(3HB) production yield. Our synthetic pathway of converting monomeric 3HB to polymeric P(3HB) is described in module one. In MAPLE system, both glucose from mixed waste and 3HB monomers are the inputs to produce P(3HB). In this section, we are going to introduce a customized metabolic model that consisted our synthetic pathway with a metabolic model built by Angela Dixon. [21] The combined metabolic model can predict the P(3HB) production with both 3HB and glucose as sources.

Introduction of the metabolic model

The metabolic model predicts the interaction of our pathway with endogenous pathways. In synthetic biology, an ideal system is fully orthogonal. However, if the synthetic pathway uses some metabolites from the metabolism pathways, the metabolic fluxes (the rate of conversion between metabolites) will be affected. If the cell is overloaded by losing too much metabolites, it would either reject the synthetic pathway or burn out. As for the optimization of PHB production in our MAPLE system, metabolic analysis need to be carried out to estimate the effects on the cell metabolism from increases in the PHB production.

Methods

We altered the metabolic model by Angela Dixon in her paper “Predictive Mathematical Model for Polyhydroxybutyrate Synthesis in Escherichia coli”. [22] The metabolic model is built in Simbiology, which is an extension in MATLAB especially for biological system. The model consists two metabolic pathways: Glycolysis and TCA cycle. The reactions in the metabolic pathways are defined by ordinary differential equations and kinetic parameters. The model also contains a synthetic pathway that produce P(3HB) from acetyl-coA. The pathway consists three enzymes which are PhaA, PhaB and PhaC. The reasons we chose this model:

1. The model is based on ODEs instead of FBA (flux balance analysis), which is a common method of metabolic analysis. However, FBA cannot determine the real time concentrations of the metabolites whereas it can be done by models with ODEs.

2. Only single gene deletion or addition can be performed in FBA. More than that, the metabolic model here can perform modification on any single or multiple genes and reactions through Simbiology.

3. The FBA model is static, it cannot analyze the dynamic interactions between our synthetic pathways with the metabolic pathways. As both our P(3HB) synthesis model and the metabolic model are in the same platform. It was really easy to integrate our pathway into the metabolic model.

4. The strain in the metabolic model is E.coli K12, the strain MG1655 we are using is also included in K12 group. The models can therefore be more easily combined.

In our project, we have a different pathway:

1. We use 3HB as one of the input.

2. ato system is present in the E.coli MG1655 strain.

3. We removed phaA in our case because we use 3HB as the major source to produce P(3HB), it will avoid too much acetyl-coA be taken from the metabolic pathways.

Therefore, we substituted the original pathway with our synthetic pathway. Moreover, we added our gene expression models for the enzymes into the metabolic model as well. In terms of the consistency of the model, we changed the units of our synthetic pathway into micro molar and seconds.

Here is the combination of our model with the metabolic model

MetaSyn.png

The pathway with the black border is our synthetic pathway.

KEYmeta.png

We also carried out a parameter estimations of the initial concentrations of several substrates in our synthetic pathway. The reason is we have single-cell scaled magnitude for them originally (fM level), but there are much higher magnitudes in any of the species in the metabolic pathways (mM level). Therefore, we scaled up the pathway in order to maintain the consistency of the model. As for the metabolites which are involved in our pathway, we cloned the blocks of them from the metabolic pathway in Simbiology and deliberately put them into our pathway. That means our pathway dynamically interact with the metabolic pathways.

Here is the table for the updated initial concentrations for the species in our synthetic pathway: All the kinetic data and initial concentrations of the metabolites in the metabolic pathways are referenced in the paper[23]:

Substrates Concentration Units Sources Assumptions
3HB 5000 uM - The initial input to the system
Acetoacetate 1000 uM - The ato system in E.coli suggest a favourable conversion from acetoacetyl-coA to acetoacetate. We assumed the concentration of acetoacetate is 10 times larger than acetoacetyl-coA[24]
Acetoacetyl-coA 100 uM [25] The ato system in E.coli suggest a favourable conversion from acetoacetyl-coA to acetoacetate. [26]
NADH 200 uM [27] -
NAD+ 1200 uM [28] -
Acetyl-coA 5000 uM [29] -
Coenzyme A 250 uM [30] -
Acetate 40000 uM [31] Assay data from the literature which the ato system efficiency was measured
3HB-coA 0 uM - One of the intermediate product of the pathway
NADPH 200 uM [32] -
NADP+ 1200 uM [33] -
PHB 0 uM - The final product of the pathway

ATP analysis:

In our optimization of P(3HB) production module, we discovered that increasing in the PhaB expression can efficiently increase the yield of P(3HB). It’s important to see whether increase in the PhaB expression would dramatically affects the metabolic pathways or not. As ATP is the essential energy carriers in the cell, we carried a parameter scan for the variation in ATP concentration under different concentrations of PhaB.

ATPPHAB.png

There is an approximately 5.3% drop in the concentration of ATP when increasing the phaB concentration from 50uM to 200uM. It suggests that increasing the flux in our pathway may reduce the energy transportation in the cell. However, it’s still unsure about the effect of 5.3% drop in ATP concentration on the cell growth.

Large scale PHB production and implementation in MAPLE

In our MAPLE system, the bioplastic P(3HB) is produced from both 3HB monomers and glucose. 3HB is the degradation product from the wasted P(3HB)plastic and glucose is the product from the degradation of mixed waste. Therefore, the overall PHB production rates need to be determined by the model.

GLU3HB.png

The result shows that the PHB yield can reach 70g/L after one day with both inputs whereas about 30g/L with glucose. It proves that using both glucose and 3HB can double the yield of P(3HB). The production rate 70g/L/Day is also relatively high for industrial scale of bioplastic production.

Assumption in the production

    • Conditions and considerations accompanying this value are:
      • volume of the system is 1 litre
      • growth medium: LB broth
      • Cultivation temperature: 37°C
      • wild-type E. coli K-12 strain MG1655 (NB: though not exactly our engineered one, it is the correct strain)
  • Assume stationary growth phase.
  • Assume that the amount of P(3HB) accumulated during the time course of simulation doesn't exceed the cell's capacity or ability to contain it.
  • Assume that all initial concentrations of the species in metabolic pathways and the synthetic pathways are in a large population level. The concentrations are maintained by approximately 1.5x1012cells in a 1L broth, calculated from value obtained from BioNumbers. A range was given: 1-2x109 cells/ml, so 1.5x109 was taken. [34]

Video: How to run a simulation in the metabolic model (please refresh if it's not here)

Metabolic model:Download

3.0 PLA Degradation Model

Introduction

The efficiency of degrading Polylactic acid (PLA) 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. It will produce useful information when degrading large amount of plastic.

2. The model also involves the gene expression model of the Proteinase K with an inducible promoter. Gene expression can be regulated by adjusting the inducer concentration.

3. A secretion module is 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 BBa_K1149002 [35]

The overall PLA degradation model is shown as below:

PKPLA.png

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.020 mM/min Please see derivation 1 below. Please see derivation 1 below.
K Activation coefficient 0.0031 mM [36] Taking the "switch point" as the activation coefficient
dmRNA mRNA degradation rate 0.10 1/min [37] 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 [38] 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) 3.0 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[39][40]
  • Average mass of a base pair = 660g/mol x 1.66x10-24 = 1.1x10-21g
  • Volume of an E.coli cell = 1µm3[41] = 1x10-15L
    • ∴Mass concentration = Imperial College Mass conc cal beta.JPG
    • ∴Molar concentration of 1 base pair in the volume of E.coli = Imperial College Molar conc cal beta.JPG = 1.66x10-6 mM
  • BioBrick assembly plasmid pSB1C3 is a high copy number plasmid (100-300 copies per cell)[42]
    • 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 proteinase K (N = 1194) in the volume of an E.coli cell is = 0.40mM
  • Transcription rate in E.coli = 80bp/s[43] = 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.40 = 0.020/min
2.Protein production rate of Proteinase K, k2
  • Average molecular weight(Mw) of an amino acid(aa)= 110g/mol[44][45]
  • Average mass of an amino acid = 110g/mol x 1.66x10-24=1.83x10-22g/L
    • ∴Mass concentration of one aa in the volume of an E.coli = Imperial College Mass conc cal beta aa.JPG = 1.83x10-6g/L
    • ∴Molar concentration of one aa = Imperial College Molar conc cal aa.JPG = 1.66x10-5mM
  • Translation rate = 20aa/s = (20 x 1.66x10-5 x 60)mM/min = 0.020mM/min
  • Proteinase K comprises of 398aa[46]
    • ∴concentration of Proteinase K's aa in the volume of an E.coli = 1.66x10-5mM x 398 = 6.61x10-3mM
  • ∴ Rate of protein production = 0.020 ÷ 4.25x10-3 = 3.0/min

Degradation

The reaction equation of the PLA degradation is:

[Polylactic acid]+[Proteinase K]= 660 [lactic acid] + [Proteinase K]

Assumptions:

1. We assumed 1 mole of polylactic acid can produce 660 moles of lactic acid.

2. The molecular weight of a single polylactic acid monomer is 90 g/mol [47]whereas the molecular weight of the solid polylactic acid is around 59500 g/mol.[48]Therefore, the short-chain polylactic acid consists approximately 660 monomers. 660 molecules of lactic acid will be produced by degrading one chain of polymer.

3. Assumed the polylactic acid only contains one length of chain, however in reality the plastic chain has a range of lengths.

4. Considerations of temperature change and PH change are not included in the model. Therefore, the assumption is that they remain constant during the degradation

5. We also assumed a simple Michaelis-Menten mechanism for proteinase K.

MMPLA.png

"[Proteinase K]" here is the concentration of the enzyme Proteinase K in mM, whereas "[Polylactic acid]" represents the concentration of polylactic acid in mM.

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 [49]

Derivation:Turnover number (Kcat) = Vmax*Mw/[proteinase K] = 2348.4

The efficiency of Secretion is assumed to be 90% secretion over 2 hours.[50] The rate of secretion in the model is therefore:

rate of secretion = 0.9[concentration of Proteinase K]/120 (mM/mins)

Simulations and Results

The initial concentration of Polylactic acid is 40mM, which can be converted to 2380g/L. The derivation is: 40*59500/1000 = 2380g/L. There is approximately 2kg of polylactic acid in the 1 litre bioreactor, which is a really large amount. The simulation result is:

PLAD.png

Results and explanations:The simulation shows a very fast degradation of large amount of polylactic acid by proteinase K. 2kg/L of the plastic can be completely degraded in 6 minutes. The reason for fast degradation is that we have relative large kinetic parameters for the enzyme. Polylactic acid was not used as the substrate in our kinetic assays. Therefore, we made an assumption that polylactic acid in the model has similar chain length as that of the substrate in the assays. In reality, the PLA plastic is in solid form which enzymes can only degrade the plastic from the surface. The efficiency of degradation will be massively reduced due the 3D geometry. A 3D model is needed for future improvements of plastic degradation. The kinetic data of degrading solid form of PLA also need to be defined.

PLA degradation model:Download

Our chassis can grow with P(3HB) and 3HB

Summary

Our chassis can grow and survive in reasonable concentrations of P(3HB) and 3HB. To increase tolerance to 3HB we have re-designed our system to include expression of Bdh2 (BBa_K1149050, 3HB dehydrogenase), which breaks down 3HB to acetoacetate. Acetoacetate toxicity was examined and these data are being utilised along with our metabolic models to determine the effects of differing levels of 3HB and acetoacetate on bioplastic (P3HB) production.

MG1655 transformed with either native phaCAB or empty vector (control) plasmid were grown for 6h in LB containing 3-HB or P-3HB to assess any potential toxicities associated with these bioplasrics.

3-hydroxybutyrate (3HB)

Toxicity of 3-hydroxybutyrate in LB with phaCAB and EV expression. 3HB was tested at 5 concentrations for toxicity - 1 μM, 100 μM, 1 mM, 10 mM and 237 mM. Between 10 mM and 237 mM, there is a clear toxicity effect observed. Data points show final time point after 6h growth for each concentration. Growth was at 37°C with shaking over 6h. Error bars are SEM, n=4. Figure made by Imperial College London 2013 iGEM.

Poly(3-hydroxybutyrate) P(3HB)

Toxicity of Poly(3-hydroxybutyrate) in LB with phaCAB and EV expression. P3HB toxicity was only tested at one concentration in LB - 31 μg/L, it was also tested at 31 mg/L however the resulting OD was over 1, so this data was not considered accurate. Growth was at 37°C with shaking over 6h. Error bars are SEM, n=4. Figure made by Imperial College London 2013 iGEM.

Conclusion:3HB is toxic above 10mM and P(3HB) had no effect on growth at 31 μg/L. These data can be utilised to determine the maximum level of 3HB which should be present in our final bioreactor design. We have shared these data with our modelling team in order to determine maximum P(3HB) production at or below 10mM (3HB). In light of these data we have modified our design to include the expression of bdh2 (BBa_K1149050)</b>, which is a 3HB dehydrogenase that can convert 3HB into acetoacetate, a common metabolite in E.coli. Expression of bdh2 should overcome the toxicity limit of 3HB, since E. coli can metabolise acetoacetate</b>To this end we tested whether an increase in acetoacetate would present any toxicities.

Acetoacetate

Toxicity of acetoacetate in LB with phaCAB and EV expression. Acetoacetate was also tested for toxicity at 4 concentrations - 1 μM, 100 μM, 1mM and 287 mM. In phaCAB and EV this manifested itself between 1 mM and 287 mM. Data points show final time point after 6h growth for each concentration. Growth was at 37°C with shaking over 6h. Error bars are SEM, n=4. Figure made by Imperial College London 2013 iGEM.

Conclusion:Acetoacetate toxicity lies somewhere between 1mM and 287 mM.

Bioplastic degradation: Phaz1 enzymatic assays

Enzyme activity of PHB depolymerase (phaz1)

It can be seen from the Western Blot results that phaz1 BBa_K1149010 was being expressed. To show that this enzyme has esterase activity, colourimetric assays were performed using the substrate analog para-Nitrophenyl butyrate. When the ester bond in this substrate is cleaved, 4-Nitrophenol is released. This is accompanied by an increase in absorbance at the wavelength 405 nm and a colour change from colourless to yellow. The concentration of 4-Nitrophenol produced could then be calculated with the Beer-Lambert Law, as the extinction coefficient of 4-NP at 405 nm is 18,000 M-1 cm-1. This experiment was performed with both the crude cell lysate and purified PHB depolymerase. Our data shows that this enzyme is definitely active.


para-Nitrophenyl butyrate is cleaved by PHB depolymerase to release 4-Nitrophenol. This is accompanied by an increase in absorbance at 405 nm. Figure adapted from TU Darmstadt 2012 iGEM
When the substrate para-Nitrophenyl butyrate is added to the reaction mixture of purified PHB depolymerase (phaz1) in phosphate buffer, there is a colour change from colourless to yellow. Figure by Imperial College London iGEM 2013.










Cell lysate assay

After 48 hours of growing and inducing the phaz1 culture, the cells were lysed by probe-sonication, spun down and resuspended in 50mM Tris-HCl buffer. A reaction mixture containing 5 µL of crude phaz1 lysate, 4 µL of para-Nitrophenyl butyrate solution and 1 mL of 50 mM Tris-HCl pH 7.4 buffer was incubated in the Eppendorf BioSpectrometer for 570 seconds, whilst the absorbance at 405 nm was automatically recorded every 30 seconds.


The increase in absorbance that accompanies the cleavage of para-Nitrophenyl butyrate by PHB depolymerase (phaz1). Figure by Imperial College London 2013 iGEM


The concentration of 4-Nitrophenol released by PHB depolymerase (phaz1) activity. Empty Vector and Substrate alone were used as negative controls. Figure by Imperial College London 2013 iGEM






















The above graph shows that greater esterase activity occurs when the phaz1 cell lysate is in the reaction mixture. The graph shows that there is also esterase activity occuring in the Empty Vector and Substrate alone reaction mixtures, but this is due to the imidazole present in the Tris-HCl buffer, which acts as a general base catalysis. (6)


Purified enzyme assay

Since the phaz1 BBa_K1149010 construct contains a His tag, the protein could be purified by metal affinity chromatography. The colourimetric assay was carried out, as with the cell lysate. 100 mM potassium phosphate (pH 7.4) buffer was used instead of the Tris-HCl buffer, to eliminate the esterase activity caused by imidazole in the buffer used previously. Every 30 seconds for 30 minutes the absorbance at 405 nm of the reaction mixture was recorded. The reaction mixture contained 1 ml of phosphate buffer, 2.5 µL of purified PHB depolyermase (stock concentration 0.230 mg/ml) and para-Nitrophenyl butyrate to make final concentrations of 50 µM, 100 µM, 150 µM, 200 µM and 250 µM of substrate.


The purified PHB depolyermase (phaz1) is definitely active. The graph shows that more 4-Nitrophenol is produced with greater concentrations of para-Nitrophenyl butyrate. Figure by Imperial College London 2013 iGEM



The above graph shows how the concentration of 4-Nitrophenol produced is greater when the para-Nitrophenyl butyrate substrate concentration is greater. This data shows that we have successfully purified active PHB depolyermase.

SDS PAGE gel shows that PHB depolymerase was purified from the crude cell extract. Figure by Imperial College London 2013 iGEM






Bioplastic degradation: Phaz1 degrades P(3HB)

Clearing Zone Assay

Clearing zone experiment was conducted to semi-quantitatively test the ability of PhaZ1 to break down P3HB. PhaZ1 cell lysate was tested on a P3HB LB agar plate. The protocol can be found here.

PhaZ1 showing P3HB degradation ability by generating a clear zone on a P3HB LB agar plate. A. No difference between empty vector cell lysate and PhaZ1 cell lysate right after they were pipetted into the wells. B. PhaZ1 in the cell lysate started to clear P3HB around the well after 1 day. C. PhaZ1 created a clear zone around the well after 3 days, whereas it was still cloudy around the empty vector cell lysate well. Figure by Imperial College London 2013 iGEM

3HB assay

In order to show that phaZ1 is functional, we have measured the amount of 3HB monomers produced due to the degradation of the polymer, when incubated with the enzyme. We have repurposed a 3HB colorimetric assay kit that is used for 3HB detection from animal tissue and bodily fluid samples. Please see the protocol and the data sections for details.

PhaZ1 treatment of P3HB sample significantly increases the amount of 3HB monomer going into sulution (p=0.0015). Both of the sample were incubated at 37°C sahking for 15 hours.

It is clear from our results that enzyme treatment produces a lot more 3HB than the control. We do expect to see some background 3HB from the untreated sample since it was in the shaking incubator at 37°C in the Buffer solution. PHB is a bioplastic and degrades spontaneously as well as enzymatically. However, the spontaneous degradation occurs slowly. We have made the process more efficient and industrially scalable in a recycling plant.


How efficient is PhaZ1 at breaking down P3HB? From the assay, we have calculated the amount of bioplastic that could be broken down by a kilogram of Phaz1. Below are the calculations and assumptions used to make these conclusions:

  • Assumptions
  • #1 Reaction is at a steady state
  • #2 Amount of 3HB detected is equivalent to the mass of P3HB broken down - i.e. 1g P3HB gives 1g 3HB


  • Concentration of purified Phaz1 - 0.23 µg / µL
  • Volume of purified Phaz1 required - 10 µL
  • Mass of Phaz1 = 10 x 0.23 - 2.3 µg


  • Concentration 3HB - 226.33 mM 3HB
  • Initial P3HB concentration - 50 mg / 800 µL ~ 6%
  • Mw (3HB) - 126.09 g/mol
  • Reaction volume - 800 µL
  • Reaction time - 15h


  • 3HB produced = (226.33/800)*126.09 = 35.67 mg / 15h
  • Actual mass of 3HB produced from 2.3 µg Phaz1 - 35.67 mg
  • 3HB produced per hour = 35.67/15 = 2.378 mg/h
  • 3HB produced per hour per microgram Phaz1 = 2.378/2.3 = 1.03 mg/h/µg


  • Scaled up to 1 tonne.
  • 1 tonne = 10^9 mg.
  • Thus: 1.03 tonnes/h/kg

Conclusion: 1kg of Phaz1 will in 1h breakdown 1 tonne of P3HB.


Our chassis can utilise 3HB

bdh2 improves growth on 3HB

We have previously observed that MG1655 cells can survive in minimal media that contains 3HB as sole carbon source. This is evidence for that some uptake mechanism and a metabolic pathway is active at a low level in the cells. We hypothesised that the 3HB dehydrogenase could improve growth since it can convert 3HB into acetoacetate, a common metabolite in E.coli.

In order to test this, we run a growth experiment with various 3HB concentrations where the growth of cells containing bdh2 (BBa_K1149050) and Empty vector was recorded on a 96 well plate. We have calculated the growth rates from this data and plot it on the graph below.

3HB sole carbon source growth in minimal media with phaCAB and EV. MG1655 phaCAB were grown in supplemented and minimal M9 media. Relative to the empty vector control, there was no significant difference in growth as t-test gave p = 0.8072 > 0.05. However, an ANOVA of the data gave (F 3,8 = 6.589, p < 0.0149), thus the null hypothesis, there is no difference between M9M and M9S must be rejected. Indeed if we look closer, we see that both M9M EV and M9S EV are significantly different from another (p = 0.045), as are M9M phaCAB and M9S phaCAB (p = 0.0284). Data points show final time point after 6h growth for each concentration. Growth was at 37°C with shaking over 6h. Error bars are SEM, n=4. Figure made by Imperial College London 2013 iGEM.
3HB sole carbon source growth in M9S with P3HB synthesis pathway genes. The 3 phaCAB promoter constructs are represented with native phaCAB, constitutive phaCAB and hybrid phaCAB. In addition bdh2 was tested to see whether its presence increased growth for the cells. This was completed with 3 3HB concentrations - 100 μM, 10 mM and 100 mM. ANOVA analysis shows that there is no significant difference between any of the constructs (F 4,10 = 0.1875, p < 0.9396). Data points show final time point after 6h growth for each concentration. Growth was at 37°C with shaking over 6h. Error bars are SEM, n=4. Figure made by Imperial College London 2013 iGEM.
Bdh2 with no pelB secretion tag growth. Bdh2 MG1655 cells were grown in M9S media to gauge growth, the growth does not differ from the control as p = 0.5543. Data points show final time point after 6h growth for each concentration. Growth was at 37°C with shaking over 6h. Error bars are SEM, n=4. Figure made by Imperial College London 2013 iGEM.
Bdh2 with no pelB secretion tag fluorescence. MG1655 cells with bdh2 were grown in M9S and then tested for fluorescence by measuring their GFP output. There is significant fluorescence as p < 0.0001. Data points show final time point after 6h growth for each concentration. Growth was at 37°C with shaking over 6h. Error bars are SEM, n=4. Figure made by Imperial College London 2013 iGEM.

Conclusion:We can conclude that there is a slight increase in the growth rate of bdh2 containing cells at 10 000 uM 3HB. This suggests that bdh2 is able to function as expected and produces acetoacetate which is used by the cell`s central metabolic pathways for growth. At 100uM or 100 000 uM, we did not observe any differences in growth. This could be be because the rate limiting step at low concentrations could be the uptake of 3HB rather than it`s conversion to acetoacetate. We hope that we could observe an increase in growth if we added the putative permease we designed to the system. At higher level, there is a drop in growth rates in both bdh2 and control cells, probably because of toxicity issues. In addition, fluorescence is highly elevated in bdh2 compared to empty vector control.

Completing the cycle

We can degrade the PHB made by our own bacteria

With this experiment, we demonstrate that the material that we purified in various experiments is Poly-Hydroxy-Butyrate bioplastc. One of the samples analysed was extracted from bacteria growing on waste, therefore we made PHB bioplastic directly from waste.

We can close the loop and degrade our own bioplastic, with our own enzyme, PhaZ1.

The bioplastic produced by our bacteria contains PHB and can be broken down into monomers by PhaZ1 enzyme (expressed by our bacteria)


The amount of 3HB monomers freed from the polymer is more in case of all PhZ1 treated samples, than the untreated sample. There is less 3HB produced from the samples that were purified from bacteria than from the pure polymer, purchased from Sigma. This could be because I put less amount in the assay to start with, because there was less available. I put about a 0.01-0.03 g of our own and 0.5 grams of Sigma-PHB.

Please see our protocols and data pages for methods and details of the data.

Bioplastic degradation: Proteinase K enzymatic assays

Proteinase K kinetic activity

According to the literature, protease K breaks down poly L-lactic acid as well as fibrous proteins, such as silk fibroin. In the active site of proteinase K is a catalytic triad containing serine, histidine and aspartic acid. L-lactic acid is analogous to L-alanine. Poly L-lactic acid, therefore, can bind to and be catalysed by the catalytic triad, as it mimics silk fibroin which is rich with L-alanine residues. Proteinase K cannot catalyse the breakdown of poly D-lactic acid, because the catalytic triad will not accommodate it. [51] [52]

To obtain parameters to use in our PLA degradation model that kinetic activity of proteinase K was measured with a substrate analogue, N-succinyl Ala-Ala-Pro-Leu p-nitroanilide (AAPL p-NA). Cleavage of AAPL p-NA releases 4-nitroaniline (4-NA), which can be recorded by an increase in absorbance at wavelength 405nm. For this assay, proteinase K and AAPL p-NA were purchased from Sigma.

Molecular structure of N-succinyl Ala-Ala-Pro-Leu p-nitroanilide. Image by Sigma
.

The extinction coefficient of 4-nitroaniline is 9,500 M-1s-1 at wavelength 405 nm, according to Sigma. Using the Beer-Lambert law, the concentrations of 4-NA were calculated and plotted against time.

The absorbance change recorded during the catalysis of different concentrations of AAPLp-NA by Proteinase K. Figure by Imperial College London 2013 iGEM
The concentration of 4-NA produced from different starting concnetrations of AAPLp-NA. Figure by Imperial College London 2013 iGEM

The initial velocities were calculated as the gradient between 0 and 30 seconds, when time = 0 seconds. These were plotted against the concentration AAPLp-NA and graphpad was used to estimate the Michaelis-Menten fit. It determined the Vmax as 46.31 µM/s and the Km as 19.28 µM.

Michaelis-Menten plot for proteinase K activity cleaving AAPLp-NA at 37ºC and pH 7.4. Figure by Imperial College London 2013 iGEM

Bioplastic degradation: ProteinaseK degrades PLA

Proteinase K degrades PLA cup

Proteinase K can degrade PLA. Please see our SEM data page for full description of expreiment and all SEM images.

Increased surface pitting and roughness indicated degradation of PLA surface after enzyme treatment. Figure by Imperial College London 2013 iGEM

References

Papers Referenced
  1. ANDERSON A, DAWES E. Occurrence, Metabolism, Metabolic Role, and Industrial Uses of Bacterial Polyhydroxyalkanoates. Microbiol Rev 1990 DEC;54(4):450-472.
  2. Harding KG, Dennis JS, von Blottnitz H, Harrison STL. Environmental analysis of plastic production processes: Comparing petroleum-based polypropylene and polyethylene with biologically-based poly-beta-hydroxybutyric acid using life cycle analysis. J Biotechnol 2007 MAY 31;130(1):57-66.
  3. Kim S, Dale BE. Energy and Greenhouse Gas Profiles of Polyhydroxybutyrates Derived from Corn Grain: A Life Cycle Perspective. Environ Sci Technol 2008 OCT 15;42(20):7690-7695.
  4. Jendrossek D, Handrick R. Microbial degradation of polyhydroxyalkanoates. Annu Rev Microbiol 2002;56:403-432.
  5. Philip S, Keshavarz T, Roy I. Polyhydroxyalkanoates: biodegradable polymers with a range of applications. Journal of Chemical Technology and Biotechnology 2007 MAR;82(3):233-247.
  6. Thomas C. Bruice , Thomas H. Fife , John J. Bruno , Patricia. BenkovicHydroxyl Group (V)1 and Imidazole (X)2 Catalysis. The General Base Catalysis of Ester Hydrolysis by Imidazole and the Influence of a Neighboring Hydroxyl Group. J. Am. Chem. Soc., 1962, 84 (15), pp 3012–3018
  7. Fusako Kawai, Kosuke Nakadai, Emiko Niskioka, Hajime Nakajima, Hitomi Ohara, Kazao Masaki, Haruyaki Lefuji.Different enantioselectivity of two types of poly (lactic acid) depolymerases toward poly (L-lactic acid) and poly (D-lactic acid).Polymer Degradation and Stability Volume 96, Issue 7, July 2011, Pages 1342–1348
  8. Jun Liao, Manfred K Warmuth, Sridhar Govindarajan, Jon E Ness, Rebecca P Wang, Claes Gustafsson, Jeremy MinshullEngineering proteinase K using machine learning and synthetic genes.BMC Biotechnology 2007, 7:16 doi:10.1186/1472-6750-7-16

Protocols

Content

MSDS

Chemical reagents

Some of the reagents we used were slightly toxic to human. In order to reduce the risks of using toxic reagents, we always wear gloves and labcoats. For some toxic reagents, we performed in fume cupboard. Click on the reagents to see the MSDS.

1. Poly(R)-3-hydroxybutyric acid.

2. Proteinase K.

3. Sodium 3-hydroxybutyrate.

4. Arabinose.

5. Acetoacetate.

We are changing the pelB secretion tag to OsmY-fusion partner for better secretion

As seen in the Western Blot results, some of the enzymes designed for secretion were not secreted. These proteins may have been stuck in the periplasmic space. Thus we are optimising enzyme secretion by building a new construct with OsmY BBa_K892008 fusion for secretion. OsmY is naturally secreted from E.coli and will carry our PUR esterase enzyme with. Once in the culture medium, we think that the enzyme is going to be functional in the fusion form but we included a TEV site in between the domains that can be cleaved by TEV protease if needed.

Future PHB PLA bricks.jpg

The linker sequence we have designed between the domains is T G S E N L Y F Q G S (ACCGGCAGCGAGAACCTGTACTTCCAAGGCAGC) and includes the 7 amino acids of the TEV site as well as 4 other amino acids for flexibility.

This construct will be a substitute for BBa_K1149010 and BBa_K1149008 biobricks which are the PHB depolymerase PhaZ1 and ProteinaseK with the pelB seretion tag.

We will construct phaCB

This construct requires removal of phaA which links central metabolism to the P(3HB) synthesis process. By removing this we may optimise plastic production by decreasing a means by which the intermediates before P(3HB) synthesis are converted into products of the central metabolism by a non-favourable equilibrium.


Proteinase K activity

The Western Blot results show that proteinase K was expressed, however, no activity was detected in a cell lysate assay and no proteinase K could be purified from the cells. The proteinase K cultures were induced overnight with 6 mM (final concentration) of arabinose. According to the literature, proteinase K self-digests (8). Perhaps the induction strength and time was too long. Next time we will reduce the time and concentration of induction for proteinase K.

Our Sponsors

TueSponsorsEppendorf.png 125px Invitrogen.jpg Geneart.jpg CSynBI.JPG