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

Welcome > Feature requests

chain continuation character in analyze

(1/2) > >>

auffinger:
Hi Xiang-Jun,

Went just through a few options.
In the analyze file, you insert the '-' character when a strand is not broken and 'x' when its broken.
Then a new chain starts. May be you could add a third '+' character for these residues ?
This is a really minor point.

Also, for the same strand P...P and C1'...C1' distances, could you add two decimals instead of one?
This could be convenient for some applications.

Best,

Pascal

xiangjun:

--- Quote ---In the analyze file, you insert the '-' character when a strand is not broken and 'x' when its broken.
Then a new chain starts. May be you could add a third '+' character for these residues ?
--- End quote ---
Could you provide an example?


--- Quote ---Also, for the same strand P...P and C1'...C1' distances, could you add two decimals instead of one?
--- End quote ---
Done -- the distribution will be updated after clarification of the above point.

Xiang-Jun

auffinger:
Well, normally after a 'x' character a new chain starts, so this is easy. Yet putting a '+' instead of a '-' at these positions seems more informative to me although it looks like a detail. A '+' should then appear for the first residue of a structure. Hence, you would have a specific marker for all 3' and 5' residues.
Of course, an isolated residue can be found and you would have to chose between '+' and 'x' (I suggest '+').

  1   (0.012) ....>A:...1_:[..C]C-          <-- a '+' here instead of a '-'
   2   (0.019) ....>A:...2_:[.DC]C-
   3   (0.038) ....>A:...3_:[.DG]G-
   4   (0.030) ....>A:...4_:[.DG]G-
   5   (0.010) ....>A:...5_:[.DC]C-
   6   (0.043) ....>A:...6_:[.DG]G-
   7   (0.012) ....>A:...7_:[.DC]C-
   8   (0.013) ....>A:...8_:[.DC]C-
   9   (0.027) ....>A:...9_:[.DG]G-
  10   (0.027) ....>A:..10_:[..G]Gx
  11   (0.015) ....>B:..11_:[..C]C-       <-- a '+' here
  12   (0.009) ....>B:..12_:[.DC]C-

Pascal

xiangjun:
Please download the updated aug06 release of 3DNA v2.1beta. Now you can specify helix begin/continuation/end and isolated bp characters through option -chain_markers. For example,


--- Code: ---find_pair 1egk.pdb stdout
    #  default as before
find_pair -chain_markers='o|x+' 1egk.pdb stdout
    #  with helix beginning character assigned to 'o'
--- End code ---

The same option can also be applied to 'analyze'.

Also, as noted previously, the same strand P...P and C1'...C1' distances are now output in two decimals.

Have a try and report back how it goes.

Xiang-Jun

auffinger:
Thanks Xiang-Jun,

Works fine, when no options for find_pairs are used. Yet, I am using following options
find_pair -s -pdbv3 -attach=off  -chain_markers='o|x+' ...
and would like to see these markers in the output of analyze

analyze -pdbv3 file.fp out

that is like

...
   3   (0.038) ....>A:...3_:[.DG]G-
   4   (0.030) ....>A:...4_:[.DG]G-
   5   (0.010) ....>A:...5_:[.DC]C-
   6   (0.043) ....>A:...6_:[.DG]G-
...

Sorry if I didn't explain myself well on this point.

Best,

Pascal

Navigation

[0] Message Index

[#] Next page

Created and maintained by Dr. Xiang-Jun Lu [律祥俊] (xiangjun@x3dna.org)
The Bussemaker Laboratory at the Department of Biological Sciences, Columbia University.

Go to full version