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

Author Topic: How is the helix axis calculated in 3DNA  (Read 18329 times)

Offline wwang02

  • with-posts
  • *
  • Posts: 6
    • View Profile
How is the helix axis calculated in 3DNA
« on: September 26, 2014, 10:12:09 am »
How does 3DNA define a helix? If I want to calculate a helix axis of a DNA, what position should I start to extend the local vector? and how should I calculate the base coordinates, by averaging all the atom coordinates of base pairs?

I need to know this, because I want to implement the analysis method of 3DNA to my simulations which are not done by the all-atom fashion. Does anyone know whether there exists any publications where I can refer to?

Thank you

Offline xiangjun

  • Administrator
  • with-posts
  • *****
  • Posts: 1646
    • View Profile
    • 3DNA homepage
Re: How is the helix axis calculated in 3DNA
« Reply #1 on: September 26, 2014, 10:29:58 am »
In 3DNA, the global linear helical axis is defined as in SCHNAaP/SCHNArP, which follows NewHelix. Basically, equivalent C1' and RN9/YN1 atom pairs along each strand are employed in a least-squares fitting procedure to derive the helix. See "Structure and conformation of helical nucleic acids: analysis program (SCHNAaP)" (http://www.ncbi.nlm.nih.gov/pubmed/9356255) for details. The source code of SCHNAaP/SCHNArP is available --- see the Section "SCHNAaP/SCHNArP" (http://forum.x3dna.org/schnaap-schnarp/).

HTH,

Xiang-Jun
 

Offline wwang02

  • with-posts
  • *
  • Posts: 6
    • View Profile
Re: How is the helix axis calculated in 3DNA
« Reply #2 on: September 26, 2014, 03:38:53 pm »
If the DNA simulation is done in a coarse grain fashion, can I calculate the center of mass of triads and build up local unit vectors, and average the all the unit vectors to define the helix axis? Do you think that would be close enough to the method you use in the paper?

Offline xiangjun

  • Administrator
  • with-posts
  • *****
  • Posts: 1646
    • View Profile
    • 3DNA homepage
Re: How is the helix axis calculated in 3DNA
« Reply #3 on: September 26, 2014, 10:48:28 pm »
It is hard to say anything specific. Why not you have a try and see what happens?

Note that in addition to the 'global linear helical axis', 3DNA also defines a local helical axis for each dinucleotide step. See the 2003 3DNA NAR paper and the doc "tech_details.pdf" distributed with 3DNA v2.1 for details. You may find them useful.

Xiang-Jun

Offline wwang02

  • with-posts
  • *
  • Posts: 6
    • View Profile
Re: How is the helix axis calculated in 3DNA
« Reply #4 on: September 29, 2014, 11:33:18 pm »
Thank you for replies!

 

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