[phenixbb] Problem with Molprobity analysis in phenix 1.6.4
Mayer, Mark (NIH/NICHD) [E]
mayerm at mail.nih.gov
Tue Oct 5 18:32:27 PDT 2010
Hi Nat,
After the long discussion, mostly on CCP4BB, about riding H,
which covered different H-bond lengths, it would be really
good if the GUI and other Phenix documentation explained
exactly how this is handled.
I think using riding H as a geometry restraint during refinement
should be standard practice, but I think 'interpreting' those H atoms
in all but the very highest resolution X-ray structures is probably
not good.
Mark
________________________________________
From: Nathaniel Echols [nechols at lbl.gov]
Sent: Tuesday, October 05, 2010 3:56 PM
To: PHENIX user mailing list
Subject: Re: [phenixbb] Problem with Molprobity analysis in phenix 1.6.4
On Tue, Oct 5, 2010 at 11:53 AM, Mark A Saper <saper at umich.edu> wrote:
> If I take my well-refined structure as output by phenix (.pdb) and use it
> directly on the Molprobity web site, the clash score is much lower: 2.75 for
> 39 total bumps. Many of these involve residues or waters with high B
> factors.
> If I first remove all hydrogens from the phenix-refined .pdb file with a
> text editor, and upload it to the Molprobity web site, then add hydrogens,
> and then calculate clash score, the score is now 8.69 or 123 bumps for 14150
> atoms. This is *identical* to what phenix calculates. Therefore, why
> doesn't phenix feed the refined coordinates directly into the Molprobity
> (probe) calculations? Has this change been made in one of the nightly
> builds?
There is a discrepancy in hydrogen bond lengths in Reduce versus the
monomer library - I believe Reduce is using the distance to the center
of the nucleus instead of the center of the electron cloud.
phenix.refine will regularize the hydrogen geometry before refining
the non-hydrogen coordinates, so it is correcting the bond lengths to
be appropriate for X-ray refinement. It looks like they're getting
stripped off and replaced when we run the clashscore calculations, so
you end up with the longer bonds and thus more clashes. I think the
consensus among the group was that these are more appropriate for
interpreting (instead of refining) the model, but I'll double-check
this with everyone. The GUI could at least be more explicit about how
the calculations are done, however.
(By the way, a clashscore of 8.69 is pretty decent anyway.)
-Nat
_______________________________________________
phenixbb mailing list
phenixbb at phenix-online.org
http://phenix-online.org/mailman/listinfo/phenixbb
More information about the phenixbb
mailing list