|
|
| (2 intermediate revisions not shown) |
| Line 1: |
Line 1: |
| - | {{:Team:TU-Eindhoven/Templates/header}}
| + | #REDIRECT [[Team:TU-Eindhoven/Modeling]] |
| - | <!-- {{:Team:TU-Eindhoven/Templates/boxes}} -->
| + | |
| - | <!-- {{:Team:TU-Eindhoven/Templates/title|colour=#386abc|title=Modeling}} -->
| + | |
| - | <!--
| + | |
| - | <html>
| + | |
| - | <img src="/wiki/images/b/b1/EindhovenBoxes.jpg" style="margin-left:4px; margin-bottom:5px; width:710px; margin-top:-2px;">
| + | |
| - | <div style="height:30px; position:relative; font-size:12px; color:#6e6e6e; line-height:15px;">
| + | |
| - | <div style="position:absolute; width:230px; left:5px; text-align:center;">Description of box 1.</div>
| + | |
| - | <div style="position:absolute; width:230px; left:240px; text-align:center;">Description of box 2.</div>
| + | |
| - | <div style="position:absolute; width:220px; left:485px; text-align:center;">Description of box 3.</div>
| + | |
| - | </div>
| + | |
| - | </html>
| + | |
| - | -->
| + | |
| - | {{:Team:TU-Eindhoven/Templates/head|image=https://static.igem.org/mediawiki/2012/7/71/Yeastmodel.jpg}}
| + | |
| - | | + | |
| - | <h3>
| + | |
| - | Global introduction
| + | |
| - | </h3>
| + | |
| - | | + | |
| - | <p>
| + | |
| - | Biological cells use <span class = "lightblue">highly regulated homeostasis systems</span> to keep a very low cytosolic Ca<sup>2+</sup> level. In normal-growing yeast the cytosolic Ca<sup>2+</sup> concentration is maintained in the range of 50-200 nM in the presence of environmental Ca<sup>2+</sup> concentrations ranging from µM to 100 mM<html><a href="#ref_miseta" name="text_miseta"><sup>[1]</sup></a></html>.
| + | |
| - | </p>
| + | |
| - | <p>
| + | |
| - | The homeostatis system in yeast cells has <span class = "lightblue">two basic characteristics</span>. The cytosolic Ca<sup>2+</sup> concentration is tightly controlled by zero steady-state error to extracellular stimuli and the system is relatively insensitive to specific kinetic parameters, due to robustness of the adaptive response<html><a href="#ref_kitano" name="text_kitano"><sup>[2]</sup></a></html>.
| + | |
| - | </p>
| + | |
| - | | + | |
| - | <p>
| + | |
| - | To achieve an accurate model, the influences of voltage-dependent calcium channels are added to a <span class = "lightblue">basic model for yeast calcium homeostasis</span>. In this model, first described by J. Cui <i>et al</i>, the main contributions of calcium transport are defined<html><a href="#ref_cui" name="text_cui"><sup>[3]</sup></a></html>.
| + | |
| - | </p>
| + | |
| - | | + | |
| - | [[File:model1.jpg|320px|left|link=]]
| + | |
| - | <br>
| + | |
| - | <h3>
| + | |
| - | Basic model
| + | |
| - | </h3>
| + | |
| - | <p>
| + | |
| - | Under normal conditions, extracellular Ca<sup>2+</sup> enters the cytosol through an unknown <span class = "lightblue">Transporter X</span>, whose encoded gene has not been identified yet. Cytosolic Ca<sup>2+</sup> can be pumped into the endoplasmic reticulum (ER) and Golgi system through <span class = "lightblue">Pmr1</span> and can be sequestered into the vacuole through <span class = "lightblue">Pmc1</span> and <span class = "lightblue">Vcx1</span>.
| + | |
| - | </p>
| + | |
| - | <p>
| + | |
| - | Both the expression and function of Pmc1, Pmr1 and Vcx1 are regulated by calcineurin, a highly conserved protein phosphatase that is activated by Ca<sup>2+</sup>-bound calmodulin<html><a href="#ref_cui" name="text_cui"><sup>[3]</sup></a></html>. Therefore, the <span class = "lightblue">transcription factor Crz1</span> can be dephosporylated by activated calcineurin. A <span class = "lightblue">conformational switch model</span> is used to simulate Crz1 translocation, as described by Okamura<i>et al</i><html> <a href="#ref_switch" name="text_switch"><sup>[4]</sup></a></html>.
| + | |
| - | </p>
| + | |
| - | <br>
| + | |
| - | | + | |
| - | [[File:model2.jpg|320px|right|link=]]
| + | |
| - | <br>
| + | |
| - | <h3>
| + | |
| - | Voltage-dependent calcium channels
| + | |
| - | </h3>
| + | |
| - | <p>
| + | |
| - | In literature, little can be found about modeling calcium channels in <i>Saccaromyces Cerevisiae</i>, most commonly known as budding yeast. Therefore we still used the model of <span class = "lightblue">sympathetic ganglion 'B' type cells of a bullfrog</span> to describe this process, since the type of voltage-dependent calcium channels is the same in both the bullfrog cells and the yeast cells.
| + | |
| - | </p>
| + | |
| - | | + | |
| - | <p>
| + | |
| - | More information on the voltage-dependent calcium channels is presented in [[Team:TU-Eindhoven/LEC/ModellingFirst|the description of the inital model]].
| + | |
| - | </p>
| + | |
| - | <br>
| + | |
| - | <br>
| + | |
| - | <br>
| + | |
| - | <br>
| + | |
| - | <br>
| + | |
| - | | + | |
| - | [[File:model3.jpg|320px|left|link=]]
| + | |
| - | <br>
| + | |
| - | <h3>
| + | |
| - | GECO-proteins
| + | |
| - | </h3>
| + | |
| - | <p>
| + | |
| - | Last, but certainly not least, we added the kinetics of calcium binding by the GECO-proteins, since we use the
| + | |
| - | fluorescent activity of CMV-X-GECO.1 proteins as light-emitting part of the cell, as a result of an increase of the cytosolic Ca<sup>2+</sup>-level.
| + | |
| - | </p>
| + | |
| - | | + | |
| - | <br>
| + | |
| - | <p>
| + | |
| - | The <span class = "lightblue">kinetic characterization</span> of the GECO-proteins:
| + | |
| - | </p>
| + | |
| - | | + | |
| - | <p>
| + | |
| - | [[File:kinchar.png|350px]]
| + | |
| - | </p>
| + | |
| - | <P ALIGN="right">
| + | |
| - | <i> (Click to enlarge) </i>
| + | |
| - | | + | |
| - | </p>
| + | |
| - | <br>
| + | |
| - | <br>
| + | |
| - | <br>
| + | |
| - | <p>
| + | |
| - | In short, we formulated <span class = "lightblue">a new basic calcium model for yeast cells</span>. This model is extended with overexpression of voltage-dependent calcium channels and addition of GECO-kinetics. The model is made in <span class = "lightblue">MATLAB</span> and the code can be found [[Media:Model_iGEM_-_TU_Eindhoven_2012.zip|here]].
| + | |
| - | </p>
| + | |
| - | <p>
| + | |
| - | For a complete explanation of the model described on this page, please view the [[Media:IGEM_TUeindhoven_-_Thesis_PJ_Alkema.pdf|bachelor thesis of Petra Alkema]], one of our team members.
| + | |
| - | </p>
| + | |
| - | <br>
| + | |
| - | | + | |
| - | <h3>
| + | |
| - | In more detail
| + | |
| - | </h3>
| + | |
| - | | + | |
| - | <p>
| + | |
| - | Since we do not want to bother the semi-interested readers with the <span class = "lightblue">complete derivation of the differential equations</span> describing this model, this can be found in the bachelor thesis of one of our team members, Petra Alkema. In order to understand the results, discussion and conclusion in more detail, we really recommend to study [[Media:IGEM_TUeindhoven_-_Calcium_Model_Yeast_Cells.pdf|this chapter]]. In an addition, a <span class = "lightblue">more detailed introduction to the conformational switch model</span>, part of the total model, can be found [[Media:IGEM_TUeindhoven_-_Conformational_Switch_Model.pdf|here]].
| + | |
| - | </p>
| + | |
| - | <br>
| + | |
| - | <p>
| + | |
| - | The main differential equations are:
| + | |
| - | <br>
| + | |
| - | [[File:DE.png|600px|link=]]
| + | |
| - | | + | |
| - | </p>
| + | |
| - | | + | |
| - | | + | |
| - | <h3>
| + | |
| - | Parameters
| + | |
| - | </h3>
| + | |
| - | | + | |
| - | <p>
| + | |
| - | The model parameters for which all results are calculated, unless otherwise stated:
| + | |
| - | </p>
| + | |
| - | <p>
| + | |
| - | [[File:symbols.jpg|700px]]
| + | |
| - | </p>
| + | |
| - | | + | |
| - | | + | |
| - | <h3>
| + | |
| - | Results
| + | |
| - | </h3>
| + | |
| - | | + | |
| - | <p>
| + | |
| - | As we could not yet accomplish a complete sensitivity analysis on the model, we only consider the <span class = "lightblue">basic characteristics</span> of the different parts of the model. Furthermore, the influences of some basic experimental setup values could be stated, i.e. the concentration of extracellular calcium, Ca<sub>{ex}</sub>, and the duration of the pulse. You can click on the figures to enlarge the results. The first two graphs present the basic model in yeast cells, with the expected oscillatory system. The next graphs show the behavior of the complete model, including the voltage pulse. More results and a more detailed description can be found in the thesis.
| + | |
| - | </p>
| + | |
| - | | + | |
| - | <P ALIGN="center">
| + | |
| - | [[File:res3.png|300px]][[File:res4.png|300px]]
| + | |
| - | [[File:res6.png|410px]][[File:res7.png|270px]]
| + | |
| - | </p>
| + | |
| - | | + | |
| - | | + | |
| - | <h3>
| + | |
| - | Discussion
| + | |
| - | </h3>
| + | |
| - | | + | |
| - | <p>
| + | |
| - | One of the most important points of discussion is the overall correctness of the calcium model. Due to the fact that the characteristics of both the voltage-dependent calcium channels and the GECO-proteins are added to the calcium homeostasis of the organism Saccaromyces Cerevisiae, a <span class = "lightblue">completely new system</span> is created. As far as we know, in literature, no experimental data can be found about this specific system. At the moment of writing, the wetwork of the iGEM team of Eindhoven University of Technology 2012 is unfortunately not able to provide applicable information. A great step forward could be made when this model can be checked with experimental data. Nevertheless this model does give some <span class = "lightblue">quantitative insight</span> in the behavior.
| + | |
| - | </p>
| + | |
| - | <p>
| + | |
| - | The simplification of the calcium homeostasis in a living yeast cell certainly leads to some imperfections upon the level of the actual physiology. For example, in the current model, we assume that the concentrations of Pmc1 and Pmr1 are directly proportional to the quantity of transcriptionally active Crz1. This is a big <span class = "lightblue">simplification</span> since in real cells, this process involves the increased gene expression through Crz1 followed by translation and transport of the proteins to the respective intracellular destinations. Moreover, the current model does not include the influence from other relevant pathways whereas in real cells, any response to given extracellular stimulus is likely to be the result of complex cross-talk between multiple pathways.
| + | |
| - | </p>
| + | |
| - | <p>
| + | |
| - | Also some other <span class = "lightblue">assumptions</span> need to be taken into account. First, we assumed that only fully dephosphorylated Crz1 molecules in the nucleus are transcriptionally active since this has been shown the case for NFAT1. Although the mechanism of Crz1 translocation in yeast cells is strikingly similar to NFAT, not all processes can be regarded the same, since there are no experimental data found to validate this assumption. Second, we assumed the behavior of the voltage-dependent calcium channels in yeast cells to be the same as in sympathetic ganglion 'B' type cells of a bullfrog.
| + | |
| - | </p>
| + | |
| - | <p>
| + | |
| - | The results, as presented in the previous section, show a large drop of the cytosolic Ca<sup>2+</sup>-level after removal of the potential difference. Although we could expect a fast recovery to the steady state, this drop is remarkably fast. This is one of the <span class = "lightblue">unexpected results</span> of this model and therefore one of the more interesting features to validate with experiments. Especially since this drop will influence the concentration of calcium bounded GECO-protein.
| + | |
| - | </p>
| + | |
| - | <p>
| + | |
| - | In order to provide more quantitative insights in the experimental setup for the total project, we have to obtain information about the fluorescent properties of the GECO-proteins. At this moment we can only predict when the CaGECO concentration reaches a set percentage of the maximum amount of Ca<sup>2+</sup>-bounded GECO-protein. An important value is the amount of calcium bounded GECO-protein that needs to be there in order to detect fluorescent light.
| + | |
| - | </p>
| + | |
| - | | + | |
| - | | + | |
| - | <h3>
| + | |
| - | Conclusion
| + | |
| - | </h3>
| + | |
| - | | + | |
| - | <p>
| + | |
| - | To examine the calcium in the cell, a dynamic calcium model was designed. As a start a <span class = "lightblue">basic model of calcium homeostasis in yeast cells</span> was presented and tested. This model was extended to include overexpression of voltage-dependent calcium channels and addition of GECO-kinetics. With this model, we can <span class = "lightblue">test and verify theoretical hypotheses</span> by comparing simulation results with corresponding experimental results and generate new hypotheses on the regulation of calcium homeostasis. On the other hand, due to the existence of unknown kinetic parameters and the lack of experimental data, this model is <span class = "lightblue">not an exact model yet</span>. However, it does can give some quantitative insight into the possible dynamics of the whole process and provide a <span class = "lightblue">general framework</span> for more elaborate investigations.
| + | |
| - | </p>
| + | |
| - | | + | |
| - | | + | |
| - | <p>
| + | |
| - | <b>We really hope that in the future of the iGEM competition one of the new iGEM teams will adapt our modeling project of calcium dynamics in yeast cells! :-)</b>
| + | |
| - | </p>
| + | |
| - |
| + | |
| - | <h3>References</h3>
| + | |
| - | <html>
| + | |
| - | <ul>
| + | |
| - | <li><a href="#text_miseta" name="ref_miseta">[1]</a> A. Miseta, L. Fu, R. Kellermayer, J. Buckley, D.M. Bedwell, “The Golgi Apparatus plays a significant role in the maintainance of Ca2+ homeostasis in the vps33 vacuolar biogenesis mutant of Saccharomyces cerevisiae”, J. Biol. Chem. 274: 5939-5947, (1999)</li>
| + | |
| - | <li><a href="#text_kitano" name="ref_kitano">[2]</a> H. Kitano, Systems biology : a brief overview, Nature 295: 1662-1664, (2002)</li>
| + | |
| - | <li><a href="#text_cui" name="ref_cui">[3]</a> J. Cui, Mathematical modeling of metal ion homeostasis and signaling systems, (2009)</li>
| + | |
| - | <li><a href="#text_switch" name="ref_switch">[4]</a> H. Okamura, J. Aramburu, C. Garca-Rodrguez, et al. Concerted dephosphorylation of the transcription factor NFAT1 induces a conformational switch that regulates transcriptional activity,
| + | |
| - | Mol. Cell 6: 539-550, (2000)</li>
| + | |
| - | </ul>
| + | |
| - | </html>
| + | |
| - | {{:Team:TU-Eindhoven/Templates/footer}}
| + | |
"
