<div dir="ltr">I just tried phenix.maps, and it yields substantially different difference maps than those obtained through my refinement process. Are there any other techniques or parameters I can modify to make residues with multiple conformations 'invisible' to phenix.refine, in addition/besides set occupancies to 0 or deleting the residue form the pdb file?<br><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Feb 17, 2016 at 3:48 PM, Pavel Afonine <span dir="ltr"><<a href="mailto:pafonine@lbl.gov" target="_blank">pafonine@lbl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">
Hi Alex,<br>
<br>
why use a cannon to kill a fly? phenix.refine is for refinement not
map calculation. I suggest to use more specific tools. phenix.maps
should do it, for example.<br>
<br>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__phenix-2Donline.org_documentation_reference_phenix-5Fmaps.html&d=AwMDaQ&c=-dg2m7zWuuDZ0MUcV7Sdqw&r=erGJlXdJeX4XOslWBGLhRT0lLKKneOe9LsS9KwzY1H0&m=yk8Yay71IL1yLpt6JP2X6qX4HpzgJjlL51q_evgsG00&s=ZgesLYvmd_lhLM6nE_t0g1iPI2GN-1LFE-2NgeXd4oQ&e=" target="_blank">http://phenix-online.org/documentation/reference/phenix_maps.html</a><br>
<br>
Pavel<br>
<br>
<blockquote type="cite">
<div dir="ltr">
<div class="gmail_quote">
<div dir="ltr">I am trying to get the difference electron
density maps for residues with multiple conformations after
I set the occupancy of their atoms to 0.
<div>My main goal is to visualize the green positive density
peaks that should appear for the atoms whose occupancy I
set to 0.</div>
<div>However, when I run the phenix.refine command (doing
all of this via cli, not gui), the output pdb file </div><span class="">
<div><br>
</div>
<div><br>
</div>
<div>I will try to illustrate my process and problem with an
example that works fine and one that doesn't, side by
side, step by step.</div>
</span><div>
<ol>
<li>I generate a pdb file that has a chosen residue set
to occupancy 0. See
'pdb_1agy_res_SER_resID_120_del.pdb' and
'pdb_1agy_res_SER_resID_129_del.pdb'. As you can see,
the former has SER 120 set to 0 occupancy and has one
conformation, and the latter has SER 129Â set to 0
occupancy for BOTH of its conformations. Both files
are correct up to this point.</li><span class="">
<li>I then execute the following command:</li>
</span></ol>
phenix.refine <b>pdb_del_file</b> <b>mtz_file</b> main.number_of_macro_cycles=5
output.write_eff_file=False output.write_geo_file=False
output.write_def_file=False
refinement.input.xray_data.r_free_flags.generate=True
--overwrite
<ul>
<li>with <b style="font-weight:bold">mtz_file </b>set
to '1agy.mtz' (from the PDB) and <b>pdb_del_file</b>Â set
to either 'pdb_1agy_res_SER_resID_120_del.pdb' and
'pdb_1agy_res_SER_resID_129_del.pdb'.</li>
</ul>
</div>
<div>
<ul>
<li>This results in 2
files pdb_1agy_res_SER_resID_XXX_del_refine_001.mtz
and pdb_1agy_res_SER_resID_XXX_del_refine_001.pdb.
Comparing these pdb files, you can see that the
occupancy for SER 120 remains set to 0 (which is what
I want), but the occupancies for the 2 conformers for
SER 129 have been filled in by phenix (not what I
want) somewhere during the refinement process. As a
result, I get the desired difference peaks for SER
120, but not for SER 129 (which seems to be averaged
out).</li>
</ul><span class="">
Other than deleting the 2nd conformation for SER 129, is
there a way to solve this problem, preferably via command
line arguments/flags?</span></div>
<div><br>
</div>
<div>Thanks,<br>
Alex</div>
</div>
</div>
</div>
<br>
<fieldset></fieldset>
<br>
<pre>_______________________________________________
phenixbb mailing list
<a href="mailto:phenixbb@phenix-online.org" target="_blank">phenixbb@phenix-online.org</a>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__phenix-2Donline.org_mailman_listinfo_phenixbb&d=AwMDaQ&c=-dg2m7zWuuDZ0MUcV7Sdqw&r=erGJlXdJeX4XOslWBGLhRT0lLKKneOe9LsS9KwzY1H0&m=yk8Yay71IL1yLpt6JP2X6qX4HpzgJjlL51q_evgsG00&s=tWFeSN7Z4_pMDE3OgVcJOC-Tl4WNFA1GwHcDRgVchmg&e=" target="_blank">http://phenix-online.org/mailman/listinfo/phenixbb</a>
Unsubscribe: <a href="mailto:phenixbb-leave@phenix-online.org" target="_blank">phenixbb-leave@phenix-online.org</a></pre>
</blockquote>
<br>
</div>
</blockquote></div><br></div></div>