<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Thank you, Edward and Oleg for the suggestions, and Pavel for the clarification. You need two selections in the constrained group and an altloc for the ligand that’s different from the protein side chain. I did not force the occupancies to add up to 1 but the
program did it anyway. Phenix refined the ligand’s occupancy to 0.68 and the PHE side chain to 0.32.
<div class=""><br class="">
</div>
<div class="">This is what I did:</div>
<div class=""><br class="">
</div>
<div class="">I did not make an alternate conformation for the PHE, so there’s no BPHE. I just added the lock to the side chain starting at carbon beta. So, I have APHE for the side chain and BLIG for all the ligand atoms and reduced their occupancies to <
1.
<div class=""><br class="">
<div class="">This is the parameter file:</div>
<div class=""><br class="">
</div>
<div class="">
<div class="">
<div class="" style="margin: 0px; line-height: normal;"><span class="" style="font-size: 12px;">refinement {</span></div>
<div class="" style="margin: 0px; line-height: normal;"><span class="" style="font-size: 12px;"> refine {</span></div>
<div class="" style="margin: 0px; line-height: normal;"><span class="" style="font-size: 12px;"> occupancies {</span></div>
<div class="" style="margin: 0px; line-height: normal;"><span class="" style="font-size: 12px;"> constrained_group {</span></div>
<div class="" style="margin: 0px; line-height: normal;"><span class="" style="font-size: 12px;"> selection = (chain B and resseq 890) and altloc A</span></div>
<div class="" style="margin: 0px; line-height: normal;"><span class="" style="font-size: 12px;"> selection = resname LIG and altloc B</span></div>
</div>
<div class="" style="margin: 0px; line-height: normal;"><span class="" style="font-size: 12px;"><span class="Apple-tab-span" style="white-space: pre;"></span> }</span></div>
<div class="" style="margin: 0px; line-height: normal;"><span class="" style="font-size: 12px;"> }<span class="Apple-tab-span" style="white-space: pre;">
</span></span></div>
<div class="" style="margin: 0px; line-height: normal;"> }</div>
<div class="" style="margin: 0px; line-height: normal;">}</div>
</div>
<div class=""><br class="">
</div>
<div class="">The rationale is in one population, the ligand is not bound and the PHE occupies the pocket. In the other, the ligand is bound the the PHE is out somewhere that I cannot observe and model, thus I did not model that PHE-out conformation. The refinement
went to completion with the occupancies refined to 0.68 and 0.32 for BLIG and APHE, respectively. R-values dropped by a few tenths of a percent. Please let me know if there’s a better way to do this refinement.</div>
<div class=""><br class="">
</div>
<div class="">@Pavel, there’s repulsion between the selections if I don’t put an altloc on the ligand as suggested by Oleg.</div>
<div class=""><br class="">
</div>
<div class="">Jason</div>
<div class=""><br class="">
<div>
<blockquote type="cite" class="">
<div class="">On Mar 1, 2018, at 4:59 PM, Pavel Afonine <<a href="mailto:pafonine@lbl.gov" class="">pafonine@lbl.gov</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div text="#000000" bgcolor="#FFFFFF" class="">Hi Jason,<br class="">
<br class="">
Oleg is correct, you can assign the ligand and PHE different non-blanc altlocs (say A and B). In this case PHE and ligand will not 'see' each other via repulsion term of restraints. Also, their occupancies will be refined independently and will be constrained
to be within 0 and 1, but it will not make sure that they add up to 1. If you want them to add up to 1, then make a parameter file with a constrained group as described in one of examples here:<br class="">
<br class="">
"13 typical occupancy refinement scenarios and available options in phenix.refine"<br class="">
<a class="moz-txt-link-freetext" href="https://na01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fphenix-online.org%2Fnewsletter%2F&data=02%7C01%7Cjason.phan%40vanderbilt.edu%7C40b308da75d345eb3cf408d57fc82098%7Cba5a7f39e3be4ab3b45067fa80faecad%7C0%7C1%7C636555419887467714&sdata=5VYe7ZcJpoBjUZVV7VTStHt2sRTUn1W4z%2FztPNHyGpE%3D&reserved=0">http://phenix-online.org/newsletter/</a><br class="">
<br class="">
You can do it in the GUI too.<br class="">
<br class="">
Let us know if you more questions or problems!<br class="">
<br class="">
Pavel<br class="">
<br class="">
<div class="moz-cite-prefix">On 3/2/18 02:57, Oleg Sobolev wrote:<br class="">
</div>
<blockquote type="cite" cite="mid:CADd3otGuVs-pmOoSsTeO6kSZM91g5ZVn8YnQkFUCNdW2KNA72w@mail.gmail.com" class="">
<div dir="ltr" class="">Hi Jason,
<div class=""><br class="">
</div>
<div class="gmail_extra">
<div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0
.8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word" class="">
<div style="margin:0px;line-height:normal" class="">
<div style="margin:0px;font-size:8px;line-height:normal" class=""><br class="">
</div>
</div>
<div class="">The ligand and APHE got pushed away from one another and off of their respective densities. Occupancies for both were reduced to 0.5 for the refinement. What am I doing wrong here? Does the ligand need to have an altloc as well? But I don’t see
an alternative conformation. </div>
</div>
</blockquote>
<div class=""><br class="">
</div>
<div class="">Yes, the ligand also needs an altloc, and it should be different from PHE altloc. This way the refinement will know that these atoms don't see each other. You don't have to put two alternative conformations for the ligand, you can have just one
with partial occupancy (<1.0) which will mean that the ligand is not always there.</div>
<div class=""> </div>
<div class="">Best regards,</div>
<div class="">Oleg Sobolev.</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px
0.8ex;border-left:1px solid
rgb(204,204,204);padding-left:1ex">
<div style="word-wrap:break-word" class=""><br class="">
<div class="">
<blockquote type="cite" class="">
<div class="">On Feb 28, 2018, at 5:35 PM, van den Bedem, Henry <<a href="mailto:vdbedem@slac.stanford.edu" target="_blank" moz-do-not-send="true" class="">vdbedem@slac.stanford.edu</a>> wrote:</div>
<br class="m_8858005827093028639Apple-interchange-newline">
<div class="">
<div class="m_8858005827093028639WordSection1" style="font-family:Helvetica;font-size:14px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
If the phe is a gatekeeper, shouldn’t the ligand be refined at partial occupancy; i.e. occupancies not necessarily have to ‘add up’ as you suggest? Maybe this link is helpful too:<span class="m_8858005827093028639Apple-converted-space"> </span><a href="https://na01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.biorxiv.org%2Fcontent%2Fearly%2F2018%2F01%2F24%2F253419&data=02%7C01%7Cjason.phan%40vanderbilt.edu%7C40b308da75d345eb3cf408d57fc82098%7Cba5a7f39e3be4ab3b45067fa80faecad%7C0%7C1%7C636555419887467714&sdata=%2BeF9t7TNZ72bkVoY8GqaBvuNd1EOp%2BaRRmhlcB8OZ1g%3D&reserved=0" style="color:rgb(149,79,114);text-decoration:underline" target="_blank" moz-do-not-send="true" class="">www.biorxiv.org/content/<wbr class="">early/2018/01/24/253419</a></div>
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
</div>
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
H</div>
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
</div>
<div style="border-style:solid none
none;border-top-width:1pt;border-top-color:rgb(181,196,223);padding:3pt
0in 0in" class="">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
<b class=""><span style="font-size:12pt" class="">From:<span class="m_8858005827093028639Apple-converted-space"> </span></span></b><span style="font-size:12pt" class=""><<a href="mailto:phenixbb-bounces@phenix-online.org" style="color:rgb(149,79,114);text-decoration:underline" target="_blank" moz-do-not-send="true" class="">phenixbb-bounces@<wbr class="">phenix-online.org</a>>
on behalf of "Phan, Jason" <<a href="mailto:jason.phan@vanderbilt.edu" style="color:rgb(149,79,114);text-decoration:underline" target="_blank" moz-do-not-send="true" class="">jason.phan@vanderbilt.edu</a>><br class="">
<b class="">Date:<span class="m_8858005827093028639Apple-converted-space"> </span></b>Wednesday, February 28, 2018 at 1:30 PM<br class="">
<b class="">To:<span class="m_8858005827093028639Apple-converted-space"> </span></b>"<a href="mailto:phenixbb@phenix-online.org" style="color:rgb(149,79,114);text-decoration:underline" target="_blank" moz-do-not-send="true" class="">phenixbb@phenix-online.<wbr class="">org</a>"
<<a href="mailto:phenixbb@phenix-online.org" style="color:rgb(149,79,114);text-decoration:underline" target="_blank" moz-do-not-send="true" class="">phenixbb@phenix-online.org</a>><br class="">
<b class="">Subject:<span class="m_8858005827093028639Apple-converted-space"> </span></b>[phenixbb] Overlapping ligand-protein side chain density</span></div>
</div>
<span class="">
<div class="">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
</div>
</div>
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
<a name="m_8858005827093028639__MailOriginalBody" moz-do-not-send="true" class="">All, <span class="m_8858005827093028639Apple-converted-space"> </span></a></div>
<div class="">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
<span class=""> </span></div>
</div>
<div class="">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
<span class="">A piece of a ligand and the side chain of a “gate-keeper” Phe occupy the same space with well-defined features observed for both (resolution is 1.7 A). It looks about 50:50. How do you refine both ligand and protein side chain in this case? A
couple of phenixbb suggestions for dealing with ligand-ligand overlapping density have been considered. With the first suggestion, the two entities still clashed and moved apart off of their respective densities. The second suggestion is not applicable in
this case since the other molecule is not a ligand but part of the protein but it was tested anyway. Although there was no bumping, the occupancies don’t add up, resulting in a big blob of negative density around the ligand piece. </span></div>
</div>
<div class="">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
<span class=""> </span></div>
</div>
<div class="">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
<span class="">------------------------------</span></div>
</div>
<div class="">
<div class="">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
<span class="">Only two comments:</span></div>
</div>
<div class="">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
<span class=""> </span></div>
</div>
<ol start="1" style="margin-bottom:0in;margin-top:0in" type="a" class="">
<li class="MsoNormal" style="margin:0in
0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">
<span class="">At that resolution, constrained group occupancy refinement should work reasonably well (provided you can model the 2 entities). Then you also do not have clashes between the molecules, because Occ(A)+Occ(B)=1, meaning when one (A) is there, the
other one (B) is not. This works with refmac (external keyword file); if you need more sophisticated occupancy re/constraints SHELXL may offer more opportunities.</span>
</li><li class="MsoNormal" style="margin:0in
0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">
<span class="">There is no necessity for the two NCS copies of the binding site to look exactly the same (non-equivalent). Maybe there is a good reason/story (accessibility, contacts etc) for one site to be occupied differently than the other one. </span>
</li></ol>
<div class="">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
<span class=""> </span></div>
</div>
<div class="">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
<span class="">Best, BR</span></div>
</div>
</div>
<div class="">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
<span class="">——————————————</span></div>
</div>
<div class="">
<div class="">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
<span class=""><span style="font-size:10.5pt;font-family:Helvetica" class="">Modeling two molecules that occupy overlapping binding sites in a structure simply involves designating them as alternate conformers, with the same chain and residue number, and an
occupancy that sums to 1.0. For example, if you have an AMP and an ADP that occupy the same binding site, you would define them as</span></span></div>
</div>
<div class="">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
<span class=""><span style="font-size:10.5pt;font-family:Helvetica" class=""> </span></span></div>
</div>
<div class="">
<div class="">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
<span class=""><span style="font-size:10.5pt;font-family:Courier" class="">AAMP B 501</span></span><span class=""><span style="font-size:10.5pt;font-family:Helvetica" class=""></span></span></div>
</div>
<div class="">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
<span class=""><span style="font-size:10.5pt;font-family:Courier" class="">BADP B 501</span></span><span class=""><span style="font-size:10.5pt;font-family:Helvetica" class=""></span></span></div>
</div>
</div>
<div class="">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
<span class=""><span style="font-size:10.5pt;font-family:Helvetica" class=""> </span></span></div>
</div>
<div class="">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
<span class=""><span style="font-size:10.5pt;font-family:Helvetica" class="">and initially set the occupancies for the atoms in each conformer to the ratio (50:50, 30:70, etc.) that you observe in the density.</span></span></div>
</div>
<div class="">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
<span class=""><span style="font-size:10.5pt;font-family:Helvetica" class=""> </span></span></div>
</div>
<div class="">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
<span class=""><span style="font-size:10.5pt;font-family:Helvetica" class="">Refinement in this manner is straightforward in PHENIX.</span></span></div>
</div>
<div class="">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
<span class=""><span style="font-size:10.5pt;font-family:Helvetica" class=""> </span></span></div>
</div>
<div class="">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
<span class=""><span style="font-size:10.5pt;font-family:Helvetica" class="">Diana</span></span></div>
</div>
<div class="">
<div style="margin:0in 0in
0.0001pt;font-size:11pt;font-family:Calibri,sans-serif" class="">
<span class=""><span style="font-size:10.5pt;font-family:Helvetica" class="">----------------------------</span></span></div>
</div>
</div>
</span></div>
</div>
</blockquote>
</div>
<br class="">
</div>
<br class="">
______________________________<wbr class="">_________________<br class="">
phenixbb mailing list<br class="">
<a href="mailto:phenixbb@phenix-online.org" moz-do-not-send="true" class="">phenixbb@phenix-online.org</a><br class="">
<a href="https://na01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fphenix-online.org%2Fmailman%2Flistinfo%2Fphenixbb&data=02%7C01%7Cjason.phan%40vanderbilt.edu%7C40b308da75d345eb3cf408d57fc82098%7Cba5a7f39e3be4ab3b45067fa80faecad%7C0%7C1%7C636555419887467714&sdata=nIHAg7svGzV23ChHuKdhLwdoEow358Mz7wadZ3TczWk%3D&reserved=0" rel="noreferrer" target="_blank" moz-do-not-send="true" class="">http://phenix-online.org/<wbr class="">mailman/listinfo/phenixbb</a><br class="">
Unsubscribe: <a href="mailto:phenixbb-leave@phenix-online.org" moz-do-not-send="true" class="">
phenixbb-leave@phenix-online.<wbr class="">org</a><br class="">
</blockquote>
</div>
<br class="">
</div>
</div>
<br class="">
<fieldset class="mimeAttachmentHeader"></fieldset> <br class="">
<pre wrap="" class="">_______________________________________________
phenixbb mailing list
<a class="moz-txt-link-abbreviated" href="mailto:phenixbb@phenix-online.org">phenixbb@phenix-online.org</a>
<a class="moz-txt-link-freetext" href="https://na01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fphenix-online.org%2Fmailman%2Flistinfo%2Fphenixbb&data=02%7C01%7Cjason.phan%40vanderbilt.edu%7C40b308da75d345eb3cf408d57fc82098%7Cba5a7f39e3be4ab3b45067fa80faecad%7C0%7C1%7C636555419887467714&sdata=nIHAg7svGzV23ChHuKdhLwdoEow358Mz7wadZ3TczWk%3D&reserved=0">http://phenix-online.org/mailman/listinfo/phenixbb</a>
Unsubscribe: <a class="moz-txt-link-abbreviated" href="mailto:phenixbb-leave@phenix-online.org">phenixbb-leave@phenix-online.org</a></pre>
</blockquote>
<br class="">
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</div>
</body>
</html>