Perhaps this should have gone a to a new thread...
I have a couple of questions about RMSD. Maybe they could lead to "requested
features":)
1) from the documentation, it seems that rmsd calc _always_ does superposition
of the molecules before calculating RMSD. From some older (2007) openbabel
thread it seams there is no way to compute RMSD without superposition. This
would be very useful, say, in ligand-protein docking calculations.
2) The computed RMSDs are "symmetry corrected". There is a global symmetry and
a local symmetry. Is the local symmetry also calculated for the purposes of
symmetry correction when doing RMSD? Actually, it's not a trivial problem; one
way to tackle it is to permute rotations around all rotatable dihedrals in a
reference molecule, and to try to obtain a nearly identical molecule ("nearly"
because of imperfections). Perhaps other algorithms exist. An example would be
a local symmetry of the phenyl ring in phenylalanine: phenylalanine is
non-symmetric, but phenyl inside it has two pairs of equivalent carbons.
Best regards,
Vis
------------------------------------------------------------------------------
RSA(R) Conference 2012
Mar 27 - Feb 2
Save $400 by Jan. 27
Register now!
http://p.sf.net/sfu/rsa-sfdev2dev2
_______________________________________________
OpenBabel-scripting mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/openbabel-scripting