In addition to VARNA, the draw program in the RNAstructure package from the Mathews Laboratory can also be used to depict DSSR-derived RNA secondary structures in connect table (.ct) format. The draw
program produces images in PostScript (or svg) format, in different styles from those generated by VARNA. Given below are a couple of examples on how to connect DSSR with draw
.
The secondary structure of the PDB entry 1msy in DSSR-derived .ct file is as below:
27 DSSR-derived secondary structure in '1msy' 1 U 0 2 0 2647 2 G 1 3 26 2648 3 C 2 4 25 2649 4 U 3 5 24 2650 5 C 4 6 23 2651 6 C 5 7 22 2652 7 U 6 8 0 2653 8 A 7 9 0 2654 9 G 8 10 0 2655 10 U 9 11 0 2656 11 A 10 12 0 2657 12 C 11 13 17 2658 13 G 12 14 0 2659 14 U 13 15 0 2660 15 A 14 16 0 2661 16 A 15 17 0 2662 17 G 16 18 12 2663 18 G 17 19 0 2664 19 A 18 20 0 2665 20 C 19 21 0 2666 21 C 20 22 0 2667 22 G 21 23 6 2668 23 G 22 24 5 2669 24 A 23 25 4 2670 25 G 24 26 3 2671 26 U 25 27 2 2672 27 G 26 0 0 2673
Let the DSSR-derived .ct file for 1msy be named 1msy.ct
, the following two draw-command runs will produce the secondary structure in PostScript (1msy.eps
) and svg (1msy.svg
) respectively.
draw 1msy.ct 1msy.eps draw 1msy.ct 1msy.svg --svg -n 1
The PDB entry 1ehz (yeast phenylalanine tRNA) has a pseudo knot, so the draw
program will create a ‘circularized’ structure as shown below:
Note the following two caveats:
- The
draw
program does not like the extra information starting with # after the 6 data columns of a .ct file produced by DSSR (by default). The program hangs at “Determining structure information…”. - For the DSSR-derived .ct file for the Dickerson DNA dodecamer,
draw
aborts with “Drawing structure 1…Segmentation fault: 11”.