Modelling
Software tool that allows comparison of three different growth models based on plate reader data
We wanted a model to extract valuable information about bacterial growth, comparing growth rates and max capacity primarily. The reason for this is that Stirling et al. (2017) indicated that kill switches always have a low level expression of the toxin, even when they are not induced. This impacts growth and means that kill switches are stable for only a finite amount of generations. Modeling these parameters can provide insight into how stable our constructs are, and if they might be outcompeted by other bacteria they would encounter in the gut.
Therefore we built a Matlab code that compared the accuracy of 3 different but widely used growth models, the exponential, logistic and Gompertz model. These models have 1, 2 and 3 free parameters, respectively. Comparing models with different numbers of free parameters is bad practice, so we made the model calculate and compare performance based on the Akaike information criterion.
After fitting it to self-obtained growth data from both Escherichia coli and Pseudomonas alcaligenes we concluded that the Gompertz model performed the best. To further improve on our efforts we looked in literature to use the Gompertz model even better. After additional literature study we found a very recent paper that detailed a parameter estimation method for the Gompertz model that was better than previous methods. We tried to incorporate this into our model, but unfortunately could not get it to work in time for the wiki freeze.
However, we still believe our growth model can contribute to future teams. It has the capacity to analyze a 96-well plate reader experiment with many repeated measurements. Based on the conditions it will group measurements together if they and compare different conditions amongst themselves with ANOVA. It will exclude wells from rows or columns that are marked as 'Empty'.
We also wanted to help other iGEM teams understand how to interpret growth models like the ones we used. So we decided to make an educational tool that allows users to change the growth rate, max capacity or lag time and visually see the effect on simulated graphs for the three models. On the wiki page Dry lab - Growth models we have posted a video where we demonstrate the tool, which we recommend checking out. This educational tool is available as a free resource.
Flux Based Analysis for Pseudomonas species extended with gene product reactions
Flux Based Analysis (FBA) is a type of modeling that maps the metabolic network and fluxes in that network. It is a constraints based model that assumes steady state. For our project we added the reactions needed for L-DOPA synthesis.
Additionally, in our project we ran into an issue that the toxin of our kill switch degrades mRNA. However, the FBA we used did not have gene product reactions (GPRs), which allow different gene expression and gene-metabolic interactions. We could therefore also not perform knock-in screening of genes that would increase or decrease yield. So we decided to add GPRs to the model.
This massively expands the capacity of the model. It can now be used to study the effects of knock-ins, differential gene expression and gene-interactions and networks. This can not only be used for L-DOPA synthesis, but for any other metabolic study of Pseudomonas species. We think this is a valuable contribution to future teams. We will provide the code as a resource.