[phenixbb] anisotropic ADP factor refinement and optimize_wxc

Pavel Afonine PAfonine at lbl.gov
Tue Jul 21 21:34:31 PDT 2009


Hi Gordon,

the target for ADP (B-factor) refinement is :

T = wxu * wxu_scale * Txray + wu * Tadp

and the target for coordinate refinement is:

T = wxc * wxc_scale * Txray + wc * Tstereochemistry .

The weights wxc (wxu) are determined automatically (similarly to CNS) as 
described in:

Adams, P.D., Pannu, N.S., Read, R.J. & Brünger, A.T. (1997). Proc. Natl. 
Acad. Sci. 94, 5018-5023;

the weight scale wxc_scale (wxu_scale) is 0.5 or 1.0 and the user can 
adjust it, wc (wu) = 1 (and the only purpose of it is to be able to turn 
off the restraints), Txray is x-ray target function (maximum-likelihood, 
least-squares, mlhl, ...), and Tadp (Tstereochemistry) is restraints term.

When you use "optimize_wxc=true" you optimize weights for coordinate 
refinement, and to do similar thing for ADP refinement you need to use 
"optimize_wxu=true".

> phenix.refine protein1.pdb 072109checks.mtz ordered_solvent=true strategy=individual_adp adp.individual.anisotropic="not element H" optimize_wxc=true --overwrite
>   

This command above refines only anisotropic B-factors of non-H atoms, 
while the coordinates are not refined at all, and therefore the 
statement "optimize_wxc=true" is not used.
I guess what you really meant is:

phenix.refine protein1.pdb 072109checks.mtz ordered_solvent=true 
adp.individual.anisotropic="not element H" 
adp.individual.isotropic="element H" optimize_wxc=true optimize_wxu=true 
--overwrite

By the way, what is the resolution of the data?

If the geometry is still not good (it's another page to write to 
describe why it is potentially possible), then instead of using 
"optimize_wxc=true" try manually set up the weight putting less weight 
on X-ray term: use wxc_scale=XXX for this, where XXX is a value less 
than the default one (0.5). Say, try 0.2, 0.1, etc.

> Is it possible to refine both anisotropic B factors and optimize ADP factors in the same run.

Yes. You can refine any set of parameters for any selected part of your 
structure in any combination.

If it is still not clear, please fell free to send me the model and data 
and I will get back to you with the optimal refinement strategy that you 
can use in future as the example.

Cheers,
Pavel.






More information about the phenixbb mailing list