Netiquette · Download · News · Gallery · G-quadruplexes · DSSR-Jmol · DSSR-PyMOL · Video Overview · DSSR v2.9.0 (DSSR Manual) · Homepage · HIRING!
# a standard B-DNA model x
ATOM 1 P DA A 1 0.621 9.421 -1.028 1.00 0.00 P
ATOM 2 O1P DA A 1 0.606 10.783 -1.583 1.00 0.00 O
ATOM 3 O2P DA A 1 1.439 9.186 0.166 1.00 0.00 O
ATOM 4 O5' DA A 1 -0.863 8.945 -0.676 1.00 0.00 O
ATOM 5 C5' DA A 1 -1.693 8.405 -1.689 1.00 0.00 C
ATOM 6 C4' DA A 1 -2.528 7.255 -1.154 1.00 0.00 C
# a mirror-image, with x-coordinates negated
ATOM 1 P DA A 1 -0.621 9.421 -1.028 1.00 0.00 P
ATOM 2 O1P DA A 1 -0.606 10.783 -1.583 1.00 0.00 O
ATOM 3 O2P DA A 1 -1.439 9.186 0.166 1.00 0.00 O
ATOM 4 O5' DA A 1 0.863 8.945 -0.676 1.00 0.00 O
ATOM 5 C5' DA A 1 1.693 8.405 -1.689 1.00 0.00 C
ATOM 6 C4' DA A 1 2.528 7.255 -1.154 1.00 0.00 C
Note: a helix is defined by base-stacking interactions, regardless of bp
type and backbone connectivity, and may contain more than one stem.
helix#number[stems-contained] bps=number-of-base-pairs in the helix
bp-type: '|' for a canonical WC/wobble pair, '.' otherwise
helix-form: classification of a dinucleotide step comprising the bp
above the given designation and the bp that follows it. Types
include 'A', 'B' or 'Z' for the common A-, B- and Z-form helices,
'.' for an unclassified step, and 'x' for a step without a
continuous backbone.
--------------------------------------------------------------------
helix#1[1] bps=40
strand-1 5'-AAAAAAAAAATTTTTTTTTTCCCCCCCCCCGGGGGGGGGG-3'
bp-type ||||||||||||||||||||||||||||||||||||||||
strand-2 3'-TTTTTTTTTTAAAAAAAAAAGGGGGGGGGGCCCCCCCCCC-5'
helix-form .......................................
1. Is my understanding correct that DSSR can identify the geometry but does not recognize the form of our L-DNA?
2. Is there any way to make DSSR recognize the form of our L-DNA?
helix#1[1] bps=40
strand-1 5'-AAAAAAAAAATTTTTTTTTTCCCCCCCCCCGGGGGGGGGG-3'
bp-type ||||||||||||||||||||||||||||||||||||||||
strand-2 3'-TTTTTTTTTTAAAAAAAAAAGGGGGGGGGGCCCCCCCCCC-5'
helix-form LLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL helix#1[3] bps=9
strand-1 5'-gcgugugug-3'
bp-type |||||||||
strand-2 3'-cgcacgcac-5'
helix-form LLxL.xLx
x3dna-dssr -i=4wb2.pdb --json | jq '.nts[] | {nt_id, is_L_sugar}'......
{
"nt_id": "D.0G35",
"is_L_sugar": true
}
{
"nt_id": "D.0C36",
"is_L_sugar": true
}
......Shear : no sign inversion
Stretch : no sign inversion
Stagger : sign inversion
Buckle : sign inversion
Propeller : sign inversion
Opening : no sign inversionShift : no sign inversion
Slide : sign inversion
Rise : no sign inversion
Tilt : sign inversion
Roll : no sign inversion
Twist : sign inversionx-displacement : no sign inversion
y-displacement : sign inversion
helical rise : no sign inversion
inclination : sign inversion
tip : no sign inversion
helical twist : sign inversionHere which makes me feel confused is that the sign inversion pattern of the base pair parameters is different from that of the step and helical parameters. I understand that if the inversion pattern of bp parameters apply to step/helical bp parameter, such as Twist and Helical Twist will be "no sign inversion" which is impossible since the most distinct difference between native and mirror-image double helix is right and left handed. I looked into X3DNA and also SCHNAaP paper which has detailed description about how to calculate the parameters which makes believe the sign inversion pattern should be same . Probably DSSR has some certain axis direction inversion when calculate step/helical bp parameters which I overlooked.# x3dna-dssr -i=1BNA_L.pdb --more
helix#1[1] bps=12
strand-1 5'-CGCGAATTCGCG-3'
bp-type ||||||||||||
strand-2 3'-GCGCTTAAGCGC-5'
helix-form LLLLLLLLLLL
helical-rise: 3.35(0.44)
helical-radius: 9.37(0.98)
helical-axis: 0.088 -0.275 -0.958
point-one: -16.795 25.817 26.180
point-two: -13.557 15.707 -9.067
1 A.DC1 B.DG24 C-G WC 19-XIX cWW cW-W
bp1-pars:* [-0.42 -0.27 -0.06 -2.76 14.18 -3.67]
step-pars: [0.36 0.15 3.52 -3.40 -6.44 -40.31]
heli-pars: [0.55 0.11 3.52 9.26 -4.88 -40.94]
bp2-pars:* [-0.02 -0.27 -0.25 4.45 10.83 -4.02]
C1'-based: rise=3.52 twist=-38.57
C1'-based: h-rise=3.52 h-twist=-39.23
...
Summary of structural features of 24 nucleotides
Note: the first five columns are: (1) serial number, (2) one-letter
shorthand name, (3) dbn, (4) id string, (5) rmsd (~zero) of base
ring atoms fitted against those in a standard base reference
frame. The sixth (last) column contains a comma-separated list of
features whose meanings are mostly self-explanatory, except for:
turn: angle C1'(i-1)--C1'(i)--C1'(i+1) < 90 degrees
break: no backbone linkage between O3'(i-1) and P(i)
1 C ( A.DC1 0.020 anti,~C3'-endo,canonical,non-pair-contact,helix-end,stem-end
2 G ( A.DG2 0.025 anti,~C3'-endo,canonical,non-pair-contact,helix,stem
3 C ( A.DC3 0.024 anti,canonical,non-pair-contact,helix,stem
4 G ( A.DG4 0.023 anti,~C3'-endo,canonical,non-pair-contact,helix,stem
5 A ( A.DA5 0.016 anti,~C3'-endo,canonical,non-pair-contact,helix,stem
................
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