Team:Manchester/Enzymetest

Results

Figure 3: Short excerpt of metabolite concentration of fatty acid biosynthesis from multiple simulations within 100 seconds at ten different time points. Colours visualise concentrations according to their amount (Dark Red: >4; Pink: > 2; Light red:> 1; White: 1-0.01; Light yellow: <0.01; Dark yellow: <0.0001)

The concentrations of the metabolites was outputted in tables, as depicted in Figure 3. Each line represents one simulation with ten different time points within 100 seconds. The whole data set of all simulations was then attributed with colours according concentration values. Another table was generated out of this chart according to the ordinal data obtained from colouring the metabolite concentrations. This was done to further improve ease of work and making the data more visual. Figure 4 shows the summary of this qualitative concentration distribution for each metabolite. Again, the brighter green a cell is in colour, the more often simulations rendered metabolite concentrations in the specified concentration interval. For example, the last metabolite in the table C18CoA is bright green, because all 41 simulations rendered between 0.01 - 1 mM. Out of this table, a clear distribution becomes obvious: Except for the first six initial replenishing reactions, all metabolite concentrations are within a small reasonable range mostly between 0.01 - 1 mM. Interestingly, in the reaction towards the end of the pathway, which are responsible for removing the metabolites from the system and therefore give rise to stearic and palmitic acid (our desired products) the range of results appears to be significantly narrower, despite the uncertainty.

The analysis of the data shows clearly, that due to a small and reasonable range of metabolite concentrations which stabilises towards the end of the model, a high validity of our functioning model can be safely assumed and demonstrates that the uncertainty is not globally deleterious. Even though the model was working with high uncertainties in data, the output is always within a valid range.

Figure 4: Comprehensive summary table of all analysed models with colour coded visualisation according to qualitative concentration distribution for each metabolite. The brighter in colour a cell is, the more often the simulations resulted in the specified concentration interval.

Upon analysing the degree of certainty in our model, and finding that it was at a level that we believe is suitable for further analysis, we were able to create a series of boxplots showing the range of values found within our simulations for species accumulation after 100 seconds. We focused on the longer chain fatty acids, which are the engineering target of our pathway. The order in which the species are shown in the box plots, Figure 5, is also the order in which they are formed. This is also shown in Figure 6, where the colour corresponds to the colour of the bar on the box plot.

Figure 5: Summarized boxplots of the last metabolite concentrations (mM) in the modelled fatty acid synthesis pathway. Indicated colours correspond to the colours and metabolites as shown in the excerpt of the KEGG pathway in Figure 6.

Figure 6: Excerpt of the simplified version of fatty acid elongation pathway. Coloured boxes and their metabolites correspond to the metabolites as indicated in Figure 5.

These results further emphasise that although we created a model based on uncertain parameters, by embracing this uncertainty we have been able to make a model that gives us useful information – and that allows us to specify for every single prediction how certain we can be of getting it right, in particular towards the end of the pathway.

Similar data analysis was carried out on the rates of the reactions, shown in Figure 7. We focused on the reactions we had labelled AAT at the end of our pathway. These are thioesterase reactions directly responsible for the formation of palmitic and stearic acid. We can see that the rates for these reactions also fall within a relatively small range.

Figure 7: Reaction rates of key reactions in the fatty acid synthesis pathway.

Conclusion

Kinetic Pathway modelling demands abundant information of the kinetic parameters. Literature research, however, showed that these were not available sufficiently or involved measurement errors. Hence this knowledge of parameter values often is uncertain. Therefore, we had to choose an approach that is able to deal with these limitations. Uncertainty modelling proved to be the most promising and useful tool for this. Even though the available data was limited, we managed to create a functioning kinetic model of the fatty acid synthesis pathway. This has not been done before and would not have been possible with any traditional approach.

A prime example of how our metabolic modelling work directly informed our experimental work is in our decision to biobrick the FabA gene (encoding β-hydroxydecanoyl-ACP dehydrase, shown by the DH_OH reactions in this model). Our uncertainty model had shown us that we would need more kinetic data on key enzymes. The least characterised reaction was catalyzed by the product of the fabA gene, therefore we wished to not only biobrick this gene, but a His-tag to purify the enzyme in order to experimental gauge its activity.

However, having taken pains to ensure our model was as realistic as possible, the idea of the insertion of a his-tag that could affect the activity of the enzyme seemed at odds to our overall goal. Therefore, we used further modelling technique to ensure the addition of this his tag would have as little overall bearing on the activity of the enzyme as possible. This can be found here

Future Applications: Potentials and Limitations

We believe that this approach to modelling could have a big impact in terms of how Synthetic Biology is modelled in the future and demonstrates a method in which, by facing the uncertainty of modelling head-on and incorporating this into our approach in a principled manner, it is possible to produce valuable models. This is particularly important in the field of Synthetic Biology, where systems, even if well characterised in one organism, are unlikely to have the same parameters when expressed in another organism.

This approach gives us the ability to model complex and poorly experimentally measured systems, where previous attempts may have produced unrepresentative models. Since the Km values can be sampled from a distribution, the model can be used to determine outcomes that may not be obvious with the use of a single Km value.

However, it is important to note that this method of modelling may not be appropriate in every case. The largest limitation of our use of this method is the inability of some of our simulations to reach steady state. This is likely to be a result of the random combination of parameter values. As the models were not fine-tuned, they will not always work. Although, we consider this as a potential strength as we can clearly highlight possible break points in the system that require further analysis. We show this in our own studies of β-hydroxydecanoyl-ACP dehydrase, described above.

Synthetic Biology operates at the cutting edge of current knowledge. Therefore, it will unavoidably face the challenge of uncertainty. Building models with incorporated acknowledgment of uncertainty will yield model predictions with specified confidence intervals, and thus will lead to more robust design strategies for a vast range of engineered cellular machines.

Appendices

Here you can download all of the spreadsheets used in the creation of this model:

The spreadsheets generated from our script can be found here:
RATES
SPECIES

Nomenclature of main metabolites