Difference between revisions of "Team:NEFU China/Modeling"

m
m
Line 93: Line 93:
 
<p>&nbsp;</p>
 
<p>&nbsp;</p>
  
<p><span style="font-family:arial,helvetica,sans-serif"><span style="font-size:16px">(1)&nbsp;Using&nbsp;the&nbsp;random&nbsp;sampling&nbsp;method&nbsp;to&nbsp;take&nbsp;parameters&nbsp;values&nbsp;out&nbsp;twice&nbsp;time, </span></span><span style="font-size:16px">、d&cent;,</span></p>
+
<p><span style="font-family:arial,helvetica,sans-serif"><span style="font-size:16px">(1)&nbsp;Using&nbsp;the&nbsp;random&nbsp;sampling&nbsp;method&nbsp;to&nbsp;take&nbsp;parameters&nbsp;values&nbsp;out&nbsp;twice&nbsp;time, <img alt="" src="https://static.igem.org/mediawiki/2015/b/b2/NEFU_China_9B9BD5F4-1F8C-41F9-91C0-39CDDED254B6.png" style="height:13px; width:10px" /></span></span><span style="font-size:16px">、<img alt="" src="https://static.igem.org/mediawiki/2015/0/0a/NEFU_China_20E9D51F-F0CC-4604-A603-0B1FCB63971C.png" style="height:14px; width:15px" />,</span></p>
  
<p><span style="font-size:16px"><span style="font-family:arial,helvetica,sans-serif">d&nbsp;=&nbsp;(h,z),&nbsp;d&cent;&nbsp;=&nbsp;(h&cent;,z&cent;)&nbsp;;</span></span></p>
+
<p><span style="font-size:16px"><span style="font-family:arial,helvetica,sans-serif"><img alt="" src="https://static.igem.org/mediawiki/2015/d/d0/NEFU_China_F26C4288-C8DB-4711-B3A6-14462BEE1336.png" style="height:14px; width:50px" />,&nbsp;<img alt="" src="https://static.igem.org/mediawiki/2015/8/87/NEFU_China_D1CB3BF5-B443-4C3D-AFA3-447EA82999BC.png" style="height:14px; width:50px" />&nbsp;;</span></span></p>
  
 
<p>&nbsp;</p>
 
<p>&nbsp;</p>
  
<p><span style="font-size:16px"><span style="font-family:arial,helvetica,sans-serif">(2)&nbsp;We&nbsp;can&nbsp;divide&nbsp;the&nbsp;parameters&nbsp;into&nbsp;two&nbsp;parts:&nbsp;the&nbsp;parameter&nbsp;K&nbsp;as&nbsp;a&nbsp;part,&nbsp;its&nbsp;value&nbsp;is&nbsp;h.&nbsp;the&nbsp;remaining&nbsp;parameters&nbsp;as&nbsp;another&nbsp;part,&nbsp;they&nbsp;correspond&nbsp;to&nbsp;z.</span></span></p>
+
<p><span style="font-size:16px"><span style="font-family:arial,helvetica,sans-serif">(2)&nbsp;We&nbsp;can&nbsp;divide&nbsp;the&nbsp;parameters&nbsp;into&nbsp;two&nbsp;parts:&nbsp;the&nbsp;parameter&nbsp;K&nbsp;as&nbsp;a&nbsp;part,&nbsp;its&nbsp;value&nbsp;is&nbsp;<img alt="" src="https://static.igem.org/mediawiki/2015/a/ac/NEFU_China_051EC49E-362E-4799-840C-21C0113BC18C.png" style="height:14px; width:10px" />.&nbsp;the&nbsp;remaining&nbsp;parameters&nbsp;as&nbsp;another&nbsp;part,&nbsp;they&nbsp;correspond&nbsp;to&nbsp;<img alt="" src="https://static.igem.org/mediawiki/2015/3/37/NEFU_China_59B7FCF0-8CD8-49D6-8A72-0559D86EC93D.png" style="height:20px; width:15px" />.</span></span></p>
  
 
<p>&nbsp;</p>
 
<p>&nbsp;</p>
  
<p><span style="font-size:16px"><span style="font-family:arial,helvetica,sans-serif">(3)&nbsp;Taking&nbsp;the&nbsp;value&nbsp;of&nbsp;the&nbsp;parameters&nbsp;in&nbsp;the&nbsp;first&nbsp;time&nbsp;into&nbsp;model&nbsp;can&nbsp;get&nbsp;a&nbsp;output&nbsp;f&nbsp;(d),&nbsp;that&nbsp;is&nbsp;the&nbsp;oscillation&nbsp;period.</span></span></p>
+
<p><span style="font-size:16px"><span style="font-family:arial,helvetica,sans-serif">(3)&nbsp;Taking&nbsp;the&nbsp;value&nbsp;of&nbsp;the&nbsp;parameters&nbsp;in&nbsp;the&nbsp;first&nbsp;time&nbsp;into&nbsp;model&nbsp;can&nbsp;get&nbsp;a&nbsp;output&nbsp;f&nbsp;(</span></span><img alt="" src="https://static.igem.org/mediawiki/2015/b/b2/NEFU_China_9B9BD5F4-1F8C-41F9-91C0-39CDDED254B6.png" style="height:13px; width:10px" /><span style="font-size:16px"><span style="font-family:arial,helvetica,sans-serif">),&nbsp;that&nbsp;is&nbsp;the&nbsp;oscillation&nbsp;period.</span></span></p>
  
 
<p>&nbsp;</p>
 
<p>&nbsp;</p>
  
<p><span style="font-size:16px"><span style="font-family:arial,helvetica,sans-serif">(4)&nbsp;We&nbsp;need&nbsp;to&nbsp;take&nbsp;the&nbsp;value&nbsp;of&nbsp;parameter&nbsp;K&nbsp;by&nbsp;the&nbsp;first&nbsp;time&nbsp;and&nbsp;the&nbsp;rest&nbsp;of&nbsp;the&nbsp;parameters&rsquo;&nbsp;numerical&nbsp;value&nbsp;by&nbsp;the&nbsp;second&nbsp;time&nbsp;into&nbsp;the&nbsp;model&nbsp;to&nbsp;get&nbsp;the&nbsp;output&nbsp;f&nbsp;=&nbsp;(h,z&cent;);Using&nbsp;the&nbsp;same&nbsp;method&nbsp;to&nbsp;get&nbsp;f&nbsp;=&nbsp;(h&cent;,z&cent;).</span></span></p>
+
<p><span style="font-size:16px"><span style="font-family:arial,helvetica,sans-serif">(4)&nbsp;We&nbsp;need&nbsp;to&nbsp;take&nbsp;the&nbsp;value&nbsp;of&nbsp;parameter&nbsp;K&nbsp;by&nbsp;the&nbsp;first&nbsp;time&nbsp;and&nbsp;the&nbsp;rest&nbsp;of&nbsp;the&nbsp;parameters&rsquo;&nbsp;numerical&nbsp;value&nbsp;by&nbsp;the&nbsp;second&nbsp;time&nbsp;into&nbsp;the&nbsp;model&nbsp;to&nbsp;get&nbsp;the&nbsp;output&nbsp;f&nbsp;=&nbsp;(</span></span><img alt="" src="https://static.igem.org/mediawiki/2015/a/ac/NEFU_China_051EC49E-362E-4799-840C-21C0113BC18C.png" style="height:14px; width:10px" /><span style="font-size:16px"><span style="font-family:arial,helvetica,sans-serif">,</span></span><img alt="" src="https://static.igem.org/mediawiki/2015/3/37/NEFU_China_59B7FCF0-8CD8-49D6-8A72-0559D86EC93D.png" style="height:20px; width:15px" /><span style="font-size:16px"><span style="font-family:arial,helvetica,sans-serif">);Using&nbsp;the&nbsp;same&nbsp;method&nbsp;to&nbsp;get&nbsp;f&nbsp;=&nbsp;(</span></span><img alt="" src="https://static.igem.org/mediawiki/2015/0/0a/NEFU_China_20E9D51F-F0CC-4604-A603-0B1FCB63971C.png" style="height:14px; width:15px" /><span style="font-size:16px"><span style="font-family:arial,helvetica,sans-serif">,</span></span><img alt="" src="https://static.igem.org/mediawiki/2015/3/37/NEFU_China_59B7FCF0-8CD8-49D6-8A72-0559D86EC93D.png" style="height:20px; width:15px" /><span style="font-size:16px"><span style="font-family:arial,helvetica,sans-serif">).</span></span></p>
  
 
<p>&nbsp;</p>
 
<p>&nbsp;</p>
Line 155: Line 155:
 
<p>&nbsp;</p>
 
<p>&nbsp;</p>
  
<p><span style="font-family:arial,helvetica,sans-serif"><span style="font-size:16px">The&nbsp;sensitivity&nbsp;coefficient&nbsp;of&nbsp;parameter&nbsp;K&nbsp;is:</span></span></p>
+
<p><span style="font-family:arial,helvetica,sans-serif"><span style="font-size:16px">The&nbsp;sensitivity&nbsp;coefficient&nbsp;of&nbsp;parameter&nbsp;K&nbsp;is:<img alt="" src="https://static.igem.org/mediawiki/2015/8/8e/NEFU_China_50E223BB-9303-4F88-9541-3A5067650F9B.png" style="height:51px; width:80px" /></span></span></p>
  
 
<p>&nbsp;</p>
 
<p>&nbsp;</p>

Revision as of 21:46, 18 September 2015

Modeling introduction:

 

For the modeling of our project, we used mathematic method to explain the progress and optimize our experimental design. We established some equations to simulate the biological 

process.

 

We have considered all of the factors that influence the blue pigment production, and created the mathematic model based on the biological dynamics. As for the future work, we can use Sobol' sensitivity analysis to predict the contribution rate of the model parameters to the results,and did a predictor among 3 different values, time, the mount of AI-2,  and  the  intensity   of pigment, using the MATLAB software to draw the 3D chart.

 

Equation constructed:

We can use mathematic models to indicate the intracellular syntheses of proteins and the AI-2 transportation and phosphorylation.

 

 

 

For these two parts, we designed a set of equations. They are listed as following:

 

The equations of intracellular synthesis of proteins:

[1]LsrACBD:         

 

[2]LsrR:     

[3]LsrK:      

The equations of AI-2 in the cells: 

[4]Intracellular AI-2:      

[5]Phosphorylated AI-2:     

*mRNA LsrACBD:      

*mRNA LsrRK:    

In those equations:

 

 

 

 

Preditor

1.Sobol' Sensitivity analysis.

Sobol’ sensitivity analysis method can predict the contribution rate of the model parameters to the results.

 

(1) Using the random sampling method to take parameters values out twice time, ,

 ;

 

(2) We can divide the parameters into two parts: the parameter K as a part, its value is . the remaining parameters as another part, they correspond to .

 

(3) Taking the value of the parameters in the first time into model can get a output f (), that is the oscillation period.

 

(4) We need to take the value of parameter K by the first time and the rest of the parameters’ numerical value by the second time into the model to get the output f = (,);Using the same method to get f = (,).

 

(5) Sobol’ sensitivity analysis uses Carlo Monte sampling method to carry out operations, the above calculation steps required to repeat N(≥104) times, the result can be closer to the 

actual situation.

 

(6) We can get these equations:

 

 

 

 

 

 

 

 

 

 

In those equations:

 

f0-----Average value of response output.

 

Du-----Variance of parameter K.

 

Dv-----The sum of the variance of the other single parameter.

 

The sensitivity coefficient of parameter K is:

 

Results analysis: 

 

In the model, the higher the sensitivity coefficient, the more important the kinetic parameters of the AI-2 transportation and phosphorylation in cells.

 

2. 3D Map

We did a predictor among 3 different values, time, the concentration of AI-2, and the intensity of pigment. Using the 3D map to find out the connection among them, that can be reference for customers to use 'yogurt gaurder'.

 

 

 

Modeling Software

MATLAB (matrix laboratory) is a numerical computing envionment and fourth-generation programming language. It is developed by MathWorks, a company in United States. MATLAB allows matrix manipulations, plotting of functions and data, implementation of algorithms, creation of user of interfaces, and interfacing with programs written in other languages, including C, C++, Java, and Fortran. Although MATLAB is intended primarily for numerical computing, an optional toolbox uses the MuPAD symbolic engine, allowing access to symbolic computing capabilities. An additional package, Simulink, adds graphical multi-domain simulation and Model-Based Design for dynamic and embedded systems.

                             

 

 

 

 

 

 

 

Back to Top