I am using the internal_coords
module from BioPython in order to compare dihedral angles of two different conformations of the same protein. The conformations are processed by me in order to make sure that their residues are aligned and they always contain exactly the same atoms, so that in the end I expect to have the same number of dihedral angles in the two structures.
However, despite the match in the residues and atoms, I obtain a different number of dihedral angles for the two conformations.
You can find a reproducible code sample here. These are the two PDB files protein1.pdb and protein2.pdb used.
As you can see, the number of chains, the number of residues, the residue names and the atom names in each residue are the same, yet there are 4 extra dihedral angles in the second protein.
My question is, am I missing something from my comparison? Is it unreasonable to expect the same dihedral angles for these two structures? (And if yes, why?). Is this the expected behaviour of the atom_to_internal_coordinates()
function?
And extra: do you have any recommendations for other tools that can easily compute all dihedral angles (including the ones in side chains) for a protein?
align
) and look at the residues by eye —PyMOL does have the commandget_dihedral
, but the laziest way would be to click the four atoms (in edit mode) and it will show the dihedral. Do the values differ? $\endgroup$