2
« on: July 27, 2015, 01:56:54 pm »
Hello Xiang-Jun and congratulations for the NAR paper!
I came into a problem using DSSR on several small structures, I will show you 3HTS as example. The input is the following:
x3dna-dssr -i=/home//3HTS.py.pdb --segid --idstr=long --u-turn --more --non-pair --po4 --nested -o=3HTS.py.dssr.out
Here the fragment of the output file which causes me troubles:
****************************************************************************
List of 3 atom-base capping interactions
dv: vertical distance of the atom above the nucleotide base
-----------------------------------------------------------
type atom nt dv
1 sugar O4'@.B000.A.DG.1. ..A.DC.12. 2.84
2 sugar O4'@..A.DG.1. .A00Z.A.DC.12. 2.84
3 sugar O4'@.B000.A.DG.1. .B00A.A.DC.12. 2.84
****************************************************************************
I attached the output file "dssr-a2bases" to show you better the issue. In fact it seems that in the capping contact between the DG 1A and DC 12A there are some troubles related to the symmetry codes. If I open the file with PyMol generating symmetry mates, I find that the residue DG 1A with symmetry A00A is stacked over DC 12A, while in the dssr output the symmetric B000 is reported. I know it's a quite tricky matter and occurs only in small structures, but I still wanted to show it to you.
Thank you in advance for your answer
Luigi