A horrible bug -
this line was somehow lost from lsqkab in the mists of time:
Lines 1029-132 - ish
should be
IF (ATNMW(I)(1:3).EQ.'CA ' .OR. ATNMW(I)(1:3).EQ.'N ' .OR.
+ ATNMW(I)(1:3).EQ.'C ' .OR. ATNMW(I)(1:3).EQ.'O ' .OR.
+ ATNMW(I)(1:3).EQ.'P ' .OR. ATNMW(I)(1:3).EQ.'O1P' .OR.
etc
NOT
IF (ATNMW(I)(1:3).EQ.'CA ' .OR. ATNMW(I)(1:3).EQ.'N ' .OR.
+ ATNMW(I)(1:3).EQ.'P ' .OR. ATNMW(I)(1:3).EQ.'O1P' .OR.
Thank you for noticing this!
Eleanor
Nathalie Colloc'h wrote:
Hello all,
I have still a question about LSQKAB
Why LSQKAB gives a rmsd (or a rms xyz I don't mind) non null for
glycine side chains ?
thanks a lot
nathalie
Eleanor Dodson a écrit :
Q1) Rms xyz and rmsd mean exactly the same . And as for what you
should report - that depends on the problem. As long as you state
clearly what has been fitted ( all atoms, CAs, loops excluded etc etc )
Q2) The algorithms for SSM and LSQKAB are different - SSM fits
secondary structure elements only, then checks the rmsd of all CAs.
This is particularly useful for fitting homolous proteins.
LSQKAB requires a list of all the atoms to be matched and that "rms
xyz" considers all those atoms and only those.
If you want a report of all differences you must click for that on
the GUI.
And as you say this is helpful for finding differences in different
copies of the same molecule.
Eleanor