<!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 Christina,<br>
<br>
if H atoms were used in refinement then they should not be removed from
the model you are going to PDB deposit. Post-refinement model
manipulation is bad. This makes the reported statistics irreproducible.
The R-factor computed for a model with and without H atoms is
different. <br>
<br>
Different programs use different libraries of ideal X-H parameters
(some put them at longer "neutron" X-H distances, and some place at the
position where the electron cloud is expected - shorter X-H distance).
The B-factors of H atoms can be inherited from atoms they are bonded to
or can be multiplied by a scale from 1 to 1.5. There is a class of X-H
geometries where this position of H atoms is not unique: for example,
C-O-H group in Tyrosine residue, where H can rotate around C-O axis.
Some programs optimize the position of H atom to accommodate potential
hydrogen bond and some doesn't. <br>
<br>
Destroying is always easier than building: removing existing hydrogens
is a matter of running one command :<br>
<br>
phenix.reduce -trim model.pdb > model_noH.pdb<br>
<br>
and takes 1 or 2 seconds. However, to build the H atoms back the
*exact* same way they were before is not that straightforward. <br>
<br>
It's just counterintuitive why you need to manipulate the refined model
for which you have obtained all the statistics numbers? If you remove
hydrogen atoms then you need to update the statistics. By stripping off
H atoms you will get somewhat smaller file, but does the file size
matter these days? <br>
<br>
Pavel.<br>
<br>
<br>
On 2/18/10 6:55 AM, Christina Bourne wrote:
<blockquote cite="mid:532838.3649.qm@web30503.mail.mud.yahoo.com"
type="cite">
<style type="text/css"><!-- DIV {margin:0px;} --></style>
<div
style="font-family: times new roman,new york,times,serif; font-size: 12pt;">
<div>Alex, Pavel,<br>
In my experience when using a riding hydrogen model they are not
included in the output pdb file. It is only when they are added
explicitly that they are given output coordinates, which should of
course be kept unchanged. Can you please verify this??<br>
-Christina<br>
</div>
<div
style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><br>
<div style="font-family: arial,helvetica,sans-serif; font-size: 13px;"><font
face="Tahoma" size="2">
<hr size="1"><b><span style="font-weight: bold;">From:</span></b>
Pavel Afonine <a class="moz-txt-link-rfc2396E" href="mailto:PAfonine@lbl.gov"><PAfonine@lbl.gov></a><br>
<b><span style="font-weight: bold;">To:</span></b> PHENIX user
mailing list <a class="moz-txt-link-rfc2396E" href="mailto:phenixbb@phenix-online.org"><phenixbb@phenix-online.org></a><br>
<b><span style="font-weight: bold;">Sent:</span></b> Wed, February
17, 2010 8:06:40 PM<br>
<b><span style="font-weight: bold;">Subject:</span></b> Re:
[phenixbb] riding model and deposition<br>
</font><br>
Hi Alex,<br>
<br>
> I have refined a 1.3A structure using individual anisotropic
B-factors (Rfactor, 15.5%; Rfree, 17.7%). Additional use of the riding
model during refinement has improvement the R values (Rf, 14.5%; Rfree,
16.6%), <br>
<br>
this is not surprising at all -- hydrogens do contribute to the
scattering.<br>
<br>
> although bond length and angle outliers did increase (anyone know
why that would be?).<br>
> <br>
<br>
Different model - slightly different Xray/Restraints weights - slightly
different statistics.<br>
<br>
> I have decided to deposit the model obtained using the riding
model, but am unsure whether to include the hydrogen coordinates in the
deposition.<br>
> <br>
<br>
If H atoms were used in refinement and final statistics calculation
(such as R-factors, etc), then the H atoms have to be included into
depositing model. Otherwise your statistics will not be reproducible.<br>
<br>
> The Phenix.refine documentation advises that they are left in, as
they are needed to reproduce the refinement statistics. However,
reviewing the high resolution structures in the PDB I found that:<br>
> <br>
<br>
Reviewing statistics in the PDB you can find many interesting things,
such as Fobs labeled as Iobs (and Iobs labeled as Fobs), atoms with
unknown scattering types labeled X (which makes R-factor impossible to
compute), O-H distances in water longer than 10A, etc. one can go on
for a long list.<br>
<br>
> 1. In the case of structures refined using REFMAC or SHELX, the
riding model does not seem to be included, although its use during
refinement is mentioned in the pdb header.<br>
> <br>
<br>
Which is, in my opinion, not good.<br>
<br>
> 2. In the case of structures refined using PHENIX, use of the
riding model is sometimes indicated in the header, but hydrogen
coordinates only appear to be included in the pdb in the case of
structures refined to a resolution better than 1.3A.<br>
> <br>
<br>
This most likely indicates a manipulation with a file contained refined
model, which is, again, in my opinion, not good.<br>
<br>
> As far as I can understand the purpose of the riding model is to
provide additional restraints for refinement, so as long as its use is
made evident in the pdb header (and maybe the methods section of the
publication) that should be enough. Inclusion of hydrogens in the
model, when you can't see them seems inappropriate.<br>
> <br>
<br>
Different programs use different libraries of ideal X-H parameters.
Depending on the program, the B-factors of H atoms can be inherited
from atoms they are bonded to or can be multiplied by a scale from 1 to
1.5. These all makes adding H program-dependent, and therefore
statistics from such models may vary.<br>
<br>
I believe that any manipulation on the refined structure is bad.<br>
<br>
Pavel.<br>
<br>
<br>
_______________________________________________<br>
phenixbb mailing list<br>
<a moz-do-not-send="true" ymailto="mailto:phenixbb@phenix-online.org"
href="mailto:phenixbb@phenix-online.org">phenixbb@phenix-online.org</a><br>
<span><a moz-do-not-send="true" target="_blank"
href="http://phenix-online.org/mailman/listinfo/phenixbb">http://phenix-online.org/mailman/listinfo/phenixbb</a></span><br>
</div>
</div>
</div>
<br>
<pre wrap="">
<hr size="4" width="90%">
_______________________________________________
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://phenix-online.org/mailman/listinfo/phenixbb">http://phenix-online.org/mailman/listinfo/phenixbb</a>
</pre>
</blockquote>
</body>
</html>