[phenixbb] Fo-Fc omit map for ligand in the protein-ligand structure
nechols at lbl.gov
Sat Sep 13 17:57:40 PDT 2014
On Sat, Sep 13, 2014 at 3:03 PM, Pavel Afonine <pafonine at lbl.gov> wrote:
> In the first case the bulk-solvent mask will be set in the ligand region
> and therefore it will mask ligand density (bulk-solvent will be filled into
> the ligand region). Depending on the strength of ligand density it may be
> masked completely or deteriorated.
If you follow the second option you will always get positive density in
> ligand area. This density may correspond to bulk-solvent, ligand or mixture
> of both. That is there will be no simple way to differentiate whether this
> density arises from the ligand or bulk-solvent.
I thought phenix.refine ignores zero-occupancy atoms when calculating the
bulk solvent mask? In fact there is an option specifically to toggle this
option, and I've noticed very different results with and without it.
Wei: It's fine to use the options you mentioned for this - however, weight
optimization is very slow and not really necessary since you are only
running SA to remove phase bias, not to get the best possible model. I
usually just use individual_sites+individual_adp plus whatever extra
restraints are appropriate (NCS, secondary structure, etc.). I wouldn't
worry if you end up with an overfit model with worse geometry - as long as
it has been sufficiently shaken (but not too much!) and the ligand
difference density is clear, you are fine. Just be careful that in the
process of running SA, you haven't actually introduced new bias into the
maps - but the mFo-DFc map is relatively safe (the 2mFo-DFc map much less
Restraining the zero-occupancy ligand (and perhaps adjacent residues) is
strongly recommended, especially if it contains heavier elements. (The
newer versions of the phenix.refine GUI have a button in the Output tab to
set this up.) Without restraints, you run the risk of refining the
surrounding model into the ligand density. As always this is worse at
lower resolution - at 2.8Å I'm not sure what to expect.
Perhaps one day we will write up something more formal about map
"validation", since there seem to be a lot of misconceptions about these
methods and they are not always used appropriately.
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the phenixbb