Take the 2-minute tour ×
Biology Stack Exchange is a question and answer site for biology researchers, academics, and students. It's 100% free, no registration required.

I have a series of protein models with DNA docked. I now want to superpose the DNA on a reference DNA molecule and extract the translational distance applied and the rotation angle used.

I can superpose DNA molecules with Chimera, but I can't extract the translational distance and the rotational angle that was applied to superpose.

I have also had a look at MODELLER but it seems MODELLER cannot not superpose DNA.

Does anyone know of software/code by which I can extract the information I want?

share|improve this question
add comment

1 Answer 1

Try this out with BioPython:

from Bio.PDB.PDBSuperimposer import PDBSuperimposer as superimposer
from Bio.PDB.PDBParser import PDBParser

parser = PDBParser()
sup = superimposer()
struct_1 = parser.get_structure("XXXX","first_pdb")
struct_2 = parser.get_structure("XXXX","second_pdb")

atoms1 = struct_1.get_atoms()
atoms2 = struct_2.get_atoms()
sup.set_atoms(atoms1,atoms2)

print sup.rotrans #print the rotation translation matrix from the SVD

#then to complete the alignment
sup.apply()  

You can look here for further details. SuperPoser

share|improve this answer
add comment

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.