Biology Stack Exchange is a question and answer site for biology researchers, academics, and students. It's 100% free, no registration required.

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

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

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()

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

#then to complete the alignment

You can look here for further details. SuperPoser

share|improve this answer

Your Answer


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.