Team:NUS Singapore-A/Modelling/Model2
EL222 is a protein that dimerizes in blue light and induces or represses a genetic circuit depending on the type of promoter used. In this study, we model an inducible and a repressible system to understand the kinetics behind the mechanism. Next, we analyze the information to understand trends observed and possible weaknesses of the study. Finally, we discuss how the information from this study impacts future work for our project.
Contents
- 1 Repressible System Modelling
- 2 Inducible System Modelling
- 3 Weakness and Future Improvement
- 3.1 All RFP matures before degradation
- 3.2 Initial mRNA concentration was assumed to be 0
- 3.3 Initial concentration of nascent RFP was assumed to be 0
- 3.4 All nascent RFP was assumed to be converted into mature RFP without degradation.
- 3.5 Assumed no reversal in the effect of EL222 concentration on the inducible promoter.
- 3.6 Changing conditions of media during the experiment
- 4 Achieving on-off Cycles
Repressible System Modelling
The first part of this study was the modeling of the repressible system by performing curve fitting on experimental data on RFP expression over time.
Experiments were performed for the following scenarios:
- 1. 8hr off
- 2. 3hr off/3hr on
- 3. 8hr on
- 4. 2hr off/4hr on
- 5. 45min off/6hr on
Curve fitting was done on scenarios 1 to 3. The resultant model was tested by simulating the system response for scenario 4 and 5 and comparing it with experimental data.
The following 10 factors were taken into account when modelling the repressible system.
- ka the rate of EL222 dimerization in blue light.
- kd the rate of degrdation of dimerized EL222.
- synmRFP the max rate of transcription.
- h the Hill coefficient for EL222 dimers binding to the promoter.
- krep the maximum amount of repression possible.
- km the concentration of EL222 dimers at which half of maximum transcription rate is reached
- degRFP the rate of mRNA degradation.
- synRFP the rate of translation.
- Kmat the rate of protein maturation.
- degRFPm the rate of degrdation of mature RFP.
The ODEs are as follows:
Unfortunately, the concentration of intermediates were not measured in light of equipment constraints. Thus, values of ka, kd, km are arbitrary and do not correspond to actual physical values. The rest of the parameters were allowed to vary within an order of magnitude in accordance literature values of similar systems.
Figure 1: Curve fitting for repressible system, light 8hr off
Figure 2: Curve fitting for repressible system, light 3hr off/3hr on
Figure 3: Curve fitting for repressible system, light 8hr on
Running the optimizer in MATLAB gave us the following values for the parameters.
Figure 4: Model testing for repressible system, light 2hr off/4hr on
Figure 5: Model testing for repressible system, light 45min off/6hr on
We found that our model was able to capture the trends in RFP concentration over time for scenario 5 quite well but less so for scenario 4. We hypothesize that more factors were present that were not accounted for and that they probably play a much smaller role when the light is on or off all the way but become important for "intermediate" scenarios. In any case, it appears that the performance of the system can be predicted with reasonable accuracy for lighting regimes close to the scenario 1 and 3. Meanwhile, more complex regimes will require a better model.
Inducible System Modelling
Our team proceeded to investigate the usefulness of our model for an blue light inducible system.
The following 10 factors were taken into account when modelling the inducible system.
- ka the rate of EL222 dimerization in blue light.
- kd the rate of degrdation of dimerized EL222.
- synmRFP the max rate of transcription.
- h the Hill coefficient for EL222 dimers binding to the promoter.
- basal the transcription rate in the absence of inducer.
- km the concentration of EL222 dimers at which half of maximum transcription rate is reached
- degRFP the rate of mRNA degradation.
- synRFP the rate of translation.
- Kmat the rate of protein maturation.
- degRFPm the rate of degrdation of mature RFP.
The ODEs are as follows:
Figure 1: Curve fitting for inducible system, light 8hr on
Figure 2: Curve fitting for inducible system, light 3hr on/4hr off
Figure 3: Curve fitting for inducible system, light 8hr off
Running the optimizer in MATLAB gave us the following values for the parameters.
Figure 4: Model testing for inducible system, light 2hr on/4hr off
Figure 5: Model testing for inducible system, light 45min on/6hr off
Reusing the same values for degmRFP, synRFP, degRFP (Control and DAS), we were able to get the following fit.
Compared to the repressible system, the fit was not as good. In particular, the model could not capture the drop in RFP concentration that appeared at the start of the experiment in many cases. In other cases, the model was unable to capture the increase in RFP concentration correctly. Unfortuantely, even when degmRFP, synRFP, degRFP were allowed to vary, we were still unable to improve the fit significantly.
Weakness and Future Improvement
We have identified the following sources of error and weaknesses in our model:
All RFP matures before degradation
Given the long maturation time for RFP, it is likely that degradation of nascent RFP could have played a significant role.
Initial mRNA concentration was assumed to be 0
This is unlikely to be true given the presence of RFP even at the start of the experiment. Future work measuring mRNA concentration will greatly reduce unvertainty during modelling.
Initial concentration of nascent RFP was assumed to be 0
This assumption is valid only if RFP matures quickly enough for the nascent form to undectable. However, our model indicated a long maturation time which implies that nascent RFP could have been present in siginificant amounts at the start.
All nascent RFP was assumed to be converted into mature RFP without degradation.
It is likely that whatever mechanisms causing degradation of mature of RFP can act on nascent RFP too even if not as efficiently. Furthermore, the long maturation time of RFP means it is possible that degradation of nascent RFP may be significant.
Assumed no reversal in the effect of EL222 concentration on the inducible promoter.
High EL222 concentration can have an inhibitory effect on the inducible promoter. The wet lab team can use a weaker promoter in the future to avoid EL222 reaching inhibitory concentrations.
Changing conditions of media during the experiment
Over the course of the experiment, nutrients are depleted from the media while metabolic waste is produced and accumulates in it. Such inconsistencies in growing conditions can affect the performance of the system. The wet-lab team attempted to minimize such effects by keeping the experiment under 8 hours. Future work starting with a lower OD, using a continuous setup or a cell-free environment can help us study the system in isolation.
Achieving on-off Cycles
We found that the slow maturation of the nascent RFP protein acted as a buffer against changes in mature RFP concentration. Using a faster maturing RFP would however significantly improve the results. This implies that the system works better on proteins that do not require additional time for steps after translation. This includes additional folding, post-translational modification and transport.