Team:TUDelft/Temperature design2

From 2008.igem.org

(Difference between revisions)
(Secondary structure stability and switching temperature)
 
(138 intermediate revisions not shown)
Line 3: Line 3:
{{Template:TUDelftiGEM2008_menu_home}}
{{Template:TUDelftiGEM2008_menu_home}}
-
>> Work in progress
+
=Design phase II: <br \><br \>Changing the switching temperature of the RNA thermometer=
-
[http://partsregistry.org/wiki/index.php?title=Part:BBa_K115016 BBa:K115016], [http://partsregistry.org/wiki/index.php?title=Part:BBa_K115017 BBa:K115017],
+
Introducing existing RNA thermometers into ''Escherichia coli'' provides us with a number of temperature sensitive switches, all with their own switching temperature. But this only provides a few switching temperatures, which are also relatively close to each other. In order to build a biothermometer, we needed more switching temperatures that have a bigger range. Note, however, that there will always be an upper bound in that the ''E. coli'' cannot survive temperatures above 40°C, and a lower bound in that the activity of the cell will drop at low temperatures.
-
[http://partsregistry.org/wiki/index.php?title=Part:BBa_K115018 BBa:K115018],
+
-
[http://partsregistry.org/wiki/index.php?title=Part:BBa_K115019 BBa:K115019]
+
-
=Changing the temperature threshold of the RNA thermometer=
+
This chapter presents a method that has been used to design RNA thermometer parts that should theoretically switch at a predefined temperature. The method used is based on the data analysis of which the results can be found on the [[Team:TUDelft/Temperature_analysis|analysis page]]. It has been used to design two new parts, one that should theoretically be switched on at 27°C ([http://partsregistry.org/wiki/index.php?title=Part:BBa_K115016 <code>BBa_K115016</code>]), and two that should theoretically be switched on at 32°C ([http://partsregistry.org/wiki/index.php?title=Part:BBa_K115017 <code>BBa_K115017</code>]).
-
The principle of the RNA thermometer is based on base pairing between the nucleotides in the Shine Dalgarno region. At a certain temperature the RNA is folded into a structure in which base-pairing nucleotides, making the Shine Dalgarno region unreachable for the ribosome.  With a rise in temperature the binding forces between the base-pairing nucleotides decrease and above a certain threshold temperature the base-pairing forces are to weak to hold the base-pairing nucleotides together, causing the RNA to (partially) unfold. The binding between the base-pairing nucleotides let loose causing the RNA to unfold and exposing the Shine-Dalgarno region enabling the ribosome to initiate translation.
+
====Secondary structure stability and switching temperature====
-
When only looking at this principle an RNA thermometer with a different temperature threshold can be designed by simply increasing or decreasing the binding forces of the base-pairing nucleotides in the Shine-Dalgarno region, shifting the temperature threshold to a higher or lower temperature respectively. Increasing the binding forces could be achieved by incorporating base-pairing C and G nucleotides, which bind relatively strong forming a stable helix. Decreasing the binding forces could be done by introducing less strong binding nucleotide base-pairs, such as A-U and G-U base pairs.
+
[[Image:HairpinStability.png|thumb|200px|right|'''Figure 1:''' '''A.''' A lot of G-C base pairs make a very stable hairpin (having a very low free energy). '''B.''' The same hairpin structure with A-U base pairs is much less stable (having a higher free energy). '''C-E''' Bulge loops and internal loops destabilize the structure even further.]]
-
But as many times is the case with biology, things are not that easy. But despite the complexity, a search to some papers indicate that the principle described above is mainly causes the RNA thermometer to be temperature sensitive. This makes it worth while trying to design an artificial RNA thermometer that operates at a different temperature.
+
The switching temperature of the RNA thermometer depends on the stability of the hairpin structure that surrounds the Shine Dalgarno region (see also [[Team:TUDelft/Temperature_analysis#Requirements_for_an_RNA_thermometer|Requirements for an RNA therometer]]). This secondary structure will melt when it is exposed to a certain amount of heat. The more stable the hairpin structure is, the more heat will be needed to melt it, and the less stable the hairpin structure is, the less heat will be needed to melt it. In case of the RNA thermometer, this melting is needed to expose the Shine Dalgarno sequence which will cause the initiation of translation.
 +
In other words, a less stable hairpin structure will melt at a lower temperature than a more stable hairpin and thus, for the less stable hairpin, the translation will initiate at a lower temperature than the more stable one. This way we could theoretically shift the switching temperature to a lower or higher temperature by making the hairpin structure surrounding the SD region less or more stable respectively.
 +
The stability of the hairpin structure depends on the base-pairing nucleotides within the secondary structure. For example, C-G base pairs are very strong, so a hairpin structure containing a lot of those base pairs will be very stable. A-U and U-G base pairs are less strong, so these can be used to destabilize the hairpin structure. The same can be achieved by incorporating bulge loops and internal loops within the hairpin structure ([https://2008.igem.org/Image:HairpinStability.png figure 1]).
 +
====Constraints====
-
Two facts that enforce this assumption...
+
There are two things that have to be taken into account. Firstly, adding mutations in order to enforce or weaken the temperature sensitive region can also cause the RNA to fold into a completely different secondary structure. This way it can lose its function as an RNA thermometer. Secondly, there are some constraints to the possible mutations. The start codon and the Shine Dalgarno sequence must off course remain unaltered and in case of a standard biobrick the scar must also be part of the temperature sensitive hairpin.
-
First, aligning the sequences of the different types of RNA thermometers produces conserved regions around the SD region.
+
-
Second, it is shown by De Smit and Van Duin (reference) that there is a strong correlation between translation efficiency and the stability of the helix containing the Shine-Dalgarno region and the initiation codon.
+
-
But will this region on itself be enough to have a working RNA thermometer? Looking at the consensus structure of the ROSE and PrfA RNA thermometer, it can be seen that there are other regions within the structure that are highly conserved. This indicates that these regions might also be of importance to the functioning of the RNA thermometer. It is still unknown if the stem-loops not containing the SD are of any help with the opening of the SD region (reference).
+
=Design starting point=
-
Pseudoknot analysis.
+
[[Image:Start.png|thumb|200px|right|'''Figure 2:''' The secondary structure of the RNA thermometer that has been used as starting point of the design process is displayed at the left. At the right the stability profile of the temperature sensitive hairpin is highlighted (light blue). Also the found stability trend (see [[Team:TUDelft/Temperature_analysis#Data_analysis|data analysis]]) and the specified boundaries (see [[Team:TUDelft/Temperature_design#Design_algorithm|design algorithm]]) are highlighted (dark blue). The grey lines are the stability profiles of the other 31 ROSE RNA thermometers from the [http://rfam.sanger.ac.uk Rfam database].]]
-
figures...(consensus structures, highly conserved SD regions)
+
Since we have a number of structures of which we know that they work as an RNA thermometer, the most logical option is to take one of those structures as the starting point for the design of an RNA thermometer with a different switching temperature. The one chosen is a ROSE RNA thermometer from ''Bradyrhizobium japonicum USDA 110'' ([http://www.ncbi.nlm.nih.gov/sites/entrez?Db=genome&Cmd=ShowDetailView&TermToSearch=272 <code>NC_004463</code>]) residing at location [http://www.ncbi.nlm.nih.gov/projects/sviewer/?id=NC_004463.1&v=5784144..5784239 <code>5784144-5784239</code>] within the genome, which is at the 5' side of the gene [http://www.ncbi.nlm.nih.gov/sites/entrez?db=gene&cmd=search&term=NC_004463%20hspb hspB]. There was no specific reason for the choice of this ROSE RNA thermometer, this is just a ROSE RNA thermometer taken from the [http://rfam.sanger.ac.uk Rfam database].
-
De Smit and Van Duin<sup id="cite_ref-1"> [[Team:TUDelft/Temperature_design#cite_note-1 |[1]]] </sup> show a correlation between the level of expression and the free energy of the secondary structure in the Shine Dalgarno region.
+
We will call this RNA thermometer a 37°C switch, while it actually reaches maximum translation several degrees higher. The reason that we call it a 37°C switch is because of the characteristic stability profile we found during the [[Team:TUDelft/Temperature_analysis#Data_analysis|analysis]], and which is calculated at 37°C. Having a 27°C switch would then mean that this is an RNA thermometer of which the stability profile at 27 degrees is similar to the trend as given in [[Team:TUDelft/Temperature_analysis#Results|figure 8 of the analysis]]. This switch should then theoretically be almost fully switched on at 27°C, as is the case for the 37°C version at 37°C.
-
<div class="toggler">
+
[https://2008.igem.org/Image:Start.png Figure 2] shows the secondary structure of the chosen RNA thermometer as predicted by ''RNAfold''. Also the stability profile that is characteristic for this RNA thermometer at 37°C is shown in light blue. The thick dark blue line is the trend as found in the [[Team:TUDelft/Temperature_analysis#Data_analysis|data analysis]] and the thin dark blue lines are design boundaries which will be discussed later in this chapter. The stability profiles of the other 31 ROSE RNA thermometers from the Rfam database are shown in grey.
-
=== test ===
+
-
</div>
+
-
<div class="element">
+
=Design constraints=
-
test test test
+
-
</div>
+
-
<div class="toggler">
+
Having a starting point, we need a way to turn this RNA thermometer into one that switches at a different temperature, and that conforms to the Standard Biological Part standard. This can be achieved by increasing or decreasing the stability of the hairpin structure, as explained at the begin of this page. But, as also mentioned, there are certain constraints that has to be laid upon the design.
-
=== test ===
+
-
</div>
+
-
<div class="element">
+
====Constraints on the sequence====
-
test test test
+
-
</div>
+
-
<div class="toggler">
+
Looking at the nucleotide sequence of the design, there is the constraint that there must be a start codon, a Shine Dalgarno sequence, and a scar present within the sequence. Note that this is the sequence after ligation of the RNA thermometer part to a protein coding part and that the scar and start codon are not part of the RNA thermometer part that we are designing. But since we are interested in the secondary structure of the end product we do have to take them into account. Also the order and distances between these regions are predefined ([https://2008.igem.org/Image:Template.png figure 3]).
-
=== test ===
+
-
</div>
+
-
<div class="element">
+
====Constraints on the secondary structure====
-
test test test
+
-
</div>
+
-
=== References ===
+
Extra constraints, which are laid upon the secondary structure of the designed RNA thermometer, are specified in order to ''1.'' simplify the design process; putting restrictions on, for example the length of the hairpin greatly reduces the solution space, making the problem easier to solve, and ''2.'' to keep the secondary structure of the design as close as possible to the secondary structure of the original RNA thermometer. The [[Team:TUDelft/Temperature_analysis|analysis]] showed that mutations influence the translation efficiency without loosing its function as RNA thermometer. It is unclear if the same would hold when for example the length of the temperature sensitive hairpin is changed<sup>1</sup>. The more we change the secondary structure, the bigger the risk that it loses its function as an RNA thermometer.
-
<ol class="references">
+
-
<li id="cite_note-1">[[Team:TUDelft/Temperature_design#cite_ref-1 | ^ ]] De Smit M H, Van Duijn J (1990). "Secondary structure of the ribosome binding site determines translation efficiency: A quantitative analysis". PNAS, 1990-10, vol.87, no.19, 7668-7672. [http://www.ncbi.nlm.nih.gov/pubmed/2217199 PMID:2217199]</li>
+
-
<li> [[Team:TUDelft/Temperature_design# | ^]] Hoe N P, Goguen J D (1993). "Temperature sensing in Yersinia pestis: Translation of the LcrF activator protein is thermally regulated". J Bacteriol, 1993 December, 175(24), 7901-7909. [http://www.ncbi.nlm.nih.gov/pubmed/7504666 PMID:7504666] </li>
+
<span class="footnote">''<sup>1</sup>During the project a paper from <span id="cite_ref_1">[[Team:TUDelft/Temperature_design2#cite_note_1|[1]]]</span> has been published in which it is shown that the length of a temperature sensitive hairpin and the size of its hairpin loop can be varied without loosing the temperature sensitivity.''</span>
-
<li> [[Team:TUDelft/Temperature_design# | ^]] Chowdhurry S, Maris C, Allain F H T, Narberhaus F (2006). "Molecular basis for temperature sensing by an RNA thermometer". The EMBO Journal, 2006, 25, 2487–2497. [http://www.ncbi.nlm.nih.gov/pubmed/16710302 PMID:16710302] </li>
+
To keep the secondary structure as close as possible to the original RNA thermometer we state that mutations may only be applied to the temperature sensitive hairpin, the rest of the structure is fixed. For the temperature sensitive hairpin it is stated that the length of the hairpin, the size of the hairpin loop and the location of the hairpin loop must remain the same.
-
<li> [[Team:TUDelft/Temperature_design# | ^]] Nocker A, Hausherr T, Balsiger S, Krstulovic N, Hennecke H, Narberhaus F (2001). "A mRNA-based thermosensor controls expression of rhizobial heat shock genes". Nucleic Acids Research, 2001 December 1, 29(23):4800-4807. [http://www.ncbi.nlm.nih.gov/pubmed/11726689 PMID:11726689] </li>
+
====List of constraints====
-
<li> [[Team:TUDelft/Temperature_design# | ^]] Balsiger S, Ragaz C, Baron C, Narberhaus F. "Replicon-specific regulation of small heat shock genes in Agrobacterium tumefaciens". Journal of Bacteriology, October 2004, p. 6824-6829, Vol.186, No.20. [http://www.ncbi.nlm.nih.gov/pubmed/15466035 PMID:15466035]</li>
+
[[Image:Template.png|thumb|200px|right|'''Figure 3:''' A template for the design of a temperature sensitive hairpin. The grey dots will not be part of the final RNA thermometer part, but they have to be taken into account because they participate in the forming of the needed hairpin structure.]]
-
<li> [[Team:TUDelft/Temperature_design# | ^]] Waldminghaus T, Heidrich N, Branti S, Narberhaus F (2007). "FourU: a novel type of RNA thermometer in Salmonella". Molecular Microbiology, Volume 65, Issue 2, 413-424. [http://www.ncbi.nlm.nih.gov/pubmed/17630972 PMID:17630972]</li>
+
Recapitulating, the list of design constraints is now as follows:
-
<li> [[Team:TUDelft/Temperature_design# | ^]] Johansson J, Mandin P, Renzoni A, Chiaruttinni C, Springer M, Cossart P. "An RNA thermosensor controls expression of virulance genes in Listeria monocytogenes". Cell , Volume 110 , Issue 5 , 551-561. [http://www.ncbi.nlm.nih.gov/pubmed/12230973 PMID:12230973] </li>
+
* The sequence before the temperature sensitive hairpin remains unaltered, so mutations are only allowed within the temperature sensitive hairpin.
 +
* The overall structure of the temperature sensitive hairpin must remain the same. This means that the loop at the end should be of the same size and at the same position and the length of the hairpin remains the same. There are no restrictions to the number and locations of internal loops and bulge loops within the double helix region of the  temperature sensitive hairpin.
 +
* The start codon, SD sequence, and scar sequence are predefined sequences and are placed at fixed location within the hairpin structure.
 +
 
 +
These constraints are also captured in [https://2008.igem.org/Image:Template.png figure 3] which can be seen as a template for the design of the temperature sensitive hairpin.
 +
 
 +
==Design algorithm==
 +
 
 +
As specified in the previous section, mutations may only be applied within the temperature sensitive hairpin region. The temperature sensitive hairpin is therefore the only sequence that has to be designed and [https://2008.igem.org/Image:Template.png figure 3] serves as a template for this design.
 +
 
 +
The data analysis of the known ROSE RNA thermometers has shown a trend within the stability profiles of the temperature sensitive hairpins ([[Team:TUDelft/Temperature_analysis#Results|figure 8 of the analysis]]). This trend has been observed for the stability profiles calculated at 37°C, and as already mentioned those RNA thermometers will therefore be called  37°C switches.
 +
 
 +
To design a switch with a different temperature switching point, for example a 27°C switch, we assume that the stability profile of this RNA thermometer, calculated at 27°C, should be similar to the trend we found for the 37°C switch (calculated at 37°C). In other words, we would like the stability of the temperature sensitive hairpin of the 27°C switch at 27°C, to be the same as the stability of the temperature sensitive hairpin of the 37°C switch at 37°C.
 +
 
 +
====Algorithm====
 +
 
 +
Such a design can be made by fitting the stability profile (calculated at the desired temperature) of the hairpin to the found stability profile trend, taking into account the mentioned constraints in the previous section. The design of a temperature sensitive hairpin that switches at a certain temperature (which is given as input to the algorithm) can be made using an algorithm that acts as follows<sup>[1]</sup>:
 +
 
 +
<span class="footnote">''Note that the description of this algorithm is purely illustrative, only showing the principle. For the implementation of such an algorithm, a different approach is needed to make it more efficient (see also [[Team:TUDelft/Temperature_software#RNA_Hairpin_Designer|RNA Hairpin Designer]]).''</span>
 +
 
 +
[[Image:Boundaries.png|thumb|200px|right|'''Figure 4:''' Boundaries are set at a distance of 2.0 ''kcal/mol'' from the trend line. A designed temperature sensitive hairpin must have a stability profile that stays within the given boundaries.]]
 +
 
 +
* Take a sequence with the predefined nucleotides (start codon, SD sequence, and the scar) present at the correct locations, the rest of the nucleotides can be filled at random.
 +
* Predict the secondary structure of the RNA thermometer (at the desired switching temperature) to check if it corresponds to the required structure. If not, go to step one.
 +
* Compare the stability profile plot (calculated for the desired switching temperature) of the temperature sensitive hairpin to the trend plot. If it stays within the boundaries that are set at a certain distance form the trend line, a good enough design is found and the algorithm is terminated. Otherwise go back to step one to give it another try.
 +
 
 +
In order to break this optimization algorithm, a conditional check is added to the last step that stops the algorithm when the stability profile of the design fits the trend line good enough. This has been done by specifying boundaries that reside at a certain distance from the trend line ([https://2008.igem.org/Image:Boundaries.png figure 4]). The design process is stopped as soon as the stability profile of the design stays completely within the boundaries. Looking at the plot in which the trend line has been observed ([[Team:TUDelft/Temperature_analysis#Results|figure 8 of the analysis]]) a boundary at a distance of 2.0 ''kcal/mol'' from the trend line seems reasonable.
 +
 
 +
====Manual approach====
 +
 
 +
Because of the restricted time of the project, there was no time to wait for a software implementation of the design algorithm. Therefore a number of parts have been designed manually. The manual design has been done using a less naive approach than the given algorithm. Taking random nucleotides for the nineteen positions that are free to choose is not efficient. The number of possible sequences is enormous while only a small fraction of these will fold into the desired secondary structure.
 +
 
 +
For the manual approach the total sequence of the RNA thermometer from [https://2008.igem.org/Image:Start.png figure 2] has been taken as starting point. Instead of randomly choosing nucleotides for all the positions that may be mutated, those positions are considered one by one to see if the double helix can be destabilized through mutations of the nucleotides at these positions. For example, at a position where there is a G-C base pair of which the C may be mutated, the C can be mutated to a U, which results in a less stable hairpin. Also a mutation to G or A would be possible in order to incorporate an internal loop, which destabilizes the double helix even further.
 +
 
 +
After each mutation, the secondary structure has been predicted, using ''RNAfold'', to check if it still folds into the desired structure. When it does, the stability profile has been plotted, using ''RNAeval'' and the [[Team:TUDelft/Temperature_software#Stability_Profile_Plotter|''Stability Profile Plotter'']], to check if it stays within the trend boundaries.
 +
 
 +
==Result==
 +
 
 +
[[Image:Design2.png|thumb|200px|right|'''Figure 5:''' The secondary structures of the RNA thermometers are given at the left. The top one is the original RNA thermometer that has been used as starting point. The bottom two are the designed RNA thermometers, one with a theoretical switching point at 32\textdegree C and one with a theoretical switching point at 27\textdegree C. The stability profiles of the temperature sensitive hairpins, calculated at three different temperatures, are given at the right. The colored stability profiles shows the stability at the theoretical switching temperature. The aim of the design process was to get this colored line within the given trend boundaries.]]
 +
[[Image:Design2 1.png|thumb|200px|right|'''Figure 6:''' The stability profiles of the original and the two designed RNA thermometers, calculated at their own theoretical switching temperature. This shows that each RNA thermometer has a similar stability at its theoretical switching temperature.]]
 +
 
 +
This manual approach has been used to design two RNA thermometers, a 27°C, and a 32°C version (<code>BBa_K115016</code> and <code>BBa_K115017</code>). The secondary structures of the resulting RNA thermometers are shown in [https://2008.igem.org/Image:Design2.png figure 5], together with the stability profiles and the trend line.
 +
 
 +
The stability profiles of the temperature sensitive hairpins off all three RNA thermometers (calculated at their switching temperature) are also displayed in one plot in [https://2008.igem.org/Image:Design2_1.png figure 6]. This shows that each of them has a similar stability at their own theoretical switching temperature.
 +
 
 +
====Temperature sensitive hairpin parts====
 +
 
 +
The only sequence that has to be designed is the sequence of the temperature sensitive hairpin. The rest of the structure is kept intact to stay as close as possible to the secondary structure of the original RNA thermometer that is used as starting point. It has been shown by Chowdhury <span id="cite_ref_2">[[Team:TUDelft/Temperature_design2#cite_note_2|[2]]]</span> that a temperature sensitive hairpin on itself, so without the rest of the ROSE structure, also works as a temperature sensitive switch.
 +
 
 +
Three temperature sensitive hairpins are added to the registry to test if these short sequences can also be used as temperature sensitive switch. These are the temperature sensitive hairpins of the ROSE RNA thermometer as designed in the first design phase ([http://partsregistry.org/wiki/index.php?title=Part:BBa_K115001 <code>BBa_K115001</code>]), and the temperature sensitive hairpins of the 27°C and 32°C switch as designed in the second design phase ([http://partsregistry.org/wiki/index.php?title=Part:BBa_K115016 <code>BBa_K115016</code>] and [http://partsregistry.org/wiki/index.php?title=Part:BBa_K115017 <code>BBa_K115017</code>]). This are the parts [http://partsregistry.org/wiki/index.php?title=Part:BBa_K115020 <code>BBa_K115020</code>], [http://partsregistry.org/wiki/index.php?title=Part:BBa_K115019 <code>BBa_K115019</code>], and [http://partsregistry.org/wiki/index.php?title=Part:BBa_K115018 <code>BBa_K115018</code>] respectively.
 +
 
 +
If these short hairpins happen to be enough to get a working RNA thermometer, this would be a nice starting point for a third design phase in which these small parts are further optimized. For us, there's not enough time left to do this, but it could be a good starting point for a next years iGEM project.
 +
 
 +
== References ==
 +
 
 +
<ol class="references">
 +
<li id="cite_note_1"> [[Team:TUDelft/Temperature_design2#cite_ref_1 | ^]] Juliane Neupert and Daniel Karcher and Ralph Bock. Design of simple synthetic RNA thermometers for temperature-controlled gene expression in Escherichia coli. ''Nucleic Acids Research'', 1-9, 2008. [http://www.ncbi.nlm.nih.gov/pubmed/18753148 PMID:18753148] </li>
 +
<li id="cite_note_2"> [[Team:TUDelft/Temperature_design2#cite_ref_2 | ^]] Saheli Chowdhury, Curdin Ragaz, Emma Kreuger, and Franz Narberhaus. Temperature-controlled Structural Alterations of an RNA Thermometer. ''The Journal of Biological Chemistry'', 278(48):47915-47921, 2003. [http://www.ncbi.nlm.nih.gov/pubmed/12963744 PMID:12963744] </li>
</ol>
</ol>
{{Template:TUDelftiGEM2008_sidebar}}
{{Template:TUDelftiGEM2008_sidebar}}

Latest revision as of 19:33, 29 October 2008

Contents

Design phase II:

Changing the switching temperature of the RNA thermometer

Introducing existing RNA thermometers into Escherichia coli provides us with a number of temperature sensitive switches, all with their own switching temperature. But this only provides a few switching temperatures, which are also relatively close to each other. In order to build a biothermometer, we needed more switching temperatures that have a bigger range. Note, however, that there will always be an upper bound in that the E. coli cannot survive temperatures above 40°C, and a lower bound in that the activity of the cell will drop at low temperatures.

This chapter presents a method that has been used to design RNA thermometer parts that should theoretically switch at a predefined temperature. The method used is based on the data analysis of which the results can be found on the analysis page. It has been used to design two new parts, one that should theoretically be switched on at 27°C ([http://partsregistry.org/wiki/index.php?title=Part:BBa_K115016 BBa_K115016]), and two that should theoretically be switched on at 32°C ([http://partsregistry.org/wiki/index.php?title=Part:BBa_K115017 BBa_K115017]).

Secondary structure stability and switching temperature

Figure 1: A. A lot of G-C base pairs make a very stable hairpin (having a very low free energy). B. The same hairpin structure with A-U base pairs is much less stable (having a higher free energy). C-E Bulge loops and internal loops destabilize the structure even further.

The switching temperature of the RNA thermometer depends on the stability of the hairpin structure that surrounds the Shine Dalgarno region (see also Requirements for an RNA therometer). This secondary structure will melt when it is exposed to a certain amount of heat. The more stable the hairpin structure is, the more heat will be needed to melt it, and the less stable the hairpin structure is, the less heat will be needed to melt it. In case of the RNA thermometer, this melting is needed to expose the Shine Dalgarno sequence which will cause the initiation of translation.

In other words, a less stable hairpin structure will melt at a lower temperature than a more stable hairpin and thus, for the less stable hairpin, the translation will initiate at a lower temperature than the more stable one. This way we could theoretically shift the switching temperature to a lower or higher temperature by making the hairpin structure surrounding the SD region less or more stable respectively.

The stability of the hairpin structure depends on the base-pairing nucleotides within the secondary structure. For example, C-G base pairs are very strong, so a hairpin structure containing a lot of those base pairs will be very stable. A-U and U-G base pairs are less strong, so these can be used to destabilize the hairpin structure. The same can be achieved by incorporating bulge loops and internal loops within the hairpin structure (figure 1).

Constraints

There are two things that have to be taken into account. Firstly, adding mutations in order to enforce or weaken the temperature sensitive region can also cause the RNA to fold into a completely different secondary structure. This way it can lose its function as an RNA thermometer. Secondly, there are some constraints to the possible mutations. The start codon and the Shine Dalgarno sequence must off course remain unaltered and in case of a standard biobrick the scar must also be part of the temperature sensitive hairpin.

Design starting point

Figure 2: The secondary structure of the RNA thermometer that has been used as starting point of the design process is displayed at the left. At the right the stability profile of the temperature sensitive hairpin is highlighted (light blue). Also the found stability trend (see data analysis) and the specified boundaries (see design algorithm) are highlighted (dark blue). The grey lines are the stability profiles of the other 31 ROSE RNA thermometers from the [http://rfam.sanger.ac.uk Rfam database].

Since we have a number of structures of which we know that they work as an RNA thermometer, the most logical option is to take one of those structures as the starting point for the design of an RNA thermometer with a different switching temperature. The one chosen is a ROSE RNA thermometer from Bradyrhizobium japonicum USDA 110 ([http://www.ncbi.nlm.nih.gov/sites/entrez?Db=genome&Cmd=ShowDetailView&TermToSearch=272 NC_004463]) residing at location [http://www.ncbi.nlm.nih.gov/projects/sviewer/?id=NC_004463.1&v=5784144..5784239 5784144-5784239] within the genome, which is at the 5' side of the gene [http://www.ncbi.nlm.nih.gov/sites/entrez?db=gene&cmd=search&term=NC_004463%20hspb hspB]. There was no specific reason for the choice of this ROSE RNA thermometer, this is just a ROSE RNA thermometer taken from the [http://rfam.sanger.ac.uk Rfam database].

We will call this RNA thermometer a 37°C switch, while it actually reaches maximum translation several degrees higher. The reason that we call it a 37°C switch is because of the characteristic stability profile we found during the analysis, and which is calculated at 37°C. Having a 27°C switch would then mean that this is an RNA thermometer of which the stability profile at 27 degrees is similar to the trend as given in figure 8 of the analysis. This switch should then theoretically be almost fully switched on at 27°C, as is the case for the 37°C version at 37°C.

Figure 2 shows the secondary structure of the chosen RNA thermometer as predicted by RNAfold. Also the stability profile that is characteristic for this RNA thermometer at 37°C is shown in light blue. The thick dark blue line is the trend as found in the data analysis and the thin dark blue lines are design boundaries which will be discussed later in this chapter. The stability profiles of the other 31 ROSE RNA thermometers from the Rfam database are shown in grey.

Design constraints

Having a starting point, we need a way to turn this RNA thermometer into one that switches at a different temperature, and that conforms to the Standard Biological Part standard. This can be achieved by increasing or decreasing the stability of the hairpin structure, as explained at the begin of this page. But, as also mentioned, there are certain constraints that has to be laid upon the design.

Constraints on the sequence

Looking at the nucleotide sequence of the design, there is the constraint that there must be a start codon, a Shine Dalgarno sequence, and a scar present within the sequence. Note that this is the sequence after ligation of the RNA thermometer part to a protein coding part and that the scar and start codon are not part of the RNA thermometer part that we are designing. But since we are interested in the secondary structure of the end product we do have to take them into account. Also the order and distances between these regions are predefined (figure 3).

Constraints on the secondary structure

Extra constraints, which are laid upon the secondary structure of the designed RNA thermometer, are specified in order to 1. simplify the design process; putting restrictions on, for example the length of the hairpin greatly reduces the solution space, making the problem easier to solve, and 2. to keep the secondary structure of the design as close as possible to the secondary structure of the original RNA thermometer. The analysis showed that mutations influence the translation efficiency without loosing its function as RNA thermometer. It is unclear if the same would hold when for example the length of the temperature sensitive hairpin is changed1. The more we change the secondary structure, the bigger the risk that it loses its function as an RNA thermometer.

1During the project a paper from [1] has been published in which it is shown that the length of a temperature sensitive hairpin and the size of its hairpin loop can be varied without loosing the temperature sensitivity.

To keep the secondary structure as close as possible to the original RNA thermometer we state that mutations may only be applied to the temperature sensitive hairpin, the rest of the structure is fixed. For the temperature sensitive hairpin it is stated that the length of the hairpin, the size of the hairpin loop and the location of the hairpin loop must remain the same.

List of constraints

Figure 3: A template for the design of a temperature sensitive hairpin. The grey dots will not be part of the final RNA thermometer part, but they have to be taken into account because they participate in the forming of the needed hairpin structure.

Recapitulating, the list of design constraints is now as follows:

  • The sequence before the temperature sensitive hairpin remains unaltered, so mutations are only allowed within the temperature sensitive hairpin.
  • The overall structure of the temperature sensitive hairpin must remain the same. This means that the loop at the end should be of the same size and at the same position and the length of the hairpin remains the same. There are no restrictions to the number and locations of internal loops and bulge loops within the double helix region of the temperature sensitive hairpin.
  • The start codon, SD sequence, and scar sequence are predefined sequences and are placed at fixed location within the hairpin structure.

These constraints are also captured in figure 3 which can be seen as a template for the design of the temperature sensitive hairpin.

Design algorithm

As specified in the previous section, mutations may only be applied within the temperature sensitive hairpin region. The temperature sensitive hairpin is therefore the only sequence that has to be designed and figure 3 serves as a template for this design.

The data analysis of the known ROSE RNA thermometers has shown a trend within the stability profiles of the temperature sensitive hairpins (figure 8 of the analysis). This trend has been observed for the stability profiles calculated at 37°C, and as already mentioned those RNA thermometers will therefore be called 37°C switches.

To design a switch with a different temperature switching point, for example a 27°C switch, we assume that the stability profile of this RNA thermometer, calculated at 27°C, should be similar to the trend we found for the 37°C switch (calculated at 37°C). In other words, we would like the stability of the temperature sensitive hairpin of the 27°C switch at 27°C, to be the same as the stability of the temperature sensitive hairpin of the 37°C switch at 37°C.

Algorithm

Such a design can be made by fitting the stability profile (calculated at the desired temperature) of the hairpin to the found stability profile trend, taking into account the mentioned constraints in the previous section. The design of a temperature sensitive hairpin that switches at a certain temperature (which is given as input to the algorithm) can be made using an algorithm that acts as follows[1]:

Note that the description of this algorithm is purely illustrative, only showing the principle. For the implementation of such an algorithm, a different approach is needed to make it more efficient (see also RNA Hairpin Designer).

Figure 4: Boundaries are set at a distance of 2.0 kcal/mol from the trend line. A designed temperature sensitive hairpin must have a stability profile that stays within the given boundaries.
  • Take a sequence with the predefined nucleotides (start codon, SD sequence, and the scar) present at the correct locations, the rest of the nucleotides can be filled at random.
  • Predict the secondary structure of the RNA thermometer (at the desired switching temperature) to check if it corresponds to the required structure. If not, go to step one.
  • Compare the stability profile plot (calculated for the desired switching temperature) of the temperature sensitive hairpin to the trend plot. If it stays within the boundaries that are set at a certain distance form the trend line, a good enough design is found and the algorithm is terminated. Otherwise go back to step one to give it another try.

In order to break this optimization algorithm, a conditional check is added to the last step that stops the algorithm when the stability profile of the design fits the trend line good enough. This has been done by specifying boundaries that reside at a certain distance from the trend line (figure 4). The design process is stopped as soon as the stability profile of the design stays completely within the boundaries. Looking at the plot in which the trend line has been observed (figure 8 of the analysis) a boundary at a distance of 2.0 kcal/mol from the trend line seems reasonable.

Manual approach

Because of the restricted time of the project, there was no time to wait for a software implementation of the design algorithm. Therefore a number of parts have been designed manually. The manual design has been done using a less naive approach than the given algorithm. Taking random nucleotides for the nineteen positions that are free to choose is not efficient. The number of possible sequences is enormous while only a small fraction of these will fold into the desired secondary structure.

For the manual approach the total sequence of the RNA thermometer from figure 2 has been taken as starting point. Instead of randomly choosing nucleotides for all the positions that may be mutated, those positions are considered one by one to see if the double helix can be destabilized through mutations of the nucleotides at these positions. For example, at a position where there is a G-C base pair of which the C may be mutated, the C can be mutated to a U, which results in a less stable hairpin. Also a mutation to G or A would be possible in order to incorporate an internal loop, which destabilizes the double helix even further.

After each mutation, the secondary structure has been predicted, using RNAfold, to check if it still folds into the desired structure. When it does, the stability profile has been plotted, using RNAeval and the Stability Profile Plotter, to check if it stays within the trend boundaries.

Result

Figure 5: The secondary structures of the RNA thermometers are given at the left. The top one is the original RNA thermometer that has been used as starting point. The bottom two are the designed RNA thermometers, one with a theoretical switching point at 32\textdegree C and one with a theoretical switching point at 27\textdegree C. The stability profiles of the temperature sensitive hairpins, calculated at three different temperatures, are given at the right. The colored stability profiles shows the stability at the theoretical switching temperature. The aim of the design process was to get this colored line within the given trend boundaries.
Figure 6: The stability profiles of the original and the two designed RNA thermometers, calculated at their own theoretical switching temperature. This shows that each RNA thermometer has a similar stability at its theoretical switching temperature.

This manual approach has been used to design two RNA thermometers, a 27°C, and a 32°C version (BBa_K115016 and BBa_K115017). The secondary structures of the resulting RNA thermometers are shown in figure 5, together with the stability profiles and the trend line.

The stability profiles of the temperature sensitive hairpins off all three RNA thermometers (calculated at their switching temperature) are also displayed in one plot in figure 6. This shows that each of them has a similar stability at their own theoretical switching temperature.

Temperature sensitive hairpin parts

The only sequence that has to be designed is the sequence of the temperature sensitive hairpin. The rest of the structure is kept intact to stay as close as possible to the secondary structure of the original RNA thermometer that is used as starting point. It has been shown by Chowdhury [2] that a temperature sensitive hairpin on itself, so without the rest of the ROSE structure, also works as a temperature sensitive switch.

Three temperature sensitive hairpins are added to the registry to test if these short sequences can also be used as temperature sensitive switch. These are the temperature sensitive hairpins of the ROSE RNA thermometer as designed in the first design phase ([http://partsregistry.org/wiki/index.php?title=Part:BBa_K115001 BBa_K115001]), and the temperature sensitive hairpins of the 27°C and 32°C switch as designed in the second design phase ([http://partsregistry.org/wiki/index.php?title=Part:BBa_K115016 BBa_K115016] and [http://partsregistry.org/wiki/index.php?title=Part:BBa_K115017 BBa_K115017]). This are the parts [http://partsregistry.org/wiki/index.php?title=Part:BBa_K115020 BBa_K115020], [http://partsregistry.org/wiki/index.php?title=Part:BBa_K115019 BBa_K115019], and [http://partsregistry.org/wiki/index.php?title=Part:BBa_K115018 BBa_K115018] respectively.

If these short hairpins happen to be enough to get a working RNA thermometer, this would be a nice starting point for a third design phase in which these small parts are further optimized. For us, there's not enough time left to do this, but it could be a good starting point for a next years iGEM project.

References

  1. ^ Juliane Neupert and Daniel Karcher and Ralph Bock. Design of simple synthetic RNA thermometers for temperature-controlled gene expression in Escherichia coli. Nucleic Acids Research, 1-9, 2008. [http://www.ncbi.nlm.nih.gov/pubmed/18753148 PMID:18753148]
  2. ^ Saheli Chowdhury, Curdin Ragaz, Emma Kreuger, and Franz Narberhaus. Temperature-controlled Structural Alterations of an RNA Thermometer. The Journal of Biological Chemistry, 278(48):47915-47921, 2003. [http://www.ncbi.nlm.nih.gov/pubmed/12963744 PMID:12963744]