<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
  <meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
Hi Carsten,<br>
<br>
they should definitely add up to 1. If the PDB file contains altloc
identifiers, then those atoms should be picked up automatically and you
don't need to specify any selections. For example, make sure you have
something like this in your PDB file (the starting occupancies can be
any):<br>
<br>
<font face="Courier New, Courier, monospace">ATOM&nbsp;&nbsp;&nbsp;&nbsp; 49&nbsp; N&nbsp;&nbsp; MET A&nbsp;
10&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3.832 -11.001&nbsp;&nbsp; 0.570&nbsp; 1.00 11.55&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; N&nbsp; <br>
ATOM&nbsp;&nbsp;&nbsp;&nbsp; 50&nbsp; CA AMET A&nbsp; 10&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2.523 -10.377&nbsp;&nbsp; 0.709&nbsp; 0.70
12.44&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; C&nbsp; <br>
ATOM&nbsp;&nbsp;&nbsp;&nbsp; 51&nbsp; CA BMET A&nbsp; 10&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2.574 -10.279&nbsp;&nbsp; 0.718&nbsp; 0.30
12.42&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; C&nbsp; <br>
ATOM&nbsp;&nbsp;&nbsp;&nbsp; 54&nbsp; CB AMET A&nbsp; 10&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2.580&nbsp; -9.129&nbsp;&nbsp; 1.601&nbsp; 0.70
14.67&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; C&nbsp; <br>
ATOM&nbsp;&nbsp;&nbsp;&nbsp; 55&nbsp; CB BMET A&nbsp; 10&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2.822&nbsp; -8.921&nbsp;&nbsp; 1.394&nbsp; 0.30
13.99&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; C&nbsp; <br>
ATOM&nbsp;&nbsp;&nbsp;&nbsp; 60&nbsp; CG AMET A&nbsp; 10&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.245&nbsp; -8.462&nbsp;&nbsp; 1.951&nbsp; 0.70
18.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; C&nbsp; <br>
ATOM&nbsp;&nbsp;&nbsp;&nbsp; 61&nbsp; CG BMET A&nbsp; 10&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.545&nbsp; -8.166&nbsp;&nbsp; 1.727&nbsp; 0.30
17.90&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; C&nbsp; <br>
ATOM&nbsp;&nbsp;&nbsp;&nbsp; 66&nbsp; SD AMET A&nbsp; 10&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.659&nbsp; -7.204&nbsp;&nbsp; 0.780&nbsp; 0.70
19.47&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; S&nbsp; <br>
ATOM&nbsp;&nbsp;&nbsp;&nbsp; 67&nbsp; SD BMET A&nbsp; 10&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.792&nbsp; -6.615&nbsp;&nbsp; 2.600&nbsp; 0.30
26.45&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; S&nbsp; <br>
ATOM&nbsp;&nbsp;&nbsp;&nbsp; 68&nbsp; CE AMET A&nbsp; 10&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.520&nbsp; -5.730&nbsp;&nbsp; 1.337&nbsp; 0.70
20.21&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; C&nbsp; <br>
ATOM&nbsp;&nbsp;&nbsp;&nbsp; 69&nbsp; CE BMET A&nbsp; 10&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2.371&nbsp; -5.566&nbsp;&nbsp; 1.271&nbsp; 0.30
26.12&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; C&nbsp; <br>
ATOM&nbsp;&nbsp;&nbsp;&nbsp; 76&nbsp; C&nbsp;&nbsp; MET A&nbsp; 10&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2.021 -10.028&nbsp; -0.671&nbsp; 1.00
10.57&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; C&nbsp; <br>
ATOM&nbsp;&nbsp;&nbsp;&nbsp; 77&nbsp; O&nbsp;&nbsp; MET A&nbsp; 10&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2.780&nbsp; -9.724&nbsp; -1.578&nbsp; 1.00
11.30&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; O&nbsp; </font><br>
<br>
and then everything should be OK.<br>
Unfortunately, I can't tell much about the selections you use without
looking at the PDB file, so in your case I would just try the above
suggestion.<br>
<br>
Pavel.<br>
<br>
<br>
On 2/4/2009 12:15 PM, Schubert, Carsten [PRDUS] wrote:
<blockquote
 cite="mid:889E20E3723D0C45A493DA5DA87D854E03413C45@JNJUSRAGMS02.na.jnj.com"
 type="cite">
  <pre wrap="">Hi,

I am trying to sort out some issues with the constrained occupancy refinement in phenix 1.4.3. The structure I am working on is fairly high resolution ~1.4A and it looks as if a couple of the alternative conformations would benefit from occupancy refinement, no hydrogens are present yet.

I defined the constrained groups like this:

    occupancies {
      individual = None
      remove_selection = None
      constrained_group {
        selection = chain A and resid 66 and altid A
        selection = chain A and resid 66 and altid B
      }
      ....

I expected that the sum of occupancies would be restrained to 1, which they should according to the manual, but the refinement gives me occupancies of 0.7 for altid A and 0.34 for altid B, starting values were both 0.5.  Similar discrepancies were also present in other groups, like 0.5 A, 0.37 B etc. The latter is the largest deviation from 1 I have seen in thus far. The alternate conformations were generated in coot, by splitting at the CA position.

Is this a rounding error I just have to contend with or is there some other parameter I need to tweak?

Cheers,

        Carsten

_______________________________________________
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="http://www.phenix-online.org/mailman/listinfo/phenixbb">http://www.phenix-online.org/mailman/listinfo/phenixbb</a>
  </pre>
</blockquote>
</body>
</html>