Question: Recalculate Atomic Coordinates After Changing Dihedral Angles
3
Daniel Standage3.9k wrote:

Disclaimer: I'm taking a class to get some experience with structural informatics and this question is related to an assignment for that class.

Given a PDB file, I have written a program that will calculate various lengths and angles from the atom coordinates provided in the file. I now need to select a certain amino acid, set its phi and psi angles to 0, recalculate the atomic coordinates for all affected atoms, and check for steric clashing.

I'm struggling conceptually with how to recalculate atomic coordinates. What types of concepts do I need to understand to tackle this type of question?

coordinates protein • 4.3k views
written 7.8 years ago by Daniel Standage3.9k
1
Daniel Standage3.9k wrote:

After some research, I ended up using Rodrigues' rotation formula. Given a vector v, an axis of rotation k, and a (right-handed) angle θ, the new rotated vector is calculated using the following formula.

v_rot = v cos(θ) + (k x v)sin(θ) + k (k · v)(1 - cos(θ))

In this case, the axis of symmetry k is the bond around which I'm rotating, the angle is θ is -phi or -psi (I'm setting the angles to 0), and the rotation is applied to every vector v subsequent to the axis of rotation.

0
dimkal730 wrote:

Daniel, i'd start looking at rotation matrices, which have to be applied to the part of the structure affected by the dihedral angle change.

Dima

0
Flow1.5k wrote:

It is much much easier than all that. You just learn some basic pymol/python scripting and just using "set_dihedral", you do not have to worry about the rest, pymol does it.