Junction definition

[lijun]

Hi,

I find that dssr treats some pseudoknots as junctions. For example in rna 4gma:
   List of 3 junctions
   1 3-way junction: nts=8; [2,0,0]; linked by [#6,#13,#7]
     summary: [3] 2 0 0 [Z.74 Z.81 Z.201 Z.193 Z.194 Z.80] 2 4 5
     nts=8 GCGAUCGC 1..Z.G.74.,1..Z.C.75.,1..Z.G.200.,1..Z.A.201.,1..Z.U.193.,1..Z.C.194.,1..Z.G.80.,1..Z.C.81.
     nts=2 CG 1..Z.C.75.,1..Z.G.200.
     nts=0
     nts=0
The stem number [#6,#13,#7] is not in the ascending order. Can I filter out these false junctions by checking the order of the stems?

Thanks,
Jun Li

[xiangjun]

Hi,

I cannot reproduce your reported result based on information provided (i.e., PDB id: 4gma). Could you please be more specific by attached the PDB file, DSSR version, and options you used?

Thanks,

Xiang-Jun

[lijun]

Hi Xiang-Jun,

I used the option -loop-exclude-icpair.  In my previous post “Stems of junction structure have only one base pair”, You told me to use this option to exclude icPairs from loops.

Thanks,

Jun Li

Hi,

I cannot reproduce your reported result based on information provided (i.e., PDB id: 4gma). Could you please be more specific by attached the PDB file, DSSR version, and options you used?

Thanks,

Xiang-Jun

[xiangjun]

OK. I got it. I'll look into the issue.

Regarding the definition of (junction) loops, the DSSR paper says:

In DSSR, a loop forms a ‘closed’ circle with any two sequential nucleotides connected either by a phosphodiester linkage or a canonical base pair, and is specified by the lengths of consecutive bridging-nucleotide segments (Figures 1F and 2C).

The situation becomes less intuitive/obvious when pseudoknots are included. You could try the --nested option to removed pseudoknots to see if the results make more sense.

Best regards,

Xiang-Jun

[lijun]

--nested option works!
Thanks,
Jun Li

OK. I got it. I'll look into the issue.

Regarding the definition of (junction) loops, the DSSR paper says:

In DSSR, a loop forms a ‘closed’ circle with any two sequential nucleotides connected either by a phosphodiester linkage or a canonical base pair, and is specified by the lengths of consecutive bridging-nucleotide segments (Figures 1F and 2C).

The situation becomes less intuitive/obvious when pseudoknots are included. You could try the --nested option to removed pseudoknots to see if the results make more sense.

Best regards,

Xiang-Jun