Netiquette · Download · News · Gallery · Homepage · DSSR Manual · G-quadruplexes · DSSR-Jmol · DSSR-PyMOL · DSSR Licensing · Video Overview· RNA Covers

Author Topic: Seeking help with DSSR (Base Pair and Stacking)  (Read 34332 times)

Offline mahfuz

  • non-commercial
  • with-posts
  • *
  • Posts: 4
    • View Profile
Seeking help with DSSR (Base Pair and Stacking)
« on: October 05, 2018, 05:49:11 pm »
Hi Dr. Xiang-Jun Lu,

I'm working on your annotation tool DSSR and trying to get the base-pair and stacking information from DSSR annotation. The base-pair part is clear enough but I'm facing problem for finding the stackings. The canonical and non-canonical stackings are not clear to me. I'm also trying to get the direction of the stackings (upword/downword/inward/outward).

Thanks

Mahfuz

Offline xiangjun

  • Administrator
  • with-posts
  • *****
  • Posts: 1650
    • View Profile
    • 3DNA homepage
Re: Seeking help with DSSR (Base Pair and Stacking)
« Reply #1 on: October 05, 2018, 10:11:08 pm »
Hi Mahfuz,

Quote
The canonical and non-canonical stackings are not clear to me.

Where did you get the concept of canonical vs non-canonical stackings from? What do you mean specifically? It is not a term used in DSSR.

In DSSR, pairwise stacking interactions are available via the --non-pair option, as documented in the -h option and the manual. Stacking interaction is detected and quantified by the overlap area of bases.

Quote
I'm also trying to get the direction of the stackings (upword/downword/inward/outward).

Check section "3.8 The --non-pair option" of the DSSR manual:

Quote
As in 3DNA [8], base-stacking is quantified as the area (in Å2) of the overlapped polygon defined by the two bases of the interacting nts, where the base atoms are projected onto the mean base plane9. In the output file, values in parentheses measure the overlap of base ring atoms only, and those outside parentheses include exocyclic atoms on the ring. Base-stacking interactions are classified into one of the following four categories: pm(>>,forward), mp(<<,backward), mm(<>,outward), and pp(><,inward). Here p and m represent the plus and minus faces of the base ring, as defined by the direction of the z-axis of the standard base reference frame. The symbols (>>, <<, <>, and ><) follow Major et al, except pm(>>) is called forward instead of upward, and mp(<<) backward instead of downward [32]. Moreover, the inter-base angle is reported; closer to zero means the two bases are nearly parallel.

Basically, each base has two faces (just as a coin). See Fig. 1 of the DSSR paper, attached below (also the 2003 3DNA paper). Each z-axis has a minus--->plus directionality, and there are 2x2=4 possible combinations of two bases in stack, as noted above.

"summary of methods to identify RNA structural components" title="summary of methods to identify RNA structural components"

HTH,

Xiang-Jun
« Last Edit: October 05, 2018, 10:13:55 pm by xiangjun »

Offline mahfuz

  • non-commercial
  • with-posts
  • *
  • Posts: 4
    • View Profile
Re: Seeking help with DSSR (Base Pair and Stacking)
« Reply #2 on: October 09, 2018, 12:16:56 pm »
Thank you for your quick response. I was just trying to compare the outputs from MC-Annotate, RNA-View and DSSR. There is something 'adjacent' and 'non-adjacent' stackings showing in MC-Annotate. I mistakenly wrote 'canonical and non-canonical' for these terms. I'm extremely sorry for my mistake.

Can you please help me to compare the 'adjacent and non-adjacent stacking' of MC-Annotate output with something in DSSR output?

Thanks

Offline xiangjun

  • Administrator
  • with-posts
  • *****
  • Posts: 1650
    • View Profile
    • 3DNA homepage
Re: Seeking help with DSSR (Base Pair and Stacking)
« Reply #3 on: October 09, 2018, 01:16:15 pm »
Quote
Can you please help me to compare the 'adjacent and non-adjacent stacking' of MC-Annotate output with something in DSSR output?

In DSSR output, the fact that stacked bases are covalently connected via a phosphodiester bond is noted by the keyword "connected". For example, for 1ehz (a tRNA), you'd see the following in DSSR output. For entry #68, the stacking between G51 and U52 is marked by connected, while the other two entries do not.

 # x3dna-dssr -i=1ehz.pdb --non-pair
List of 93 non-pairing interactions
......
36 A.A21    A.C48    stacking: 5.9(2.9)--mm(<>,outward) interBase-angle=5 min_baseDist=3.28
68 A.G51    A.U52    stacking: 6.8(4.0)--pm(>>,forward) interBase-angle=6 connected min_baseDist=3.22
71 A.G53    A.A62    stacking: 4.2(2.0)--mm(<>,outward) interBase-angle=8 min_baseDist=3.15


It is up to you to compare the DSSR output with that from MC-Annotate. You're welcome to post your finding here for the benefit of other viewers of the thread.

HTH,

Xiang-Jun

Offline mahfuz

  • non-commercial
  • with-posts
  • *
  • Posts: 4
    • View Profile
Re: Seeking help with DSSR (Base Pair and Stacking)
« Reply #4 on: October 12, 2018, 05:48:04 pm »
Thank you so much. I've another query. In the base-pair list I found the following symbols which I understand: cWW, tWW, cSH, tSH, cHS, tHS

But I cannot be able to get the meaning of c.H, t.H, c.S, t.S, c.W, t.W - what is the . (dot) here? Can you please explain me.

Sorry for asking so many questions.

Thanks.

Offline xiangjun

  • Administrator
  • with-posts
  • *****
  • Posts: 1650
    • View Profile
    • 3DNA homepage
Re: Seeking help with DSSR (Base Pair and Stacking)
« Reply #5 on: October 12, 2018, 06:40:11 pm »
You're welcome to ask any DSSR-related questions on the Forum. On the other hand, please respond to my follow-up questions, for the benefit of all parties involved.

Quote
But I cannot be able to get the meaning of c.H, t.H, c.S, t.S, c.W, t.W - what is the . (dot) here? Can you please explain me.

Please read the DSSR User Manual (specifically section 3.2.2 "Base pairs"). Let me know how the corresponding description can be improved.

HTH,

Xiang-Jun



Offline mahfuz

  • non-commercial
  • with-posts
  • *
  • Posts: 4
    • View Profile
Re: Seeking help with DSSR (Base Pair and Stacking)
« Reply #6 on: October 15, 2018, 12:29:15 pm »
I'm extremely sorry that I haven't seen it before. I think
Quote
The dot symbol denotes cases where edges or orientations cannot be defined (see bp #7 in the list of base pairs on Page 14, and Figure 5(D)).
is my answer.

 

Funded by the NIH R24GM153869 grant on X3DNA-DSSR, an NIGMS National Resource for Structural Bioinformatics of Nucleic Acids

Created and maintained by Dr. Xiang-Jun Lu, Department of Biological Sciences, Columbia University