Pavel Afonine wrote: . .
The underlying procedure would do the following: - extract two sets of coordinates of atoms corresponding to two provided atom selections; - draw two optimal lines (LS fit) passing through the above sets of coordinates; - compute and report angle between those two lines?
This could be innacurate for very short helices (admittedly not the case one usually would be looking for angles), or determining the axis of a short portion of a curved helix. A more accurate way to determine the axis- have a long canonical duplex constructed with its axis along Z (0,0,1). Superimpose as many residues of that as required on the duplex being tested, using only backbone atoms or even only phosphates. Operate on (0,0,1) with the resulting operator (i.e. take the third column of the rotation matrix) and use that as a vector parallel to the axis of the duplex being tested.