[phenixbb] measuring the angle between two DNA duplexes
tg at shelx.uni-ac.gwdg.de
Tue Jan 21 01:25:15 PST 2014
that's the method described in
http://journals.iucr.org/a/issues/2011/01/00/sc5036/index.html ;-) based
on the moments of inertia (a computer scientist might name it
differently). I am not sure, though, you would get the desired result
for short helices. E.g. a helix defined by three atoms the eigenvalue
would point roughly in the direction of the external phosphates, which
is far from parallel with the helix axis.
On 01/21/2014 04:20 AM, Pavel Afonine wrote:
> Hi Ed,
> interesting idea! Although I was thinking to have a tool that is a
> little more general and a little less context dependent. Say you have
> two clouds of points that are (thinking in terms of macromolecules) two
> alpha helices (for instance), and you want to know the angle between the
> axes of the two helices. How would I approach this?..
> First, for each helix I would compute a symmetric 3x3 matrix like this:
> sum(xn-xc)**2 sum(xn-xc)*(yn-xc) sum(xn-xc)*(zn-zc)
> sum(xn-xc)*(yn-xc) sum(yn-yc)**2 sum(yn-yc)*(yz-zc)
> sum(xn-xc)*(zn-zc) sum(yn-yc)*(yz-zc) sum(zn-zc)**2
> where (xn,yn,zn) is the coordinate of nth atom, the sum is taken over
> all atoms, and (xc,yc,zc) is the coordinate of the center of mass.
> Second, for each of the two matrices I would find its eigen-values and
> eigen-vectors, and select eigen-vectors corresponding to largest
> Finally, the desired angle is the angle between the two eigen-vectors
> found above, which is computed trivially.
> I think this a little simpler than finding the best fit for a 3D line.
> What you think?
> On 1/20/14, 2:14 PM, Edward A. Berry wrote:
>> 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.
>> phenixbb mailing list
>> phenixbb at phenix-online.org
> phenixbb mailing list
> phenixbb at phenix-online.org
Dr Tim Gruene
Institut fuer anorganische Chemie
GPG Key ID = A46BEE1A
-------------- next part --------------
A non-text attachment was scrubbed...
Size: 234 bytes
Desc: OpenPGP digital signature
More information about the phenixbb