Valencia/Controller

From 2007.igem.org

(Difference between revisions)
Line 30: Line 30:
-
par={αx=20,αy=20,αz=20,αv=20,αu=20,α2=20,n1=2,ny=2,n2=2,nx=2,n3=2,nu=2,nv=2,nz=2,k1=40,ky=40,k2=40,kx=40,k3=40,ku=40,kv=40,kz=40,βx=0.069,βy=0.069,βz=0.069 10-3,βv=0.069,βu=0.069,β2=0.069 10-3,β1=0.069 10-5,β3=0.069 10-5,γx=0.01,γy=0.01,γz=0.01,γv=0.01,γu=0.01,γp2=0.000001,x(0)=0,y(0)=0,z(0)=0,v(0)=0,u(0)=0,p2(0)=3,p1(0)=5,p3(0)=1};
+
par={αx=20,αy=20,αz=20,αv=20,αu=20,α2=20,n1=2,ny=2,n2=2,nx=2,n3=2,nu=2,nv=2,nz=2,k1=40,ky=40,k2=40,kx=40,k3=40,ku=40,kv=40,kz=40,
 +
βx=0.069,βy=0.069,βz=0.069 10-3,βv=0.069,βu=0.069,β2=0.069 10-3,β1=0.069 10-5,β3=0.069 10-5,γx=0.01,γy=0.01,γz=0.01,γv=0.01,γu=0.01,
 +
γp2=0.000001,x(0)=0,y(0)=0,z(0)=0,v(0)=0,u(0)=0,p2(0)=3,p1(0)=5,p3(0)=1};
Line 36: Line 38:
-
[[Image:Vlalala copia.jpg|500 px]]
+
[[Image:VLalala copia.jpg|500 px]]
In this figure we can appreciate the typical behaviour of a P-Controller with the associated steady state error. Depending on the desired precision in the control of the variable, this controller could be enough or we need to add another device which performs the integer operation of a PI-controller.   
In this figure we can appreciate the typical behaviour of a P-Controller with the associated steady state error. Depending on the desired precision in the control of the variable, this controller could be enough or we need to add another device which performs the integer operation of a PI-controller.   
Line 65: Line 67:
Some representative results can be seen in the following figures:
Some representative results can be seen in the following figures:
{|
{|
-
|[[Image:V1-controller.jpg|thumb|left|1050 px|Change of initial concentrations of (x,y), (u,v) and (x,y,u,v) ]]
+
|[[Image:V1-controller.jpg|thumb|left|900 px|Change of initial concentrations of (x,y), (u,v) and (x,y,u,v) ]]
|-
|-
-
|[[Image:V2-controller.jpg|thumb|left|1050 px|Study of the z basal state dependency of the system. a) z~0 b) z>> basal state of x, y, u and v c) z ~ basal state of x, y, u and v ]]
+
|[[Image:V2-controller.jpg|thumb|left|900 px|Study of the z basal state dependency of the system. a) z~0 b) z>> basal state of x, y, u and v c) z ~ basal state of x, y, u and v ]]
|-
|-
-
|[[Image:V3-controller.jpg|thumb|left|1050 px|Study of the p2 basal state. a) p2*100  b) p2*1000 c) p2*10000  ]]
+
|[[Image:V3-controller.jpg|thumb|left|900 px|Study of the p2 basal state. a) p2*100  b) p2*1000 c) p2*10000  ]]
|-
|-
-
|[[Image:V4-controller.jpg|thumb|left|1050 px|A progressive reduction of the degradation constant of x shows the behaviour observed in the first two figs, the reduction of the degradation constant of y is show the last one ]]
+
|[[Image:V4-controller.jpg|thumb|left|900 px|A progressive reduction of the degradation constant of x shows the behaviour observed in the first two figs, the reduction of the degradation constant of y is show the last one ]]
|-
|-
-
|[[Image:V5-controller.jpg|thumb|center|1050 px|A progressive reduction of the degradation constant of u shows the behaviour observed in the first two figs, the reduction of the degradation constant of v is show the last one]]
+
|[[Image:V5-controller.jpg|thumb|center|900 px|A progressive reduction of the degradation constant of u shows the behaviour observed in the first two figs, the reduction of the degradation constant of v is show the last one]]
|-
|-
|[[Image:V6-controller.jpg|thumb|left|700 px|Z degradation constant dependency study: a) 10 times higher b) 100 times higher
|[[Image:V6-controller.jpg|thumb|left|700 px|Z degradation constant dependency study: a) 10 times higher b) 100 times higher

Revision as of 07:13, 26 September 2007

VBacteriaAnimada.gif
Summary The idea Promoter calibrator Biological comparator Electronic comparator Biological controller Lab Work Silico Work Results Adversities

Contents

An alternative use of our comparator: The P controller

As it is explained in the idea section, a possible application of the comparator could be to couple two of these devices to build a controller that controls the concentration of certain substance in order to maintain it between certain values.

Electronic Proportional Controller

The controller which we are going to describe is very similar to the what is call in engineering as a proportional controler. Therefore, it could be suitable to give first a brief introduction to this devices in the context of electronics.

These controllers are characterized by the fact that they have a proportional response to the error value of the controlled magnitud.The proportional response can be adjusted by multiplying the error by a constant Kp, called the proportional gain. A typical response of this devices is shown in the following figure:

VControlador proporcional.jpg

A high proportional gain results in a large change in the output for a given change in the error. If the proportional gain is too high, the system can become unstable (See the section on Loop Tuning). In contrast, a small gain results in a small output response to a large input error, and a less responsive (or sensitive) controller. If the proportional gain is too low, the control action may be too small when responding to system disturbances.

Moreover, from that figure we can see an importat characteristic of these systems: In the absence of disturbances pure proportional control will not settle at its target value, but will retain a steady state error that is a function of the proportional gain and the process gain. Depending on the control accuracy desired either we can avoid to consider this error or we should complement this device with another devices which complement it. This devices could be an integer controller, in that case we have a PI controller or also a derivative controller and we obtain a PID controller quite extended in automatic systems.

Biological Proportional Controller

To build this device we should couple two comparators with the proper promoters (by proper we mean, properly calibrated to produce an action at certain level) in such a way that if the concentration is higher or lower than the reference value they activate the required metabolic circuits to reduce or increase the amount of the controlled protein. A basic scheme of the P-Controller is shown in the following figure:

VEstabilizador.jpg

According to the explanation of the last sections, the system of equations which describe that system should be:

VEcuacioneswikicontroller.jpg

Solving this system of equations with the following parameters value:


par={αx=20,αy=20,αz=20,αv=20,αu=20,α2=20,n1=2,ny=2,n2=2,nx=2,n3=2,nu=2,nv=2,nz=2,k1=40,ky=40,k2=40,kx=40,k3=40,ku=40,kv=40,kz=40, βx=0.069,βy=0.069,βz=0.069 10-3,βv=0.069,βu=0.069,β2=0.069 10-3,β1=0.069 10-5,β3=0.069 10-5,γx=0.01,γy=0.01,γz=0.01,γv=0.01,γu=0.01, γp2=0.000001,x(0)=0,y(0)=0,z(0)=0,v(0)=0,u(0)=0,p2(0)=3,p1(0)=5,p3(0)=1};


Where the parameters values were selected according the ones given by ref[1], and performing the simulation for t=100000.We observe the system behaviour shown in the following figure


VLalala copia.jpg

In this figure we can appreciate the typical behaviour of a P-Controller with the associated steady state error. Depending on the desired precision in the control of the variable, this controller could be enough or we need to add another device which performs the integer operation of a PI-controller.


Sensibility of the system to the different parameters

Later on we are going to describe the sensibility of the device to different parameters of the model in order to understand in more detail the capabilities of the system and in which conditions could operate:

In this analysis, the dependence of the system with the following parameters was studied:

  • Initial concentrations of “x”, “y”, “u”, “v”
  • Basal constants of those variables
  • Basal constant of “z”.
  • Basal constant of “p2”.
  • Degradation constants of “x”, “y”, “u”, “v”
  • Degradation constant of “z”
  • Degradation constant of “p2”
  • Degradation constants of “p1” and “p3”

Some representative results can be seen in the following figures:

Change of initial concentrations of (x,y), (u,v) and (x,y,u,v)
Study of the z basal state dependency of the system. a) z~0 b) z>> basal state of x, y, u and v c) z ~ basal state of x, y, u and v
Study of the p2 basal state. a) p2*100 b) p2*1000 c) p2*10000
A progressive reduction of the degradation constant of x shows the behaviour observed in the first two figs, the reduction of the degradation constant of y is show the last one
A progressive reduction of the degradation constant of u shows the behaviour observed in the first two figs, the reduction of the degradation constant of v is show the last one
Z degradation constant dependency study: a) 10 times higher b) 100 times higher
P2 degradation tendency study:a) representative behaviour of an increase in the constant b) representative behaviour of a decrease in the constant
Change in the basal constant of x, y, u and v
Study of the influence of the degradation constant of p1 and p3 in the performance of the controller.


The simulations of that system can be found in the mathematica file: [http://www.igem.upv.es/igem07/images/3/37/Stabilizator.rar Mathematica model file]

The most important conclusions of that work were:

  • The initial concentrations of x,y,u and v do not have too much influence in the results
  • The Z basal state reduce under the consigned levels the concentration of the controlled magnitude (p2)
  • The basal concentration of p2 should be in the same range as the basal concentration of the controllers because if we increase it too much the controller does not work
  • The basal state of x,y,u and v do not have any influence in the final controller result
  • The degradation constant of x and u reduces significantly the efficiency of the controller, but the degradation constants of y and v do not have any influence on the system
  • If the degradation constant of z increases the system lose its hability to operate as a controller
  • If the controlled variable (p2) suffers a too large or too small degradation the controller do not work
  • Although p1 and p3 becomes degradated the system works


It's a pity that it is not so easy to manipulate the value of the different parameter of a biological system (like the strenght of a promoter) than to manipulate the parameters of an electronic device :-( !!

Nevertheless, this theoretical development constitutes a very good example of how modelling could be combined with biology in order to implement a system.

Regarding the stability of the system, developing a similar study to the one developed for the promoter calibrator, we see that the local lyapunov coefficients are also lower than 0, therefore this system is also stable as it can be expected from the previous figures.