Team:BGU Israel/Modeling/scRNA

From 2014.igem.org

(Difference between revisions)
Line 4: Line 4:
-
      <section>
+
<section>
       <div style="height:12px"> </div>
       <div style="height:12px"> </div>
     <div style="margin-bottom:10px"> <img src="https://static.igem.org/mediawiki/2014/5/5f/BGU14scrna.png" width="1079" height="197" alt=""/></div>
     <div style="margin-bottom:10px"> <img src="https://static.igem.org/mediawiki/2014/5/5f/BGU14scrna.png" width="1079" height="197" alt=""/></div>
Line 53: Line 53:
-
   <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><img  src="https://static.igem.org/mediawiki/2014/9/99/BGU14model1fig1.png" style="height:300px;border:solid"/><br>Figure 1 – Need to put content</p>
+
   <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/9/99/BGU14model1fig1.png"><img  src="https://static.igem.org/mediawiki/2014/9/99/BGU14model1fig1.png" style="height:300px;border:solid"/></a><br>Figure 1 – ‘Intelligent Medication’ scRNA mechanism</p>
</p>
</p>
Line 59: Line 59:
         </div>
         </div>
         <div class="col2">
         <div class="col2">
-
   <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><img  src="https://static.igem.org/mediawiki/2014/a/a5/BGU14fig2.png" style="height:130px;"/>
+
   <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/a/a5/BGU14fig2.png"><img  src="https://static.igem.org/mediawiki/2014/a/a5/BGU14fig2.png" style="height:130px;"/></a>
-
   <br>  Figure 2</p>
+
   <br>  Figure 2- Lengths of subsequences</p>
    
    
<p>We wrote a python program which can create the required scRNA, basing on online and offline calculations. The program gets the trigger gene and the target gene as inputs. The output of the program is a list of scRNAs (A•B) and their detection targets (X), that predicted to be the most effective for this mechanism, considering the trigger and target genes.</p>
<p>We wrote a python program which can create the required scRNA, basing on online and offline calculations. The program gets the trigger gene and the target gene as inputs. The output of the program is a list of scRNAs (A•B) and their detection targets (X), that predicted to be the most effective for this mechanism, considering the trigger and target genes.</p>
Line 104: Line 104:
         We can see it in the following example:</li>
         We can see it in the following example:</li>
</ul>
</ul>
-
<p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><img  src="https://static.igem.org/mediawiki/2014/2/27/BGU14fig3.png" style=" width::250px;"/>
+
<p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/2/27/BGU14fig3.png"><img  src="https://static.igem.org/mediawiki/2014/2/27/BGU14fig3.png" style=" width::250px;"/></a>
-
  <br>
+
    
-
   Figure 3</p>
+
</p>
    
    
   <p>The idea of this score is that we  don&rsquo;t want a situation which a high score component compensate the other  component, this is due to the nature of the mechanism (both states are equally  important, with interdependency). Therefore, we chose multiplication as the  scoring function:<br>
   <p>The idea of this score is that we  don&rsquo;t want a situation which a high score component compensate the other  component, this is due to the nature of the mechanism (both states are equally  important, with interdependency). Therefore, we chose multiplication as the  scoring function:<br>
Line 139: Line 139:
           </ul>
           </ul>
          
          
-
           <p  style="margin-left:-40px"><img src="https://static.igem.org/mediawiki/2014/e/ec/BGU14cmdDir.png" style="height:190px"  /></p>
+
           <p  style="margin-left:-40px"><a href="https://static.igem.org/mediawiki/2014/e/ec/BGU14cmdDir.png"><img src="https://static.igem.org/mediawiki/2014/e/ec/BGU14cmdDir.png" style="height:190px"  /></a></p>
         </ol>
         </ol>
       </div>
       </div>
Line 146: Line 146:
           <li><span dir="LTR"> </span>Open the command window and  run:<br>
           <li><span dir="LTR"> </span>Open the command window and  run:<br>
             &gt; <strong>python  scRNAalgo.py trigger_file target_file</strong></li>
             &gt; <strong>python  scRNAalgo.py trigger_file target_file</strong></li>
-
           <p align="center" style="margin-left:-40px"><img src="https://static.igem.org/mediawiki/2014/0/0c/BGU14cmdRun.png" /></p>
+
           <p align="center" style="margin-left:-40px"><a href="https://static.igem.org/mediawiki/2014/0/0c/BGU14cmdRun.png"><img src="https://static.igem.org/mediawiki/2014/0/0c/BGU14cmdRun.png" /></a></p>
           <li><span dir="LTR"> </span>During the running of the  program, it creates a set of all the results from the servers in the directory,  for calculation purposes. After the program ends, an appropriate message will  appear on the screen and the results will be written into the <strong>output</strong> file.</li>
           <li><span dir="LTR"> </span>During the running of the  program, it creates a set of all the results from the servers in the directory,  for calculation purposes. After the program ends, an appropriate message will  appear on the screen and the results will be written into the <strong>output</strong> file.</li>
         </ol>
         </ol>
Line 156: Line 156:
           Score: %%%%%%%%<br>
           Score: %%%%%%%%<br>
           A: %%%%%%%%%<br>
           A: %%%%%%%%%<br>
-
           B: %%%%%%%%%<br><div id="test3"></div>
+
           B: %%%%%%%%%<br>
           X: %%%%%%%%%</p>
           X: %%%%%%%%%</p>
 +
          <div id="test3"></div>
         </div>
         </div>
         </div>
         </div>
Line 163: Line 164:
          
          
          
          
-
         <div class="textCont" style="height: 1300px; margin-bottom: 12px;">
+
         <div class="textCont" style="height: 1500px; margin-bottom: 12px;">
       <h3 style="border-bottom:dashed;border-color:#000000">How We Used It</h3>
       <h3 style="border-bottom:dashed;border-color:#000000">How We Used It</h3>
         <p>We ran the program on the  following sequences, located in the two input files:</p>
         <p>We ran the program on the  following sequences, located in the two input files:</p>
         <p style="margin-left:180px"><u>Trigger gene sequence - <strong>Apolipoprotein3</strong></u></p>
         <p style="margin-left:180px"><u>Trigger gene sequence - <strong>Apolipoprotein3</strong></u></p>
-
         <p align="center"><img src="https://static.igem.org/mediawiki/2014/8/8f/BGU14fig4.png"  /></p>
+
         <p style="font-size:10px;margin-left:180px">
 +
        TGCTCAGTTCATCCCTAGAGGCAGCTGCTCCAGGAACAGAGGTGCCATGCAGCCCCGGGTACTCCTTGTTGTTGCCCTCCTGGCGCTCCTGGCCTCTGCCCGAGCTTCAGAGGCCGA<br>
 +
        GGATGCCTCGCTTCATGCAGGGTTACATGAAGCACGCCACCAAGACCGCCAAGGATGCACTGAGCAGCGTGCAGGAGTCCCAGGTGGCCCAGCAGGCCAGGGGCTGGGTGACCGATG<br>
 +
        GCTTCAGTTCCCTGAAAGACTACTGGAGCACCGTTAAGGACAAGTTCTCTGAGTTCTGGGATTTGGACCCTGAGGTCAGACCAACTTCAGCCGTGGCTGCCTGAGACCTCAATACCC<br>
 +
        CAAGTCCACCTGCCTATCCATCCTGCGAGCTCCTTGGGTCCTGCAATCTCCAGGGCTGCCCCTGTAGGTTGCTTAAAAGGGACAGTATTCTCAGTGCTCTCCTACCCCACCTCATGC<br>
 +
        CTGGCCCCCCTCCAGGCATGCTGGCCTCCCAATAAAGCTGGACAAGAAGCTGCTATG</p>
 +
       
         <p style="margin-left:180px"><u>Target gene sequence - <strong>PTP1B</strong></u></p>
         <p style="margin-left:180px"><u>Target gene sequence - <strong>PTP1B</strong></u></p>
-
         <p align="center"><img src="https://static.igem.org/mediawiki/2014/8/81/BGU14fig5.png"  /></p>
+
         <p style="font-size:10px;margin-left:180px">
-
        <div style="margin-left:180px">
+
        CCTTCTGCTTCAGGGGCGGAGCCCCTGGCAGGCGTGATGCGTAGTTCCGGCTGCCGGTTGACATGAAGAAGCAGCAGCGGCTAGGGCGGCGGTAGCTGCAGGGGTCGGGGATTGCAG<br>
-
        <p><strong>Results:</strong><br><br>
+
        CGGGCCTCGGGGCTAAGAGCGCGACGCGGCCTAGAGCGGCAGACGGCGCAGTGGGCCGAGAAGGAGGCGCAGCAGCCGCCCTGGCCCGTCATGGAGATGGAAAAGGAGTTCGAGCAG<br>
 +
        ATCGACAAGTCCGGGAGCTGGGCGGCCATTTACCAGGATATCCGACATGAAGCCAGTGACTTCCCATGTAGAGTGGCCAAGCTTCCTAAGAACAAAAACCGAAATAGGTACAGAGAC<br>
 +
        GTCAGTCCCTTTGACCATAGTCGGATTAAACTACATCAAGAAGATAATGACTATATCAACGCTAGTTTGATAAAAATGGAAGAAGCCCAAAGGAGTTACATTCTTACCCAGGGCCCT<br>
 +
        TTGCCTAACACATGCGGTCACTTTTGGGAGATGGTGTGGGAGCAGAAAAGCAGGGGTGTCGTCATGCTCAACAGAGTGATGGAGAAAGGTTCGTTAAAATGCGCACAATACTGGCCA<br>
 +
        CAAAAAGAAGAAAAAGAGATGATCTTTGAAGACACAAATTTGAAATTAACATTGATCTCTGAAGATATCAAGTCATATTATACAGTGCGACAGCTAGAATTGGAAAACCTTACAACC<br>
 +
        CAAGAAACTCGAGAGATCTTACATTTCCACTATACCACATGGCCTGACTTTGGAGTCCCTGAATCACCAGCCTCATTCTTGAACTTTCTTTTCAAAGTCCGAGAGTCAGGGTCACTC<br>
 +
        AGCCCGGAGCACGGGCCCGTTGTGGTGCACTGCAGTGCAGGCATCGGCAGGTCTGGAACCTTCTGTCTGGCTGATACCTGCCTCTTGCTGATGGACAAGAGGAAAGACCCTTCTTCC<br>
 +
        GTTGATATCAAGAAAGTGCTGTTAGAAATGAGGAAGTTTCGGATGGGGCTGATCCAGACAGCCGACCAGCTGCGCTTCTCCTACCTGGCTGTGATCGAAGGTGCCAAATTCATCATG<br>
 +
        GGGGACTCTTCCGTGCAGGATCAGTGGAAGGAGCTTTCCCACGAGGACCTGGAGCCCCCACCCGAGCATATCCCCCCACCTCCCCGGCCACCCAAACGAATCCTGGAGCCACACAAT<br>
 +
        GGGAAATGCAGGGAGTTCTTCCCAAATCACCAGTGGGTGAAGGAAGAGACCCAGGAGGATAAAGACTGCCCCATCAAGGAAGAAAAAGGAAGCCCCTTAAATGCCGCACCCTACGGC<br>
 +
        ATCGAAAGCATGAGTCAAGACACTGAAGTTAGAAGTCGGGTCGTGGGGGGAAGTCTTCGAGGTGCCCAGGCTGCCTCCCCAGCCAAAGGGGAGCCGTCACTGCCCGAGAAGGACGAG<br>
 +
        GACCATGCACTGAGTTACTGGAAGCCCTTCCTGGTCAACATGTGCGTGGCTACGGTCCTCACGGCCGGCGCTTACCTCTGCTACAGGTTCCTGTTCAACAGCAACACATAGCCTGAC<br>
 +
        CCTCCTCCACTCCACCTCCACCCACTGTCCGCCTCTGCCCGCAGAGCCCACGCCCGACTAGCAGGCATGCCGCGGTAGGTAAGGGCCGCCGGACCGCGTAGAGAGCCGGGCCCCGGA<br>
 +
        CGGACGTTGGTTCTGCACTAAAACCCATCTTCCCCGGATGTGTGTCTCACCCCTCATCCTTTTACTTTTTGCCCCTTCCACTTTGAGTACCAAATCCACAAGCCATTTTTTGAGGAG<br>
 +
        AGTGAAAGAGAGTACCATGCTGGCGGCGCAGAGGGAAGGGGCCTACACCCGTCTTGGGGCTCGCCCCACCCAGGGCTCCCTCCTGGAGCATCCCAGGCGGGCGGCACGCCAACAGCC<br>
 +
        CCCCCCTTGAATCTGCAGGGAGCAACTCTCCACTCCATATTTATTTAAACAATTTTTTCCCCAAAGGCATCCATAGTGCACTAGCATTTTCTTGAACCAATAATGTATTAAAATTTT<br>
 +
        TTGATGTCAGCCTTGCATCAAGGGCTTTATCAAAAAGTACAATAATAAATCCTCAGGTAGTACTGGGAATGGAAGGCTTTGCCATGGGCCTGCTGCGTCAGACCAGTACTGGGAAGG<br>
 +
        AGGACGGTTGTAAGCAGTTGTTATTTAGTGATATTGTGGGTAACGTGAGAAGATAGAACAATGCTATAATATATAATGAACACGTGGGTATTTAATAAGAAACATGATGTGAGATTA<br>
 +
        CTTTGTCCCGCTTATTCTCCTCCCTGTTATCTGCTAGATCTAGTTCTCAATCACTGCTCCCCCGTGTGTATTAGAATGCATGTAAGGTCTTCTTGTGTCCTGATGAAAAATATGTGC<br>
 +
        TTGAAATGAGAAACTTTGATCTCTGCTTACTAATGTGCCCCATGTCCAAGTCCAACCTGCCTGTGCATGACCTGATCATTACATGGCTGTGGTTCCTAAGCCTGTTGCTGAAGTCAT<br>
 +
        TGTCGCTCAGCAATAGGGTGCAGTTTTCCAGGAATAGGCATTTGCCTAATTCCTGGCATGACACTCTAGTGACTTCCTGGTGAGGCCCAGCCTGTCCTGGTACAGCAGGGTCTTGCT<br>
 +
        GTAACTCAGACATTCCAAGGGTATGGGAAGCCATATTCACACCTCACGCTCTGGACATGATTTAGGGAAGCAGGGACACCCCCCGCCCCCCACCTTTGGGATCAGCCTCCGCCATTC<br>
 +
        CAAGTCAACACTCTTCTTGAGCAGACCGTGATTTGGAAGAGAGGCACCTGCTGGAAACCACACTTCTTGAAACAGCCTGGGTGACGGTCCTTTAGGCAGCCTGCCGCCGTCTCTGTC<br>
 +
        CCGGTTCACCTTGCCGAGAGAGGCGCGTCTGCCCCACCCTCAAACCCTGTGGGGCCTGATGGTGCTCACGACTCTTCCTGCAAAGGGAACTGAAGACCTCCACATTAAGTGGCTTTT<br>
 +
        TAACATGAAAAACACGGCAGCTGTAGCTCCCGAGCTACTCTCTTGCCAGCATTTTCACATTTTGCCTTTCTCGTGGTAGAAGCCAGTACAGAGAAATTCTGTGGTGGGAACATTCGA<br>
 +
        GGTGTCACCCTGCAGAGCTATGGTGAGGTGTGGATAAGGCTTAGGTGCCAGGCTGTAAGCATTCTGAGCTGGGCTTGTTGTTTTTAAGTCCTGTATATGTATGTAGTAGTTTGGGTG<br>
 +
        TGTATATATAGTAGCATTTCAAAATGGACGTACTGGTTTAACCTCCTATCCTTGGAGAGCAGCTGGCTCTCCACCTTGTTACACATTATGTTAGAGAGGTAGCGAGCTGCTCTGCTA<br>
 +
        TATGCCTTAAGCCAATATTTACTCATCAGGTCATTATTTTTTACAATGGCCATGGAATAAACCATTTTTACAAAAATAAAAACAAAAAAAGCAAGGTGTTTTGGTATAATACCTTTT<br>
 +
        CAGGTGTGTGTGGATACGTGGCTGCATGACCGGGTGGGTGGGGGGGAGTGTCTCAGGGTCTTCTGTGACCTCACAGAACTGTCAGACTGTACAGTTTTCCAACTTGCCATATTCATG<br>
 +
        ATGGGTTTGCATTTTAGCTGCAACAATAAAATTTTTTTCTAAAGAACATGAAAAAAAAAAAAA</p>
 +
 
 +
<div>
 +
      <p><strong><u>Results:</u></strong><br><br>
           The beginning of the output file:<br>
           The beginning of the output file:<br>
         <p  style="font-size:14px">Score: 9.999665000418599e-09<br>
         <p  style="font-size:14px">Score: 9.999665000418599e-09<br>
Line 192: Line 231:
<p>
<p>
-
This scRNA targets X which starting at the 215th basis of the target gene up to the 243th:</p> <div id="test4"></div>
+
<div id="test4"></div>
 +
This scRNA targets X which starting at the 215th basis of the target gene up to the 243th:</p>
<p style="font-size:14px"><strong> X:</strong> CAGCAGGCCAGGGGCUGGGUGACCGAUGG (215 - 243)
<p style="font-size:14px"><strong> X:</strong> CAGCAGGCCAGGGGCUGGGUGACCGAUGG (215 - 243)
</p></div>
</p></div>
Line 208: Line 248:
The duplex A·B is the only complex (100%) that formed from A and B:  </p>
The duplex A·B is the only complex (100%) that formed from A and B:  </p>
-
<p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><img src="https://static.igem.org/mediawiki/2014/4/4a/BGU14noTriggerAB.png" height="350"  style="border:solid;"/><br>   
+
<p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/4/4a/BGU14noTriggerAB.png"><img src="https://static.igem.org/mediawiki/2014/4/4a/BGU14noTriggerAB.png" height="350"  style="border:solid;"/></a><br>   
-
   Figure 3</p>
+
   Figure 3- The secondary structure of duplex A·B</p>
         <p><u>ON state (with trigger)</u></p>
         <p><u>ON state (with trigger)</u></p>
<p>There is a concentration of 50% of B and 50% of X·A.  That is to say, all of the B forms will transform into hairpins, and will look like  in the figure below: </p>
<p>There is a concentration of 50% of B and 50% of X·A.  That is to say, all of the B forms will transform into hairpins, and will look like  in the figure below: </p>
-
  <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><img src="https://static.igem.org/mediawiki/2014/3/3c/BGU14withTriggerB.png" height="350"  style="border:solid;"/><br> Figure 4</p>
+
  <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/3/3c/BGU14withTriggerB.png"><img src="https://static.igem.org/mediawiki/2014/3/3c/BGU14withTriggerB.png" height="350"  style="border:solid;"/></a><br> Figure 4- The secondary structure of B</p>
   
   
   
   
Line 219: Line 259:
  <p>The A which has  separated from B, assembles, along with X, the complex X·A:</p>
  <p>The A which has  separated from B, assembles, along with X, the complex X·A:</p>
   
   
-
  <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><img src="https://static.igem.org/mediawiki/2014/9/98/BGU14withTriggerATrigger.png" height="350"  style="border:solid;"/><br>Figure 5</p>
+
  <p align="center" style="font-size:14px; font-weight:bold; line-height:normal"><a href="https://static.igem.org/mediawiki/2014/9/98/BGU14withTriggerATrigger.png"><img src="https://static.igem.org/mediawiki/2014/9/98/BGU14withTriggerATrigger.png" height="350"  style="border:solid;"/></a><br>Figure 5- The secondary structure of complex X·A</p>
  <p><strong>Conclusions</strong>:<br>
  <p><strong>Conclusions</strong>:<br>
   We used the  program we wrote in the design and test of our scRNA mechanism. It helped us to  choose the scRNA we wanted to test in our lab, among a very large number of  potential scRNAs, and is also even spared us the manual analysis of a single  scRNA. We offer the code to of the program as open source, from the notion that  it can help anyone who want to design the same mechanism we used, for other  purposes. Of course the program is adjustable to the needs of any objective,  with a different trigger gene and a different target gene.</p>
   We used the  program we wrote in the design and test of our scRNA mechanism. It helped us to  choose the scRNA we wanted to test in our lab, among a very large number of  potential scRNAs, and is also even spared us the manual analysis of a single  scRNA. We offer the code to of the program as open source, from the notion that  it can help anyone who want to design the same mechanism we used, for other  purposes. Of course the program is adjustable to the needs of any objective,  with a different trigger gene and a different target gene.</p>

Revision as of 16:52, 17 October 2014

Overview

A python program we wrote in order to simplify the design stage of the scRNA in the mechanism of ‘Intelligent Medication’. Our program constructs the most suitable scRNAs for this mechanism, out of various automated calculations.

How It Works


User Manual


How We Used It


Construction of our scRNA

Jump To Results Analysis

Overview


Purpose: Construct the most suitable scRNA, for the following mechanism.


Figure 1 – ‘Intelligent Medication’ scRNA mechanism


Figure 2- Lengths of subsequences

We wrote a python program which can create the required scRNA, basing on online and offline calculations. The program gets the trigger gene and the target gene as inputs. The output of the program is a list of scRNAs (A•B) and their detection targets (X), that predicted to be the most effective for this mechanism, considering the trigger and target genes.



How It Works

Main Goal:

Finding X, A and B that will make the optimal distinction between OFF and ON states:

  • OFF state - minimal production of Dicer substrate B in the absence of detection target X, the presence of mRNA silencing target Y, or the presence of mRNA off-target Z.
  • ON state - strong production of B in the presence of short RNA detection target Xs (‘a-b-c’) or full-length mRNA detection target X.

Steps:

  1. From the target gene (given as an input), find an appropriate sub-sequence to be the silencing target with a length of |y|+|z|, and consequently define y and z.
  2. For each optional sub-sequence with the length of X (|a|+|b|+|c|), of the trigger gene (given as an input), do:
    1. Define this sub-sequence as X.
    2. Split X into three parts, and consequently define a, b and c.
    3. Define A = z · c* · b* · a*.
    4. Define B = z · c* · b · c · z* · y*.
    5. Scoring the complex of A·B (scRNA) and X. **
  3. Create a list of all the scRNAs (and their mRNA detection target X), sorted by descending score, and set it as the output. The scRNA with the best score is the first.
* - Reverse complement.
** - Using an external online service.

Scoring a scRNA A scRNA’s score consists of two components:

  • ScoreOFF (OFF state) – the relative concentration of duplex A•B in the absence of X.
  • ScoreON (ON state) – the relative concentration of hairpin B in the presence of X. We tried to find a simple scoring function which can describe the suitability of a scRNA in a proper way. There are two simple functions which are candidates for the scoring function, summation and multiplication. We can see it in the following example:

The idea of this score is that we don’t want a situation which a high score component compensate the other component, this is due to the nature of the mechanism (both states are equally important, with interdependency). Therefore, we chose multiplication as the scoring function:
The scRNA’s total score = ScoreOFF * ScoreON

External online services usage:

  • InvivoGen - for finding an appropriate sub-sequence for mRNA detection target.
  • NUPACK Server - for analyzing secondary structures and finding the concentration of complexes.

User Manual

System Requirements:

  1. Python 3.x installed.
  2. Internet connection.
  3. Free space of a few gigabytes (depends on the trigger gene length).

Running the program:

  1. The program consist of two python files that can be downloaded from here:
  2. [link 1], [link 2].
  3. Put the two ‘.py’ files (“serverHandlers.py”, “scRNA.py”) in the same directory.
  4. Create two files in this directory:
    • Trigger file (e.g. “trigger.txt”), which consist the trigger gene sequence.
    • Target file (e.g. “target.txt”), which consist the target gene sequence.

  1. Open the command window and run:
    > python scRNAalgo.py trigger_file target_file
  2. During the running of the program, it creates a set of all the results from the servers in the directory, for calculation purposes. After the program ends, an appropriate message will appear on the screen and the results will be written into the output file.

Output File (Results):
During the running, the program creates a directory which called “All Results” (if it not exists). In this directory the program will create a sub-directory for every run, with the time and date that make up its name. After the run, this sub-directory will contain all the data that the program has processed of every potential scRNA.
The output file is located in ‘All Results/000000000000 run/000000000000 run - Output.txt’
The output consists all of the scRNAs the program analyzed, in descending order by the total score. A single scRNA will showed like that:
Score: %%%%%%%%
A: %%%%%%%%%
B: %%%%%%%%%
X: %%%%%%%%%

How We Used It

We ran the program on the following sequences, located in the two input files:

Trigger gene sequence - Apolipoprotein3

TGCTCAGTTCATCCCTAGAGGCAGCTGCTCCAGGAACAGAGGTGCCATGCAGCCCCGGGTACTCCTTGTTGTTGCCCTCCTGGCGCTCCTGGCCTCTGCCCGAGCTTCAGAGGCCGA
GGATGCCTCGCTTCATGCAGGGTTACATGAAGCACGCCACCAAGACCGCCAAGGATGCACTGAGCAGCGTGCAGGAGTCCCAGGTGGCCCAGCAGGCCAGGGGCTGGGTGACCGATG
GCTTCAGTTCCCTGAAAGACTACTGGAGCACCGTTAAGGACAAGTTCTCTGAGTTCTGGGATTTGGACCCTGAGGTCAGACCAACTTCAGCCGTGGCTGCCTGAGACCTCAATACCC
CAAGTCCACCTGCCTATCCATCCTGCGAGCTCCTTGGGTCCTGCAATCTCCAGGGCTGCCCCTGTAGGTTGCTTAAAAGGGACAGTATTCTCAGTGCTCTCCTACCCCACCTCATGC
CTGGCCCCCCTCCAGGCATGCTGGCCTCCCAATAAAGCTGGACAAGAAGCTGCTATG

Target gene sequence - PTP1B

CCTTCTGCTTCAGGGGCGGAGCCCCTGGCAGGCGTGATGCGTAGTTCCGGCTGCCGGTTGACATGAAGAAGCAGCAGCGGCTAGGGCGGCGGTAGCTGCAGGGGTCGGGGATTGCAG
CGGGCCTCGGGGCTAAGAGCGCGACGCGGCCTAGAGCGGCAGACGGCGCAGTGGGCCGAGAAGGAGGCGCAGCAGCCGCCCTGGCCCGTCATGGAGATGGAAAAGGAGTTCGAGCAG
ATCGACAAGTCCGGGAGCTGGGCGGCCATTTACCAGGATATCCGACATGAAGCCAGTGACTTCCCATGTAGAGTGGCCAAGCTTCCTAAGAACAAAAACCGAAATAGGTACAGAGAC
GTCAGTCCCTTTGACCATAGTCGGATTAAACTACATCAAGAAGATAATGACTATATCAACGCTAGTTTGATAAAAATGGAAGAAGCCCAAAGGAGTTACATTCTTACCCAGGGCCCT
TTGCCTAACACATGCGGTCACTTTTGGGAGATGGTGTGGGAGCAGAAAAGCAGGGGTGTCGTCATGCTCAACAGAGTGATGGAGAAAGGTTCGTTAAAATGCGCACAATACTGGCCA
CAAAAAGAAGAAAAAGAGATGATCTTTGAAGACACAAATTTGAAATTAACATTGATCTCTGAAGATATCAAGTCATATTATACAGTGCGACAGCTAGAATTGGAAAACCTTACAACC
CAAGAAACTCGAGAGATCTTACATTTCCACTATACCACATGGCCTGACTTTGGAGTCCCTGAATCACCAGCCTCATTCTTGAACTTTCTTTTCAAAGTCCGAGAGTCAGGGTCACTC
AGCCCGGAGCACGGGCCCGTTGTGGTGCACTGCAGTGCAGGCATCGGCAGGTCTGGAACCTTCTGTCTGGCTGATACCTGCCTCTTGCTGATGGACAAGAGGAAAGACCCTTCTTCC
GTTGATATCAAGAAAGTGCTGTTAGAAATGAGGAAGTTTCGGATGGGGCTGATCCAGACAGCCGACCAGCTGCGCTTCTCCTACCTGGCTGTGATCGAAGGTGCCAAATTCATCATG
GGGGACTCTTCCGTGCAGGATCAGTGGAAGGAGCTTTCCCACGAGGACCTGGAGCCCCCACCCGAGCATATCCCCCCACCTCCCCGGCCACCCAAACGAATCCTGGAGCCACACAAT
GGGAAATGCAGGGAGTTCTTCCCAAATCACCAGTGGGTGAAGGAAGAGACCCAGGAGGATAAAGACTGCCCCATCAAGGAAGAAAAAGGAAGCCCCTTAAATGCCGCACCCTACGGC
ATCGAAAGCATGAGTCAAGACACTGAAGTTAGAAGTCGGGTCGTGGGGGGAAGTCTTCGAGGTGCCCAGGCTGCCTCCCCAGCCAAAGGGGAGCCGTCACTGCCCGAGAAGGACGAG
GACCATGCACTGAGTTACTGGAAGCCCTTCCTGGTCAACATGTGCGTGGCTACGGTCCTCACGGCCGGCGCTTACCTCTGCTACAGGTTCCTGTTCAACAGCAACACATAGCCTGAC
CCTCCTCCACTCCACCTCCACCCACTGTCCGCCTCTGCCCGCAGAGCCCACGCCCGACTAGCAGGCATGCCGCGGTAGGTAAGGGCCGCCGGACCGCGTAGAGAGCCGGGCCCCGGA
CGGACGTTGGTTCTGCACTAAAACCCATCTTCCCCGGATGTGTGTCTCACCCCTCATCCTTTTACTTTTTGCCCCTTCCACTTTGAGTACCAAATCCACAAGCCATTTTTTGAGGAG
AGTGAAAGAGAGTACCATGCTGGCGGCGCAGAGGGAAGGGGCCTACACCCGTCTTGGGGCTCGCCCCACCCAGGGCTCCCTCCTGGAGCATCCCAGGCGGGCGGCACGCCAACAGCC
CCCCCCTTGAATCTGCAGGGAGCAACTCTCCACTCCATATTTATTTAAACAATTTTTTCCCCAAAGGCATCCATAGTGCACTAGCATTTTCTTGAACCAATAATGTATTAAAATTTT
TTGATGTCAGCCTTGCATCAAGGGCTTTATCAAAAAGTACAATAATAAATCCTCAGGTAGTACTGGGAATGGAAGGCTTTGCCATGGGCCTGCTGCGTCAGACCAGTACTGGGAAGG
AGGACGGTTGTAAGCAGTTGTTATTTAGTGATATTGTGGGTAACGTGAGAAGATAGAACAATGCTATAATATATAATGAACACGTGGGTATTTAATAAGAAACATGATGTGAGATTA
CTTTGTCCCGCTTATTCTCCTCCCTGTTATCTGCTAGATCTAGTTCTCAATCACTGCTCCCCCGTGTGTATTAGAATGCATGTAAGGTCTTCTTGTGTCCTGATGAAAAATATGTGC
TTGAAATGAGAAACTTTGATCTCTGCTTACTAATGTGCCCCATGTCCAAGTCCAACCTGCCTGTGCATGACCTGATCATTACATGGCTGTGGTTCCTAAGCCTGTTGCTGAAGTCAT
TGTCGCTCAGCAATAGGGTGCAGTTTTCCAGGAATAGGCATTTGCCTAATTCCTGGCATGACACTCTAGTGACTTCCTGGTGAGGCCCAGCCTGTCCTGGTACAGCAGGGTCTTGCT
GTAACTCAGACATTCCAAGGGTATGGGAAGCCATATTCACACCTCACGCTCTGGACATGATTTAGGGAAGCAGGGACACCCCCCGCCCCCCACCTTTGGGATCAGCCTCCGCCATTC
CAAGTCAACACTCTTCTTGAGCAGACCGTGATTTGGAAGAGAGGCACCTGCTGGAAACCACACTTCTTGAAACAGCCTGGGTGACGGTCCTTTAGGCAGCCTGCCGCCGTCTCTGTC
CCGGTTCACCTTGCCGAGAGAGGCGCGTCTGCCCCACCCTCAAACCCTGTGGGGCCTGATGGTGCTCACGACTCTTCCTGCAAAGGGAACTGAAGACCTCCACATTAAGTGGCTTTT
TAACATGAAAAACACGGCAGCTGTAGCTCCCGAGCTACTCTCTTGCCAGCATTTTCACATTTTGCCTTTCTCGTGGTAGAAGCCAGTACAGAGAAATTCTGTGGTGGGAACATTCGA
GGTGTCACCCTGCAGAGCTATGGTGAGGTGTGGATAAGGCTTAGGTGCCAGGCTGTAAGCATTCTGAGCTGGGCTTGTTGTTTTTAAGTCCTGTATATGTATGTAGTAGTTTGGGTG
TGTATATATAGTAGCATTTCAAAATGGACGTACTGGTTTAACCTCCTATCCTTGGAGAGCAGCTGGCTCTCCACCTTGTTACACATTATGTTAGAGAGGTAGCGAGCTGCTCTGCTA
TATGCCTTAAGCCAATATTTACTCATCAGGTCATTATTTTTTACAATGGCCATGGAATAAACCATTTTTACAAAAATAAAAACAAAAAAAGCAAGGTGTTTTGGTATAATACCTTTT
CAGGTGTGTGTGGATACGTGGCTGCATGACCGGGTGGGTGGGGGGGAGTGTCTCAGGGTCTTCTGTGACCTCACAGAACTGTCAGACTGTACAGTTTTCCAACTTGCCATATTCATG
ATGGGTTTGCATTTTAGCTGCAACAATAAAATTTTTTTCTAAAGAACATGAAAAAAAAAAAAA

Results:

The beginning of the output file:

Score: 9.999665000418599e-09
A: UGGAGAAAGGUUCGUUAAACCAUCGGUCACCCAGCCCCUGGCCUGCUG
B: UGGAGAAAGGUUCGUUAAACCAGGCUGGGUGACCGAUGGUUUAACGAACCUUUCUCCAUC
 X: CAGCAGGCCAGGGGCUGGGUGACCGAUGG (215 - 243)

Score: 9.999613001825999e-09
A:UGGAGAAAGGUUCGUUAAAUGAGGUCUCAGGCAGCCACGGCUGAAGUU
B:UGGAGAAAGGUUCGUUAAAUGAGGCUGCCUGAGACCUCAUUUAACGAACCUUUCUCCAUC
X:AACUUCAGCCGUGGCUGCCUGAGACCUCA (325 - 353)

Score: 9.9993870000612e-09
A:UGGAGAAAGGUUCGUUAAAGCUCGCAGGAUGGAUAGGCAGGUGGACUU
B:UGGAGAAAGGUUCGUUAAAGCUCUAUCCAUCCUGCGAGCUUUAACGAACCUUUCUCCAUC
X:AAGUCCACCUGCCUAUCCAUCCUGCGAGC (361 - 389)

The scRNA with the best score (highest) is the first:

A: UGGAGAAAGGUUCGUUAAACCAUCGGUCACCCAGCCCCUGGCCUGCUG
B: UGGAGAAAGGUUCGUUAAACCAGGCUGGGUGACCGAUGGUUUAACGAACCUUUCUCCAUC

This scRNA targets X which starting at the 215th basis of the target gene up to the 243th:

X: CAGCAGGCCAGGGGCUGGGUGACCGAUGG (215 - 243)

Results Analysis

Analysis of the chosen scRNA

We chose the first scRNA, and verified the results by putting A and B together in an identical concentration, with and without the trigger gene.

OFF state (without trigger)

No concentration of B (0%) - this is a very good OFF state.
The duplex A·B is the only complex (100%) that formed from A and B:


Figure 3- The secondary structure of duplex A·B

ON state (with trigger)

There is a concentration of 50% of B and 50% of X·A. That is to say, all of the B forms will transform into hairpins, and will look like in the figure below:


Figure 4- The secondary structure of B

The A which has separated from B, assembles, along with X, the complex X·A:


Figure 5- The secondary structure of complex X·A

Conclusions:
We used the program we wrote in the design and test of our scRNA mechanism. It helped us to choose the scRNA we wanted to test in our lab, among a very large number of potential scRNAs, and is also even spared us the manual analysis of a single scRNA. We offer the code to of the program as open source, from the notion that it can help anyone who want to design the same mechanism we used, for other purposes. Of course the program is adjustable to the needs of any objective, with a different trigger gene and a different target gene.