Difference between revisions of "Team:UFSCar-Brasil/part1.html"
Celio.dias (Talk | contribs) |
Celio.dias (Talk | contribs) |
||
Line 30: | Line 30: | ||
<h5 class="ui center aligned header"><b>z = \(z_0\)+ax+by+cx²+dy²+fxy</b>.</h5> | <h5 class="ui center aligned header"><b>z = \(z_0\)+ax+by+cx²+dy²+fxy</b>.</h5> | ||
− | <p>This curve obtained r² = 0.95457, | + | <p>This curve obtained r² = 0.95457, which means that more than 95% of data are modeled by this function. Using the <b>MatLab</b> software to function fitting. We also find the value of equation constants (Appendix Table 2).</p> |
<p>The surface equation is:</p> | <p>The surface equation is:</p> | ||
Line 59: | Line 59: | ||
<h5 class="ui center aligned header"> (Eq. 6) $$ H[z(x,y)]=(-0,08774)(0,00178)-(0,00319)^2=-0,0001663533 $$</h5> | <h5 class="ui center aligned header"> (Eq. 6) $$ H[z(x,y)]=(-0,08774)(0,00178)-(0,00319)^2=-0,0001663533 $$</h5> | ||
<p>As H[f(x,y)] | <p>As H[f(x,y)] | ||
− | <0, we know that | + | <0, we know that this function posses a saddle point and there are no minimum or maximum values (Guidorizzi, 1997). Saddle point is understood as the point where the slope is zero, but it is not a maximum or minimum value. In it there |
is the highest elevation in one direction and lowest in the perpendicular direction.</p> | is the highest elevation in one direction and lowest in the perpendicular direction.</p> | ||
<p>Now, known the point we want to find, we take the partial derivatives of the first order function and made them equal to zero. To know:</p> | <p>Now, known the point we want to find, we take the partial derivatives of the first order function and made them equal to zero. To know:</p> | ||
Line 71: | Line 71: | ||
<img class="ui centered 400x400 image" src="https://static.igem.org/mediawiki/2015/f/f4/UFSCar-Brasil_fitted.jpg"> | <img class="ui centered 400x400 image" src="https://static.igem.org/mediawiki/2015/f/f4/UFSCar-Brasil_fitted.jpg"> | ||
<h5 class="ui big center aligned header"><b>Figure 3:</b> Extrapolated fitted surface for extreme values of PEG concentrations and time.</h5> | <h5 class="ui big center aligned header"><b>Figure 3:</b> Extrapolated fitted surface for extreme values of PEG concentrations and time.</h5> | ||
− | <p>It is possible to observe that for very long time values the absorbance decrease to zero very fast and the time do not | + | <p>It is possible to observe that for very long time values the absorbance decrease to zero very fast and the time do not surpass 8 weeks. The PEG concentrations also increase the absorbance values for both extreme concentrations, when approximate |
of zero and pass of 30% of PEG 6000 the absorbance value explodes. Despite mathematically valid, biologically this means that near to zero there will be no control of cell division, and for values near to 30% of PEG 6000, culture reach a plateau | of zero and pass of 30% of PEG 6000 the absorbance value explodes. Despite mathematically valid, biologically this means that near to zero there will be no control of cell division, and for values near to 30% of PEG 6000, culture reach a plateau | ||
of minimum cell division, followed by a decreasing function.</p> | of minimum cell division, followed by a decreasing function.</p> |
Revision as of 02:39, 17 September 2015
Plasmolysis
How long Bug Shoo could be stored?
Fitting Data
Starting from the experimental data presented in Table 1 in the appendix, a 3D graph (Figure 1) was constructed. It was used to better system understanding, which is the relationship between absorbance, time and PEG 6000 concentration.
![](https://static.igem.org/mediawiki/2015/8/87/UFSCar-Brasil_abs_t_peg.jpg)
Figure 1: Experimental plots of absorbance versus time versus PEG concentration.
This graph was used to find best-fitted surface to experimental points. To this, parameter r² was used, because it shows the percentage of points that fit in a certain curve. Thus, we observed the following equation to explain this behaviour:
z = \(z_0\)+ax+by+cx²+dy²+fxy.
This curve obtained r² = 0.95457, which means that more than 95% of data are modeled by this function. Using the MatLab software to function fitting. We also find the value of equation constants (Appendix Table 2).
The surface equation is:
z = 1.07+0.06x-0.044y-0.04x²+0.0009y²+0.003xy (Eq. 1)
Where, X means KKK and Y means KKKK. From known curve was possible to simulate a surface that models most likely behavior of the three variables in question.
![](https://static.igem.org/mediawiki/2015/1/16/UFSCar-Brasil_Fitted_surface.png)
Figure 2: Fitted surface obtained using the Eq. 1.
Results and Discussion
Our goal to find main points of system, meaning those points where plane tangent in surface has a null gradient value. First, it was analyzed (using a Hessian matrix) on these important points to make possible to know if these are the maximum points, minimum or saddle points of the function.
The calculation of the Hessian matrix is given by:
(Eq. 2) $$ H[z(x,y)] = \begin{bmatrix} \frac{\partial^2 z}{\partial x^2} & \frac{\partial^2 z}{\partial z,\partial y} \\ \frac{\partial^2 z}{\partial z,\partial y} & \frac{\partial^2 z}{\partial y^2} \end{bmatrix} $$
thus to find the value of the Hessian of z(x,y) we obtained the partial derivatives of the second order, whose values are below:
(Eq. 3) $$ \frac{\partial^2 z}{\partial x^2} = 2c = - 0,08774 $$
(Eq. 4) $$ \frac{\partial^2 z}{\partial y^2} = 2d = - 0,00178 $$
(Eq. 5) $$ \frac{\partial^2 z}{\partial x,\partial y} = \frac{\partial^2 z}{\partial y,\partial x} = f = 0,00319 $$
with the known partial derivatives, we have the Hessian:
(Eq. 6) $$ H[z(x,y)]=(-0,08774)(0,00178)-(0,00319)^2=-0,0001663533 $$
As H[f(x,y)] <0, we know that this function posses a saddle point and there are no minimum or maximum values (Guidorizzi, 1997). Saddle point is understood as the point where the slope is zero, but it is not a maximum or minimum value. In it there is the highest elevation in one direction and lowest in the perpendicular direction.
Now, known the point we want to find, we take the partial derivatives of the first order function and made them equal to zero. To know:
(Eq. 7) $$ \frac{dz}{dx}=a+2cx+fy=0 $$
(Eq. 8) $$ \frac{dz}{dy}=b+2dy+fx=0 $$
Therefore:
(Eq. 9) $$ -\frac{a+2cx}{f}=y=\frac{-b+fx}{2d} $$
replacing the values that are found in Table 2 of the appendix, we have that $$x=\frac{(2{da}-bf)} {(f²-4cd)}=1,5$$
And thus y=22.2 and z=0.63 . Then it is possible to concluded that for optimal absorbance point (z = 0.63), it is necessary 22.2% PEG 6000 where it gets 1.5 weeks of plasmolysis (Figure 3).
![](https://static.igem.org/mediawiki/2015/f/f4/UFSCar-Brasil_fitted.jpg)
Figure 3: Extrapolated fitted surface for extreme values of PEG concentrations and time.
It is possible to observe that for very long time values the absorbance decrease to zero very fast and the time do not surpass 8 weeks. The PEG concentrations also increase the absorbance values for both extreme concentrations, when approximate of zero and pass of 30% of PEG 6000 the absorbance value explodes. Despite mathematically valid, biologically this means that near to zero there will be no control of cell division, and for values near to 30% of PEG 6000, culture reach a plateau of minimum cell division, followed by a decreasing function.
Conclusions
Through these analyzes we could first find one function that suited in r² = 95%, which is very satisfactory, and validates our whole discussion about the obtained surface. Then, using Hessian matrix, we found a saddle point, which is particularly biologically relevant. The variation of absorbance along time indicates two things to us: first, cell numbers, and second, the metabolic activity of cells, since they are active, they should divide increasing absorbance. Due to its plasmolysed state, variation rates of absorbance values should be low, but not zero, because if so, they would be dead and would not serve for our propose. As it is a saddle point, we know that it is not null. Finally, the point was found through gradient of the function and we concluded that for the minimum value of absorbance the best PEG concentration is 22,2% with a maximum time of 8 weeks of storage.
Appendix
List of tables
Table 1: Experimental dates obtained in microbiology lab
Our amazing sponsors!
![](https://static.igem.org/mediawiki/2015/2/24/UFSCar-Brasil_logos1.png)