![](https://secure.gravatar.com/avatar/7984547b53fca6ae837ba286f9c92e29.jpg?s=120&d=mm&r=g)
Hello Ivan,
I want to see a ligand in the structure I have. I am using Phenix for refinement. When I do a molecular replacement with the model or part of the model, it looks like there is continues density of the ligand on 2Fo-Fc and Fo-Fc map at lower counter level.
I presume at this point you have not placed the ligand yet, and therefore the ligand map you look at is a "discovery map" in Dale Tronrud's terms (see my previous email where I mention this).
1st QUESTION: Once we got the phases from the model which serves for MR, what is the counter level limit for a 2Fo-Fc and Fo-Fc maps we can look at ??
It is a good question with no universal answer. First, it depends how you scale the map; for recent sort of review see: Acta Cryst. (2014). D70, 2593-2606 Metrics for comparison of crystallographic maps A. Urzhumtsev, P. V. Afonine, V. Y. Lunin, T. C. Terwilliger and P. D. Adams If map is "sigma scaled" then the rule of thumb is 1sigma for 2mFo-DFc and 3sigma for mFo-DFc. Now, this applies for fully occupied typical atoms. For mobile partially partially occupied atoms (like those of ligands) you can go well below these 1 or 3 "sigma". "1sigma" rule applies to composite residual OMIT map computed with phenix.composite_omit_map (see previous email). If you use FEM map: http://phenix-online.org/presentations/fem_06MAY2014.pdf then contouring doesn't really matter, and anything below or equal to "1sigma" is suitable.
Further going on the structure refinement, at some point I add the compound I poorly see and I do 3 different refinements !! One with compound occupancy 0.8, another one with 0.7 and another one with 0.6.
Well, I presume you do not just set occupancy to these numbers, but you actually refine group occupancy starting from these values. If this is correct, then I hope you are getting the same result. If you just set these occupancies but you do not refine them, then I don't think it is a good idea.
Then the density for the compound after refinement appears. At that point I am on the second option in Pavel's suggestions below - I get probably a biased positive density in the map. BUT, I would like to notice that my Fwork and Ffree go significantly down. 2nd QUESTION: Does it mean the compound is there if statistics do down??
Usually, a ligand generates a small contribution to the total scattering so it is next to "invisible" by global metrics such as Rwork and Rfree. Therefore using local map and model quality metrics is a better way to do assessments. However, if adding it improves both Rs then I would say it is a good sign. All the best, Pavel