Team:BGU Israel/Modeling/Aspiration Shift
From 2014.igem.org
(4 intermediate revisions not shown) | |||
Line 47: | Line 47: | ||
<div class="col2"> | <div class="col2"> | ||
<p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/9/9c/BGU14MODEL2glucose_utilisation_vs_pgc1a.png"><img src="https://static.igem.org/mediawiki/2014/9/9c/BGU14MODEL2glucose_utilisation_vs_pgc1a.png" style="height:350px;border:solid"/></a><br> | <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/9/9c/BGU14MODEL2glucose_utilisation_vs_pgc1a.png"><img src="https://static.igem.org/mediawiki/2014/9/9c/BGU14MODEL2glucose_utilisation_vs_pgc1a.png" style="height:350px;border:solid"/></a><br> | ||
- | Figure 1 | + | Figure 1 – Comparison between relative changes in PGC-1α overexpression (%)<br> |
and changes in insulin-stimulated glucose utilization (%)</p> | and changes in insulin-stimulated glucose utilization (%)</p> | ||
Line 68: | Line 68: | ||
<p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/d/d9/BGU14MODEL2modelDiagram.png"><img src="https://static.igem.org/mediawiki/2014/d/d9/BGU14MODEL2modelDiagram.png" style="height:215px;border:solid"/></a><br> | <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/d/d9/BGU14MODEL2modelDiagram.png"><img src="https://static.igem.org/mediawiki/2014/d/d9/BGU14MODEL2modelDiagram.png" style="height:215px;border:solid"/></a><br> | ||
- | Figure 2 | + | Figure 2 – The diagram of the model: The blue objects represent the variables and the yellow circles represent the interrelationship between the variables. The treatment comes in the form of the two plasmids, Dp and Dr, that we secrete into the body<br> |
</p> | </p> | ||
<p> A reaction of translation, for example, assumed to have a Michaelis-Menten behavior, when only the DNA transcription factor complex is able to make mRNA. On equation 1.1, the first part demonstrates the synthesis of the PGC-1α mRNA while the repressor acts as an inhibitor.</p> | <p> A reaction of translation, for example, assumed to have a Michaelis-Menten behavior, when only the DNA transcription factor complex is able to make mRNA. On equation 1.1, the first part demonstrates the synthesis of the PGC-1α mRNA while the repressor acts as an inhibitor.</p> | ||
Line 86: | Line 86: | ||
<p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/e/ee/BGU14MODEL2fig4.png"><img src="https://static.igem.org/mediawiki/2014/e/ee/BGU14MODEL2fig4.png" style=" width::400px;"/></a> | <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/e/ee/BGU14MODEL2fig4.png"><img src="https://static.igem.org/mediawiki/2014/e/ee/BGU14MODEL2fig4.png" style=" width::400px;"/></a> | ||
<br> | <br> | ||
- | + | </p> | |
<p style="font-size:14px">(*) will not change during treatment</p> | <p style="font-size:14px">(*) will not change during treatment</p> | ||
<p><strong>Parameters of reactions</strong></p> | <p><strong>Parameters of reactions</strong></p> | ||
<p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/7/72/BGU14MODEL2fig5.png"><img src="https://static.igem.org/mediawiki/2014/7/72/BGU14MODEL2fig5.png" style=" height:230px;"/></a> | <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/7/72/BGU14MODEL2fig5.png"><img src="https://static.igem.org/mediawiki/2014/7/72/BGU14MODEL2fig5.png" style=" height:230px;"/></a> | ||
<br> | <br> | ||
- | + | </p> | |
<p><strong>Doses of the treatment</strong></p> | <p><strong>Doses of the treatment</strong></p> | ||
<p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/e/e1/BGU14MODEL2fig6.png"><img src="https://static.igem.org/mediawiki/2014/e/e1/BGU14MODEL2fig6.png" style=" height:62px;"/></a> | <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/e/e1/BGU14MODEL2fig6.png"><img src="https://static.igem.org/mediawiki/2014/e/e1/BGU14MODEL2fig6.png" style=" height:62px;"/></a> | ||
<br> | <br> | ||
- | + | </p> | |
<br> | <br> | ||
<br> | <br> | ||
Line 104: | Line 104: | ||
<br> | <br> | ||
<br> | <br> | ||
- | + | ||
</div> | </div> | ||
</div> | </div> | ||
- | + | <div id="test2" style="position:relative; bottom:70px"></div> | |
<div class="textCont" style="height: 2100px; margin-bottom: 12px;"> | <div class="textCont" style="height: 2100px; margin-bottom: 12px;"> | ||
<h3 style="border-bottom:dashed;border-color:#000000">Simulations Results & Conclusion</h3> | <h3 style="border-bottom:dashed;border-color:#000000">Simulations Results & Conclusion</h3> | ||
Line 118: | Line 118: | ||
<a href="https://static.igem.org/mediawiki/2014/5/53/BGU14MODEL2fig8.png"><img src="https://static.igem.org/mediawiki/2014/5/53/BGU14MODEL2fig8.png" style=" height:280px; border:solid; "/></a> | <a href="https://static.igem.org/mediawiki/2014/5/53/BGU14MODEL2fig8.png"><img src="https://static.igem.org/mediawiki/2014/5/53/BGU14MODEL2fig8.png" style=" height:280px; border:solid; "/></a> | ||
<br> | <br> | ||
- | Figure | + | Figure 3, 4 - The general description of the influence of the two plasmids on the system, while the two plasmids (DP and DR) are added at time=5. The elements of PGC-1α, the active PPARγ and the repressor are rising until a specific value and stay constant at steady state</p> |
<p>Using Wolfram Mathematica we studied the behavior of the dynamic system. The system showed similar behavior in various ranges of parameters values. A steady state is always achieved, stabilizes the levels of the variables in the system into values that depend in the given parameters.</p> | <p>Using Wolfram Mathematica we studied the behavior of the dynamic system. The system showed similar behavior in various ranges of parameters values. A steady state is always achieved, stabilizes the levels of the variables in the system into values that depend in the given parameters.</p> | ||
<p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/7/73/BGU14MODEL2mathematica1.png"><img src="https://static.igem.org/mediawiki/2014/7/73/BGU14MODEL2mathematica1.png" style=" height:225px; border:solid"/></a> | <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/7/73/BGU14MODEL2mathematica1.png"><img src="https://static.igem.org/mediawiki/2014/7/73/BGU14MODEL2mathematica1.png" style=" height:225px; border:solid"/></a> | ||
<br> | <br> | ||
- | Figure | + | Figure 5 - Parameters manipulation in Wolfram Mathematica</p> |
<p>Another method for investigation the steady state behavior is the eigenvalues of the corresponding Jacobian matrix of the dynamic system. By defining our system variable as the vector (PGC1αmRNA, PGC1α, REPmRNA, REP), the Jacobian matrix is:</p> | <p>Another method for investigation the steady state behavior is the eigenvalues of the corresponding Jacobian matrix of the dynamic system. By defining our system variable as the vector (PGC1αmRNA, PGC1α, REPmRNA, REP), the Jacobian matrix is:</p> | ||
<p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/6/6e/BGU14MODEL2fig12.png"><img src="https://static.igem.org/mediawiki/2014/6/6e/BGU14MODEL2fig12.png" style=" height:96px;border:solid;"/></a> | <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/6/6e/BGU14MODEL2fig12.png"><img src="https://static.igem.org/mediawiki/2014/6/6e/BGU14MODEL2fig12.png" style=" height:96px;border:solid;"/></a> | ||
<br> | <br> | ||
- | + | </p> | |
<p>In every set of parameters chosen, we could calculate the eigenvalues of J and found the behavior of the system at steady state: | <p>In every set of parameters chosen, we could calculate the eigenvalues of J and found the behavior of the system at steady state: | ||
Line 137: | Line 137: | ||
<p>Assignment of the parameters values we showed above shows that the system will be in a stable stationary point after the introduction of the treatment. </p> | <p>Assignment of the parameters values we showed above shows that the system will be in a stable stationary point after the introduction of the treatment. </p> | ||
<p><strong>Regulated vs. Unregulated treatment</strong><br> | <p><strong>Regulated vs. Unregulated treatment</strong><br> | ||
- | The modest expression mechanism has to meet two challenges. First to prevent too much synthesis of PGC-1α, that would make an increment over 80% in its concentration. The second challenge is the need to be very potent in cells with low PGC-1α activity, but to have minimum effect in cells with regular PGC-1α activity. The synthesis of the repressor supposes to meet the first challenge. As we can see in Figure | + | The modest expression mechanism has to meet two challenges. First to prevent too much synthesis of PGC-1α, that would make an increment over 80% in its concentration. The second challenge is the need to be very potent in cells with low PGC-1α activity, but to have minimum effect in cells with regular PGC-1α activity. The synthesis of the repressor supposes to meet the first challenge. As we can see in Figure 6, an unregulated mechanism without a repressor would lead to much more production of PGC-1α and could cross to a risk increment. To prove our mechanism produce around 80% of production we should test it in the lab, the model here just illustrate the potential. </p> |
Line 147: | Line 147: | ||
<p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/9/91/BGU14MODEL2fig9.png"><img src="https://static.igem.org/mediawiki/2014/9/91/BGU14MODEL2fig9.png" style=" height:270px; border:solid"/></a> | <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/9/91/BGU14MODEL2fig9.png"><img src="https://static.igem.org/mediawiki/2014/9/91/BGU14MODEL2fig9.png" style=" height:270px; border:solid"/></a> | ||
<br> | <br> | ||
- | Figure | + | Figure 6 - The threshold located at 80% increment of PGC-1α. The PGC-1α level does not pass the threshold as a result of the regulated treatment, unlike the unregulated treatment</p> |
<p>The input level of PGC-1α, in a steady state before the introduction of the plasmids, can be simply determined like this:</p> | <p>The input level of PGC-1α, in a steady state before the introduction of the plasmids, can be simply determined like this:</p> | ||
<p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/f/fe/BGU14MODEL2fig13.PNG"><img src="https://static.igem.org/mediawiki/2014/f/fe/BGU14MODEL2fig13.PNG" style=" height:60px;border:solid;"/></a> | <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/f/fe/BGU14MODEL2fig13.PNG"><img src="https://static.igem.org/mediawiki/2014/f/fe/BGU14MODEL2fig13.PNG" style=" height:60px;border:solid;"/></a> | ||
<br> | <br> | ||
</p> | </p> | ||
- | <p>To face the second challenge we examined different values of natural rate of production of PGC-1α. We ran the model with different input rates and measured the increment in PGC-1α from its initial value to its final value on steady state (Figure | + | <p>To face the second challenge we examined different values of natural rate of production of PGC-1α. We ran the model with different input rates and measured the increment in PGC-1α from its initial value to its final value on steady state (Figure 7).</p> |
<p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/7/77/BGU14MODEL2fig10.png"><img src="https://static.igem.org/mediawiki/2014/7/77/BGU14MODEL2fig10.png" style=" height:275px;border:solid;"/></a><br> | <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/7/77/BGU14MODEL2fig10.png"><img src="https://static.igem.org/mediawiki/2014/7/77/BGU14MODEL2fig10.png" style=" height:275px;border:solid;"/></a><br> | ||
- | Figure | + | Figure 7 - For different initial levels of PGC-1α, the treatment increased it in different relative amounts</p> |
<p>PGC-1α level stabilizes at the steady state of the system. If we give our regulated treatment to people with different natural production rate of PGC-1α from each other, the system stabilizes on different increment. This important feature provides the personalized medicine properties to our regulated treatment we designed.</p> | <p>PGC-1α level stabilizes at the steady state of the system. If we give our regulated treatment to people with different natural production rate of PGC-1α from each other, the system stabilizes on different increment. This important feature provides the personalized medicine properties to our regulated treatment we designed.</p> | ||
Line 161: | Line 161: | ||
<p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/1/11/BGU14MODEL2fig11.png"><img src="https://static.igem.org/mediawiki/2014/1/11/BGU14MODEL2fig11.png" style=" height:400px; border:solid"/></a> | <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/1/11/BGU14MODEL2fig11.png"><img src="https://static.igem.org/mediawiki/2014/1/11/BGU14MODEL2fig11.png" style=" height:400px; border:solid"/></a> | ||
<br> | <br> | ||
- | Figure | + | Figure 8 - The increment of PGC-1α decreases with the growth of the natural production rate</p> |
- | <p>At low input of PGC-1α we can see a much more drastic raise in its concentration and a fixation on higher increment than high input values. On Figure | + | <p>At low input of PGC-1α we can see a much more drastic raise in its concentration and a fixation on higher increment than high input values. On Figure 8 we can only see the final increment on the steady state. We clearly see a decline which shows that in cells with natural high levels of PGC-1α the influence of the plasmids will be insignificant. </p> |
<p style="font-size:20px"><strong><u>Conclusion:</u></strong></p> | <p style="font-size:20px"><strong><u>Conclusion:</u></strong></p> | ||
<p>The model shows that our treatment is personalized and self-regulated. In other words, for those who have inefficient internal production rate, the increment of PGC-1α, as a result of our treatment, will be the more significant than healthy and lean people. This concludes and confirms our method of design in order to achieve the main goal for this treatment. </p> | <p>The model shows that our treatment is personalized and self-regulated. In other words, for those who have inefficient internal production rate, the increment of PGC-1α, as a result of our treatment, will be the more significant than healthy and lean people. This concludes and confirms our method of design in order to achieve the main goal for this treatment. </p> |
Latest revision as of 00:18, 18 October 2014
![](https://static.igem.org/mediawiki/2014/6/6a/BGU14topD.png)
![](https://static.igem.org/mediawiki/2014/b/b2/BGU14Model2banner.png)
Overview
We designed a mathematical model of our ‘Aspiration Shift’ mechanism, and analyzed the treatment effects on the dynamic system of the genetic circuit.
The Mathematical Model
We described the dynamic circuits with differential equations, and used MATLAB SimBiology and Wolfram Mathemitica to analyze their numeric and analytical solution.
Results & Conclusion
The model describes the regulation advantages and the flexibility nature of the mechanism.
Overview
Our objective in ‘Aspiration Shift’ was to create a self-regulated mechanism leading to a limited expression of PGC-1a. To ensure this important feature, we designed a mathematical model that describes this mechanism. We wanted to explore the effects of the treatment, which are the introduction of two plasmids, on PGC-1a levels.
Lower and upper bound of the treatment effects:
Lower Bound
The model is based on the fact that the internal production of PGC-1a in a person cannot change by the treatment. In other words, we wanted our treatment to only make an improvement in the PGC-1a production and the production rate will not decrease below the natural rate.
Upper Bound
The main goal of our model is to verify that the effect of the treatment will not be the same for lean person and for obese person. In lean healthy people, the mean PGC-1α is higher than in obese people [1]. So regardless to the different starting points, the increment of PGC-1α due to our treatment will not reach to a dangerous level, no matter whether it is lean or obese person.
Figure 1 – Comparison between relative changes in PGC-1α overexpression (%)
and changes in insulin-stimulated glucose utilization (%)
We know that an excessive increment in PGC-1α could harm the glucose utilization [2]. We set this threshold of increment to 80%. An increment to any person that will make it over 80%, will not be beneficial, and will even decrease the glucose utilization.
The Mathematical Model
Methods & Tools
While working on the model we used MATLAB SimBiology for the design of the model diagram and for running simulations on various parameters and variables. The solving of the differential equation system were conducted using Simbiology, as well as the creation of most of the plots that we show here. Another powerful tool we used for analyzing our model was Wolfram Mathematica. We used Mathematica for investigating the characteristics of the equations in a more analytical way.
The Model Design
We conducted the design of our model of the dynamic system using MATLAB SimBiology. We defined all the variables and the dynamic interactions, which are relevant to the genetic circuit. The behavior of the circuit depends on a group of differential equations that we constructed from the diagram, which describe chemical reactions and the change of the dynamic variables with time.
Figure 2 – The diagram of the model: The blue objects represent the variables and the yellow circles represent the interrelationship between the variables. The treatment comes in the form of the two plasmids, Dp and Dr, that we secrete into the body
A reaction of translation, for example, assumed to have a Michaelis-Menten behavior, when only the DNA transcription factor complex is able to make mRNA. On equation 1.1, the first part demonstrates the synthesis of the PGC-1α mRNA while the repressor acts as an inhibitor.
The system of the differential equations
Because the repressor and the transcription factor bind to different sites on the plasmid, the inhibition is assumed to be noncompetitive. The 1.2 equation describes the synthesis of PGC-1α in translation from its mRNA, and its natural synthesis in the cell. Equation 2.3 describes equilibrium between PPARγ and its PGC-1α Ligand. The RXR and other factors, which have an influence on the active complex, were neglected because we assumed they can join into the constant value, so the PPARγactive element is acting as an active transcription factor. Although degradation of mRNA and translation of a protein are dependent on the amount of base pairs, we assumed the same values of the PGC-1α gene to the repressor gene.
The initial values of the system variables
(*) will not change during treatment
Parameters of reactions
Doses of the treatment
Simulations Results & Conclusion
General behavior
In order to simulate the effects of our treatment, we added the two plasmids, DP and DR, and created a plot of the concentration levels of all the relevant variables.
Figure 3, 4 - The general description of the influence of the two plasmids on the system, while the two plasmids (DP and DR) are added at time=5. The elements of PGC-1α, the active PPARγ and the repressor are rising until a specific value and stay constant at steady state
Using Wolfram Mathematica we studied the behavior of the dynamic system. The system showed similar behavior in various ranges of parameters values. A steady state is always achieved, stabilizes the levels of the variables in the system into values that depend in the given parameters.
Figure 5 - Parameters manipulation in Wolfram Mathematica
Another method for investigation the steady state behavior is the eigenvalues of the corresponding Jacobian matrix of the dynamic system. By defining our system variable as the vector (PGC1αmRNA, PGC1α, REPmRNA, REP), the Jacobian matrix is:
In every set of parameters chosen, we could calculate the eigenvalues of J and found the behavior of the system at steady state:
- All the eigenvalues have negative real parts: Stationary point (stable)
- One of the eigenvalues has a positive real part: Stationary point (unstable)
- If the largest real part of the eigenvalues is zero: No stable point
Assignment of the parameters values we showed above shows that the system will be in a stable stationary point after the introduction of the treatment.
Regulated vs. Unregulated treatment
The modest expression mechanism has to meet two challenges. First to prevent too much synthesis of PGC-1α, that would make an increment over 80% in its concentration. The second challenge is the need to be very potent in cells with low PGC-1α activity, but to have minimum effect in cells with regular PGC-1α activity. The synthesis of the repressor supposes to meet the first challenge. As we can see in Figure 6, an unregulated mechanism without a repressor would lead to much more production of PGC-1α and could cross to a risk increment. To prove our mechanism produce around 80% of production we should test it in the lab, the model here just illustrate the potential.
Figure 6 - The threshold located at 80% increment of PGC-1α. The PGC-1α level does not pass the threshold as a result of the regulated treatment, unlike the unregulated treatment
The input level of PGC-1α, in a steady state before the introduction of the plasmids, can be simply determined like this:
To face the second challenge we examined different values of natural rate of production of PGC-1α. We ran the model with different input rates and measured the increment in PGC-1α from its initial value to its final value on steady state (Figure 7).
Figure 7 - For different initial levels of PGC-1α, the treatment increased it in different relative amounts
PGC-1α level stabilizes at the steady state of the system. If we give our regulated treatment to people with different natural production rate of PGC-1α from each other, the system stabilizes on different increment. This important feature provides the personalized medicine properties to our regulated treatment we designed.
Figure 8 - The increment of PGC-1α decreases with the growth of the natural production rate
At low input of PGC-1α we can see a much more drastic raise in its concentration and a fixation on higher increment than high input values. On Figure 8 we can only see the final increment on the steady state. We clearly see a decline which shows that in cells with natural high levels of PGC-1α the influence of the plasmids will be insignificant.
Conclusion:
The model shows that our treatment is personalized and self-regulated. In other words, for those who have inefficient internal production rate, the increment of PGC-1α, as a result of our treatment, will be the more significant than healthy and lean people. This concludes and confirms our method of design in order to achieve the main goal for this treatment.
References
- RK Semple, VC Crowley, CP Sewter, M Laudes, C Christodoulides, RV Considine, A Vidal-Puig,, S O’Rahilly. Expression of the thermogenic nuclear hormone receptor coactivator PGC-1a is reduced in the adipose tissue of morbidly obese subjects (2003). doi:10.1038/sj.ijo.0802482.
- C. R. Benton, G. P. Holloway, X.-X. Han, Y. Yoshida, L. A. Snook, J. Lally, J. F. C. Glatz, J. J. F. P. Luiken, A. Chabowski, A. Bonen. Increased levels of peroxisome proliferator-activated receptor gamma, coactivator 1 alpha (PGC-1α) improve lipid utilisation, insulin signalling and glucose transport in skeletal muscle of lean and insulin-resistant obese Zucker rats (2010). doi:10.1007/s00125-010-1773-1.