Difference between revisions of "Team:Uppsala/CrocinPathway"

Summary:We wanted to extend the saffron pathway from FFP to zeaxanthin to crocin, but the pathway was poorly characterizes. We identified and created sequence verified BioBricks out of three enzymes that can perform the three step conversion from zeaxanthin to crocin: CaCCD2, CsADH2946 and UGTCs2. We have also characterized these enzymes with experiments and simulations. Above all, we are the first to purify and confirm activity of CsADH2946 as well as measuring the kinetic parameters of the enzyme (K_M = 20.7842 µM ± 3.5264. In addition, we performed steered molecular dynamics (pulling) with CsADH2946 and the substrate crocetin dialdehyde, which showed that CsADH2946 has a high affinity towards the substrate. Our experimental data and modeling results show that CsADH2946 is a very good enzyme for this crocin pathway reaction.

We are the first to express and characterize CsADH2946 (Crocus Sativus aldehyde dehydrogenase 2946)! This aldehyde dehydrogenase gene from Crocus Sativus has previously only been identified as a candidate gene through proteome analysis, and has thus never been isolated or characterized before (3). We successfully made a sequence verified BioBrick of CsADH2946 with his-tag (BBa_K2423007). The BioBrick was also combined with the other steps in the pathway and inserted into the zeaxanthin producing E. coli strain for a complete pathway from FPP to crocin. See the result here! In summary, our experimental data and modeling results show that CsADH2946 is a very good enzyme for this reaction.

CsADH2946 was transformed and expressed in E. coli strain BL21(DE3*) and purified using IMAC on an ÄKTA protein purification system. We used a gradient of imidazole concentration from 20–500 mM, in order to get our enzyme as separated as possible from other proteins that ends up in the fractions. The peak pointed at by the arrow in the chromatogram (figure 2) indicates protein that elutes at high imidazole concentration, i.e our desired his-tagged CsADH2946. The purification was followed by SDS-PAGE to analyse the fractions, control purity and verify the protein product. In figure 3a the band at around 60 kDa in the crude pellet indicate an overexpression of a protein in that size range. In the SDS gel of fractions 16–26 collected between 115–145 mL elution volume (figure 3b) there is a strong band at 60 kDa corresponding to the molecular weight of CsADH2946, indicating that our protein was successfully overexpressed and well-separated.

As can be seen in figure 4, the absorbance of the product crocetin increases over time in well 2 containing enzyme and the substrate crocetin dialdehyde. After 9 hours of reaction, the blue curve corresponding to the enzyme + substrate mixture has increased its absorbance in the exact range of the product. The negative and positive control curves look similar to time point zero, apart from some precipitation of product and substrate indicated by the decreased curves. A definite evidence that we succeeded to produce a functional CsADH2946 enzyme. Using this data, we could estimate K_M = 20.7842 µM ± 3.5264.

In addition, in figure 5 we can see that well 2 containing enzyme and crocetin dialdehyde has changed color compared to the negative control, to become more yellow like the product crocetin in well 8. This also shows that CsADH2946 was produced and that it converts crocetin dialdehyde into crocetin.

Since the enzyme is poorly characterized, we created a homology model and performed stability simulations to verify that our model was reasonable. The homology modeling revealed that CsADH2946 is in fact tetrameric, which helped us in the purification and characterization process. We performed a pulling simulation between the enzyme and its substrate in order to estimate binding energy and calculate a theoretical K_d (=4.9321 µM). The resulting structure of the homology modeling and a plot of the pulling simulation can be seen in figure 6. Using the results from the activity measurement, the earlier unknown Michaelis-Menten kinetic parameters of the reaction could also be estimated using a Bayesian inference algorithm. With this method we got K_M (=20.7842 µM). Read more about the homology modeling, dynamics modeling and the kinetic parameter estimation in the Modeling section.

We successfully made a sequence verified BioBrick of UGTCs2 (Crocus Sativus UDP-glucuronosyltransferase 2) with his-tag (BBa_K2423008). The BioBrick was also combined with the other steps in the pathway and inserted into the zeaxanthin producing E. coli strain for a complete pathway from FPP to crocin. See the result here!

Since the enzyme is poorly characterized, we created a homology model and performed stability simulations to verify that our model was reasonable. The resulting structure of the homology modeling and its corresponding RMSD plot from simulations can be seen in figure 7. The RMSD plot indicated that the model conforms to a stable structure. Read more about the homology modelling and dynamics modelling in the Modeling section.