Netiquette · Download · News · Gallery · Homepage · DSSR Manual · G-quadruplexes · DSSR-Jmol · DSSR-PyMOL · DSSR Licensing · Video Overview· RNA Covers
Questions and answers > DSSR-NAR paper
Figure 4 -- analysis of the env22 twister ribozyme (4rge)
(1/1)
xiangjun:
--- Quote ---Figure 4: DSSR discloses complexity in the folding of the env22 twister ribozyme not apparent in the two-armed tertiary structure (chain A, PDB id: 4rge (43)). (A) The software automatically detects the long helical arm with five coaxially stacked stems and the short single-stemmed arm of the molecule. Failing to account for the pseudoknots within the structure leads to a characterization of the molecule very different from its real organization. When pseudoknots are omitted, the RNA appears to form a simplified [2,1,3] three-way junction as shown in both planar (B) and linear (C) secondary structure diagrams. In reality, the DSSR-derived dot-bracket notation points to a double-pseudoknotted structure (D) with two types of brackets distinguishing the pseudoknotted pairs (matched [] and {}), and uncovers a novel [4,2,2,0,1,3,0,0,1,1] ten-way junction loop (D,E). The junction, which can be traced by following the arrows along the red arcs and bases (starting from U3, marked with *) in D, contains both ends of four of the six stems and follows a supercoiled pathway in 3D (Supplementary Figure S5). In contrast, without consideration of pseudoknots (F), the junction forms a simple relaxed circle (Supplementary Figure S5). DSSR also detects three previously ignored base pairs that help to anchor the consecutive A-minor motifs reported in the literature (43) (G). U41 pairs with A42 and A43 through bifurcated hydrogen bonding, as well as with A26 (Supplementary Figure S4C,D). Moreover, U41 and A42 constitute a UpA dinucleotide platform, and in combination with G25 and A26, create a unique network of eight interacting nucleotides (G). All eight nucleotides are involved in the ten-way junction loop (labeled red in (E)).
--- End quote ---
Here is the tarball (fig4-twister-ribozyme-4rge.tar.gz) with the script and all related data files.
The content of the full script (named tasks) is shown below. Please see also notes for "Figure 2 -- analysis of the yeast phenylalanine tRNA (1ehz)".
--- Code: Bash ---# Step #1 -- reorient the twister ribozyme verticallypdb_frag A 1:56 4rge.pdb 4rge-A.pdbx3dna-dssr -i=4rge-A.pdb -o=4rge-A.out --more --prefix=4rge-A # To get the result illustrated in panel D, load '4rge-A-2ndstrs.ct'# or '4rge-A-2ndstrs.dbn' into VARNA to draw the linear secondary# structure diagram, exported as .svg for annotation in Inkscape. # Extract the two helical axes from 4rge-A.out to file: 4rge-A.rot1# then reorient the structure vertically: 4rge-A.rot2rotate_mol -t=4rge-A.rot1 4rge-A.pdb 4rge-A-rot1.pdbrotate_mol -r=4rge-A.rot2 4rge-A-rot1.pdb 4rge-A-ok.pdb # Step #2 -- get the cartoon-block representation with the two# ls-fitted helical axes.x3dna-dssr -i=4rge-A-ok.pdb --helical-axis -o=temp\mv dssr-helicalAxes.pdb 4rge-A-ok-helices.pdbx3dna-dssr -i=4rge-A-ok.pdb --block-file -o=4rge-A-ok-blocks.r3d # Step #3 -- simplified representation of the [4,2,2,0,1,3,0,0,1,1]# 10-way junction in 3D -- panel E# note the '--raw-xyz' option: it keeps the original coordinatesx3dna-dssr -i=4rge-A-ok.pdb --raw-xyz --simple-junction -o=temp\mv dssr-simplifiedJcts.pdb 4rge-A-ok-jct.pdb# see file: 4rge-A-ok-jct.pmlpymol -qkc 4rge-A-ok-jct.pmlconvert -trim +repage -border 10 -bordercolor white 4rge-A-ok-jct-pymol.png 4rge-A-ok-jct.png # see file: 4rge-A-ok-full.pml (cartoon-block with the schematic# junction overlaid) -- panel Apymol -qkc 4rge-A-ok-full.pmlconvert -trim +repage -border 10 -bordercolor white 4rge-A-ok-full-pymol.png 4rge-A-ok-full.png # Step #4 -- remove pseudoknots to get a fully nested structure. It now# has only a [2,1,3] 3-way junction -- panels B, C, and F\cp 4rge-A-ok.pdb 4rge-nested.pdbx3dna-dssr -i=4rge-nested.pdb --nested --raw-xyz --simple-junction --prefix=4rge-nested -o=4rge-nested.out # The planar (panel B) and linear (panel C) secondary structure# diagrams are produced by loading '4rge-nested-2ndstrs.ct' or# '4rge-nested-2ndstrs.dbn' into VARNA, exported as .svg, and# annotated with Inkscape. # see file: 4rge-nested-jct.pml -- panel Fpymol -qkc 4rge-nested-jct.pmlconvert -trim +repage -border 10 -bordercolor white 4rge-nested-jct-pymol.png 4rge-nested-jct.png # Step #5 -- bifurcated U-A pairs in a network of 8 nucleotidespdb_frag A 13:14 A 25:26 A 36 A 41:43 4rge-A-ok.pdb 4rge-bifurcated.pdbx3dna-dssr -i=4rge-bifurcated.pdb --block-file -o=4rge-bifurcated.r3d# see file: 4rge-bifurcated.pml -- panel Gpymol -qkc 4rge-bifurcated.pmlconvert -trim +repage -border 10 -bordercolor white 4rge-bifurcated-pymol.png 4rge-bifurcated.png
Here are the images generated from the above script:
Navigation
[0] Message Index
Created and maintained by Dr. Xiang-Jun Lu [律祥俊] (xiangjun@x3dna.org)
The Bussemaker Laboratory at the Department of Biological Sciences, Columbia University.
