Team:Uppsala University/Modelling

From 2012.igem.org

(Difference between revisions)
 
(30 intermediate revisions not shown)
Line 9: Line 9:
</div>
</div>
</div>
</div>
-
 
-
<br>
 
<table id="ttable" align="center">
<table id="ttable" align="center">
<tr>
<tr>
Line 21: Line 19:
<td colspan="2">
<td colspan="2">
<table>
<table>
-
<tr>
+
<tr><td class="subtext"><h2>Modelling of small RNA interacting with the AAC(6)’5-UTR</h2></td>
-
 
+
-
 
+
-
<td class="subtext"><h2>Modelling of small RNA interacting with the AAC(6)’5-UTR</h2></td>
+
<td valign="bottom"><a id="top" href="#top">Back to top</a></td></tr>
<td valign="bottom"><a id="top" href="#top">Back to top</a></td></tr>
</table>
</table>
</td>
</td>
</tr>
</tr>
-
 
<tr>
<tr>
Line 36: Line 30:
<td valign="top">
<td valign="top">
<p>
<p>
-
The modelling of our constructed sRNA is a vital part for evaluating hypotheses of how sRNA molecules interact with their target. One canonical approach for sRNA stability is to evaluate secondary structures by minimum free energy (MFE) approaches. This gives a pointer for approximating the probability of different kind of interactions.
+
The modelling of our synthetic sRNA is a vital part for evaluating hypotheses of how sRNA molecules interact with their target mRNA. One canonical approach for sRNA stability is to evaluate secondary structures by minimum free energy (MFE) approaches. This gives a pointer for approximating the probability of different kind of interactions.
</p>
</p>
 +
<br>
 +
<p>
 +
<b>Stabile structures prevents translation</b>
 +
A thermodynamic scheme of RNA interactions at different stages were made to show how an interaction between mRNA and sRNA could possibly work. So called toeholds are belived to play an important role in the interaction between sRNA and mRNA. [1]
 +
Understanding the interaction between RNA with different structures is a key to learn how to design a sRNA. Often the hybridization reaction between the sRNA and the mRNA starts with the unpaired nucleotides, the so called toeholds. To start a basepairing reaction, the sRNA needs a few basepairs to start hybridize to. The sRNA and the mRNA can then create a more stable secondary complex, hindering translation. </p>
 +
<br>
 +
<table id="ptable">
 +
<td id="ptd1">
 +
<tr>
 +
  <a href="https://static.igem.org/mediawiki/2012/a/a6/SRNAstructure_legend_low.png"><img src="https://static.igem.org/mediawiki/2012/a/a6/SRNAstructure_legend_low.png" width="250"></a>
 +
  </td>
 +
</tr>
 +
</table>
 +
<br>
 +
<table id="ptable">
 +
<b>Reaction diagram sRNA</b>
 +
<tr>
 +
  <td id="ptd1">
 +
  <a href="https://static.igem.org/mediawiki/2012/7/78/Reaction_diagrame_sRNA_low.png"><img src="https://static.igem.org/mediawiki/2012/7/78/Reaction_diagrame_sRNA_low.png" width="580"></a>
 +
  </td>
 +
</tr>
 +
<tr style="border-top:thin solid;border-color:black">
 +
  <td>
 +
  <p>The sRNA 17 and mRNA for the AAC(6)-UTR were first calculated in CLC main workbench separately, after that we calculated the interaction structure between them .
 +
<br>
 +
This is a model of how a secondary complex can be created between the sRNA and mRNA. </p>
 +
  </td>
 +
</tr>
 +
</table>
 +
<br>
<p>
<p>
<br>
<br>
Line 46: Line 70:
<br>
<br>
<p>
<p>
-
Spot42 has like many other smallRNA, two distinctive parts. One that binds to a mRNA sequence, and another sequence that interacts with the Hfq-rna binding protein that is belived to play an important role in the function of smallRNA. (Holmquist, 2012)
+
The sRNA Spot42 has like many other sRNA two distinctive functional parts, a recognition and a scaffold sequence. The recognition sequence binds the mRNA, while the scaffold sequence interacts with Hfq, an RNA binding protein. The sRNA-Hfq complex play an important role in the function of sRNA translational regulation. (Holmquist, 2012)
</p>
</p>
<br>
<br>
-
<p style="margin-right:10px;font-size:10px;padding:0 0 0 0;float:left"><a href="https://static.igem.org/mediawiki/2012/b/b5/Natural_spot42_structure_medium.png"><img src="https://static.igem.org/mediawiki/2012/b/b5/Natural_spot42_structure_medium.png" width=300></a><br>
+
<p style="margin-right:10px;font-size:10px;padding:0 0 0 0;float:left"><a href="https://static.igem.org/mediawiki/2012/f/f6/Spot42.2d.3d.1.png"><img src="https://static.igem.org/mediawiki/2012/f/f6/Spot42.2d.3d.1.png" width="350"></a><br>
-
Here you can see the native spot42 with the scaffold marked in blue.
+
Figure above. Panel A shows the secondary structure of scaffold<br>
 +
Spot42 and panel B shows a MMB generated three dimensional of <br>
 +
the Spot42 scaffold. HP1 and HP2 shows the positions of the two hairpins.
</p>
</p>
<p>
<p>
-
<b>Picture of native spot42</b><br>
+
The secondary structure of the scaffold of Spot42 (figA) is known. However, no 3D structure of spot42 has been solved to date.  We have generated a 3D structure from the information in the secondary structure with molecular dynamics software MacroMoleculeBuilder (MMB) [1] and the 3D structure of Spot42 scaffold is presented in figure B. MMB allows the user to have full control of the biomolecule (DNA, RNA and/or protein) by applying forces between atoms, constraints and flexibility to different parts of the biomolecule.
-
To gain better understanding of the mechanisms of smallRNA translation inhibition we have made 2D models of all our found smallRNAs. The 2D models were calculated in CLC main workbench.<br>
+
 
-
IntaRNA, an RNA-interaction prediction software adapted for sRNA and ncRNA interactions (Smith et. al., 2010) was used  to predict the sRNA-mRNA interactions of the candidate sRNAs.
+
[3]
 +
<br>
 +
<br>
 +
<br>
 +
To gain better understanding of the mechanisms of small RNA translation inhibition we have made 2D models of all our synthetic small RNAs. The 2D models were calculated in CLC main workbench using adapted versions of maximum free energy algorithms developed by [2] and thermodynamical parameters of Mfold version 3.
 +
IntaRNA, an RNA-interaction prediction software adapted for sRNA and ncRNA interactions [5] was used  to predict the sRNA-mRNA interactions of the candidate sRNAs.
</p>
</p>
<br>
<br>
<p>
<p>
-
The sRNA were isolated, sequenced and analyzed to find the hybridizing base pairs. The sequencing of the sRNAs that showed downregulation of YFP also had a matching sequence in the 5’UTR of our target mRNA.<br>
+
The sRNAs were isolated, sequenced and analyzed to find the hybridizing base pairs. The sequences of the sRNAs that showed downregulation of YFP mainly had matching sequences in the 5’UTR of our target mRNA, although some of the sRNAs was shown to hybridize at the SYFP2 coding mRNA region. Two of these were further studied and modelled.
-
Altough, a few of the sRNA that downregulated YFP was shown to hybridize at the YFP mRNA region. Two of these were further studied and modelled. At last, a prediction of the structure between the sRNA UU17 and AAC(6’)UTR mRNA were modelled.
+
At last, a prediction of the structure between the sRNA UU17 and AAC(6’)UTR mRNA was modelled.
</p>
</p>
-
<br>
 
Line 68: Line 98:
<td colspan="2">
<td colspan="2">
<table>
<table>
 +
<a name="candsrnas">
<tr><td class="subtext"><h2>Candidate sRNAs</h2></td>
<tr><td class="subtext"><h2>Candidate sRNAs</h2></td>
<td valign="bottom"><a id="top" href="#top">Back to top</a></td></tr>
<td valign="bottom"><a id="top" href="#top">Back to top</a></td></tr>
Line 80: Line 111:
<table id="ptable">
<table id="ptable">
<p>
<p>
-
As you can see below, our candidate sRNA that downregulate both the fluorescence and the antibiotic resistance gene also shows to hybridize at the region close to the RBS and the start codon. This supports the idea that many sRNA prevents the ribosome from binding to the RBS, and thereby preventing translation. (Erik Holmquist, 2012)<br>
+
As you can see below, our sRNAs that downregulate the SYFP2 and the antibiotic resistance gene also seem to hybridize at the region close to the RBS and the start codon. This supports the idea that many sRNAs prevent the ribosome from binding to the RBS, thereby preventing translation. (Erik Holmquist, 2012)
-
Also our predictions data from CLC shows that the hybridizing area have a strong secondary structure, with small hairpin loops.
+
Our structure prediction data from CLC also shows that the hybridizing region of our silencing synthetic sRNAs often have strong secondary structures, with small hairpin loops.  
</p>
</p>
<br>
<br>
Line 95: Line 126:
  <tr>
  <tr>
   <td id="ptd1">
   <td id="ptd1">
-
  ∆G =  -31.2 kcal/mol<br>
+
<p>sRNA UU17 was found to be one of the two most effective sRNAs against the resistance gene AAC(6’).
 +
  <br><br> ∆G =  -31.2 kcal/mol<br></p>
   </td>
   </td>
   <td id="ptd2">
   <td id="ptd2">
Line 101: Line 133:
   Number of hybridizing base pairs = 17<br>
   Number of hybridizing base pairs = 17<br>
   Maximum number of hybridizing bases in a row = 13<br>
   Maximum number of hybridizing bases in a row = 13<br>
-
   Result on Etest > 256 µg/ml<br>
+
   Result on Etest = 48 µg/ml<br>
   ∆G = -14.3kcal/mol<br>
   ∆G = -14.3kcal/mol<br>
 +
SYFP2 downregulation = 78%
 +
   </td>
   </td>
  </tr>
  </tr>
Line 121: Line 155:
  <tr>
  <tr>
   <td id="ptd1">
   <td id="ptd1">
-
   ∆G = -31.9 kcal/mol
+
<p>sRNA UU37 was the the other of the two most effective sRNAs against the resistance gene AAC(6’).
 +
   <br><br>∆G = -31.9 kcal/mol
   </td>
   </td>
   <td id="ptd2">
   <td id="ptd2">
Line 129: Line 164:
   Result on Etest = 64 µg/ml<br>
   Result on Etest = 64 µg/ml<br>
   ∆G = -15.7 kcal/mol<br>
   ∆G = -15.7 kcal/mol<br>
 +
SYFP2 downregulation: 82%
   </td>
   </td>
  </tr>
  </tr>
Line 142: Line 178:
   </td>
   </td>
   <td id="ptd1">
   <td id="ptd1">
-
   <a href="https://static.igem.org/mediawiki/2012/c/c4/SRNA_46-mRNA_medium.png"><img src="https://static.igem.org/mediawiki/2012/c/c4/SRNA_46-mRNA_medium.png" id="intpic"></a>
+
   <a href="https://static.igem.org/mediawiki/2012/c/c4/SRNA_46-mRNA_medium.png"><img src="https://static.igem.org/mediawiki/2012/c/c4/SRNA_46-mRNA_medium.png" width="350"></a>
   </td>
   </td>
  </tr>
  </tr>
  <tr>
  <tr>
   <td id="ptd1">
   <td id="ptd1">
-
   ∆G =  -28.7 kcal/mol<br>
+
<p>Complimentary to AAC(6) UTR<br><br>
 +
   ∆G =  -28.7 kcal/mol<br></p>
   </td>
   </td>
   <td id="ptd2">
   <td id="ptd2">
Line 154: Line 191:
   Maximum number of hybridizing bases in a row = 8 bp<br>
   Maximum number of hybridizing bases in a row = 8 bp<br>
   Result on Etest > 256µg/ml <br>
   Result on Etest > 256µg/ml <br>
 +
SYFP2 downregulation: 83%
   ∆G = -11.2kcal/mol<br>
   ∆G = -11.2kcal/mol<br>
   </td>
   </td>
Line 168: Line 206:
   </td>
   </td>
   <td id="ptd1">
   <td id="ptd1">
-
   <a href="https://static.igem.org/mediawiki/2012/3/36/SRNA55-mRNA_medium.png"><img src="https://static.igem.org/mediawiki/2012/3/36/SRNA55-mRNA_medium.png" id="intpic"></a>
+
   <a href="https://static.igem.org/mediawiki/2012/3/36/SRNA55-mRNA_medium.png"><img src="https://static.igem.org/mediawiki/2012/3/36/SRNA55-mRNA_medium.png" width="300"></a>
   </td>
   </td>
  </tr>
  </tr>
  <tr>
  <tr>
-
   <td id="ptd1">
+
   <td id="ptd1"><p>
-
   ∆G =  -31.2 kcal/mol<br>
+
Complimentary to AAC(6) UTR<br><br>
 +
   ∆G =  -31.2 kcal/mol<br></p>
   </td>
   </td>
   <td id="ptd2">
   <td id="ptd2">
Line 180: Line 219:
   Maximum number of hybridizing bases in a row = 8 <br>
   Maximum number of hybridizing bases in a row = 8 <br>
   Result on Etest = 124 µg/ml <br>
   Result on Etest = 124 µg/ml <br>
 +
SYFP2 downregulation= 58%
   ∆G = -10.2kcal/mol<br>
   ∆G = -10.2kcal/mol<br>
   </td>
   </td>
Line 190: Line 230:
<td colspan="2">
<td colspan="2">
<table>
<table>
 +
<a name="srnayfp">
<tr><td class="subtext"><h2>sRNAs matching SYFP2</h2></td>
<tr><td class="subtext"><h2>sRNAs matching SYFP2</h2></td>
<td valign="bottom"><a id="top" href="#top">Back to top</a></td></tr>
<td valign="bottom"><a id="top" href="#top">Back to top</a></td></tr>
Line 205: Line 246:
</b>
</b>
<br>
<br>
-
Below are two examples of structures that were found to match regions in the coding sequence of the SYFP2. It is apparent that false positives will appear when using a reporter system that itself can be repressed, but by modelling interactions with the target mRNA, we save much time knowing which sRNAs are not of direct interest.
+
Below are two examples of structures that were found to match regions in the coding sequence of the SYFP2. It is apparent that false positives will appear when using a reporter system that itself can be repressed, but by modelling interactions with the target mRNA, we save much time knowing which sRNAs are not of direct interest for downregulation of the actual target gene.
</p>
</p>
<br>
<br>
Line 216: Line 257:
   </td>
   </td>
   <td id="ptd1">
   <td id="ptd1">
-
   <a href="https://static.igem.org/mediawiki/2012/6/63/SRNA52-mRNA_medium.png"><img src="https://static.igem.org/mediawiki/2012/6/63/SRNA52-mRNA_medium.png" id="intpic"></a>
+
   <a href="https://static.igem.org/mediawiki/2012/6/63/SRNA52-mRNA_medium.png"><img src="https://static.igem.org/mediawiki/2012/6/63/SRNA52-mRNA_medium.png" width="300"></a>
   </td>
   </td>
  </tr>
  </tr>
  <tr>
  <tr>
-
   <td id="ptd1">
+
   <td id="ptd1"><p>
-
   ∆G =  -28.3kcal/mol<br>
+
sRNA basepairing to reporter gene SYFP2<br><br>
 +
   ∆G =  -28.3kcal/mol<br></p>
   </td>
   </td>
   <td id="ptd2">
   <td id="ptd2">
Line 227: Line 269:
   Number of hybridizing base pairs =10 <br>
   Number of hybridizing base pairs =10 <br>
   Maximum number of hybridizing bases in a row = 6 <br>
   Maximum number of hybridizing bases in a row = 6 <br>
-
   Gibbs free energy = -5.6 kcal/mol<br>
+
   ∆G = -5.6 kcal/mol<br>
   </td>
   </td>
  </tr>
  </tr>
Line 248: Line 290:
  </tr>
  </tr>
  <tr>
  <tr>
-
   <td id="ptd1">
+
   <td id="ptd1"><p>
-
   ∆G =  -34.5 kcal/mol<br>
+
sRNA basepairing to reporter gene SYFP2<br><br>
 +
   ∆G =  -34.5 kcal/mol<br></p>
   </td>
   </td>
   <td id="ptd2">
   <td id="ptd2">
Line 255: Line 298:
   Number of hybridizing base pairs =21 <br>
   Number of hybridizing base pairs =21 <br>
   Maximum number of hybridizing bases in a row =9  <br>
   Maximum number of hybridizing bases in a row =9  <br>
-
   Gibbs free energy = -10.9kcal/mol<br>
+
   ∆G = -10.9kcal/mol<br>
   </td>
   </td>
  </tr>
  </tr>
Line 262: Line 305:
</table>
</table>
-
 
+
<br>
<table id="ptable">
<table id="ptable">
  <tr>
  <tr>
 +
<p>
 +
Although a randomized library was used to find an optimized sRNA for the downregulation of sRNA, we asked ourselves how an sRNA with perfect complementarity to the RBS and UTR region would work. This complementary sRNA was synthesized and we tested its function.
 +
 +
The results were not simple to decipher, because it seems that this sRNA targets more than the AAC(6’) UTR. The construct confered a large fitness cost on the bacteria, with very slow growth and all sorts of different colony morphologies and measurements of SYFP2 activity on the reporter strain. Some results showed a significant decrease in SYFP2 activity, while measurements of a different clone showed only moderate to non existing downregulation of SYFP2.<br>
 +
In addition to this, the sRNA selected from the random library seemed to have secondary structure on the hybridizing part, while the optimal secondary structure of our perfectly matching sRNA is without any base pairing within the antisense region. This makes us speculate that this sRNA might have unspecific base pair matching with many different bacterial mRNAs.
 +
<br> <br>The approach to construct a perfect complementary is maybe not the best way to make a working sRNA, because there are so many other mechanisms involved that are difficult to predict.
 +
</p>
 +
<br>
<b>sRNA UU01</b>
<b>sRNA UU01</b>
   <td id="ptd1">
   <td id="ptd1">
Line 270: Line 321:
   </td>
   </td>
   <td id="ptd1">
   <td id="ptd1">
-
   <a href="https://static.igem.org/mediawiki/2012/6/6a/SRNA01_medium.png"><img src="https://static.igem.org/mediawiki/2012/6/6a/SRNA01_medium.png" id="intpic"></a>
+
   <a href="https://static.igem.org/mediawiki/2012/6/6a/SRNA01_medium.png"><img src="https://static.igem.org/mediawiki/2012/6/6a/SRNA01_medium.png" width="400"></a>
   </td>
   </td>
  </tr>
  </tr>
  <tr>
  <tr>
-
   <td id="ptd1">
+
   <td id="ptd1"><p>
 +
Perfect complimentary to AAC(6)<br><br></p>
   ∆G =  -31.2 kcal/mol<br>
   ∆G =  -31.2 kcal/mol<br>
   </td>
   </td>
Line 282: Line 334:
   Maximum number of hybridizing bases in a row = 8 <br>
   Maximum number of hybridizing bases in a row = 8 <br>
   Result on Etest = 124 µg/ml <br>
   Result on Etest = 124 µg/ml <br>
-
   Gibbs free energy = -10.2kcal/mol<br>
+
   ∆G = -10.2kcal/mol<br>
   </td>
   </td>
  </tr>
  </tr>
Line 288: Line 340:
<br>
<br>
<br>
<br>
-
<br>
 
-
<br>
 
-
 
 +
<br>
 +
<b>AAC(6')-SYFP2 and sRNA UU17 interaction </b>
<table id="ptable">
<table id="ptable">
-
<b>sRNA UU17</b>
+
<tr>
-
<tr>
+
-
  <td id="ptd1">
+
-
  <a href="https://static.igem.org/mediawiki/2012/f/f8/SRNA17structure_medium.png"><img src="https://static.igem.org/mediawiki/2012/f/f8/SRNA17structure_medium.png" id="rnapic"></a>
+
-
  </td>
+
<td id="ptd1">
<td id="ptd1">
-
   <a href="https://static.igem.org/mediawiki/2012/3/3d/AAC%286%29%27UTR_Structure_medium2.png"><img src="https://static.igem.org/mediawiki/2012/3/3d/AAC%286%29%27UTR_Structure_medium2.png" id="rnapic"></a>
+
   <a href="https://static.igem.org/mediawiki/2012/3/3d/AAC%286%29%27UTR_Structure_medium2.png"><img src="https://static.igem.org/mediawiki/2012/3/3d/AAC%286%29%27UTR_Structure_medium2.png" width="300" ></a>
   </td>
   </td>
<td id="ptd1">
<td id="ptd1">
   <a href="https://static.igem.org/mediawiki/2012/a/aa/AAC%286%29%27UTR%2BsRNA17_medium.png"><img src="https://static.igem.org/mediawiki/2012/a/aa/AAC%286%29%27UTR%2BsRNA17_medium.png" id="rnapic"></a>
   <a href="https://static.igem.org/mediawiki/2012/a/aa/AAC%286%29%27UTR%2BsRNA17_medium.png"><img src="https://static.igem.org/mediawiki/2012/a/aa/AAC%286%29%27UTR%2BsRNA17_medium.png" id="rnapic"></a>
   </td>
   </td>
-
 
-
 
-
 
  </tr>
  </tr>
  <tr>
  <tr>
   <td id="ptd1">
   <td id="ptd1">
-
   ∆G =  -31.2 kcal/mol<br>
+
<p>The mRNA of AAC(6') 5-UTR and a part of the coding strand for SYFP2 was included.<br><br>
 +
   ∆G =  -23,4 kcal/mol<br></pr>
   </td>
   </td>
   <td id="ptd2">
   <td id="ptd2">
 +
<p>The mRNA and sRNA UU17 in a combined structure<br><br>
 +
  ∆G =  -60,8 kcal/mol kcal/mol<br></p>
   </td>
   </td>
  </tr>
  </tr>
</table>
</table>
 +
<br>
 +
<tr>
 +
    <td colspan="2">
 +
        <table>
 +
        <tr>
 +
            <td class="subtext">
 +
            <h2>
 +
                References
 +
            </h2>
 +
            </td>
 +
            <td valign="bottom">
 +
            <a id="top" href="#top" name="top">Back to top</a>
 +
            </td>
 +
        </tr>
 +
        </table>
 +
    </td>
 +
    </tr>
 +
    <tr>
 +
    <td valign="top">
 +
        <p>
 +
[1] Rodrigo, G., Landrain, T.E., Jaramillo, A., 2012. De novo automated design of small RNA circuits for engineering synthetic riboregulation in living cells. Proc. Natl. Acad. Sci. U.S.A. 109, 15271–15276.
 +
<br>
 +
[2] Zuker, M., 1989b. The use of dynamic programming algorithms in  RNA secondary structure prediction, 159–184. In
 +
Waterman, M., ed. Mathematical Methods for DNA Sequences, CRC Press, Boca Raton, FL<br>
 +
[3]. Samuel C Flores and Russ B Altman. 2010 Turning limited experimental information intio 3D models of RNA. RNA 16(9):1769-78.
 +
<br>Notes about CLC workbench: Uses thermodynamical parameters of Mfold, version 3, found on http://www.bioinfo.rpi.edu/zukerm/rna/energy/
 +
<br>[4] Erik Holmquist
 +
Macromolecular Matchmaking
 +
Mechanisms and Biology of Bacterial Small RNAs, 2012
 +
<br>
 +
[5]
 +
Smith, C., Heyne, S., Richter, A.S., Will, S., Backofen, R., 2010. Freiburg RNA Tools: a web server integrating INTARNA, EXPARNA and LOCARNA.
 +
Nucleic Acids Research 38, W373–W377.
 +
        </p>
 +
    </td>
-
 
+
    </tr>
-
 
+
-
 
+
-
 
+
-
 
+
-
 
+
-
<br>
+
</table></table>
</table></table>
-
 
-
 
-
</div>
 
-
<br><br>
 
</html>
</html>
{{Uppsala_footer}}
{{Uppsala_footer}}

Latest revision as of 03:09, 27 October 2012

Team Uppsala University – iGEM 2012

Modelling of small RNA interacting with the AAC(6)’5-UTR

Back to top


The modelling of our synthetic sRNA is a vital part for evaluating hypotheses of how sRNA molecules interact with their target mRNA. One canonical approach for sRNA stability is to evaluate secondary structures by minimum free energy (MFE) approaches. This gives a pointer for approximating the probability of different kind of interactions.


Stabile structures prevents translation A thermodynamic scheme of RNA interactions at different stages were made to show how an interaction between mRNA and sRNA could possibly work. So called toeholds are belived to play an important role in the interaction between sRNA and mRNA. [1] Understanding the interaction between RNA with different structures is a key to learn how to design a sRNA. Often the hybridization reaction between the sRNA and the mRNA starts with the unpaired nucleotides, the so called toeholds. To start a basepairing reaction, the sRNA needs a few basepairs to start hybridize to. The sRNA and the mRNA can then create a more stable secondary complex, hindering translation.



Reaction diagram sRNA

The sRNA 17 and mRNA for the AAC(6)-UTR were first calculated in CLC main workbench separately, after that we calculated the interaction structure between them .
This is a model of how a secondary complex can be created between the sRNA and mRNA.



Using a native scaffold to stabilize the interactions
The wildtype small RNA of E.coli K-12 MG1655 that we used as a template for constructing our own smallRNA is named Spot42 and has been shown to be interacting with the Hfq protein.


The sRNA Spot42 has like many other sRNA two distinctive functional parts, a recognition and a scaffold sequence. The recognition sequence binds the mRNA, while the scaffold sequence interacts with Hfq, an RNA binding protein. The sRNA-Hfq complex play an important role in the function of sRNA translational regulation. (Holmquist, 2012)



Figure above. Panel A shows the secondary structure of scaffold
Spot42 and panel B shows a MMB generated three dimensional of
the Spot42 scaffold. HP1 and HP2 shows the positions of the two hairpins.

The secondary structure of the scaffold of Spot42 (figA) is known. However, no 3D structure of spot42 has been solved to date. We have generated a 3D structure from the information in the secondary structure with molecular dynamics software MacroMoleculeBuilder (MMB) [1] and the 3D structure of Spot42 scaffold is presented in figure B. MMB allows the user to have full control of the biomolecule (DNA, RNA and/or protein) by applying forces between atoms, constraints and flexibility to different parts of the biomolecule. [3]


To gain better understanding of the mechanisms of small RNA translation inhibition we have made 2D models of all our synthetic small RNAs. The 2D models were calculated in CLC main workbench using adapted versions of maximum free energy algorithms developed by [2] and thermodynamical parameters of Mfold version 3. IntaRNA, an RNA-interaction prediction software adapted for sRNA and ncRNA interactions [5] was used to predict the sRNA-mRNA interactions of the candidate sRNAs.


The sRNAs were isolated, sequenced and analyzed to find the hybridizing base pairs. The sequences of the sRNAs that showed downregulation of YFP mainly had matching sequences in the 5’UTR of our target mRNA, although some of the sRNAs was shown to hybridize at the SYFP2 coding mRNA region. Two of these were further studied and modelled. At last, a prediction of the structure between the sRNA UU17 and AAC(6’)UTR mRNA was modelled.

Candidate sRNAs

Back to top

As you can see below, our sRNAs that downregulate the SYFP2 and the antibiotic resistance gene also seem to hybridize at the region close to the RBS and the start codon. This supports the idea that many sRNAs prevent the ribosome from binding to the RBS, thereby preventing translation. (Erik Holmquist, 2012) Our structure prediction data from CLC also shows that the hybridizing region of our silencing synthetic sRNAs often have strong secondary structures, with small hairpin loops.


sRNA UU17

sRNA UU17 was found to be one of the two most effective sRNAs against the resistance gene AAC(6’).

∆G = -31.2 kcal/mol

Length in base pairs = 82 bp
Number of hybridizing base pairs = 17
Maximum number of hybridizing bases in a row = 13
Result on Etest = 48 µg/ml
∆G = -14.3kcal/mol
SYFP2 downregulation = 78%

sRNA UU37

sRNA UU37 was the the other of the two most effective sRNAs against the resistance gene AAC(6’).

∆G = -31.9 kcal/mol

Length in base pairs = 83 bp
Number of hybridizing base pairs = 24
Maximum number of hybridizing bases in a row = 7
Result on Etest = 64 µg/ml
∆G = -15.7 kcal/mol
SYFP2 downregulation: 82%

sRNA UU46

Complimentary to AAC(6) UTR

∆G = -28.7 kcal/mol

Length in base pairs = 65 bp
Number of hybridizing base pairs = 8 bp
Maximum number of hybridizing bases in a row = 8 bp
Result on Etest > 256µg/ml
SYFP2 downregulation: 83% ∆G = -11.2kcal/mol

sRNA UU55

Complimentary to AAC(6) UTR

∆G = -31.2 kcal/mol

Length in base pairs = 83
Number of hybridizing base pairs = 8
Maximum number of hybridizing bases in a row = 8
Result on Etest = 124 µg/ml
SYFP2 downregulation= 58% ∆G = -10.2kcal/mol

sRNAs matching SYFP2

Back to top


Below are two examples of structures that were found to match regions in the coding sequence of the SYFP2. It is apparent that false positives will appear when using a reporter system that itself can be repressed, but by modelling interactions with the target mRNA, we save much time knowing which sRNAs are not of direct interest for downregulation of the actual target gene.


sRNA UU52

sRNA basepairing to reporter gene SYFP2

∆G = -28.3kcal/mol

Length in base pairs =82
Number of hybridizing base pairs =10
Maximum number of hybridizing bases in a row = 6
∆G = -5.6 kcal/mol

sRNA UU53

sRNA basepairing to reporter gene SYFP2

∆G = -34.5 kcal/mol

Length in base pairs =85
Number of hybridizing base pairs =21
Maximum number of hybridizing bases in a row =9
∆G = -10.9kcal/mol

Although a randomized library was used to find an optimized sRNA for the downregulation of sRNA, we asked ourselves how an sRNA with perfect complementarity to the RBS and UTR region would work. This complementary sRNA was synthesized and we tested its function. The results were not simple to decipher, because it seems that this sRNA targets more than the AAC(6’) UTR. The construct confered a large fitness cost on the bacteria, with very slow growth and all sorts of different colony morphologies and measurements of SYFP2 activity on the reporter strain. Some results showed a significant decrease in SYFP2 activity, while measurements of a different clone showed only moderate to non existing downregulation of SYFP2.
In addition to this, the sRNA selected from the random library seemed to have secondary structure on the hybridizing part, while the optimal secondary structure of our perfectly matching sRNA is without any base pairing within the antisense region. This makes us speculate that this sRNA might have unspecific base pair matching with many different bacterial mRNAs.

The approach to construct a perfect complementary is maybe not the best way to make a working sRNA, because there are so many other mechanisms involved that are difficult to predict.


sRNA UU01

Perfect complimentary to AAC(6)

∆G = -31.2 kcal/mol
Length in base pairs = 83
Number of hybridizing base pairs = 8
Maximum number of hybridizing bases in a row = 8
Result on Etest = 124 µg/ml
∆G = -10.2kcal/mol



AAC(6')-SYFP2 and sRNA UU17 interaction

The mRNA of AAC(6') 5-UTR and a part of the coding strand for SYFP2 was included.

∆G = -23,4 kcal/mol

The mRNA and sRNA UU17 in a combined structure

∆G = -60,8 kcal/mol kcal/mol


References

Back to top

[1] Rodrigo, G., Landrain, T.E., Jaramillo, A., 2012. De novo automated design of small RNA circuits for engineering synthetic riboregulation in living cells. Proc. Natl. Acad. Sci. U.S.A. 109, 15271–15276.
[2] Zuker, M., 1989b. The use of dynamic programming algorithms in RNA secondary structure prediction, 159–184. In Waterman, M., ed. Mathematical Methods for DNA Sequences, CRC Press, Boca Raton, FL
[3]. Samuel C Flores and Russ B Altman. 2010 Turning limited experimental information intio 3D models of RNA. RNA 16(9):1769-78.
Notes about CLC workbench: Uses thermodynamical parameters of Mfold, version 3, found on http://www.bioinfo.rpi.edu/zukerm/rna/energy/
[4] Erik Holmquist Macromolecular Matchmaking Mechanisms and Biology of Bacterial Small RNAs, 2012
[5] Smith, C., Heyne, S., Richter, A.S., Will, S., Backofen, R., 2010. Freiburg RNA Tools: a web server integrating INTARNA, EXPARNA and LOCARNA. Nucleic Acids Research 38, W373–W377.

Sponsors