Team:TU Munich/Modeling

From 2010.igem.org

Revision as of 20:51, 26 October 2010 by Sgude (Talk | contribs)

Navigation:

Home →  Modeling

iGEM MainPage

Contents

Overview

We simulated the termination and antitermination properties of our signal-terminator constructs with the Kinefold web server[1] and used some standard estimations for diffusive terms. Our main goal was to prove that our constructs work and that termination is stopped efficiently, that is, that the signal molecule binds and antiterminations occurs before the RNA polymerase falls off.
The Kinefold webserver provides a web interface for stochastic folding simulations of nucleic acids and offers the choice of renaturation or co-transcriptional folding. The folding paths are simulated at the level of helix formation and dissociation as these stochastic formation and the removal of individual helices are known to be the limiting steps of RNA folding kinetics.
For our purposes, co-transcriptional folding was the appropriate choice: Folding proceeds while the sequence is being synthesized from its 5' to 3' ends at a speed of 3 ms per newly added base (for RNA polymerase T7 phage). Thus, the transcript starts to fold before the whole sequence is fully available.
Kinefold offers the possibility to include additional bases (X) which do not pair to model hybridization dynamics between two sequences. In order to simulate how the binding of the signal molecule prevents termination we linked the signal via a linker sequence consisting of 'X' bases to the sequence of the terminator.
For each signal-terminator pair we did batch simulations with various random seeds in order to guarantee accuracy.
We also varied signal length form two base pairs to full signal length which provides insight in how long the signal needs to be in order to bind to the terminator and how long this process takes at least.

Results

Diffusion

The question whether antitermination occurs is not only guided by the folding process of the signal-terminator pair, but also by how long the signal takes to diffuse to the terminator sequence.

Read more

To account for the diffusion time, we estimated the hit rate τ (following 6.), which is the time until the signal meets the terminator sequence for the first time:

TU Munich iGEM2010 sstraub tau.png


where D is the diffusion constant, a the radius of gyration of the signal molecule and r the radius of the cell.
For E. coli r is 1 μm. The radius of gyration a can be estimated using the worm-like-chain model by

TU Munich iGEM2010 sstraub a.png


where n is the length of the signal which is 0,3 nm/monomer, l is the persistency length which is following (5.) 2nm for single-stranded RNA. Thus, for a signal of length 32 nt, a = 6,4 nm.

The diffusion constant D was obtained by

TU Munich iGEM2010 sstraub D.png


where kB is the Boltzmann constant and T is the absolute temperature.
Using the a we calculated above, we get D = 3,4318 m2/s.

Thus, for a cell containing 100 signal molecules, the signal needs 0,1518 sec until it first hits the terminator sequence.

Close

As the folding time (at least 1,7 sec for full signal length) is significantly larger than the estimated diffusion time (0,1518 sec), diffusion does not play an important role and thus can be neglected in the modeling our signal-terminator constructs.

Switch

His-terminator

Terminator sequence:
5' UCGGCUUCAACGUGCUCCACGAAAGCCCCCGGAAGAUGCAUCUUCCGGGGGCUUUUUUUUU 3'
Signal sequence:
5' CCGGGGGCUUUCGUGGAGCACGUUGAAGCCGA 3'

Both terminator and signal sequence comprise two parts:

  • His-terminator and a part which partially binds to the His-terminator respectively (marked green),
  • random signaling part and antisense to the random signaling part respectively (marked red).

We modeled the signal-terminator interaction for 2, 4-32 nt signal length and proved that the signal sequence binds to the terminator sequence and impedes the terminator folding, thus the polymerase doesn't fall off and the output signal is produced.

The minimal free-energy structure of the folded terminator including the random signaling part is shown in the figure "His-terminator" below:

His-terminator

Our full signal sequence of length 32 nt completly binds to the terminator and antitermination is successful:

His-terminator plus signal

On the other hand the picture for a signal length of only 2 nt shows that such a small signal length is not sufficient for antitermination. But already a signal of length 16 nt binds to the terminator so that it cannot fold completely anymore.

signal length 2 nt
signal length 16 nt

As we observed a dependence of the folding time and minimal energy respectively from the linker length we used linkers of length 12, 16, 20, 24, 30, 34, 40, 50, 60, 120 'X' bases. For the minimal energy there is only a very small linker length dependence as one can see in the figure below. The figure also shows that the minimal free-energy structure for the signal sequence bound to the terminator sequence is much more stable, around 35 kcal/mol, than the minimal free-energy structure for the terminator sequence alone.

Linker length dependence: minimal energy

However, for the folding time linker length dependence is significant as one can see in the figure below.

Linker length dependence: folding time

But, as, for us, the fact whether the terminator is folding or not is crucial and the videos below show that for appropriate signal length the terminator is not folding. Concerning this fact, for all linker lengths we tested there was only a dependence on the signal length but the results were independent from the linker length. Hence, the linker length dependence can be neglected.

RNA folding path videos

The video on the left for 16 nt signal length shows that the terminator is already hindered from folding completely and the video on the right, which is done for full signal length of 32 nt, displays that the terminator does not fold at all as the signal immediately binds to the terminator sequence as it is synthesized.

TU Munich iGEM2010 sstraub 16.gif TU Munich iGEM2010 sstraub 32.gif

Optimization of the signal sequence

We also tried to optimize our signal sequence to find the minimal sequence length which still enables antitermination. For this optimization kept the random signaling part constant and varied the part which partially binds to the His-terminator and vice versa to see which length is sufficient for antitermination for each case. Then we combined the results to find a optimal signal.

When varying the part which partially binds to the His-terminator we found that having only the 2 nt of this part together with the random signaling part suffices for antitermination as the picture and the corresponding video below show.

Signal consists of the random signaling part (20 nt) and only 2 nt of the part which partially binds to the His-terminator.
TU Munich iGEM2010 sstraub ter2.gif

For the full 12 nt of the part which partially binds to the His-terminator, the random signaling part has to have at least 2 nt in order to prevent termination as one can see in the picture below.

Signal consists of the part which partially binds to the His-terminator (12 nt) and only 2 nt of the random signaling part.
TU Munich iGEM2010 sstraub ran2.gif

We then combined these two results and found that a signal with overall 8 nt, 4 nt for the random signaling part and 4 nt for the part which partially bind to the His-terminator, 5' UUUCGUGG 3', suffice to enable antitermination as the picture and the corresponding video below display.

Signal consists of 4 nt for the part which partially binds to the His-terminator and 4 nt of the random signaling part.
TU Munich iGEM2010 sstraub signalopt.gif

Trp-terminator

Terminator sequence:
5' ACGCGACAUAUCAGCUUAGAAUACCCAGCCCGCCUAAUGAGCGGGCUUUUUUUU 3'
Signal sequence:
5' GGCGGGCUGGGUAUUCUAAGCUGAUAUGUCGCGU 3'

Both terminator and signal sequence comprise three parts:

  • Trp-terminator and a part which partially binds to the Trp-terminator respectively (marked green),
  • a natural sequence part upstream of the Trp-terminator binding to the terminator in antitermination state and a part which binds to the upstream sequence of the terminator respectively (marked yellow),
  • random signaling part and antisense to the random signaling part respectively (marked red).


We did similar simulations as we did for the His-terminator for the Trp-terminator and it turned out to be less stable than the His-terminator. The folded terminator is not the minimal free-energy structure for the given the RNA sequence, it's only the third minimal one at -12.9 kcal/mol compared with the His-terminator which has a minimal free-energy structure at -31.3 kcal/mol. The folded terminator including the natural sequence part upstream of the Trp-terminator binding to the terminator in antitermination state and the random signaling part is shown in the figure "Trp-terminator" below:

Trp-terminator

Our full signal sequence of length 34 nt completly binds to the terminator and antitermination is successful:

Trp-terminator plus signal
TU Munich iGEM2010 sstraub trp.gif

Optimization of the signal sequence

As for the His-terminator, we tried to optimize our signal sequence to find the minimal sequence length which still enables antitermination.

It turned out that a signal with overall 10 nt, 2 nt for the random signaling part, all 4 nt for the part which binds to the upstream sequence of the terminator and 2 nt for the part which partially bind to the Trp-terminator, 5' CUGGGUAUUC 3', suffice to enable antitermination as the picture and the corresponding video below display.

Minimal signal length that still enables antitermination.
TU Munich iGEM2010 sstraub trpopt.gif


T7-terminator

Terminator sequence:
5' CUCUGCCGCACGUAU AGCAUAACCCCUUGGGGCCUCUAAACGGGUCUUGAGGGGUUUUUUG 3'
Signal sequence:
5' AAGGGGUUAUGCUAUACGUGCGGCAGAG 3'

Both terminator and signal sequence comprise two parts:

  • T7-terminator and a part which partially binds to the T7-terminator respectively (marked green),
  • random signaling part and antisense to the random signaling part respectively (marked red).


We also modeled the T7-terminator which is also less stable than the His-terminator. Its minimal free-energy structure is at -23.1 kcal/mol. The folded terminator including the random signaling part is shown in the figure "T7-terminator" below:

T7-terminator

Our full signal sequence of length 28 nt completly binds to the terminator and anti-termination is successful:

T7-terminator plus signal
TU Munich iGEM2010 sstraub t7.gif

Optimization of the signal sequence

As for the His-terminator, we tried to optimize our signal sequence to find the minimal sequence length which still enables antitermination.

We finaly found that a signal with overall 8 nt, 3 nt for the random signaling part and 5 nt for the part which partially bind to the T7-terminator, 5' AUGCUAUA 3', suffice to enable anti-termination as the picture and the corresponding video below display.

Minimal signal length that still enables anti-termination.
TU Munich iGEM2010 sstraub t7opt.gif


Network

Outlook

In future our switch could be used to interconnect BioBricks as demonstrated on our [TU Munich/Software| Software] page.

Furthermore we have discovered that the currently available RNA folding programs are not fully suitable for describing complex RNA interactions between multiple RNA sequences. This obstacle has to be overcome to allow reliable in silcio modeling of RNA devices in the future.

References

[1] A. Xayaphoummine, T. Bucher and H. Isambert, Kinefold web server for RNA/DNA folding path and structure prediction including pseudoknots and knots, 2005, Nucleic Acids Research, Vol. 33, Web Server issue W605–W610. [2] A. Xayaphoummine, T. Bucher, F. Thalmann, and H. Isambert, Prediction and statistics of pseudoknots in RNA structures using exactly clustered stochastic simulations, December 23, 2003, 15310–15315 PNAS vol. 100 no. 26. [3] Christoph Flamm, Walter Fontana, Ivo L. Hofacker, and Peter Schuster, RNA folding at elementary step resolution, RNA (2000), 6:325–338 Cambridge University Press. [4] David H. Mathews, Jeffrey Sabina, Michael Zuker and Douglas H. Turner, Expanded Sequence Dependence of Thermodynamic Parameters Improves Prediction of RNA Secondary Structure, 1999, J. Mol. Biol. (1999) 288, 911-940. [5] Rippe, Making contacts on a nucleic acid polymer, 2001, TRENDS in Biochemical Sciences. [6] Berg, v. Hippel, 1985.