Team:TU Munich/Modeling

From 2010.igem.org

Revision as of 16:40, 20 November 2010 by Sstraub (Talk | contribs)

Navigation:

Home →  Modeling

iGEM MainPage

Contents

Overview

We probed the static and the dynamical properties our transmitter-switch (also signal-terminator) constructs with two subsequent approaches: 1) finding a terminator and a corresponding transmitter using thermodynamically NUPACK simulations 2) investigating kinetic parameters regarding diffusion and transmitter-switch binding dynamics.

As described on the project page, one key aspect of our switches is the idea, that a RNA transmitter molecule is capable to shift the state of a switch only if its trigger site is present and its identity site corresponds to the recognition site of the switch.
We simulated the termination and antitermination properties of our transmitter-switch constructs with the NUPACK and the kinetically parameters with [http://kinefold.curie.fr/| Kinefold web server] 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 transmitter molecule binds and antitermination occurs before the RNA polymerase falls off.

Results

We successfully constructed several switches and their corresponding transmitter RNA in silico on a thermodynamical basis. We modified different transcriptional terminators in such a way, that the formation of the terminator was prevented by a transmitter molecule. As desired, this only occurred if the transmitter molecule contained both, a trigger and an identity site. Analogously, we were able to design and verify a NOT gate using the same thermodynamical approach.
Our estimations for the switches obtained by the in silico approach indicated that diffusion is negligible. Our RNA dynamics simulations showed for each switch that for appropriate transmitter length the terminator is not folding, hence antitermination is successful. We also found that much smaller transmitter length is suitable for antitermination. Finally, we showed using the example of an AND gate that logic gates based on our transmitter-switch constructs work well.

in silico design based on thermodynamic calculations

In the following section we will describe how we theoretically designed our switches and their respective transmitter molecules. Our goal was to develop synthetic modified terminators consisting of a target and a recognition site in the described manner (compare Implementation). As explained, only the complete transmitter RNA may shift the switches´ state, but neither the trigger site nor the identity site alone. This is the prerequisite to fulfill all mentioned requirements for a functional switch allowing to construct logical gates as discussed previously. We implemented this by optimizing the thermodynamically parameters for pairing of each subsegment (identity site, trigger site, recognition site, target site).

Read more

For all NUPACK simulations, following parameters were used if not mentioned explicitly. A temperature of 37°C was used in combination with the default Serra and Turner (1995) energy parameters.Default 1.0 M Na+ salt concentrations were used. For simulating several strand species, 3 nM switch concentration and 50 nM transmitter concentration were used, which equals conditions in a e.coli cell with a low copy number plasmid coding the switch and accumulation of transmitter RNA.

switches based on attenuation principle

On a first approach, we designed switches based on terminators performing antitermination in the context of their natural occurrence within the scope of attenuation, the His and Trp terminator.[1][2]. A 20 bp random sequence was derived from [http://www.faculty.ucr.edu/~mmaduro/random.htm| random sequence generator] and put upstream the terminator loop. A complementary sequence, reaching within the terminator´s stem loop (the transmitter) was stepwise shortened to find the length, where the formation of the terminator is thermodynamically favored compared to the strand displacement by this sequence. The final transmitter sequence was then defined by selecting this shortest strand which is still able to "destroy" the stem loop. Subsequently, the complementary random sequence (identity site) was removed an the remaining trigger site was tested in regard of not being able to resolve the stem loop on its on. As illustrated for the His-based switch below, the terminator is thermodynamically favored toward the trigger unit, but in combination with the identity site, destroying the terminator´s stem loop becomes possible.


element free energy [kcal/mol] ratio of pairing [%]
switch (terminator secondary structure) -30.10 -
transmitter + switch -72.88 100
trigger site binding to single target site -24.21 0
identity site binding to single recognition site -38.78 0
switch based on His terminator[1], using a 20 bp recognition site: Comparison of free energy for secondary structure formation at 37°C in the presence and absence of a respective transmitter, simulated by using NUPACK analysis.



Pairing probabilities of the His-Terminator based switch simulated by NUPACK secondary structure analysis tool
Illustration of the His-Terminator based switch secondary structure, generated by NUPACK utility tools
Pairing probabilities of the His-Terminator based switch in presence of the respective transmitter RNA, simulated by NUPACK secondary structure analysis tool
Illustration of the His-Terminator based switch secondary structure in presence of the respective transmitter RNA, generated by NUPACK utility tools
Pairing probabilities of the His-Terminator based switch in presence of the transmitter RNA´s trigger site only, simulated by NUPACK secondary structure analysis tool
Illustration of the His-Terminator based switch secondary structure in presence of the transmitter RNA´s trigger site only, generated by NUPACK utility tools
Pairing probabilities of the His-Terminator based NOT-gate in presence of the transmitter RNA´s identity site only, simulated by NUPACK secondary structure analysis tool
Illustration of the His-Terminator based NOT-gate secondary structure in presence of the transmitter RNA´s identity site only, generated by NUPACK utility tools



NOT Gate based on His terminator switch
The basic composition of the NOT-gate was described in the implementation section. Putting this principle into practice based on the His terminator switch by further NUPACK simulation was successful. The terminator is only built in presence of a transmitter molecule and the trigger or identity site are not able to shift this state on their on. Results are presented in the following:

element free energy [kcal/mol] ratio of paaring [%]
switch (NOT-gate secondary structure) -68.18 -
transmitter + switch -98,18 100
trigger unit + target site -24.01 0
identity site + recognition site -38.78 0
NOT-gate based on modified His terminator[1] based switch:Comparison of free energy for secondary structure formation at 37°C, simulated by using NUPACK analysis.


Pairing probabilities of the His-Terminator based NOT-gate simulated by NUPACK secondary structure analysis tool
Illustration of the His-Terminator based NOT-gate secondary structure, generated by NUPACK utility tools
Pairing probabilities of the His-Terminator based NOT-gate in presence of the respective transmitter RNA, simulated by NUPACK secondary structure analysis tool
Illustration of the His-Terminator based NOT-gate secondary structure in presence of the respective transmitter RNA, generated by NUPACK utility tools
Pairing probabilities of the His-Terminator based NOT-gate in presence of the transmitter RNA´s trigger site only, simulated by NUPACK secondary structure analysis tool
Illustration of the His-Terminator based NOT-gate secondary structure in presence of the transmitter RNA´s trigger site only, generated by NUPACK utility tools
Pairing probabilities of the His-Terminator based NOT-gate in presence of the transmitter RNA´s identity site only, simulated by NUPACK secondary structure analysis tool
Illustration of the His-Terminator based NOT-gate secondary structure in presence of the transmitter RNA´s identity site only, generated by NUPACK utility tools


Furthermore, we successfully designed a switch based on the Trp terminator[2], using the same approach as pointed out above. In the following, energies for pairing only trigger site and target site / identity site and recognition site are not mentioned explicitly. Nevertheless, for each switch it was proofed that trigger or identity site on its one are not able to shift the switches´ state, respectively.


element free energy [kcal/mol] ratio of pairing [%]
switch (terminator secondary structure) -73,88 -
transmitter + switch -73,88 100
switch based on Trp terminator[2], using a 20 bp recognition site:Comparison of free energy for secondary structure formation at 37°C in the presence and absence of a respective transmitter, simulated by using NUPACK analysis.


switches using ubiquitous terminators

We applied the same principle to terminators which do not exhibit a natural occurring antitermination and were able to establish switches based on the synthetic terminator BBa_B1006 [3] and the T500 terminators[4]. We could show the length of the identity site can be varied. Nevertheless, How this will influence kinetic parameters discussed in the Kinetic simulations using Kinefold section

element free energy [kcal/mol] ratio of pairing [%]
switch (terminator secondary structure) -19.70 -
transmitter + switch -70.98 100
switch based on BBa_B1006 synthetic terminator[3], using a 20 bp recognition site:Comparison of free energy for secondary structure formation at 37°C in the presence and absence of a respective transmitter, simulated by using NUPACK analysis.


element free energy [kcal/mol] ratio of pairing [%]
switch (terminator secondary structure) -29.80 -
transmitter + switch -81.08 100
switch based on T500 terminator[4], using a 30 bp recognition site:Comparison of free energy for secondary structure formation at 37°C in the presence and absence of a respective transmitter, simulated by using NUPACK analysis.
element free energy [kcal/mol] ratio of pairing [%]
switch (terminator secondary structure) -29.80 -
transmitter + switch -81.08 100
switch based on T500 terminator[4], using a 10 bp recognition site:Comparison of free energy for secondary structure formation at 37°C in the presence and absence of a respective transmitter, simulated by using NUPACK analysis.

switches relying on tiny abortive RNAs

We also implemented switches based on a the phenomenon of tiny abortive RNAs[5]. This principle is in contradiction to NUPACK-simulations, suggesting the T7-terminator to be less stable than simulated in NUPACK. Nevertheless, we tried to implement switches integrating tiny abortive RNAs in different length as trigger site. Results are illustrated in the in vitro Screening section.

Close

Kinetic calculations using Kinefold

In order to prove if the antitermination process is not only possible on a thermodynamically base, but also regarding kinetics, we performed the following modeling.

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. The folding of the signal-terminator pair has to be kinetically on the scale of the folding of the terminator structure on its on, to be a valid competitive reaction. This means, the antitermination has to occur before the polymerase will fall of the DNA strand.

Read more

To account for the diffusion time, we estimated the hit rate τ[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 compartment.
For E. coli r is approximately 1 μm. The radius of gyration a can be estimated using the worm-like-chain model[7b 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[8] which is 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 and T = 300 K, we get D = 3,4318 m2/s.

Thus, for a cell containing 100 signal molecules, the signal needs approximately 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 modeling our signal-terminator constructs. Thus the folding process is the rate limiting factor.

Shifting the switch

To understand the RNA-folding dynamics of our switches, we performed Kinefold simulations for the His-terminator, the Trp-terminator and for the T7-terminator, each with the corresponding signal. For each terminator-signal construct, there are folding path videos available below and we optimized the signal sequences.

The Kinefold web server[9] 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 transmitter molecule prevents termination we linked the transmitter via a linker sequence consisting of 'X' bases to the sequence of the switch.
For each transmitter-switch pair we did batch simulations with various random seeds in order to guarantee accuracy.
We also varied transmitter length from two base pairs to full transmitter length which provides insight in how long the transmitter needs to be in order to bind to the switch and how long this process takes at least.
We also examined whether diffusion plays a role in the dynamics of our constructs or not.


We checked that our switch works well for all three terminators.

I. His-terminator

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

Both terminator and signal sequence comprise two parts:

  • His-terminator, containing the target site, and a part which partially binds to the His-terminator, called trigger site, respectively (marked green),
  • a random signaling part, also called identity site, and an antisense to the random signaling part, also called recognition site, respectively (marked red).


For further details on the switch design see the Project page.

We modeled the signal-terminator interaction for 2, 4-32 nt signal length and find that the signal sequence binds to the terminator sequence and impedes the terminator folding, on a timescale, that the polymerase doesn't fall off and the output signal can be produced.

Read more

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 completely binds to the terminator leading to antitermination.

His-terminator plus signal 32 nt

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

signal length 2 nt
signal length 16 nt

Furthermore, we explored whether the artificial linker, introduced to allow simulations in kinefold, which can only model intramolecular folding, has any influence on the results.

Read more

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 bases, indicated as 'X'. 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 which -31,3 kcal/mol (and thus close to the minimal free-energy for the terminator connected to a signal of length only 2 nt).

Linker length dependence: minimal energy

However, for the folding time the linker length dependence is significant as one can see in the figure below, but only because the production of this linker is also included in the calculations. However, since our system uses intermolecular RNA/RNA interaction, the time for the linker production can be neglected.

Linker length dependence: folding time

Read more

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

Finally we optimized the signal sequence and obtained a signal sequence of only 8 nt length.

Read more

Optimization of the signal sequence

We tried to optimize our signal sequence to find the minimal sequence trigger sequence length, which still enables antitermination. For this optimization we kept the random signaling part constant and varied the part which partially binds to the His-terminator to see which length is sufficient for antitermination for each case. Subsequently, we optimized the identity site lenght using the detected optimal trigger site lenght. Then we combined the results to find a optimal overall signal lenght.

When varying the trigger site, 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. This is in contradiction to the NUPACK results, and is probably due to the fact Kinefold takes RNA-production into account. Thus, the identity site can bind before any part of the terminator it self is produced. It seems like the kinetically slow folding of the terminator can not compete with the signal, which is already bound to the nascent switch strand. If this simulation depicts the real situation, than our trigger site could be even modified and would be much more efficient as expected by only thermodynamically considerations.

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

In contrast, using a 12 nt trigger site only requires 2 nt identity site 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: Above one always had the full length of either one of the two parts of the signal, as we now wanted to modify both parts and only 2 nt of each part were not sufficient for antitermination, we finally discovered that a signal with overall 8 nt, 4 nt for the random signaling part and 4 nt for the part which partially binds to the His-terminator, 5' UUUCGUGG 3', suffices to enable antitermination as the picture and the corresponding video below display.
This highlights how strongly the length of identity site and trigger site are both responsible for antitermination, since both of them influence the overall free energy of the competitive reaction to termination. Thus, this proves, that a meaningful selection and optimization of trigger and identity site allows antitermination by RNA/RNA interaction and performing a specific addressing of single switches by only varying the identity site. This manifests the fact, that our idea could be implemented and the specific addressing and switching could be separated, as for example showed for the His-switch in the NUPACK-section.

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

Close

Close

Close

Close

Results

We showed that in theoretical simulations the His-terminator works perfectly as the terminator itself is relatively stable whereas for appropriate signal length the terminator is not folding at all, that is antitermination occurs. We also demonstrated that a quite short signal length of 8 nt is sufficient for antitermination to happen. Nevertheless longer signals provide a larger number of implementable switches. 8 nt is only the minimal length for our switch to function, and the simulations showed that longer identity sites do marginally increase the antitermination process. This is because for RNA/RNA interaction, the rate limiting step is the binding of the first 3-5 bases.

II. Trp-terminator

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

Both terminator and signal sequence comprise two parts:

  • Target site (marked green),
  • recognition site (marked red).

We modeled the Trp-terminator-signal interaction, generated folding path videos that give more insight in its dynamics and optimized the signal sequence.


Read more

The Trp-terminator 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, 8 nt for the random signaling part, and 2 nt for the part which partially binds 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

Close

Results

We showed that the Trp-terminator should function properly. Although the terminator itself is less stable than the His-terminator for appropriate signal length the terminator is not folding, thus antitermination occurs. We also demonstrated that a signal length of 10 nt is sufficient for successful antitermination.

III. 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),
  • a random signaling part and an antisense to the random signaling part respectively (marked red).


There are folding path videos available for the T7-terminator as well as an optimization of the signal sequence.

Read more

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 finally found that a signal with overall 8 nt, 3 nt for the random signaling part and 5 nt for the part which partially binds 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

Close

Results

We demonstrated that in theoretical simulations the T7-terminator work well. It's less stable than the His-terminator (and rather a lot of bases are unpaired), but its more stable than the Trp-terminator. For appropriate signal length the terminator is not folding at all, that is antitermination occurs. We also demonstrated that a signal length of 8 nt is sufficient for antitermination to happen.

Logic gate

After establishing that our switch works as designed, we combine two switches to build an AND gate an investigate it with further simulations. We used the His-terminator, the Trp-terminator and two corresponding signals with different random signaling parts. Having different signaling parts is necessary, as if one uses the same one in both signals, the AND gate will be transformed into a OR gate and will have the output 'true' even in the case that only one signal is present. We linked all the sequences with a 34 X linker. As above both terminator and signal sequences comprise two and three parts respectively:

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

His-terminator sequence 1:
UCGGCUUCAACGUGCUCCACGAAAGCCCCCGGAAGAUGCAUCUUCCGGGGGCUUUUUUUUU 3'
Signal sequence 1:
5' CCGGGGGCUUUCGUGGAGCACGUUGAAGCCGA 3'
Trp-terminator sequence 2:
5' ACGCGACAUAUCAGCUUAGAAUACCCAGCCCGCCUAAUGAGCGGGCUUUUUUUU 3'
Signal sequence 2:
5' GGCGGGCUGGGUAUUCUAAGCUGAUAUGUCGCGU 3'

We simulated the binding kinetics of the AND gate, created folding path videos which nicely show how the gate is working and applied the results from our signal sequence optimization.

Read more

First we considered the case that only one signal sequence is present, that is at most one terminator gets destroyed, which one can see in the picture below.

AND gate, only one signal present.

As the picture below displays both of the signals successfully bound to the terminators, thus as both input signals are present (that it both inputs are '1') the AND gate also returns the value '1' as both terminators are destroyed.

AND gate, both signals present.

We also tried to model an AND gate with the optimized signals which were sufficient to destroy a single terminator above, but didn't work out for the AND gate as the picture below shows. The signal sequences were to short to bind to the terminator probably due to hindering by other structures, e.g. by the linker.

AND gate and the minimal signal length

RNA folding path videos

We created folding path videos for the AND gate. For only one signal present the video shows how the Trp-terminator is not destroyed.

TU Munich iGEM2010 sstraub AND only one signal.gif

On the other hand, when both signals are present, one can see how both signal sequences bind to the terminator sequences as those are synthesized and none of the terminators can fold.

TU Munich iGEM2010 sstraub AND.gif

Close

Results

We were able to model a logic gate consisting of two of our switches we introduced above. We proved this logic AND gate to work well, that is if the two input sequences are present the terminators are destroyed and if only one or none signal is present termination occurs.

Outlook

So far there is no RNA folding program that is specialized on modeling the interaction of two sequences not to speak of a whole network consisting of RNA based switches and signaling sequences. Combining the concepts of Kinefold[10] and Vienna[11] and using sets of secondary structures that are adjusted to having two sequences instead of only one sequence one could get rid of the bulky linkers. The linkers we used did not play a significant role for our simulation as we were just handling small RNA devices, but for larger networks they are likely to be sterically hindering. This obstacle has to be overcome to allow reliable in silcio modeling of RNA devices.

With appropriate software RNA devices like our switches could in future be used to interconnect BioBricks, as demonstrated on our Software page, and to simulate those interconnections at elementary step resolution.

References

[1] Gollnick P, Babitzke P.,Biochim Biophys Acta., Sep 2002 13;1577(2):240-50 pp. 901-1062 [2] Yanofski et al., Nucleic Acids Research, 1981, Volume 9 Number 24 [3] http://partsregistry.org/Part:BBa_B1006:Design [4] Larson et al., March 2008, Cell, Volume 132, Issue 6 [5] Lee, S., H.M. Nguyen, and C. Kang, Tiny abortive initiation transcripts exert antitermination activity on an RNA hairpin-dependent intrinsic terminator. Nucleic Acids Research. [6] Berg, v. Hippel, 1985. [7] Doi and Edwards, 1999. The Theory of Polymer Dynamics. [8] Rippe, Making contacts on a nucleic acid polymer, 2001, TRENDS in Biochemical Sciences. [9] 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. [10] 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. [11] Christoph Flamm, Walter Fontana, Ivo L. Hofacker, and Peter Schuster, RNA folding at elementary step resolution, RNA (2000), 6:325–338 Cambridge University Press.