[phenixbb] some general refine questions

Pavel Afonine pafonine at lbl.gov
Wed Jul 4 22:26:25 PDT 2007


Hi Bill and all interested in this,

> Why is the weight multiplied by 0.5 by default?  Since changing that  
> to 1.0 slightly increases freeR, I assume it is the default for some  
> reason.
>   

the total target for refinement of coordinates in phenix.refine is 
defined as:
Etotal = wc_scale * wxc * Exray + wc * Egeom

Similar for refinement of isotropic or anisotropic ADP.

Here:
- Exray is least-squares (LS), maximum-likelihood (ML) or any other kind 
of x-ray target;
- Egeom is geometry restrains target (Ebonds + Eangles + Eplanarities + 
Enonboned +...);
- wc = 1.0 and used to exclude geometry restraints by setting wc = 0.0;
- wxc ~ ratio of gradient's norms for geometry and x-ray targets as 
defined in (Adams et al, 1997, PNAS, Vol. 94, p. 5018);
- wc_scale is 'ad hoc' scale found empirically to be ok for most of the 
cases at medium resolution. It should be optimized at higher and lower 
resolutions. The default values are most likely NOT optimal at 
resolutions higher then 1.5A. Currently we are working on making it 
resolution dependent rather than one single value for all resolutions.

It is advised to play with wc_scale manually if:
- Rfree or Rfree - Rwork look not realistic;
- geometry deviations look not good (too tight or too loose at given 
resolution).
However, once again, the default values for most of cases at "normal" 
resolution are OK.

In future we are planning to work out more robust and sophisticated 
methodology for automatic weights determination.

Pavel.



More information about the phenixbb mailing list