That is definitely some strange behavior going on there. The position of the N-H atom is in part determined by the position of carbonyl carbon in the preceding residue, but usually it would have the matching altID associated with it. In the example you've sent, it looks like Reduce may not be adding this hydrogen correctly.<div>
<br></div><div>Would you be willing to send me your entire PDB file directly (NOT to the list) so that I could see why Reduce is behaving that way?</div><div><br></div><div>I'm not sure about the occupancy problem, though. That should only happen if there were multiple positions for that hydrogen.</div>
<div><br></div><div>Jeff<br><div><br><div class="gmail_quote">On Fri, Oct 26, 2012 at 4:07 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
I did not look closely yet but my understanding is that its position is dictated by the next to it residue that has non-blank altlocs, and therefore this hydrogen inherits the conformers. So this behavior is expected. I hoped Jeff will comment on this as he was the first who explained it to me.<span class="HOEnZb"><font color="#888888"><br>
<br>
Pavel</font></span><div class="HOEnZb"><div class="h5"><br>
<br>
On 10/26/12 11:39 AM, Nathaniel Echols wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
On Thu, Oct 25, 2012 at 8:05 PM, Alexander Scouras <<a href="mailto:scouras@berkeley.edu" target="_blank">scouras@berkeley.edu</a>> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
After a series of residues (1+) with alternate conformations, the H atom of the next residue seems to lock its occupancy to that of the final conformation of the previous residue. �Hence I'm getting several H's with 0.1-0.5 occupancy on single conformation residues, where the entire rest of the residue has occupancy 1. �This happens after every altconf series, and showed up in two separate structures I've worked on this week.<br>
<br>
Within a series of 2+ residues with altconfs, the H's refine correctly with their altconf groups. �So it isn't just that all subsequent residue's H's are getting overridden.<br>
<br>
This behavior seems to have started with 1.8.1-1168, or at least previous refinements of the same structure didn't do this. �64 bit build on Retina MBP running OS X 10.8.2.<br>
<br>
I couldn't google any reference to this behavior, so assuming it's a new bug I can send along source files.<br>
</blockquote>
Very interesting - this is also reproducible with a simple test case<br>
and the latest code. �Judging from the revision history, it looks like<br>
this probably would have happened with version 1.8 too. �I am not<br>
entirely sure what the code is doing, but I think it is partly based<br>
on the assumption that if you have an alternate conformation for the C<br>
and/or CA atoms of a residue, the H atom for the next residue will<br>
also have alternate conformations. �This is in fact what<br>
phenix.ready_set will do - if you re-run it on your input PDB file,<br>
you will have A/B/C conformations for Met30 H. �So I'm inclined to<br>
suggest simply doing this (it is the most chemically sensible<br>
approach), but the behavior of phenix.refine when no alternate<br>
conformer is present for H does appear to be a bug.<br>
<br>
-Nat<br>
______________________________<u></u>_________________<br>
phenixbb mailing list<br>
<a href="mailto:phenixbb@phenix-online.org" target="_blank">phenixbb@phenix-online.org</a><br>
<a href="http://phenix-online.org/mailman/listinfo/phenixbb" target="_blank">http://phenix-online.org/<u></u>mailman/listinfo/phenixbb</a><br>
</blockquote>
<br>
______________________________<u></u>_________________<br>
phenixbb mailing list<br>
<a href="mailto:phenixbb@phenix-online.org" target="_blank">phenixbb@phenix-online.org</a><br>
<a href="http://phenix-online.org/mailman/listinfo/phenixbb" target="_blank">http://phenix-online.org/<u></u>mailman/listinfo/phenixbb</a><br>
</div></div></blockquote></div><br></div></div>