Thanks, Ian, this is excellent. It appears that depending on the sequence the "ideal target rmsd" may vary from 0.018 for a poly-H to 0.024 for a poly-P. Except for some really short sequences, in PDB the variation is generally between 0.021-0.022, indeed undetectable.
On Fri, 2010-10-15 at 10:26 +0100, Ian Tickle wrote: > Ed, > > > On a practical side, I now have a misleading bond_rmsd value. The > > correct one can be calculated as described, but this may make geometry > > weight optimization cumbersome. Do I understand correctly that an > > alternative is to monitor Zbonds, with a rule that it should be around > > 1.0? > > No, the average value of Zbonds over the whole PDB is almost exactly > 0.5 (there's no clear resolution-dependence), so that would be a good > starting point. Re-refinement of ~98% of the PDB (PDB-REDO project) > using weight optimisation showed a linear dependence of Zbonds on > resolution in the resolution range 1 to 3 Ang.: > > Zbonds = 0.85 - 0.146*dmin (SD = 0.2) > > so e.g. <Zbonds> = 0.56 @ 2 Ang. > > Cheers > > -- Ian > > > And more generally, shouldn't we not look at rmsd_bonds at all and > > only use Zbonds instead (which is, I assume, an average bond length > > deviation to the target value ratio?) > > Correct. > > > I suspect that "acceptable" > > bond_rmsd value is slightly affected by sequence. > > Not really, the effect of sequence on the bond lengths is > insignificant in comparison with the uncertainty. > > Cheers > > -- Ian -- "I'd jump in myself, if I weren't so good at whistling." Julian, King of Lemurs