Messages

1126

RNA structures (DSSR) / Re: list nucleotide/nucleotide contacts involving a phosphate group.

« on: April 21, 2013, 06:23:55 pm »

Hi Pascal,

I've just released DSSR beta-r09-on-20130421 which contains a new option -non-pair (or --non-pair) to detect/output non-pairing interactions, including H-bonds and base stacking. As an example, see below for PDB entry 1msy which contains a GNRA tetra-loop.

x3dna-dssr --non-pair -i=1msy -o=1msy.out

The output file '1msy.out' contains the following:

****************************************************************************
List of 12 non-pairing interaction(s)
   1 A.G2648          A.G2673         base-overlap-area=2.0   
       H-bonds[0]: ""
   2 A.U2650          A.G2671         base-overlap-area=0.1   
       H-bonds[0]: ""
   3 A.C2652          A.G2669         base-overlap-area=0.2   
       H-bonds[0]: ""
   4 A.A2654          A.U2656         base-overlap-area=3.7   
       H-bonds[1]: "O4'*O4'[3.05]"
   5 A.G2655          A.G2664         base-overlap-area=4.4   
       H-bonds[1]: "O2'(hydroxyl)-O6(carbonyl)[3.09]"
   6 A.G2655          A.A2665         base-overlap-area=0.0   
       H-bonds[3]: "N1(imino)-OP2[2.77]; N2(amino)-OP2[3.34]; N2(amino)-O5'[2.89]"
   7 A.U2656          A.G2664         base-overlap-area=0.0   
       H-bonds[2]: "OP2-N1(imino)[3.04]; OP2-N2(amino)[2.94]"
   8 A.A2657          A.A2665         base-overlap-area=3.7   
       H-bonds[0]: ""
   9 A.G2659          A.A2661         base-overlap-area=0.0   
       H-bonds[1]: "O2'(hydroxyl)-N7[2.60]"
  10 A.G2659          A.G2663         base-overlap-area=3.9   
       H-bonds[0]: ""
  11 A.U2660          A.A2661         base-overlap-area=7.5   
       H-bonds[0]: ""
  12 A.A2661          A.A2662         base-overlap-area=6.3   
       H-bonds[0]: ""
****************************************************************************

Please let me know how you'd like to revise the content/format. As always, concrete examples work the best.

Xiang-Jun

1127

General discussions (Q&As) / Re: 3DNA installation problems on Mac OS X

« on: April 17, 2013, 04:40:38 pm »

Glad to hear that you've fixed the installation problems. Could you share what went wrong, and how you solved the issues? Some instructions are likely to be unclear in the post "How to install 3DNA on Linux and Windows?" that need to be improved.

Thanks for using 3DNA and posting on the forum.

Xiang-Jun

1128

General discussions (Q&As) / Re: BI and BII conformations

« on: April 12, 2013, 12:34:45 pm »

I am glad to see the happy . This thread help to illustrate clearly that a concrete example is the most effective and unambiguous way to get a (technical) point across. So in the future, whenever you have a 3DNA-related question, please be specific, and do not hesitate to post on the forum.

Xiang-Jun

1129

General discussions (Q&As) / Re: BI and BII conformations

« on: April 12, 2013, 11:44:05 am »

if epsilon = 0 to +360 and zeta= 0 to +360 e-z = -360     0        +360, but how can get negative values in this way?

Well, I am a bit confused by your reply too. If epsilon=60, zeta=160, then e-z=60-160=-100, which is NEGATIVE, right?

That's why I asked you to provide some concrete examples to walk through. Okay, let's use PDB  id 355d as an example, if you run the command:

Code: [Select]

analyze -tor=355d.tor 355d.pdb, you will have the following in file 355d.tor:

Code: [Select]

              base      chi A/S     alpha    beta   gamma   delta  epsilon   zeta     e-z BI/BII
------------------------------------------------------------------------------------------------
  10 A:..10_:[.DG]G   -83.6 anti    -60.3   163.2    39.5   143.2  -100.0   146.3   113.6  BII
  11 A:..11_:[.DC]C  -112.8 anti    -73.1   144.3    50.8   143.5  -164.4  -126.1   -38.3  BI
For nt C11, epsilon=-164.4, corresponding to -164.4+360=195.6; whilst zeta=-126.1, corresponding to -126.1+360=233.9; 195.6-233.9=-38.3 which is the value reported in 3DNA (see above). Since -38.3 < 0, 360 is added to it: -38.3+360=321.7, which is out of [20, 200], so it is assigned BI. Please work out the numbers for G10, and report back here.

Does this clarify your confusion?

Xiang-Jun

1130

General discussions (Q&As) / Re: BI and BII conformations

« on: April 11, 2013, 08:54:36 am »

Thanks for your followup. After checking the source code carefully and a sample 3DNA output, I noticed that my previous reply is incomplete (inaccurate)  . Thus the discrepancy: you found that e-z can be negative in analyze -tor output, while I said that it's in range [0, 360] which should always be positive.

The missing piece is that the 3DNA reported difference (d=e-z) is the raw value, which can be positive or negative. Only for classifying a BI/BII conformation, negative d value is added 360 to make it positive; this simplify the code since the continuous range 20, 200 is checked for BII conformation.

Does this solve your puzzle? If still not, please post some concrete examples so we can walk them through to get to the bottom of the issue.

Xiang-Jun

1131

General discussions (Q&As) / Re: BI and BII conformations

« on: April 10, 2013, 05:16:41 pm »

It is a good question; the short answer is that the range of e-z is in [0, 360].

In 3DNA, the backbone BI/BII classification is based on e-z and works as below:

• It follows the definition as given by IMB Jena, "Nucleic acid backbone parameters". Specifically,

  BI:  (epsilon-zeta) = -160° ... +20°
  BII: (epsilon-zeta) = +20° ... +200°

• In 3DNA, all torsion angles are given in the range of [-180, +180]. To calculate e-z, epsilon and zeta are first converted to [0, 360]. If the difference (d=e-z) is < 0, then d=d+360. So d should be in [0, 360]. If d is in [20, 200], it is classified as BII, otherwise, BI.

Does that make sense? If you know of any more elaborate definition of the BI/BII conformations, please let me know.

Xiang-Jun

1132

RNA structures (DSSR) / Re: list nucleotide/nucleotide contacts involving a phosphate group.

« on: April 05, 2013, 11:12:40 am »

Hi Pascal,

Thanks for your request for adding non-pairing nucleotide-nucleotide contacts (nt-nt) in 3DNA output. Currently, there is no such an option in 3DNA, but I do have this topic in mind for quite a while. It won't be hard to implement this functionality in DSSR/3DNA, so I will probably be able to get something done next week.

While we are at it, in addition to non-pairing nt-nt interactions via H-bond, how about those with only base-stacking? What output format would be prefer?

Xiang-Jun

1133

MD simulations / Re: Base stacking from x3dna_ensemble

« on: April 01, 2013, 12:06:03 pm »

theta  = arccos (dot(n1,n2)) where n1 and n2 are normal vectors from the b1n file

I wonder whether the angle I am getting is in radians?

You are getting the angle in radians instead of degrees. To verify, let n1 = [1, 0, 0]; n2 = [0, 0, 1], the angle should be 90 degrees.

Xiang-Jun

1134

General discussions (Q&As) / Re: zero or negative helical rise?

« on: March 28, 2013, 03:36:17 pm »

DSSR does not calculate helical/step parameters (at least not yet), as shown clearly in the enclosed file for 1xvk in my previous reply. See also my reply to "large deviations in output values".

In 3DNA distribution, there is also a program called 'cehs' (as in SCHNAaP) which uses RC8--YC6 as the bp long axis, and mean bp normal as the z-axis. In dinucleotide steps with non-canonical bps, 'cehs' normally provides a more 'sensible' rise/twist values. Please have a try on 1xvk and report back how it goes.

Xiang-Jun

1135

General discussions (Q&As) / Re: large deviations in output values

« on: March 27, 2013, 12:54:49 pm »

There are large deviations from average values in some of the parameters. Are the values real? What is the best way to view them in structures if so ?

The short answer is yes, they are 'real', as would be expected from 3DNA. The large deviations of 3DNA parameters for 3sj2 (and many other PDB entries) are due to non-canonical base pairs (bp). In the case of 3sj2, they are the three G+G pairs.

3DNA adopts the standard base reference frame, which is based on what would be a perfectly planar Watson-Crick (WC) bp geometry. Thus, by definition, WC pairs have the six bp parameters (Shear, Stretch, Stagger, Buckle, Propeller, Opening) close to zeros, with certain variations, mostly in Stagger, Buckle,  and Propeller. Non-canonical bps can have numerous ways to deviate from a WC bp; however, whatever the case, they can be rigorously quantified by the six rigid-body bp parameters. See the 3DNA NAR03 paper for details.

A non-cannonical bp not only has a set of characteristic bp parameters (see below for G+G bps in 3sj2), it could also greatly effect the (middle) bp reference frame used to derive the bp step and helical parameters. That's why you see large deviations of those parameters from normal WC steps (see below).

     bp        Shear    Stretch   Stagger    Buckle  Propeller  Opening
    3 G-U      -2.33     -0.47      0.13      2.30    -12.27      1.61
    4 C-G       0.22     -0.06     -0.00      3.70    -18.31      2.91
    5 C-G       0.22     -0.10      0.10      0.15     -4.59     -0.12
    6 G+G      -1.48     -3.61     -0.15     11.94      4.31     87.72

    step       Shift     Slide      Rise      Tilt      Roll     Twist
   4 CC/GG     -0.61     -1.97      3.29     -1.87      9.30     30.50
   5 CG/GG      0.04     -3.25     -1.34   -170.79     31.71    160.42
   6 GG/CG     -0.47     -3.69     -3.19    128.85   -110.61     97.02

3DNA derives a complete set of bp parameters to rigorously characterize the relative base geometry. For example, run the following 3DNA commands to see how 'analyze' and 'rebuild' complement each other to illustrate the point:

Code: [Select]

find_pair 3sj2.pdb stdout | analyze stdin
rebuild -atomic bp_step.par 3sj2-3dna.pdb
# superimpose '3sj2-3dna.pdb' onto '3sj2.pdb' using only base atoms, the rmsd would be ~0
find_pair 3sj2-3dna.pdb stdout | analyze stdin
# compare '3sj2-3dna.out' and '3sj2.out', the bp parameters are virtually identical

I see it as an advantage for 3DNA to report those 'weird' parameters, as it would (should) draw a user's attention. These bps and steps should be treated separately from the normal variations in structures (fragments) consisting of only WC bps. There could be ad hoc ways to made 'weird' values look normal, but they lack rigor and consistency. Again, see the 3DNA NAR03 paper for more info.

Note also that for a relatively straight duplex structure like 3sj2, 3DNA also output a set of "Global parameters based on C1'-C1' vectors:" as shown below, which you may find useful:

Code: [Select]

disp.: displacement of the middle C1'-C1' point from the helix
angle: inclination between C1'-C1' vector and helix (subtracted from 90)
twist: helical twist angle between consecutive C1'-C1' vectors
rise:  helical rise by projection of the vector connecting consecutive
       C1'-C1' middle points onto the helical axis

     bp       disp.    angle     twist      rise
   1 G-C      8.82     12.28     28.16      3.19
   2 G-C      8.09     11.46     32.47      2.62
   3 G-U      6.75     11.40     28.94      2.40
   4 C-G      7.06      8.02     31.16      2.87
   5 C-G      7.45      8.60     33.16      2.88
   6 G+G      7.07      7.95     29.95      2.72
   7 G-C      7.20      8.14     30.17      3.24
   8 C-G      6.86      9.94     34.81      2.83
   9 G+G      6.13     10.17     28.18      2.62
  10 G-C      6.56      6.75     33.11      2.91
  11 C-G      6.71      6.30     30.30      2.70
  12 G+G      6.05      9.99     33.37      2.80
  13 G-C      6.48     10.27     29.88      3.13
  14 G-C      6.82      7.48     27.46      2.97
  15 U-G      7.18      5.34     38.30      2.79
  16 C-G      7.16     10.13     26.27      2.90
  17 C-G      7.17     12.43      ---       ---
HTH,

Xiang-Jun

1136

RNA structures (DSSR) / Re: DSSR output - Base pair characteristics

« on: March 26, 2013, 01:33:06 pm »

Hi Jose,

Thanks for trying out DSSR and for your kind comment about the program. User feedback like yours is a great incentive for me to make DSSR a better tool to serve the RNA structure community.

To start with, I have tried hard to made DSSR easy to set up and play with. Based on my reading of literature in structural biology, I've made the DSSR output more intuitive (compared to previous 3DNA programs). For example, A.U2647 means U2647 on chain A. So far, I have not heard of any installation problem yet, and I am glad that you can make sense of most items in the DSSR output.

Your question regarding the meaning of the last column is well expected. It represents my own notation to specify a base pair, as elaborated below:

• Each base has three edges: W for the Watson-Crick edge, M for the major groove edge, and m for the minor groove edge. M corresponds to the Hoogsteen (or C-H) edge of the Leontis-Westhof nomenclature, and for the majority of cases (where the glycosidic bond is anti) m agrees with the 'sugar' edge. Note that in DSSR, the edges are defined purely on the geometry of the base plane as would be in a Watson-Crick base pair, and it is not related to sugar. See my post "The chi (χ) torsion angle characterizes base/sugar relative orientation". The DSSR definition applies to RNA as well as DNA, with either syn or anti glycosidic bond.
• In some boundary cases, the two bases in a pair may not be directly interacting edge-to-edge, where it is not straightforward to clearly designate which edge is involved. This is where the '.' comes in.
• The DSSR notation contains 4 characters of the pattern: [ct][WMm.][+-][WMm.]. The third position is either '+' or '-', and it designates the relative orientation of the two bases (flipped or normal) as has been consistently used in 3DNA. For example, see the difference in A+U Hoogsteen pair vs. A-U Watson-Crick pair.
• The first position is either 'c' for cis and 't' for trans of the two glycosidic bonds. It is defined by the 'virtual' torsion angle tor(N1-C1'-C1'-N9) reported in the DSSR output.

HTH,

Xiang-Jun
PS. Does "PDBID 3OXO" correspond to an RNA structure?

1137

General discussions (Q&As) / Re: zero or negative helical rise?

« on: March 25, 2013, 04:20:39 pm »

The issue you noticed with negative rise in 1xvk etc is due to the extensive non-canonical base pairs. Under such circumstance, as in more commonly seen in RNA structures, the meaning of base-pair step and helical parameters may not make (much) intuitive sense. Yet, these parameters are required to rigorously characterize the structure. Note that the base-pair parameters (Shear, Stretch, Stagger, Buckle, Propeller and Opening) still have their normal interpretation.

I'd recommend you use DSSR, which shows clearly the two Hoogsteen pairs, among other things, as show below. Check also files dssr-torsions.dat, dssr-pairs.pdb etc.

HTH,

Xiang-Jun

****************************************************************************
    DSSR: Software for Defining the (Secondary) Structures of RNA
      by Xiang-Jun Lu (xiangjun@x3dna.org), beta-r08-on-20130323

   The program is currently under active development. As always, we
   greatly appreciate your feedback! Please report all DSSR-related
   issues on the 3DNA Forum (http://forum.x3dna.org/), and I strive
   to promptly respond to any questions posted there.
****************************************************************************
Date and time: Mon Mar 25 16:15:39 2013
File name: 1xvk.pdb1
    no. of DNA/RNA chains: 1 [A=16]
    no. of nucleotides:    16
    no. of waters:         112
    no. of metals:         2 [Mg=2]
****************************************************************************
List of 8 base pair(s)
   1 1:A.DG1          2:A.DC8          [G+C]              00-n/a    cHW cM+W
       74.1(syn) C2'-endo lambda=47.5; -100.4(anti) C2'-endo lambda=61.6
       d(C1'-C1')=8.31 d(N1-N9)=6.64 d(C6-C8)=6.16 tor(N1-C1'-C1'-N9)=-0.1
       H-bonds[2]: "N7*N3[2.69]; O6(carbonyl)-N4(amino)[2.85]"
       bp_pars: [0.46    -3.41   -0.35   3.71    -5.92   67.37]
   2 1:A.DC2          2:A.DG7          [C-G] WC           19-XIX    cWW cW-W
       -103.9(anti) C4'-exo  lambda=55.5; -105.3(anti) C1'-exo  lambda=52.5
       d(C1'-C1')=10.53 d(N1-N9)=8.80 d(C6-C8)=9.65 tor(N1-C1'-C1'-N9)=0.0
       H-bonds[3]: "O2(carbonyl)-N2(amino)[2.76]; N3-N1(imino)[2.89]; N4(amino)-O6(carbonyl)[2.74]"
       bp_pars: [0.27    -0.18   0.35    -22.28  3.73    -2.75]
   3 1:A.DG3          2:A.DC6          [G-C] WC           19-XIX    cWW cW-W
       -108.4(anti) C1'-exo  lambda=54.1; -107.8(anti) C4'-exo  lambda=54.2
       d(C1'-C1')=10.47 d(N1-N9)=8.77 d(C6-C8)=9.65 tor(N1-C1'-C1'-N9)=-2.6
       H-bonds[3]: "O6(carbonyl)-N4(amino)[2.87]; N1(imino)-N3[2.88]; N2(amino)-O2(carbonyl)[2.81]"
       bp_pars: [-0.38   -0.18   0.41    22.84   2.54    -2.67]
   4 1:A.DT4          2:A.DA5          [T+A] Hoogsteen    23-XXIII  cWH cW+M
       -95.4(anti) C2'-endo lambda=59.7; 68.3(syn) C1'-exo  lambda=53.5
       d(C1'-C1')=8.37 d(N1-N9)=6.77 d(C6-C8)=6.17 tor(N1-C1'-C1'-N9)=0.5
       H-bonds[2]: "N3(imino)-N7[2.86]; O4(carbonyl)-N6(amino)[2.80]"
       bp_pars: [-0.69   3.57    0.31    -3.48   7.41    -70.55]
   5 1:A.DA5          2:A.DT4          [A+T] Hoogsteen    23-XXIII  cHW cM+W
       68.3(syn) C1'-exo  lambda=53.5; -95.4(anti) C2'-endo lambda=59.7
       d(C1'-C1')=8.37 d(N1-N9)=6.77 d(C6-C8)=6.17 tor(N1-C1'-C1'-N9)=0.5
       H-bonds[2]: "N7-N3(imino)[2.86]; N6(amino)-O4(carbonyl)[2.80]"
       bp_pars: [0.69    -3.57   -0.31   3.48    -7.41   70.54]
   6 1:A.DC6          2:A.DG3          [C-G] WC           19-XIX    cWW cW-W
       -107.8(anti) C4'-exo  lambda=54.2; -108.3(anti) C1'-exo  lambda=54.1
       d(C1'-C1')=10.47 d(N1-N9)=8.77 d(C6-C8)=9.65 tor(N1-C1'-C1'-N9)=-2.6
       H-bonds[3]: "O2(carbonyl)-N2(amino)[2.81]; N3-N1(imino)[2.88]; N4(amino)-O6(carbonyl)[2.87]"
       bp_pars: [0.38    -0.18   0.41    -22.84  2.54    -2.67]
   7 1:A.DG7          2:A.DC2          [G-C] WC           19-XIX    cWW cW-W
       -105.3(anti) C1'-exo  lambda=52.5; -103.9(anti) C4'-exo  lambda=55.5
       d(C1'-C1')=10.53 d(N1-N9)=8.80 d(C6-C8)=9.65 tor(N1-C1'-C1'-N9)=0.0
       H-bonds[3]: "O6(carbonyl)-N4(amino)[2.74]; N1(imino)-N3[2.89]; N2(amino)-O2(carbonyl)[2.76]"
       bp_pars: [-0.27   -0.18   0.35    22.28   3.73    -2.75]
   8 1:A.DC8          2:A.DG1          [C+G]              00-n/a    cWH cW+M
       -100.4(anti) C2'-endo lambda=61.5; 74.1(syn) C2'-endo lambda=47.6
       d(C1'-C1')=8.31 d(N1-N9)=6.64 d(C6-C8)=6.16 tor(N1-C1'-C1'-N9)=-0.1
       H-bonds[2]: "N3*N7[2.69]; N4(amino)-O6(carbonyl)[2.86]"
       bp_pars: [-0.46   3.41    0.35    -3.71   5.93    -67.37]
****************************************************************************
List of 1 helix
  helix=1[2] bps=8
   1 1:A.DG1          2:A.DC8          [G+C]              00-n/a    cHW cM+W
   2 1:A.DC2          2:A.DG7          [C-G] WC           19-XIX    cWW cW-W
   3 1:A.DG3          2:A.DC6          [G-C] WC           19-XIX    cWW cW-W
   4 1:A.DT4          2:A.DA5          [T+A] Hoogsteen    23-XXIII  cWH cW+M
   5 1:A.DA5          2:A.DT4          [A+T] Hoogsteen    23-XXIII  cHW cM+W
   6 1:A.DC6          2:A.DG3          [C-G] WC           19-XIX    cWW cW-W
   7 1:A.DG7          2:A.DC2          [G-C] WC           19-XIX    cWW cW-W
   8 1:A.DC8          2:A.DG1          [C+G]              00-n/a    cWH cW+M
****************************************************************************
List of 2 stems
  stem=1[#1] bps=2
   1 1:A.DC2          2:A.DG7          [C-G] WC           19-XIX    cWW cW-W
   2 1:A.DG3          2:A.DC6          [G-C] WC           19-XIX    cWW cW-W

  stem=2[#1] bps=2
   1 1:A.DC6          2:A.DG3          [C-G] WC           19-XIX    cWW cW-W
   2 1:A.DG7          2:A.DC2          [G-C] WC           19-XIX    cWW cW-W
****************************************************************************
List of 1 coaxial stack(s)
   1 Helix#1 contains 2 stems: [#1, #2]
****************************************************************************
List of 1 internal loop(s)
   1 symmetric internal loop: 8 nts; [2x2]; linked by [#1, #2]
       1:A.DG3+1:A.DT4+1:A.DA5+1:A.DC6+2:A.DG3+2:A.DT4+2:A.DA5+2:A.DC6 [GTACGTAC]
****************************************************************************
>chain-A #1 DNA* with 16 nts
GCGTACGCGCGTACGC
.((..((..))..)).

1138

General discussions (Q&As) / Re: Helical parametrs for 2′,5′-Linked RNA or DNA

« on: March 22, 2013, 11:25:30 pm »

Thanks for pointing out the PDB id (2kd4) of the structure you are interested in, and attaching the corresponding 3DNA output files. As always, such information is useful by making our discussions concrete.

Browsing through the output file (2kd4.out) and looking at the structure with Jmol, it appears to me 3DNA has no problem to analyze this structure.

First, the twist angles associate with intercalated steps are smaller, while those for the flanking steps are larger, than normal A-DNA twist values.

Local base-pair step parameters
    step       Shift     Slide      Rise      Tilt      Roll     Twist
   1 GC/GC      1.21      0.19      3.77      5.32    -16.06      9.14
   2 CC/GG      0.42     -1.82      4.09    -13.22      0.06     37.79
   3 CG/CG     -0.08     -1.17      6.68     -7.87    -11.95     20.28
   4 GC/GC      0.07     -0.19      3.40      0.22     -0.57     49.80
   5 CG/CG     -0.14     -1.18      6.62      3.50     -6.88     21.72
   6 GG/CC      0.16     -1.82      3.69     14.10     12.69     39.30
   7 GC/GC     -1.33      0.15      3.44     -6.43    -13.80      6.98

Second, the output for the sugar torsion angles are also as expected (see below for the output for strand I). Because of the 2'--5' backbone linkage, angles alpha [O3'(i-1)-P-O5'-C5'], epsilon [C4'-C3'-O3'-P(i+1)] and zeta [C3'-O3'-P(i+1)-O5'(i+1)] are not defined in the conventional sense.

Strand I
  base    alpha    beta   gamma   delta  epsilon   zeta    chi
   1 G     ---     ---   -135.8   105.5    ---     ---     18.3
   2 C     ---    145.3    62.2    86.8    ---     ---   -126.7
   3 C     ---    176.8   -35.2   134.4    ---     ---   -103.9
   4 G     ---    115.9    64.8    61.7    ---     ---   -137.3
   5 C     ---    163.8    48.4   113.5    ---     ---   -114.8
   6 G     ---    178.6    50.6    83.4    ---     ---   -120.9
   7 G     ---   -149.2    33.4   139.0    ---     ---    -83.2
   8 C     ---    154.9    20.2   120.1    ---     ---   -132.2

Note the following output from find_pair:

Code: [Select]

^^vv opposite bp direction: 1(8) 1(1)-2(2)
^^vv opposite bp direction: 1(8) 7(7)-8(8)
I have attached the stacking diagram of the first step, where one can see clearly the two base pairs have opposite orientation. The same is true for the last step. Such base flapping occurs around B-Z junctions, as in 2acj. I do not know how reliable this part of the structure is, or its relevance.

You may also want to analyze this structure (or related ones) using Curves+. I noticed that Horowitz et al. used Curves to analyze 2kd4 in their JACS publication.

HTH,

Xiang-Jun

1139

General discussions (Q&As) / Re: Helical parametrs for 2′,5′-Linked RNA or DNA

« on: March 20, 2013, 01:53:17 pm »

Could you try 3DNA on the structure you linked to ("Solution Structure and Thermodynamics of 2′,5′ RNA Intercalation"), and report back any problem you have?

Xiang-Jun

PS: In posting a question, it'd be very helpful to attach a structure file, or provide a PDB id.

1140

RNA structures (DSSR) / Re: Bug report of DSSR beta

« on: March 19, 2013, 12:36:23 am »

Hi Marc,

I've just released DSSR beta-r06-on-20130319 which should have fixed the "segmentation fault" bug for PDB entry 2a64. I've also taken this opportunity to update the -h (--help) message.

Have a try and report back how it goes!

Xiang-Jun

1141

RNA structures (DSSR) / Re: Bug report of DSSR beta

« on: March 18, 2013, 01:33:30 pm »

Hi Marc,

Thanks for trying out DSSR and reporting back a bug with the 64-bits unix version of beta-r05-on-20130316 on PDB entry 2a64. I've verified the bug and will get it fixed ASAP.

I've also tried the same 2a64 entry with the Mac OS X version without the reported problem. Since the "Segmentation fault (core dumped)" occurs after the output of "total number of junctions: 4", it is likely due to the ribose zipper detection function recently added into DSSR.

In any event, keep testing and report back any further issues you encounter. Also stay tuned for the next release!

Xiang-Jun

1142

RNA structures (DSSR) / Re: DSSR - List of bases involved in hairpins?

« on: March 14, 2013, 03:29:27 pm »

I've just updated DSSR to beta-r04-on-20130314. Among other things, now all nucleotides in hairpin loops are explicitly listed in addition to information reported before. So for 1msy, the new output looks like below:

Code: [Select]

List of 1 hairpin loop(s)
   1 nts=4 GUAA closed by pair {A.C2658+A.G2663 [CG], #-1}
       A.C2658+A.G2659+A.U2660+A.A2661+A.A2662+A.G2663 [CGUAAG]
I have some reasons to report hairpin loops differently from other types of loops (bulges, internal loops, junctions), one being to follow the convention. For example, at first glance, one would immediately see that 1msy contains a GUAA tetra-loop (of the most common GNRA type).

I am not aware of a consistent way to name other loops in the literature of RNA structures, so I've come up with my own convention. I'd like to hear what the community has to say when DSSR gains more popularity in the RNA structural world, and make adjustments accordingly.

HTH,

Xiang-Jun

1143

RNA structures (DSSR) / Re: DSSR - List of bases involved in hairpins?

« on: March 12, 2013, 06:40:11 pm »

Thanks for your feedback. I will get DSSR updated in a couple of days, where the requested info will be added to hairpin loops. In the meantime, please test DSSR more thoroughly and report back your thoughts that would make it better tool from a user's perpective.

Best regards,

Xiang-Jun

1144

RNA structures (DSSR) / Re: What are the definitions of helix, stem, step etc in DSSR?

« on: March 10, 2013, 10:44:05 am »

Thanks for downloading the beta testing version of DSSR. I am working on a manuscript, and more details will be made available in due course. In the meantime, here is a short answer to your questions:

• DSSR defines a helix from purely a base stacking perspective, regardless of base-pair type or backbone connection. So the tRNA example (1ehz) has two helices, corresponding to the two arms of the stereotypical L-shaped tertiary structure.
• The stem are further restricted by canonical base pairs (Watson-Crick and G-U wobble), and backbone connections (each strand treated separately). So the same tRNA (1ehz) has four stems, corresponding to the famous tRNA cloverleaf secondary structure.
• The shortest helix or stem is composed of two base pairs. DSSR does not define a separate class of 'step'.
• As you may have observed in the two example output files, DSSR does treat lone (isolated) Watsoc-Crick (or G-U wobble) base pair separately.

The default settings currently adopted in DSSR are based on my understanding of the conventions in the RNA structural world. I'd like to hear what users have to say and will make refinements accordingly. More functionality will be added the DSSR; with the current beta release, I've just kicked the ball rolling. Based on experience of supporting 3DNA and using other bioinformatics tools, I've created DSSR from the ground up to be trivial to set up and easy to use. Just play with it, and report back any issues you have.

HTH,

Xiang-Jun

1145

RNA structures (DSSR) / Re: Bug report of DSSR beta

« on: March 09, 2013, 04:06:55 pm »

Following my previous post, I've just released DSSR beta-r03-on-20130309. I tested this version against all nucleic-acid-containing structures in the PDB, so all the bugs you reported yesterday should have been fixed. I was actually approaching the end of such an extensive tests when I say your post, and your reported cases were already handled at that time.

It's worthy noting that some of the DSSR buggy cases are due to problematic data quality of the corresponding PDB entries. For example, 1bdn apparently have two structures overlapped (with 4 chains instead of 2) -- check the structure with PyMOL or Jmol, and you will seem my point.

While more bugs are bound to be uncovered, and further refinements are necessary, I am confident to say that DSSR (beta-r03-on-20130309) is robust enough for real-world applications. As always, do not be shy to report back any issues you encounter.

Xiang-Jun

1146

RNA structures (DSSR) / Re: Bug report of DSSR beta

« on: March 08, 2013, 04:42:50 pm »

I wholeheartedly welcome bug reports, and this latest is of no exception . I'll get them fixed shortly -- hopefully by tomorrow. With the next release, the program would be quite robust, I believe.

Have a great weekend.

Xiang-Jun

1147

General discussions (Q&As) / Re: Install 3DNA v2.1 on Windows

« on: March 07, 2013, 09:08:04 pm »

Did you read the post "How to install 3DNA on Linux and Windows?". The last part is about installing 3DNA "On MinGW/MSYS--based Windows".

HTH,

Xiang-Jun

1148

General discussions (Q&As) / Re: plot to show DNA conformation

« on: March 07, 2013, 09:21:25 am »

I'd suggest that you pay attention to the Zp parameter, along with some other more commonly used ones.

Good luck with you research!

Xiang-Jun

1149

General discussions (Q&As) / Re: plot to show DNA conformation

« on: March 06, 2013, 12:28:44 pm »

Hi Harry,

Thanks for posting on the Forum. 3DNA per se does not generate Ramachandran-like plot for DNA or RNA, or any plot of that kind. It calculates various structural parameters that can be plotted using third-party tools. Regarding Ramachandran-like plot for RNA (based on pseudo torsion angles), you may check papers from the Pyle lab.

3DNA v2.1 has a new option "analyze -torsion" to calculate commonly used torsion angles, including the following section:

Code: [Select]

Pseudo (virtual) eta/theta torsion angles:

Note: eta:    C4'(i-1)-P(i)-C4'(i)-P(i+1)
      theta:  P(i)-C4'(i)-P(i+1)-C4'(i+1)

      eta':   C1'(i-1)-P(i)-C1'(i)-P(i+1)
      theta': P(i)-C1'(i)-P(i+1)-C1'(i+1)

      eta":   Borg(i-1)-P(i)-Borg(i)-P(i+1)
      theta": P(i)-Borg(i)-P(i+1)-Borg(i+1)

To follow DNA conformational changes, you may find our 2000 JMB paper titled "A-form conformational motifs in ligand-bound DNA structures" relevant. The work examined a variety of parameters, and found Zp to be the most discriminative between A- and B-form DNAs.

Furthermore, 3DNA v2.1 has a Ruby script x3dna_ensemble that makes analysis of MD simulation trajectories straightforward. We have a new 3DNA JoVE paper in press that provides protocols for such MD analysis.

HTH,

Xiang-Jun

 

1150

RNA structures (DSSR) / Re: Bug report of DSSR beta

« on: March 06, 2013, 12:05:45 pm »

Hi,

All the bugs you reported should have been fixed, and DSSR beta-r02-on-20130306 released. Please have a try and report back how it goes.

As you may notice, I have split the thread from the original DSSR announcement post. I also put the long list of buggy cases reported in your second post into the code tag for easy browsing.

Do not be shy in bringing up more bugs --- the more, the merrier! 

Xiang-Jun