3DNA Forum

Questions and answers => General discussions (Q&As) => Topic started by: wolson on February 22, 2012, 12:23:05 pm

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

Title: These distinctly different structures give the same base-pair step parameters.
Post by: wolson on February 22, 2012, 12:23:05 pm
According to their 3DNA output files, all four of the attached structures have the same positive rise (+3.4 A), the PP and MM files have the same positive twist (+36 deg), and the PM and MP files have the same negative twist (-36 deg). The structures should have all four possible combinations of the signs for rise and twist, i.e., (+3.4, +36) for PP, (–3.4, -36) for MM, (+3.4, -36) for PM, and (-3.4, +36) for MP.

The error is related to the implementation of the Cambridge Convention in the computation of base-pair step parameters in 3DNA. Users interested in structures like these are interested in rigid-body parameters that reproduce the input structures rather than adhering to the Convention.

A useful option for such users would be the option NOT to follow the Cambridge Convention.
Title: Re: These distinctly different structures give the same base-pair step parameters.
Post by: xiangjun on February 22, 2012, 11:16:30 pm
Hi Wilma,

Thanks for bringing up this "issue". As always, the four concrete PDB files helped clarify everything. In short, 3DNA is behaving properly: because PP and MM are IDENTICAL (RMSD=0 Å), so are PM and MP (RMSD=0.0174364 Å). There are only two structures; the PP/MM pair is right-handed with Rise=+3.4Å, and Twist=+36°, whilst the PM/MP pair is left-handed Rise=+3.4Å, and Twist=–36°.

Now let's get into details to see why PP=MM, and PM=MP.

I do not quite understand how Srini and you come to the conclusion that 3DNA is in error here. Unless I am missing something obvious, it is hard for me to imagine that simply rotate a DNA structure by 180 degrees about the y-axis should reverse its Rise and Twist. Maybe Srini can shed more light on his thought?

Xiang-Jun
Title: Re: These distinctly different structures give the same base-pair step parameters.
Post by: wolson on February 23, 2012, 09:54:41 am
Hi, Xiang-Jun.

While the structures may appear to be identical at the level of successive base pairs, the backbone connections/directions differ. The direction is clear from the color-coding of bases and strands in the attached images. Look, for example, at the different pathways of the red, blue, yellow, green base sequence on the red strand.

Best,
w
 
Title: Re: These distinctly different structures give the same base-pair step parameters.
Post by: xiangjun on February 23, 2012, 10:13:02 am
Then we may have different understanding as to what it means to be of the same structure: to me, since PP and MM have an RMSD of 0, they are IDENTICAL. Naturally, 3DNA should output the SAME parameters, as it does. As mentioned in my previous post, the case of PM vs MP follows the same argument.

Ask Sriri to show here in details most possible, how and why PP and MM are different.

Xiang-Jun
Title: Re: These distinctly different structures give the same base-pair step parameters.
Post by: xiangjun on February 23, 2012, 04:30:38 pm
Hurrying for a meeting this morning, I did not notice the new set of 4 structures you attached when I posted my previous reply. Now things are becoming quite interesting. My arguments with regard to the first set of four structures you provided (Srini_PP.pdb, Srini_MM.pdb, Srini_PM.pdb and Srini_MP.pdb) still hold, i.e., PP=MM, and PM≈MP. There is no such things as "four distinctly different structures" there.

The new set of four structures (Srini_pp1.pdb, Srini_mm1.pdb, Srini_pm1.pdb, Srini_mp1.pdb), hereafter referred as PP1, MM1, PM1 and MP1, are completely different from the first set. 3DNA has no problem in identifying the four distinct forms, based on exactly the same algorithm as described in the 1997 JMB SCHNAaP paper (http://www.sciencedirect.com/science/article/pii/S0022283697913462). For example, for MM1, the output from 3DNA is as below:

Code: [Select]
****************************************************************************
Structure classification:

This is a right-handed unknown R-form structure
****************************************************************************

More specifically, Figure 7 of the SCHNAaP paper answers this question:

(http://ars.sciencedirect.com/content/image/1-s2.0-S0022283697913462-gr7.jpg)
Quote
Figure 7. A representation of four possible arrangements for antiparallel nucleic acid duplexes. Left-handed W and Z-DNA are shown on the left (the characteristic zig-zag backbone pattern is not represented for simplicity). Right-handed A/B and hypothetical R-DNA are shown on the right. The Twist free ladder forms are shown in the middle column. In the top row, the minor groove faces the viewer, while in the bottom row, the major groove faces the viewer. The SCHNAaP coordinate system is also shown. These structures were generated using SCHNArP (see accompanying paper) with Twist= ±36° (0° for the ladder forms), Rise=3.34 Å, and all other step parameters are set to zero. Color scheme: the minor groove side, dark green; the major groove side, light green; and the backbone, red.

In connection with the new set of 4 structures, they correspond to the four forms classified in 3DNA (SCHNAaP) as below:
Also as in SCHNAaP, right-handed structures (PP1/MM1) have positive Twist, and left-handed structures (PM1/MP1) have negative Twist. Moreover, they all have positive Rise.

I am really pleased to see the model structures representing the 4 possible distinct forms of double helices. I will consider to include them in future releases of the 3DNA distribution.

Xiang-Jun

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