[phenixbb] SA omit map
PAfonine at lbl.gov
Thu Sep 10 21:00:41 PDT 2009
> I have a theoretical/practical question about omit maps and refinement.
> I am completing the refinement (in Phenix) of a protein-ligand complex
> at 1.3A resolution. I solved it by MR and automatic rebuilding of the
> protein alone first then built in the ligand. Rfree and Rfac are
> 19.4%/18.2% after TLS and water-picking in Phenix. The model includes
> everything protein, water, ligand and some ions.
Just a few remarks:
- At this resolution you should refine individual anisotropic ADPs,
rather than constrained ADPs (using TLS model), for all macromolecular
atoms. You may try to do so for water as well, but I would try both and
see which gives more reasonable Rwork, Rfree and Rfree-Rwork;
- Riding H atoms should be included as well;
- As a final step, you may want to run the refinement with
"optimize_wxc=true optimize_wxu=true" which may take a while to run, but
will optimize relative X-ray/Restraints weights.
> However I have some slight doubts one region in my ligand.
> What would be the best map, less biased, to look at this "very late"
> stage of the refinement. Are composite or systematic SA omit map
> useful options at this stage ?
As a possible and easy to do option you can use Average Kick Maps, which
we recently implemented in PHENIX. For details, here is the reference:
Acta Cryst. (2009). D65, 921-931. "Averaged kick maps: less noise, more
signal...and probably less bias".
I just noticed Tom's reply, but since I went this far, I will continue
and send my version of the reply as well...
In summary, the idea is:
An average kick map (AK map) is computed as following (Gunc(ar et al.,
2000; Turk, 2007; Pranikar et al., 2009): a large ensemble of
structures (several hundreds) is created where the coordinates of each
structure from the ensemble are all randomly shaken. The shake amount
(rmsd distortion introduced to coordinates) varies from 0 to 1.0 Å. Then
for each structure a map is computed ((mFobs-DFmodel)exp(i?model) or
(2mFobs-DFmodel)exp(i?model) or any other map, for example a ligand-omit
map). Finally, all maps are averaged out to produce one averaged kick
map. An AK map is expected to have less or no bias, less noise, enhance
existing signal and potentially can clear up some initially bad densities.
So, to compute it, omit the ligand from your structure and follow the
example in Tom's email.
More information about the phenixbb