$ RNA2Dfold -p < 2dfold.inp > 2dfold.out
The outputfile 2dfold.out should look like below, check it out using less.
CGUCAGCUGGGAUGCCAGCCUGCCCCGAAAGGGGCUUGGCGUUUUGGUUGUUGAUUCAACGAUCAC ((((((((((....)))))..(((((....))))).)))))...(((((((((...))))))))). (-30.60) ((((((((((....)))))..(((((....))))).)))))...(((((((((...))))))))). (-30.60) <ref 1> .................................................................. ( 0.00) <ref 2> free energy of ensemble = -31.38 kcal/mol k l P(neighborhood) P(MFE in neighborhood) P(MFE in ensemble) MFE E_gibbs MFE-structure 0 24 0.28019139 1.00000000 0.28019156 -30.60 -30.60 ((((((((((....)))))..(((((....))))).)))))...(((((((((...))))))))). 1 23 0.10926548 0.43037770 0.04702541 -29.50 -30.02 ((((((((((....)))))..(((((....))))).)))))....((((((((...)))))))).. 2 22 0.01833047 0.36607549 0.00671033 -28.30 -28.92 ((((.(((((....)))))..(((((....)))))..))))....((((((((...)))))))).. 2 24 0.00472936 0.53598082 0.00253484 -27.70 -28.08 ((((((((((....))))...(((((....)))))))))))...(((((((((...))))))))). 3 21 0.00214640 0.10357506 0.00022231 -26.20 -27.60 ((((.(((((....)))))..(((((....)))))..)))).....(((((((...)))))))... 3 23 0.00147662 0.28811146 0.00042543 -26.60 -27.37 ((((((((((....))))...(((((....)))))))))))....((((((((...)))))))).. [...]
For visualizing the output the ViennaRNA Package includes two scripts 2Dlandscape_pf.gri, 2Dlandscape_mfe.gri located in VRP/share/ViennaRNA/. gri (a language for scientific graphics programing) is needed to create a colored postscript plot. We use the partition function script to show the free energies of the distance classes (graph below, left):
$ gri ../Progs/VRP/share/ViennaRNA/2Dlandscape_pf.gri 2dfold.out
Compare the output file with the colored plot and determine the MFE minima with corresponding distance classes. For easier comparision the outputfile of RNA2Dfold can be sorted by a simple sort command. For further information regarding sort use the -help option.
$ sort -k6 -n 2dfold.out > sort.outNow we choose the structure with the lowest energy besides our startstructure, replace the open chain structure from our old input with that structure and repeat the steps above with our new values
The new projection (right graph) shows the two major local minima which are separated by 39 bp (red dots in figure below)
and both are likely to be populated with high probability. The landscape gives an estimate of
the energy barrier separating the two minima (about -20 kcal/mol).
The red dots mark the distance from open chain to the MFE structure respectively the
distance from the 2nd best structure to the MFE.
Sven Findeiss 2013-11-22