Difference between revisions of "Team:Valencia UPV/Modeling/Simulations"
Line 244: | Line 244: | ||
</header> | </header> | ||
<p>We knew from literature, that changes induced by light are instantaneous. However, there are several differences between switching on and off these switches. For example, if we had implemented a red toggle switch, we would be able to ensure the stop of production using far red light, whereas dark periods induce a slower decay. </p> | <p>We knew from literature, that changes induced by light are instantaneous. However, there are several differences between switching on and off these switches. For example, if we had implemented a red toggle switch, we would be able to ensure the stop of production using far red light, whereas dark periods induce a slower decay. </p> | ||
+ | <p>This means that although light irradiation is stopped, the presence of photosensible elements may will decrease slowly, and production will not occur only when light is given. In other words: production will not be totally over light control.</p> | ||
+ | <p>In order to ensure this coordination between light and presence of active elements, we suggested our wet lab mates tagging PhyB and ePDZ, they may will be degraded fast enough to minimize their presence when light has been eliminated. This is the reason why B and D have higher degradation rates than other proteins.</p> | ||
+ | <p>So, if dark periods mean loss of activators, pauses between first and second inputs should be as short as they can be, in order to maximize production. The shorter pause possible, is a no-pause, i.e continuous irradiation.</p> | ||
+ | <p>Firstly, we will compare production of alpha with different pauses between both red inputs, in order to analyze the influence of the dark period length. We implemented in Matlab a for loop to facilitate the comparison:</p> | ||
+ | <p><font face="Courier New,Courier, monospace">pause_dur=[0,2,30,60,80,100]; % vector containing different duration of pauses.</p> | ||
+ | <p>for i=1:length(pause_dur) | ||
+ | <p>light_dur=length(input_t)-351-pause_dur(i);% duration of irradiation</p> | ||
+ | <p>input_red=[zeros(251,1);red_intensity*ones(light_dur/2,1);</p> | ||
+ | <p>zeros(pause_dur(i),1);red_intensity*ones(light_dur/2,1);</p> | ||
+ | <p>zeros(100,1)];</p> | ||
+ | <p>input_blue=[zeros(3101,1)];</p> | ||
+ | <p>[t,x] = ode23s(@model_switch_v2, time_span, ic, opt,input_t,input_red,input_blue,param);</p> | ||
+ | <p>figure (1)</p> | ||
+ | <p>subplot(221);plot(t,[x(:,19)],'g');xlabel('t');ylabel('alpha[···]');grid on;hold on</p></font> | ||
+ | <p>Here we plot the results. All time axis are expressed in minutes.</p> | ||
+ | |||
+ | |||
+ | <p>Fixing our attention in alpha:</p> | ||
+ | |||
+ | <p>It is clear that production decays with the increase of the pause between inputs. In other words, production is increased with more irradiation. Then, combinations of the same color repeated (red_red and blue_blue), will provide better results if they are one single pulse instead of two.</p> | ||
<header id="scrollsect2" class="major"> | <header id="scrollsect2" class="major"> | ||
<h3 style="text-align:left">Inputs combination: Sequences<br /> | <h3 style="text-align:left">Inputs combination: Sequences<br /> |
Revision as of 11:29, 18 September 2015
Before studying saturation when light pulses were given, we must take in advice that first level (producing A, B, C and D), starts expressing from minute 0. Activation of next level is directly proportional to the amount of those proteins that were produced firstly, as they are binding domains of secondly activated constructions. Thus, more expression will be obtained as more binding domains are synthetized. Inputs used in our circuit are determined by the value of two variables: color and time. Distinguish different colors is easy due to their different wavelength, resulting in the stimulation of different photoreceptors that regulate each switch. But how can our genetic constructions assign an order to inputs of the same colors? Recombinases are the biological elements which grant memory to AladDNA. These enzymes work as excisionases and integrases, cutting the gen which is surrounded by recognition sites specific for those recombinases. In our model this excision has been represented by the change from g+ to g-, referring to the capability of the gen for being translated. In our device, recombinases used are BxB1 and θC31. They are produced with the first input: Both second levels have the recognition sites for recombinases produced in the analogous level. Thus, enzymes firstly produced recombinate the option which has not been activated. When the second input is introduced, production in the second level will not be possible, as genetic codes have been flipped. In other words, if the second light is a different color, it will only mean production regulated by switches of the final level, resulting in a circuit that reminds which was the first input. As we explain in the biochemical reactions section, when recombinases are produced they instantaneously dimerize, and two dimmers will be needed to recombinate each gene. The duration of this process will determine the moment when it will be save to give the second input. Firstly, we will plot the results of a single red input in order to appreciate the curve described by g+ in the analogous level. Gene copy number used for this performance is of a single copy. t_final=3100; % with input_t=(0:Ts:t_final)' (minutes) input_red =[zeros(251,1);red_intensity*ones(length (input_t)-351,1) ;zeros(100,1)]; input_blue=[zeros(251,1); zeros(length (input_t)-351,1) ;zeros(100,1)]; On the other hand, the moment when recombinases switch off gH+, gI+ and gJ+, is approx. the minute 280, as we can see in the graphic below: In order to prove recombinases efficiency, we will simulate two followed inputs of different colors. If they work, the expected result is that there is no expression of that genes later, when the other light is given. t_final=3100; % with input_t=(0:Ts:t_final)' (minutes) t_recomb=280; %time needed for recombinases work input_red =[zeros(251,1); red_intensity*ones(t_recomb),1); zeros(length(intput_t)-t_recomb-251,1)]; As we guessed, no expression of the second level is produced after a previous input of 280 minutes, which is the time that ensure us the end of recombinases action. Different vectors can be used to introduce our genetic constructions in the cell. The amount of gene’s copies varies in each kind of plasmid, and this is a critical parameter. There will be a trade-off between the amount of genes expressing enzymes, and the number of genes that those will have to recombinate. We performed simulations with different values of gene copy numbers relying on some vectors typically used to transform E.coli. Values considered for this parameter were: gin=[1,5,10,17,400,600]; %Initial amount of genes From 1 to 17 copies: With more genes, there is more basal expression of recombinases, but also of those genes that BxB1 and θc31 must recombinate. Moreover, levels of basal expression of the analogous recombinases will also increase. This is the reson of the decay in gE+, which is more abrupted as the gene copy numer rises. That basal recombination is undesired, because the amount of binding domain produced also decays, as we can appreciate in the second group of curves (E). However, with 5 copies the difference is positive, but this arise does not persist when the number of copies keeps on increasing. Graphic above suggests that recombinases efficiencyis not affected by the numer of copies, because they can get the same result (g+ 0) in the same time. Then, if increasing the number of copies means less amount of expression and it makes no difference in recombinases efficiency, these results suggest that low numbers of copies are more suitable to get the desired behavior. Results with 400 and 600 copies support that higher amounts of copies induce more basal expression of recombinases, which leads in an undesired behaviour of our circuit: No protein production is appreciated. Looking closer to the first 100 minutes of simulation: Basal expression of both recombinases is so high that even when light is still off, they are produced in enough quantities to totally switch off production in the second level. Hello world We knew from literature, that changes induced by light are instantaneous. However, there are several differences between switching on and off these switches. For example, if we had implemented a red toggle switch, we would be able to ensure the stop of production using far red light, whereas dark periods induce a slower decay. This means that although light irradiation is stopped, the presence of photosensible elements may will decrease slowly, and production will not occur only when light is given. In other words: production will not be totally over light control. In order to ensure this coordination between light and presence of active elements, we suggested our wet lab mates tagging PhyB and ePDZ, they may will be degraded fast enough to minimize their presence when light has been eliminated. This is the reason why B and D have higher degradation rates than other proteins. So, if dark periods mean loss of activators, pauses between first and second inputs should be as short as they can be, in order to maximize production. The shorter pause possible, is a no-pause, i.e continuous irradiation. Firstly, we will compare production of alpha with different pauses between both red inputs, in order to analyze the influence of the dark period length. We implemented in Matlab a for loop to facilitate the comparison: pause_dur=[0,2,30,60,80,100]; % vector containing different duration of pauses. for i=1:length(pause_dur)
light_dur=length(input_t)-351-pause_dur(i);% duration of irradiation input_red=[zeros(251,1);red_intensity*ones(light_dur/2,1); zeros(pause_dur(i),1);red_intensity*ones(light_dur/2,1); zeros(100,1)]; input_blue=[zeros(3101,1)]; [t,x] = ode23s(@model_switch_v2, time_span, ic, opt,input_t,input_red,input_blue,param); figure (1) subplot(221);plot(t,[x(:,19)],'g');xlabel('t');ylabel('alpha[···]');grid on;hold on Here we plot the results. All time axis are expressed in minutes. Fixing our attention in alpha: It is clear that production decays with the increase of the pause between inputs. In other words, production is increased with more irradiation. Then, combinations of the same color repeated (red_red and blue_blue), will provide better results if they are one single pulse instead of two.Saturation
Figure 1.
Figure 2.
Figure 3.
Figure 4.
Figure 5.
Figure 6.
Figure 7.
Figure 8.
Recombinases action
How do they behave?
Figure 8.
Figure 8.
Figure 8.
Figure 8.
Figure 8.
Figure 8.
Figure 8.
Influence of gene copy number
Figure 8.
Figure 8.
Figure 8.
Figure 8.
Figure 8.
Light stimuli
Pause between pulses
Inputs combination: Sequences
Inputs choice for each output