Hi all--
I'm trying to make a 2Fo-Fc map with density wrapped around my asymmetric unit to make a figure showing density for my sub-8Å structure. I am not building into this map or doing MR. My model was refined in Phenix.refine (R's in the high 30's), and so I have the output model .pdb and the output map coefficients .mtz. When I input these into cut_out_density, the only "map labels" input options it gives me are 2FOFCWT, PH2FOFCWT, which seem fine. Using either box, sphere, or model for the cutout type, I get a VERY nice looking map that looks extremely solvent flattened and tight around my model, albeit still to 8Å.
Suspicious of anything that looks good, I went to try the same thing in FFT (CCP4). There I generated a simple map (using the same .mtz) to cover all atoms in a pdb file- the same model I put into cut_out_density. For the input I selected 2FOFCWT and PH2FOFCWT, leaving Sigma and Weight unassigned. The map that comes out of this looks terrible, with a lot of noise and little density corresponding to the model aside from what one might expect from a very poor resolution experimental map. I've also tried using FFT in Phenix to do the same thing, and it gives like better maps than CCP4 FFT but worse than cut_out_density. It almost appears as if there is some additional solvent flattening going on, which I know happens at the beginning of refinement, but that should be intrinsic to the map coefficients file which is common to everything I've tried.
My questions are (1) what should I use for what I am trying to do (just display 2Fo-Fc density around asymmetric unit for a figure) and (2) what is the difference between what these three programs (cut_out, FFT-CCP4, FFT-Phenix) are calculating? The inputs are exactly the same -- mtz (2FOFCWT and PH2FOFCWT) and pdb output from refinement.
Thanks for the help,
Yunji
