Team:Sheffield/Chimeric protein system
From 2010.igem.org
(Difference between revisions)
(New page: {{Template:Team:Sheffield/Navigation2}} <div style="text-align: left;"> == '''Model of the chimeric protein system''' == '''Explanation of the system''' The diagram below explains the p...) |
|||
Line 77: | Line 77: | ||
d1 = transcription degradation | d1 = transcription degradation | ||
+ | After the reporter DNA genes have been transcribed to mRNAs, the next step in the gene expression is the mRNAs being translated to proteins. In the presence of tRNAs, the mRNAs get translated into proteins and therefore to construct a mathematical model for translation for our chimeric system, the following assumptions were also made; | ||
+ | |||
+ | 1. The reporter DNA gene expression is limited at the beginning of the translation process. [3] | ||
+ | |||
+ | 2. Also, tRNAs are abundant in the E. Coli and that there are over – expressing rare tRNAs and so all the mRNAs are translated to proteins. [4] | ||
+ | |||
+ | 3. The degradation rate of the protein is very negligible and therefore the protein produced have longer half – life. | ||
+ | |||
+ | With these assumptions we create following differential equation for the translation. | ||
+ | |||
+ | [[Image:Equation4.JPG]] | ||
+ | |||
+ | where : | ||
+ | |||
+ | t2 = the rate of translation | ||
+ | |||
+ | d2 = translation degradation | ||
+ | |||
+ | |||
+ | '''Simulation of the model''' | ||
+ | |||
+ | We first analysed the behaviour of the different layers of proteins in the phosphorylation pathway by simulating the model for the signal transduction using ode solver in matlab to solve the Equations (6) – (9). | ||
+ | |||
+ | For the simulations we set the total amount of each protein layer to 10 and all the rate constants to 1 and the steady state responses of each of the layers were investigated for different levels of input. Graph 1 below shows the results from this simulation. | ||
+ | |||
+ | [[Image:Graph1.JPG]] | ||
+ | |||
+ | Graph 1 | ||
+ | |||
+ | From this graph it could be seen that the amount of phosphorylated proteins in each layer increases as the input increase and they saturate after a certain input level with the 4th layer being the fastest. The saturation rate of each protein depends on the saturation of the succeeding layers and these levels could be changed by altering the rate constants for each of the reactions [2]. | ||
+ | |||
+ | In the complete system these would not always tend to a steady state as the response regulator at the final layer binds to the reporter DNA to facilitate transcription of the reporter gene. | ||
+ | |||
+ | In simulating the mathematical model for the complete system consisting of the six differential equations (Equations (6) – (11)) we again input the equations into matlab and used an ode solver to solve the equations. (The code can be found in the Appendix) | ||
+ | |||
+ | Again we set the total amount of each protein layer to 10 and all the rate constants to 1 except for transcription/translation degradation which are at 0.01. | ||
+ | |||
+ | We first investigated the expression of the reporter gene over 1000 seconds to determine the production of the colour. | ||
+ | |||
+ | [[Image:Graph2.JPG]] | ||
+ | |||
+ | Graph 2 | ||
+ | |||
+ | It can be seen that as in Graph 1 the colour increases and tend to a steady state at around 900 seconds. From this we could estimate that it will take about 15 minutes for the colour to be produced. The input was kept constant at 1 over this period. | ||
+ | |||
+ | The curve shows that the initial rate of translation was low for sometime before becoming high as the translation curves were both sigmoidal. From this, it therefore could be concluded that fluorescence light produced from the reporter protein in the system may not be visible to the naked eye initially when the E.coli detects the CAI-1. | ||
+ | |||
+ | As the real system parameters will not be 1, we also created a stochastic model where the system was simulated with parameters selected at random. The result for this simulation is given in Graph 3 below. | ||
+ | |||
+ | [[Image:Graph3.JPG]] | ||
+ | |||
+ | |||
+ | Graph 3 | ||
+ | |||
+ | This graph shows that the system still follows the same shape as before with only a little amount of noise. Therefore it can be deduced that the chimeric protein system is highly robust for variations in system parameters. | ||
+ | |||
+ | We then investigated how the steady state response for the reporter protein would change with different levels of inputs. This is shown by Graph 2 below. | ||
+ | |||
+ | [[Image:Graph4.JPG]] | ||
+ | |||
+ | |||
+ | Graph 4 | ||
+ | |||
+ | It can be seen that from 0 to 1.2 moles of CAI-1 (estimated values with systems parameters set to 1) the steady state levels increase and at these levels there would be a difference in colour intensity seen. After 1.2 the concentration of the protein tends to a steady state and any more increase in the input would not change the intensity of the colour. Therefore we can deduce that we can use the biosensor to measure the amount of cholera present only up to a certain level of autoinducer. | ||
+ | |||
+ | |||
+ | Different parameter and input values could change the output values of the system this can be analysed using the interface we have created in matlab. More information can be found at [[Team:Sheffield/Graphical User Interface|Graphical User Interface]]. | ||
+ | |||
+ | '''References''' | ||
+ | [1] Kim JR, Cho KH ,2006 ,The multi-step phosphorelay mechanism of unorthodox two-component systems in E-coli realizes ultrasensitivity to stimuli while maintaining robustness to noises, Computational Biology and Chemistry , [online] Available at: <http://www.sciencedirect.com/science?_ob=ArticleURL&_udi=B73G2-4PNJ1K7-7&_user=10&_coverDate=12%2F31%2F2006&_rdoc=1&_fmt=high&_orig=search&_origin=search&_sort=d&_docanchor=&view=c&_acct=C000050221&_version=1&_urlVersion=0&_userid=10&md5=4ed1aecda98446af7276385b09dff8b2&searchtype=a> [Accessed 24 October 2010] pp. 438 - 439 | ||
+ | |||
+ | [2] Csikasz-Nagy A. ,Cardelli, L. Soyer OS.(2010) ,Response Dynamics of phosphorelays suggest their potential utility in cell signalling, [online] Available at: http://rsif.royalsocietypublishing.org/content/early/2010/08/06/rsif.2010.0336.short?rss=1 [Accessed 24 October 2010] pp. 1 – 3 ,7-8 | ||
+ | [3]Kudla G, Murray AW, Tollervery D Plotkin JB (2009); “Coding sequence determinants of gene expression in Escherichia coli” | ||
+ | |||
+ | [4] Kane JF (1995); Effects of rare codon clusters on high - levels expression of the heterologous proteins in Escherichia coli. Curr Opin Biotechnol 6: 494 - 500 | ||
Revision as of 00:15, 25 October 2010