Fundamentals of Sequence Analysis, 1998-1999
Problem set 9:  RNA folding.

If you get stuck, refer to the OpenVMS and GCG resources in the 
class home page.

References:
 
  See the GCG documentation.

Problem group 1.  Stem-loop structures

1A.  Are there any stem/loop structures in the tet gene of pBR322 that
     have at least 16 hydrogen bonds in the stem?


Problem group 2.  RNA folding

2A.  Does the tet RNA from pBR322 have a defined secondary structure?


2B.  Does the 3' end of the Drosophila bicoid mRNA have a defined
     structure?  (gb_in:X14458, 1550-2456)

Problem group 3.  RNA folding with pseudoknots
It is a good idea to rerun MFold and PlotFold for
sequences which may include pseudoknots.  However, doing so is
actually a little tricky.  The problem is that if you try to force
MFold to have pseudoknots, it will ignore one of the two
helices specified, unfortunately, without giving any warning that it
is doing so.  If you try to exclude these regions with /Closedexcise,
MFold will ignore that directive in some instances.  This is a
bug of some sort.  This leaves only the /prevent qualifier,
which may be used to force the pseudoknot regions into a large,
noninteracting loops. 

3A.  Use RNABOB to analyze gendocdata:alucons.seq, then run
     MFold and PlotFold to see how removing the
     pseudoknot regions affects the rest of the fold. 

3B.  Use PlotFold to output a connect file, manually edit
     it so that the pseudoknot links are present, then use 
     Squiggles to display it.  What happens?