Team:BGU Israel/Modeling/Aspiration Shift

From 2014.igem.org

(Difference between revisions)
 
(11 intermediate revisions not shown)
Line 6: Line 6:
-
 
 
       <section>
       <section>
       <div style="height:12px"> </div>
       <div style="height:12px"> </div>
Line 19: Line 18:
         <h3 style="border-bottom:dashed;border-color:#000000">The Mathematical Model<a href="#" onclick="goToByScroll('test'); return false;" class="right"></a></h3>
         <h3 style="border-bottom:dashed;border-color:#000000">The Mathematical Model<a href="#" onclick="goToByScroll('test'); return false;" class="right"></a></h3>
         <br>
         <br>
 +
We described the dynamic circuits with differential equations, and used MATLAB SimBiology and Wolfram Mathemitica to analyze their numeric and analytical solution.
         </div>
         </div>
          
          
Line 24: Line 24:
         <h3 style="border-bottom:dashed;border-color:#000000">Results & Conclusion<a href="#" onclick="goToByScroll('test2'); return false;" class="right"></a></h3>
         <h3 style="border-bottom:dashed;border-color:#000000">Results & Conclusion<a href="#" onclick="goToByScroll('test2'); return false;" class="right"></a></h3>
         <br>
         <br>
-
          
+
         The model describes the regulation advantages and the flexibility nature of the mechanism.
         </div>
         </div>
         <div class="clear"></div>
         <div class="clear"></div>
Line 35: Line 35:
       <h3 style="border-bottom:dashed;border-color:#000000">Overview </h3>
       <h3 style="border-bottom:dashed;border-color:#000000">Overview </h3>
         <div class="col2">
         <div class="col2">
-
           <p>Our objective in &lsquo;Aspiration  Shift&rsquo; 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.</p>
+
           <p>Our objective in &lsquo;Aspiration  Shift&rsquo; 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.</p>
         <p><strong>Lower and upper bound of the treatment effects:</strong><br>
         <p><strong>Lower and upper bound of the treatment effects:</strong><br>
           <u>Lower Bound</u><br>
           <u>Lower Bound</u><br>
-
           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.<br>
+
           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.<br>
   <u>Upper Bound</u></p>
   <u>Upper Bound</u></p>
-
<p> The main goal of our model is to verify that the  effect of the treatment will <u>not</u> be the same for lean person and for  obese person. In lean healthy people, the mean PGC-1a is higher than in obese people [1].So regardless to the  different starting points, the increment of PGC-1a due to our treatment will not reach to a dangerous  level, no matter whether it is lean or obese person.<div id="test"></div>
+
<p> The main goal of our model is to verify that the  effect of the treatment will <u>not</u> be the same for lean person and for  obese person. In lean healthy people, the mean PGC-is higher than in obese people [1]. So regardless to the  different starting points, the increment of PGC-due to our treatment will not reach to a dangerous  level, no matter whether it is lean or obese person.<div id="test"></div>
        
        
</p>
</p>
Line 46: Line 46:
         </div>
         </div>
         <div class="col2">
         <div class="col2">
-
   <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><img  src="https://static.igem.org/mediawiki/2014/9/9c/BGU14MODEL2glucose_utilisation_vs_pgc1a.png" style="height:350px;border:solid"/><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 –Comparison between relative changes in PGC-1a overexpression (%)<br>
+
           Figure 1 – Comparison between relative changes in PGC-overexpression (%)<br>
 and changes in insulin-stimulated  glucose utilization (%)</p>
 and changes in insulin-stimulated  glucose utilization (%)</p>
    
    
-
<p>We know that an excessive increment in PGC-1a 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.</p>
+
<p>We know that an excessive increment in PGC-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.</p>
Line 61: Line 61:
         <div class="col2" style="margin-left:10px">
         <div class="col2" style="margin-left:10px">
           <p><strong>Methods &amp; Tools</strong><br>
           <p><strong>Methods &amp; Tools</strong><br>
-
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.</p>
+
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.</p>
           <p></p>
           <p></p>
-
           <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><img  src="https://static.igem.org/mediawiki/2014/8/8a/BGU14MODEL2matlab-logo.jpg" style="height:150px;"/><img src="https://static.igem.org/mediawiki/2014/0/00/BGU14MODEL2wolfram-mathematica-logo.jpg" style="height:150px;"/></p>
+
           <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/8/8a/BGU14MODEL2matlab-logo.jpg"><img  src="https://static.igem.org/mediawiki/2014/8/8a/BGU14MODEL2matlab-logo.jpg" style="height:150px;"/></a><a href="https://static.igem.org/mediawiki/2014/0/00/BGU14MODEL2wolfram-mathematica-logo.jpg"><img src="https://static.igem.org/mediawiki/2014/0/00/BGU14MODEL2wolfram-mathematica-logo.jpg" style="height:150px;"/></a></p>
            
            
-
<p><strong>The Model Design</strong><br>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.</p>
+
<p><strong>The Model Design</strong><br>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.</p>
          
          
-
<p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><img  src="https://static.igem.org/mediawiki/2014/d/d9/BGU14MODEL2modelDiagram.png" style="height:215px;border:solid"/><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 –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>
+
           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-1a 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-mRNA while the repressor acts as an inhibitor.</p>
<p><strong>The system of the  differential equations</strong></p>
<p><strong>The system of the  differential equations</strong></p>
-
<p align="center" style="font-size:14px; font-weight:bold; line-height:normal">
+
<p align="center" style="font-size:14px; font-weight:bold; line-height:normal">
 +
<a href="https://static.igem.org/mediawiki/2014/b/b5/BGU14MODEL2fig3.png">
<img  src="https://static.igem.org/mediawiki/2014/b/b5/BGU14MODEL2fig3.png" style="height:215px;  
<img  src="https://static.igem.org/mediawiki/2014/b/b5/BGU14MODEL2fig3.png" style="height:215px;  
-
border:solid"/>
+
border:solid"/></a>
Line 81: Line 82:
          
          
         <div class="col2" style="margin-left:30px">
         <div class="col2" style="margin-left:30px">
-
<p>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-1a in translation from its mRNA, and its natural synthesis in the cell. Equation 2.3 describes equilibrium between PPAR? and its PGC-1a 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-1a gene to the repressor gene.</p><br>
+
<p>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-in translation from its mRNA, and its natural synthesis in the cell. Equation 2.3 describes equilibrium between PPARγ and its PGC-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-gene to the repressor gene.</p><br>
  <p><strong>The initial values of the system variables</strong></p>
  <p><strong>The initial values of the system variables</strong></p>
-
<p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><img  src="https://static.igem.org/mediawiki/2014/e/ee/BGU14MODEL2fig4.png" style=" width::400px;"/>
+
<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>
-
   Figure 3</p>
+
   </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"><img  src="https://static.igem.org/mediawiki/2014/7/72/BGU14MODEL2fig5.png" style=" height:230px;"/>
+
<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>
-
   Figure 4</p>
+
   </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"><img  src="https://static.igem.org/mediawiki/2014/e/e1/BGU14MODEL2fig6.png" style=" height:62px;"/>
+
<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>
-
   Figure 5</p>
+
   </p>
   <br>
   <br>
   <br>
   <br>
Line 103: Line 104:
   <br>
   <br>
   <br>
   <br>
-
  <div id="test2"></div>
+
 
         </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 114: Line 115:
            
            
-
  <p align="center" style="font-size:14px; font-weight:bold;line-height:normal"><img  src="https://static.igem.org/mediawiki/2014/d/d3/BGU14MODEL2fig7.png" style=" height:280px;border:solid;"/>
+
  <p align="center" style="font-size:14px; font-weight:bold;line-height:normal"><a href="https://static.igem.org/mediawiki/2014/d/d3/BGU14MODEL2fig7.png"><img  src="https://static.igem.org/mediawiki/2014/d/d3/BGU14MODEL2fig7.png" style=" height:280px;border:solid;"/></a>
-
  <img  src="https://static.igem.org/mediawiki/2014/5/53/BGU14MODEL2fig8.png" style=" height:280px; border:solid; "/>
+
  <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 6,7- 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-1a, the active PPAR? and the repressor are rising until a specific value and stay constant at steady state</p>
+
   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-, 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"><img src="https://static.igem.org/mediawiki/2014/7/73/BGU14MODEL2mathematica1.png"  style=" height:225px; border:solid"/>
+
   <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 8- Parameters manipulation in Wolfram Mathematica</p>
+
   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 (PGC1amRNA, PGC1a, 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"><img  src="https://static.igem.org/mediawiki/2014/6/6e/BGU14MODEL2fig12.png" style=" height:96px;border:solid;"/>
+
  <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>
-
   Figure 9</p>
+
   </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 136: 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-1a, 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-1a activity, but to have minimum effect in cells with regular PGC-1a activity. The synthesis of the repressor supposes to meet the first challenge.  As we can see in Figure 5, an unregulated mechanism without a repressor would  lead to much more production of PGC-1a 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>
+
     The modest  expression mechanism has to meet two challenges. First to prevent too much  synthesis of PGC-, 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-activity, but to have minimum effect in cells with regular PGC- 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-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 144: Line 145:
   <br>
   <br>
-
   <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><img  src="https://static.igem.org/mediawiki/2014/9/91/BGU14MODEL2fig9.png" style=" height:270px; border:solid"/>
+
   <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 10-The threshold located at 80% increment of PGC-1a. The PGC-1a level does not pass the threshold as a result of the regulated treatment, unlike the unregulated treatment</p>
+
   Figure 6 - The threshold located at 80% increment of PGC-. The PGC-level does not pass the threshold as a result of the regulated treatment, unlike the unregulated treatment</p>
-
   <p>The input level of PGC-1a, in a steady state before the introduction of the plasmids, can be simply determined like this:</p>
+
   <p>The input level of PGC-, 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"><img  src="https://static.igem.org/mediawiki/2014/f/fe/BGU14MODEL2fig13.PNG" style=" height:60px;border:solid;"/>
+
<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>
-
   Figure 11</p>
+
   </p>
-
   <p>To face the second challenge we examined different values of natural rate of production of PGC-1a. We ran the model with different input rates and measured the increment in PGC-1a from its initial value to its final value on steady state (Figure 11).</p>
+
   <p>To face the second challenge we examined different values of natural rate of production of PGC-. We ran the model with different input rates and measured the increment in PGC-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"><img  src="https://static.igem.org/mediawiki/2014/7/77/BGU14MODEL2fig10.png" style=" height:275px;border:solid;"/><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 12-For different initial levels of PGC-1a, the treatment increased it in different relative amounts</p>
+
   Figure 7 - For different initial levels of PGC-, the treatment increased it in different relative amounts</p>
-
<p>PGC-1a level stabilizes at the steady state of the system. If we give our regulated treatment to people with different natural production rate of PGC-1a 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-level stabilizes at the steady state of the system. If we give our regulated treatment to people with different natural production rate of PGC-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 align="center" style="font-size:14px; font-weight:bold; line-height:normal"><img  src="https://static.igem.org/mediawiki/2014/1/11/BGU14MODEL2fig11.png" style=" height:400px; border:solid"/>
+
   <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 13-The increment of PGC-1a decreases with the growth of the natural production rate</p>  
+
   Figure 8 - The increment of PGC-decreases with the growth of the natural production rate</p>  
-
   <p>At low input of PGC-1a we can see (Figure 6) a much more drastic raise in its concentration and a fixation on higher increment than high input values. On Figure 7 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-1a the influence of the plasmids will be insignificant. </p>
+
   <p>At low input of PGC-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-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-, 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>
    
    
    
    
Line 175: Line 176:
       <ol>
       <ol>
         <li><span dir="LTR"> </span>RK Semple, VC  Crowley, CP Sewter, M Laudes, C Christodoulides, RV Considine, A Vidal-Puig,, S  O&rsquo;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.</li>
         <li><span dir="LTR"> </span>RK Semple, VC  Crowley, CP Sewter, M Laudes, C Christodoulides, RV Considine, A Vidal-Puig,, S  O&rsquo;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.</li>
-
         <li><span dir="LTR"> </span>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-1a) 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.</li>
+
         <li><span dir="LTR"> </span>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-) 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.</li>
       </ol>
       </ol>
         </div>
         </div>
Line 185: Line 186:
          
          
       </section>
       </section>
-
 
-
 

Latest revision as of 00:18, 18 October 2014

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

  1. 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.
  2. 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.