<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Thanks - we could have an argument about what is "natural" here, but that's philosophical. <div><br></div><div>The only option that works for me so far is to move the other model to the neighboring unit cell and disable symmetry clashes. As I said, it's an ugly hack, but it seems to work. I do run into further issues because omit map calculation fails complaining about clashes, but that's a problem for another day.</div><div><br></div><div>Problem with alternate conformer based approach is that it's pain to merge the two models - by hand it would take too long, but I do have reasonable coding abilities so I wrote a little tool that merges two pdb files into one and marks models with altlocs. However, phenix.refine fails when trying to read my NCS definitions even though it appears to have validated them just fine (see log file excerpt below). Any suggestions - am I maybe defining NCS groups wrong?</div><div><br></div><div><div>Validating user-supplied NCS groups...</div><div> Validating:</div><div>ncs_group {</div><div> reference = "altloc A"</div><div> selection = "altloc B"</div><div>}</div><div> OK. All atoms were included in validated selection.</div></div><div> ...</div><div><br></div><div><div>==================== Process input NCS or/and find new NCS ====================</div><div><br></div><div>Number of NCS groups: 1</div><div><br></div><div>Traceback (most recent call last):</div><div> File "/ibbr/phenix-1.14-3228/build/../modules/phenix/phenix/command_line/refine.py", line 11, in <module></div><div> command_line.run(command_name="phenix.refine", args=sys.argv[1:])</div><div> File "/ibbr/phenix-1.14-3228/modules/phenix/phenix/refinement/command_line.py", line 106, in run</div><div> log=log)</div><div> File "/ibbr/phenix-1.14-3228/modules/phenix/phenix/refinement/driver.py", line 450, in __init__</div><div> ncs_groups_as_phil = iotbx.phil.parse(ncs_groups_str)</div><div> File "/ibbr/phenix-1.14-3228/modules/cctbx_project/iotbx/phil.py", line 81, in parse</div><div> process_includes=process_includes)</div><div> File "/ibbr/phenix-1.14-3228/modules/cctbx_project/libtbx/phil/__init__.py", line 2178, in parse</div><div> primary_parent_scope=result)</div><div> File "/ibbr/phenix-1.14-3228/modules/cctbx_project/libtbx/phil/parser.py", line 120, in collect_objects</div><div> scope_extract_call_proxy_cache=scope_extract_call_proxy_cache)</div><div> File "/ibbr/phenix-1.14-3228/modules/cctbx_project/libtbx/phil/parser.py", line 134, in collect_objects</div><div> words=collect_assigned_words(word_iterator, lead_word),</div><div> File "/ibbr/phenix-1.14-3228/modules/cctbx_project/libtbx/phil/parser.py", line 31, in collect_assigned_words</div><div> str(lead_word), lead_word.where_str()))</div><div>RuntimeError: Missing value for selection (input line 3)</div></div><div><br></div></div></div></div></div><br><div class="gmail_quote"><div dir="ltr">On Thu, Nov 1, 2018 at 7:50 PM Oleg Sobolev <<a href="mailto:osobolev@lbl.gov">osobolev@lbl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear Edwin,<div> <br></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 dir="ltr"><div><div></div><div>The alternative is probably to define the second orientation as an alternative conformation. Aside from aesthetic objections, I don't want to lose the ability to restrain NCS - I need two orientations to be identical structures, and I am not sure yet whether I can use NCS groups defined via alternative conformers.</div></div></div></blockquote><div><br></div><div>I think this would be the easiest way forward and it looks like a very natural use of alternative conformations. You can definitely define NCS groups and even automatic search works. I'm attaching a small dummy example I experimented on. Two overlapping helices: chain A with altloc A and occupancy 0.5, chain B with altloc B and occupancy 0.5. Nothing prevents you from defining NCS groups using "altloc A" and it should work, but I imagine there should always be equivalent selection in terms of "chain A".</div><div><br></div><div>Best regards,</div><div>Oleg Sobolev.</div><div><br></div></div></div></div>
</blockquote></div>