Dear all,
I am refining a protein structure at 1.74 A resolution with phenix.refine
(version 1.10.1_2155) and I have noticed that some side chains don't remain
in the electronic density even after careful hand adjustment in which they
are almost perfect. It seems pretty clear where the side chain should
(must) be located, moreover it is consistent with two rotamers (either
rotamer tp60 - #4 in coot - or rotamer tp-100 - #9 in coot -, which are
approx. at 180 deg for chi3). In spite of that, after reciprocal space
refinement, I observe that the end of side chain rotates chi3 angle and
both N and O gets far from where they started (there is this difference in
the chi3 angle).
My concerns are about the proper interpretation of geometrical restraints
listed in .geo file, specially what the ideal value is, how weights are
applied to ideal value and why phenix doesn't reach convergence in
reciprocal space refinement. Below are the values from .geo for the
referred Gln:
dihedral pdb=" CB GLN B 21 "
pdb=" CG GLN B 21 "
pdb=" CD GLN B 21 "
pdb=" OE1 GLN B 21 "
ideal model delta sinusoidal sigma weight residual
0.00 23.41 -23.41 2 3.00e+01 1.11e-03 8.42e-01
is there a description for us to interpret what these numbers are? To my
understanding, the ideal values are established as being either 0 or 180
degrees (yet sinusoidal 2). I cannot say about how heavy the weights are on
these (in fact, I left my refinement to optimize weights), but anyway, if I
understand correctly, the ideal value (for chi3) should be rotamer
dependent. Taking Gln as an example, I see several rotamers do have values
close to either 0 or 180 degree, but what about rotamers 4 and 9? I might
imagine a similar situation for residues like Asn and His, though I have
not searched deep for these yet.
Thanks,
Renato
