Bernhard, I think for the loop that is about to close taking the average from flanking residues is good enough. Not sure what you mean by the second example ("last residue has a low B").
I think what causes problems in refmac is a sharp discontinuity, as it contributes greatly to the B-factor restraints. Refinement will adjust B-factors properly if they are higher/lower, but not if they jump from 90 to 20 from residue to residue. The way I adjust B manually is by looking at the B-factor of the next residue, which is not too sophisticated approach and should be easy to implement. Well, in the end it's a minor problem, and easy to fix manually. Cheers, Ed. On Thu, 2010-11-04 at 16:11 +0100, Bernhard Lohkamp wrote: > Dear Ed, > > I think I get your idea, however this is not a straight forward issue I > believe. What do you imagine happens for a loop then? You may end up > connecting a high B residue to a low one when closing it? Does that make > refmac happy? Or your last residue (since it is still well ordered) has > a low B, then you dont gain anything? More appropriate may be a > weighting and adjustment of B-factors based on e.g. real-space CC. But I > guess until we have tools in Coot to refine (or in the first instance > guess) B-factors you better adjust the Bs yourself. Your guess (knowing > the structure and seeing the map) is better than any approximation we > can easily come up with (at least for now). > > My 50 Oere (just fazed out BTW), > > Bernhard > > > When adding residues to the N- or C-terminus, it may be desirable that > > the B-factors of the new residue match that of the previous terminal > > residue. Right now it uses defaults (20 for main chain, 30 for side > > chain) which are sensible fir a well-ordered model, but this causes > > problems with the subsequent refmac run for somewhat disordered regions > > where<B> may be significantly higher. > > > > The current workaround is to modify residue properties manually (or keep > > resetting the defaults). Naturally, the wisdom of adding extra residues > > with<B>~100 is questionable :) > > > > Cheers, > > > > Ed. > > > > -- Coot verendus est