Team:Calgary/Project/OSCAR/Desulfurization

From 2012.igem.org

Revision as of 08:28, 3 October 2012 by Rmmayall (Talk | contribs)

Hello! iGEM Calgary's wiki functions best with Javascript enabled, especially for mobile devices. We recommend that you enable Javascript on your device for the best wiki-viewing experience. Thanks!

Desulfurization

Why Remove Sulfur?

Sulfur is the third most abundant element in crude oil (Ma, 2010), and when sulfur containing hydrocarbons are burned they release S02 and S03 gasses into the atmosphere. Not only does this reduce the efficiency and value of our product, but it also contributes to global warming, acid rain, and various health issues due to the pollution (Reichmuth et al., 2000). Strict regulation on sulfur in fuels are now in place and low-sulfur gasoline is mandated across all of Canada (Source: Environment Canada). To upgrade the quality of our fuel we need to remove the sulfur but keep the hydrocarbon backbone for combustion.

Our Vision

Though a few pathways for biodesulfurization exist in the microbial world, most involve the destruction of part of the carbon skeleton (an example would be the Kodama pathway). This would effectively reduce the quality of our product. With this in mind the pathway we have chosen is the 4S pathway found in Rhodococcus spp. It has been characterized and shown to remove sulfur from the model substrate dibenzothiophene (DBT) and convert it to 2-hydroxybiphenyl (2-HBP) in a non-destructive manner. DBT and its derivatives make up 70% of the organic sulfur compounds found in crude oil (Ma 2010), and are also some of the most difficult to remove through chemical means. By using the 4S pathway we will be able to upgrade our fuel and remove recalcitrant compounds at the same time.

Figure 1: The 4S Desulfurization Pathway, showing the desulfurization of the model compound DBT by DszA, DszB, DszC, and DszD.

4S pathway

Four enzymes are involved in the 4S pathway, 3 of which are directly involved in the conversion of DBT to 2-HBP. Dibenzothiophene monooxygenase (DszC) is responsible for the first two steps of the pathway, converting DBT to DBT-sulfoxide and finally to DBT-sulfone (DBTO2) through the addition of 2 oxygen atoms to the sulfur atom. DBT-sulfone monooxygenase (DszA) then carries out the next step in the pathway, producing 2-hydroxybiphenyl-2-sulfinic acid (HBPS) through addition of a final oxygen to the heteroatom. This causes cleavage of the chemical bonds at the sulfur, breaking the ring and converting the compound from a 3-ring structure to a 2-ring structure. HBPS is then converted to the final product of the 4S pathway by HBPS desulfinase (DszB), producing 2-HBP. At this point, the sulfur has been released from the hydrocarbon in the form of sulfite.

The first three steps of the 4S pathway require FMNH2 and subsequently reduces the reductive power of the cell. WIn order to regain this power an oxidoreductase (DszD) uses NADH to recycle the FMNH2, allowing the reaction to proceed. Without DszD the desulfurization pathway would grind to a halt.

The dszA,B, and C genes form an operon on the pSOX plasmid of R. erythropolis, while dszD is found in the chromosome. Naturally this pathway is slow, however using synthetic biology approaches this process can be optimized.

Our Approach

1) Find the genes!

We isolated the plasmid containing the dsz genes from a desulfurising environmental isolate of Rhodococcus using a modified miniprep pocedure. As the native promoter has been shown to be repressed by various sulfur-containing compounds (Li et al., 1996), we designed primers for just the coding sequences of the A, B, and C genes. As these genes all have some illegal cutsites in them we constructed them into the PSB1C3 vector and started our mutagenesis protocol.

Figure 2: Rhodococcus cells were grown in a modified M9 media (refer to protocols for details). Samples were taken out at different time points and were run through GCMS to detect the amount of DBT.

Figure 3: The peak in this mass spectra demonstrates presence of DBT based on its molecular weight of 184 g/mol. This peak is based on the average of our samples at retention time of 13.9 minute (refer to previous graph).

Figure 4: Rhodococcus cells were grown in a modified M9 media (refer to protocols for details). Samples were taken out at different time points and were run through GCMS to detect the amount of DBT.

2) Mutagenesis: Biobrick Compatability and Increasing DszB Activity

In total the dszABC genes had 7 PstI sites and 1 NotI site that needed to be mutated for the biobrick standard. The primers were designed such that the site was removed without the amino acid being changed. It was also shown that a point mutation changing dszB's 63rd amino acid from Y to F increases the activity of the protein (Oshiro et al., 2007). This mutation was also included in the mass mutagenesis we undertook.

3) Replacing DszD with HpaC & Introducing Catalase

In the 4S pathway, FMNH2 is consumed by the first three steps in the pathway carried out by DszA and DszC. When this compound is consumed, the speed of oxidation of DBT is reduced as the pathway begins to halt. A high concentration of FMNH2 on the other hand produces hydrogen peroxide (H2O2) that is lethal to cell in high concentrations (Gala´n et al. 2000).

To overcome these problems in our system, two plans of optimization were implemented. First, the hpaC gene coding for another flavin reductase in E. coli W is used in place of the dszD gene which is native to the pathway. It has been shown that HpaC improves the rate of desulfurization by 7-10 times when compared to DszD (Gala´n et al., 2000). Second, a catalase gene was included in the final construct in order to take care of any excess hydrogen peroxide formed during the recycling of FMNH2 and prevent stress on the cell ( pLacI-katG-LAA).

Figure 5: Diagrammatic representation of the full "optimization circuit", consisting of the oxidoreductase HpaC and a catalase (KatG).

Results

In order to determine if the KatG-LAA protein was functional, catalase activity was tested using a turbitity comparison between E. coli expressing native catalase levels (carrying only the pLacI biobrick) and cultures of pLacI-katG-LAA. Cultures were grown overnight, induced with 100 uM IPTG, and subsequently subcultured into various concentrations of hydrogen peroxide (0 mM, 1 mM, 5 mM, and 10 mM). These cultures were then grown overnight. It was found that the negative control exhibited no growth after 12h at 1 mM peroxide, however cultures with induced expression of catalase were turbid after 12 h of growth at this concentration (Fig. 3). This demonstrated the ability of the catalase to protect the cells from excessive peroxide concentrations.

Figure 6: Catalase Assay. Overnight cultures of pLacI and pLacI-KatGLAA were innoculated into 0 mM, 1 mM, 5 mM, and 10 mM peroxide. Cultures were grown overnight and turbidity was observed. It was found that at 1 mM of peroxide, cultures with just the lacI promotor perished, however when KatG-LAA was expressed, the cells survived.

HpaC activity was tested by comparing cell lysate ability to consume NADH between cultures over-expressing HpaC to a control culture that is not. This is done by monitoring the decrease in absorbance at 340 nm (Kamali et al., 2010). In order to do this, cultures of pLacI-hpaC and pLacI-dszB were grown up overnight in LB with appropriate antibiotics. Following this, protein expression was induced with IPTG, after which the assay was carried out as described in the following figure and on the protocols page.

Figure 7: HpaC Assay. Cultures of pLacI-hpaC and pLacI-dszB were grown up overnight in LB with appropriate antibiotics. The following morning, cells were subcultured 1/4 into LB with 200 µM IPTG and allowed to grow for 2h in order to induce protein expression. 1 mL samples of cells were then transferred to 2 mL tubes, washed twice in 50 mM Tris-HCl (pH 7.5) and resuspended in this buffer. Samples were then subjected to 5 freeze-thaw cycles in order to lyse cells. After spinning down samples, 1 mL of cell lysate was transferred to a cuvette, and a spectrophotometer was blanked at 340 nm with this sample. 140 µM NADH and 20 µM FMN was then added, the cuvette was quickly inverted, and readings were taken at 340 nm. pLacI-dszB was used as a control to measure native amounts of oxidoreductase activity, whereas the pLacI-hpaC cultures were used to measure activity when HpaC was expressed. The control was just Tris-HCl buffer with the NADH and FMN compounds added. Decrease in absorbance at 340 nm corresponds to the loss of NADH as it is converted to NAD+.

When the assay was run, it was found that NADH does not convert readily to NAD+ on its own. When cell lysate containing the naturally expressed amounts of oxidoreductase was added, a decrease in absorbance could quickly be observed as the NADH was converted to NAD+. When cultures over-expressing HpaC were tested, the absorbance levels were found to start much lower than the control. We believe that this is because with the amount of cell lysate tested, when the HpaC protein is overexpressed the NADH is consumed almost immediately and therefore the data reflecting the drop in absorbance is missed. Further tests will use differing amounts of cell lysate in order to try to capture data that shows the drop in absorbance for HpaC cultures.

4) Optimizing Gene Order

After the genetic elements we intended to use in the pathway were selected, more research was performed into further ways to optimize the system for maximum output and efficiency in order to overcome the issues present in the natural system. Rearranging the order of the genes in the operon has also shown to be beneficial in improving the efficiency of the 4S pathway. The natural order of genes in the operon is dszABC, however the relative activity of the enzymes does not reflect this pattern, as the catalytic activity of DszA:DszB:DszC is 25:1:5 respectively. Because of this, rearranging the operon to dszBCA produces more of the enzyme with lower activity and less of the enzyme with high activity. This is accomplished simply by having more of the less active enzyme transcribed, which translates into more enzyme being present to compensate for its lower activity (Li et al., 2008). In addition, the initiation codon of dszB overlaps with the termination codon of dszA in the natural operon, which leads to translational coupling and further reduces the amount of DszB, leading to a bottleneck in the pathway. By seperately PCR'ing out each of the genes and arranging them in this fashion in a biobrick plasmid, with each separated by an RBS site part (BBa_B0034), this issue is also eliminated.

Figure 8: Method of optimizing gene order. The top circuit represents that found natively in the organism, with the bottom circuit representing our modified version.

Final Constructs

After all of the above, we hope to develop four final constructs for our system which will allow us to test the desulfurization under a few different conditions.

Figure 9: First set of final constructs for the desulfurization operon, with constitutive Dsz expression and inducible expression of the optimization proteins; either HpaC on its own or coexpressed with KatG

The first two constructs have the modified dsz operon (dszB, dszC, dszA) under the control of a constitutive TetR promotor (BBa_J13002) This is to allow for the testing of the optimization circuit, which is under the control of a lacI promotor inducible by IPTG (BBa_J04500). The set-up of these two constructs will therefore allow for the expression of the dsz genes with the ability to test and compare their desulfurization rates
A) On their own
B) With the addition of hpaC
C) With the addition of both hpaC and katG-LAA

This will allow us to determine what the optimal construct and expression levels of the additional genes must be in order to have the most effective sulfur removal system.

Figure 10: Second set of final constructs for the desulfurization operon, with all genes under an IPTG inducible promotor.

Because there are a large number of proteins being expressed in this system, we were concerned about the possibility of inclusion bodies forming which would essentially result in an unusable system, in addition to possible other toxic effects due to protein expression. Because of this a backup system was built, where both the optimization circuit and the dsz operon were under the control of the lacI promoter. This system would allow us to tune the expression of the genes, and determine which expression level is optimal. Once tuning was performed, a constitutive promotor of that strength could then be implemented to maintain optimal gene expression.

Currently, assembly of these final constructs is underway, with only a couple more construction steps before functionality tests can begin.