Team:SDU-Denmark/project-bc
From 2010.igem.org
(Difference between revisions)
(→The real system) |
(→Model results) |
||
(29 intermediate revisions not shown) | |||
Line 13: | Line 13: | ||
=== The real system === | === The real system === | ||
- | Several chemosensory systems exist in bacteria. The one most studied is the chemotaxis of ''E. coli''. Polar clusters of membrane-spanning methyl accepting chemotaxis proteins (MCPs) are positioned in the ends of the rod shaped ''E. coli''[[https://2010.igem.org/Team:SDU-Denmark/project-bc#References 1]].<br> | + | Several chemosensory systems exist in bacteria. The one most studied is the chemotaxis of ''E. coli''. Polar clusters of membrane-spanning methyl accepting chemotaxis proteins (MCPs) are positioned in the ends of the rod shaped ''E. coli'' [[https://2010.igem.org/Team:SDU-Denmark/project-bc#References 1]].<br> |
- | Tar (responding to aspartate), Tsr (responding to serine), Tap (responding todipetides), Trg (responding to galactose) ang Aer (responding to oxygen) are the most common of these receptors seen in ''E. coli'' strains[[https://2010.igem.org/Team:SDU-Denmark/project-bc#References | + | Tar (responding to aspartate), Tsr (responding to serine), Tap (responding todipetides), Trg (responding to galactose) ang Aer (responding to oxygen) are the most common of these receptors seen in ''E. coli'' strains[[https://2010.igem.org/Team:SDU-Denmark/project-bc#References 2]]. Transduction of the signal between the MCPs and the flagella motor is managed by a particular fine-tuned signalling cascade, which is governed by several intracellular proteins.<br> |
- | [[Image: Team-SDU-Denmark-Chemotaxis.png |thumb|center|550px| '''Figure 1:''' Schematics of the bacterial chemotaxis pathways [[https://2010.igem.org/Team:SDU-Denmark/project-bc#References | + | [[Image: Team-SDU-Denmark-Chemotaxis.png |thumb|center|550px| '''Figure 1:''' Schematics of the bacterial chemotaxis pathways [[https://2010.igem.org/Team:SDU-Denmark/project-bc#References 2]]]] |
Chemotaxis consists of three important phases: reaction, adaptation and relaxation. In the reaction phase, one of the MCPs sense the specific chemical gradient it responds to and the membrane protein activates the CheW enzyme. The active CheW enzyme suppresses the auto-phosphorylation of the CheA enzyme and when no or very little CheA enzyme is phosphorylated, no phosphorylation of the CheY and CheB enzymes occurs. In the end, this increases run, because phosphorylated CheY binds to the flagella motors and induce tumbling. When the bacteria have reached the area with a higher concentration of the chemical in question adaptation sets in.<br><br> | Chemotaxis consists of three important phases: reaction, adaptation and relaxation. In the reaction phase, one of the MCPs sense the specific chemical gradient it responds to and the membrane protein activates the CheW enzyme. The active CheW enzyme suppresses the auto-phosphorylation of the CheA enzyme and when no or very little CheA enzyme is phosphorylated, no phosphorylation of the CheY and CheB enzymes occurs. In the end, this increases run, because phosphorylated CheY binds to the flagella motors and induce tumbling. When the bacteria have reached the area with a higher concentration of the chemical in question adaptation sets in.<br><br> | ||
Because no CheB have been phosphorylated during the reaction period CheR is “allowed” to methylate the MCPs, which increases the auto-phosphorylation of the CheA enzyme. This increase in auto-phosphorylation of CheA results in phosphorylation of CheY that leads to tumbling of the bacteria. In the last phase, relaxation, CheB has again been phosphorylated by CheA and is now able to demethylate the MCPs, returning the bacteria to its normal run/tumble frequency with a 0,1 second tumbling period every 1 second.<br><br> | Because no CheB have been phosphorylated during the reaction period CheR is “allowed” to methylate the MCPs, which increases the auto-phosphorylation of the CheA enzyme. This increase in auto-phosphorylation of CheA results in phosphorylation of CheY that leads to tumbling of the bacteria. In the last phase, relaxation, CheB has again been phosphorylated by CheA and is now able to demethylate the MCPs, returning the bacteria to its normal run/tumble frequency with a 0,1 second tumbling period every 1 second.<br><br> | ||
- | One other thing to consider in the model is the distance from the CheY phosphorylation site to the flagella motors that are randomly located on the entire cell body[[https://2010.igem.org/Team:SDU-Denmark/project-bc#References 1]]. Several modelling methods have been proposed and nearly all of them include these considerations[[https://2010.igem.org/Team:SDU-Denmark/project-bc#References 1]]:<br> | + | One other thing to consider in the model is the distance from the CheY phosphorylation site to the flagella motors that are randomly located on the entire cell body [[https://2010.igem.org/Team:SDU-Denmark/project-bc#References 1]]. Several modelling methods have been proposed and nearly all of them include these considerations [[https://2010.igem.org/Team:SDU-Denmark/project-bc#References 1]]:<br> |
"Attractant is assumed to be in excess and binding is assumed to be rapid". One of the most common receptor types to be modelled is the Tar receptor. This is not unexpected given that a large degree of the experimental literature has focused on the response of ''E. coli'' to aspartate which is detected by this receptor. Tar is also one of the most abundant receptors within the cytoplasmic membrane.<br> | "Attractant is assumed to be in excess and binding is assumed to be rapid". One of the most common receptor types to be modelled is the Tar receptor. This is not unexpected given that a large degree of the experimental literature has focused on the response of ''E. coli'' to aspartate which is detected by this receptor. Tar is also one of the most abundant receptors within the cytoplasmic membrane.<br> | ||
In models which have focused on the adaptation and/or the phosphorylation cascade, receptors are commonly modelled as complexes consisting of MCP, CheW and CheA. This is justified by the tight association of CheW and CheA to the MCP receptors.<br> | In models which have focused on the adaptation and/or the phosphorylation cascade, receptors are commonly modelled as complexes consisting of MCP, CheW and CheA. This is justified by the tight association of CheW and CheA to the MCP receptors.<br> | ||
Line 25: | Line 25: | ||
=== Description of the model === | === Description of the model === | ||
The model we will be using to predict biochemical behaviour of the system is a simplification of the real system.<br> | The model we will be using to predict biochemical behaviour of the system is a simplification of the real system.<br> | ||
- | The simplified model does not consider the normal metabolism in the cell. It only includes the most important enzymes in the normal ''E. coli'' chemotaxis pathway and the proteins responsible for coupling | + | The simplified model does not consider the normal metabolism in the cell. It only includes the most important enzymes in the normal ''E. coli'' chemotaxis pathway and the proteins responsible for coupling phototaxis to chemotaxis.<br> |
- | Below, you will | + | Below, you will find a schematic representation of the system our biochemical modelling was performed on. The crossed-over enzymes and proteins have not been modelled as a differential equation but considered as constants in the system.<br> |
- | [[Image: Team-SDU-Denmark-biomodelling.png |thumb|550px |center| '''Figure 2''': Schematics of the system we want to make. The enzymes that are crossed-over are considered constant in our system]] | + | [[Image: Team-SDU-Denmark-biomodelling.png |thumb|550px |center| '''Figure 2''': Schematics of the system we want to make. The enzymes that are crossed-over are considered constant in our system.]] |
The important enzymes in the system are mathematically represented by a differential equation under the assumption of steady state in the system and modelled in a demo version of the program [http://www.berkeleymadonna.com/ Berkeley Madonna]. | The important enzymes in the system are mathematically represented by a differential equation under the assumption of steady state in the system and modelled in a demo version of the program [http://www.berkeleymadonna.com/ Berkeley Madonna]. | ||
<br> | <br> | ||
- | [[Image: Team-SDU-Denmark-Madonnapic.png |thumb|center|550px|'''Figure | + | [[Image: Team-SDU-Denmark-Madonnapic.png |thumb|center|550px|'''Figure 3''': Picture from the modelling program Berkeley Madonna, showing a model example.]] |
<br> | <br> | ||
- | When light hits the photo receptor | + | When light hits the photo receptor SopII, the signal is transduced through the HtrII and Tar and onto the enzyme CheW. CheW is then activated because it is a straight-forward signal transduction. The proteins SopII and HtrII are not considered in the model and the equation for the activation of CheW is: |
<br> | <br> | ||
[[Image:Team-SDU-Denmark-Chewa.png|center]] | [[Image:Team-SDU-Denmark-Chewa.png|center]] | ||
Line 44: | Line 44: | ||
[[Image:Team-SDU-Denmark-CheYp.png|center]] | [[Image:Team-SDU-Denmark-CheYp.png|center]] | ||
<br> | <br> | ||
- | Phosphorylated CheB is dephosphorylated when methylated | + | Phosphorylated CheB is dephosphorylated when methylated Tar is demethylated, this can be described mathematically by this equation: |
<br> | <br> | ||
[[Image:Team-SDU-Denmark-CheBp.png|center]] | [[Image:Team-SDU-Denmark-CheBp.png|center]] | ||
<br> | <br> | ||
- | Methylated | + | Methylated Tar is constantly produced by the CheR enzyme which, in this model, is also considered a constant. The demethylation is performed by CheB-p. The concentration of methylated Tar can be described by this equation: |
<br> | <br> | ||
[[Image:Team-SDU-Denmark-TSRm.png|center]] | [[Image:Team-SDU-Denmark-TSRm.png|center]] | ||
<br> | <br> | ||
- | In Berkeley Madonna it is possible to define different parameters as sliders which enable the user to change parameters continually while the model is running. The sliders used in the model are: | + | In Berkeley Madonna it is possible to define different parameters as sliders which enable the user to change parameters continually while the model is running. The sliders used in the model are:<br> |
+ | [[Image: Team-SDU-Denmark-Sliders.png|thumb|center|550px|'''Figure 4''': Picture from the modelling program Berkeley Madonna, showing the sliders used in the model.]] | ||
<br><br> | <br><br> | ||
=== Model results === | === Model results === | ||
- | With the Madonna program it | + | With the Madonna program it is possible to model the system parameters, the most important of them are:<br> |
TAR/TARm<br> | TAR/TARm<br> | ||
CheA/CheAp<br> | CheA/CheAp<br> | ||
Line 62: | Line 63: | ||
CheY/CheYp<br> | CheY/CheYp<br> | ||
Flagella moters CCW/CW<br> | Flagella moters CCW/CW<br> | ||
- | The time scale of the model might be to large, because it’s not possible to define units to the | + | The time scale of the model might be to large, because it’s not possible to define units to the different parameters.<br> |
- | Focusing on the behaviour of TAR/TARm in the model, it demonstrates that the concentration of methylated TAR rises if no phosphorylated CheB is present in the system.<br> When the concentration of phosphorylated CheB rises the concentration of methylated TAR drops.<br> | + | In the program defining time steps as seconds, minutes and so on, is not possible this problem also involves defining concentrations of the enzymes.<br> |
+ | Some of the initial reaction rate constants and enzyme concentrations was from articles, and afterwards optimized for this model [[https://2010.igem.org/Team:SDU-Denmark/project-bc#References 3]].<br> | ||
+ | Focusing on the behaviour of TAR/TARm in the model (Figure 5), it demonstrates that the concentration of methylated TAR rises if no phosphorylated CheB is present in the system.<br> When the concentration of phosphorylated CheB rises the concentration of methylated TAR drops.<br> | ||
This corresponded well with the known theory for the real system. | This corresponded well with the known theory for the real system. | ||
- | |||
- | |||
<br> | <br> | ||
- | Focusing on the CheA/CheAp behaviour in the model, it demonstrates that CheA auto-phophorylation rate rises when no TAR | + | [[Image: Team-SDU-Denmark-Madonnapic.pngTSRmodel.png |thumb|center|550px|'''Figure 5''': graph showing the concentration of phosphorylated CheB affecting the concentration of methylated TAR, making them drop.]] |
- | The dephosphorylation of CheA is performed by both CheB and CheY | + | <br> |
+ | Focusing on the CheA/CheAp behaviour in the model (The top graph in figure 6), it demonstrates that the CheA auto-phophorylation rate rises when no TAR is methylated, thus producing phosphorylated CheA.<br> | ||
+ | The dephosphorylation of CheA, is performed by both CheB and CheY (the bottom graph in figure 6), rise as the concentrations of both phosphorylated CheB and CheY rise and thereby decreases the concentration of phosphorylated CheA.<br> | ||
This corresponded well with the known theory for the real system. | This corresponded well with the known theory for the real system. | ||
<br> | <br> | ||
- | [[Image:Team-SDU-Denmark-CheAmodel.png |thumb|center|550px|'''Figure | + | [[Image:Team-SDU-Denmark-CheAmodel.png |thumb|center|550px|'''Figure 6''': '''The top graph''' is showing the concentration of phosphorylated CheA rising when the concentration of methylated TSR drops. '''The bottom graph''' is showing the concentrations of CheB and CheY rising when the concentration of phosphorylated CheA rises.]] |
<br> | <br> | ||
- | Focusing on the CheB/CheBp behaviour in the model which illustrate that concentrations of phosphorylated CheB rises when phosphorylated CheA is present, as seen in figure | + | Focusing on the CheB/CheBp behaviour in the model which illustrate that concentrations of phosphorylated CheB rises when phosphorylated CheA is present, as seen in figure 6 and its dephosphorylation occurs during the demethylation of TAR as seen in figure 5.<br> |
- | + | ||
This corresponded well with the known theory for the real system.<br> | This corresponded well with the known theory for the real system.<br> | ||
- | Focusing on the CheY/CheYp and the flagella motor rotational direction the model illustrates that when concentrations of phosphorylated CheA rises, the concentrations of phosphorylated CheY rises and the flagellar activates clockwise rotation.<br> | + | Focusing on the CheY/CheYp (Figure 6) and the flagella motor rotational direction the model illustrates that when concentrations of phosphorylated CheA rises, the concentrations of phosphorylated CheY rises and the flagellar activates clockwise rotation.<br> |
The concentration of phosphylated CheY is proportional to the flagella motor rotational direction, because phosphylated CheY activates Clockwise rotation in the flagellar motors which without phosphylated CheY turns Counter Clockwise.<br> | The concentration of phosphylated CheY is proportional to the flagella motor rotational direction, because phosphylated CheY activates Clockwise rotation in the flagellar motors which without phosphylated CheY turns Counter Clockwise.<br> | ||
The CheY er constantly dephosphorylated by CheZ, which in this case is modelled as the constant k.<br> | The CheY er constantly dephosphorylated by CheZ, which in this case is modelled as the constant k.<br> | ||
This corresponded well with the known theory for the real system. | This corresponded well with the known theory for the real system. | ||
<br> | <br> | ||
- | [[Image: Team-SDU-Denmark-CheYmodel.png |thumb|center|550px|'''Figure | + | [[Image: Team-SDU-Denmark-CheYmodel.png |thumb|center|550px|'''Figure 7''': graph showing the concentration of CheY rising when the concentration of phosphorylated CheA rises.]] |
<br> | <br> | ||
- | This is a simplified model and | + | This is a simplified model and only describes a small part of the real system, this makes the data produced by the model incorrect prediction for actual bacterial behaviour.<br> |
+ | Also the too large time scale is making it impossible too describe what happens when the photosensor in the cell is activated by light. <br> | ||
+ | When trying too introduce the photosensor and the effect of light in the model, the increase in light doesn’t change the system behaviour to a degree that affects the output of the program.<br> | ||
+ | This might be because the concentration of the enzymes, the rate constants in the model and the other model parameters have been set to high or to low.<br> | ||
+ | Although several problems exist in the model it is still good enough to model the system behaviour and enzyme interaction but on a larger timescale than seen in nature.<br> | ||
=== References === | === References === | ||
- | [1] M.J. Tindall, S.L. Porter, P.K. Maini, G. Gaglia, J.P. Armitage [http://www.springerlink.com/content/f602r72767124602/ Overview of Mathematical Approaches Used to Model Bacterial Chemotaxis I: The Single Cell] Bulletin of Mathematical Biology (2008) 70: 1525–1569<br> | + | [1] Suzanne R. Lybarger, Janine R. Maddock [http://jb.asm.org/cgi/reprint/183/11/3261 Polarity in Action: Asymmetric Protein Localization in Bacteria] journal of bacteriology (2001) 3261-3267<br> |
- | [ | + | [2] M.J. Tindall, S.L. Porter, P.K. Maini, G. Gaglia, J.P. Armitage [http://www.springerlink.com/content/f602r72767124602/ Overview of Mathematical Approaches Used to Model Bacterial Chemotaxis I: The Single Cell] Bulletin of Mathematical Biology (2008) 70: 1525–1569<br> |
+ | [3] Nikita Vladimirov, Linda Løvdok, Dirk Lebiedz, Victor Sourjik [http://www.ploscompbiol.org/article/info%3Adoi%2F10.1371%2Fjournal.pcbi.1000242 Dependence of Bacterial Chemotaxis on Gradient Shape and Adaptation Rate] Computational Biology (2008) 4(12): e1000242 | ||
</div> | </div> |
Latest revision as of 23:16, 27 October 2010