Hi

1) refmac's default option has been "generate all" starting from version 5 (it is around 14 years). The reason is for most tests generating hydrogens gave good R or geometry. And if you look at the original paper by Ramachandran (I think it was pulished in 1963) you see using hard ball approximation and including all atoms generated plot is very similar to what we know as Ramachandran plot. A little bit more work on included hydrogen plus phi/psi torsion angle restraints (in free form, not from Ramachandran plot) should generate faithfully Ramchandran plot. It needs to be tested properly and carefully. refmac version 4 did not have hydrogen generation option.
2) ccp4i's option has always been use hydrogens if present in file (I personally do not like this option: this option should be used with care)
3) If there is no non-bondedn ... others then it is very likely that hydrogens were not used. 
4) I am not sure that there is a fullproof way to know if hydrogens were added or not.

Regards
Garib


On 7 Mar 2012, at 15:00, Paul Smith wrote:

Hello CCP4bb,

Firstly, thanks to all for your comments.  However, I'm still unsure how to sort all of this riding hydrogen business out.

Robbie's comments seem particularly apt:

"Because there were some reporting errors in the past (http://proteincrystallography.org/ccp4bb/message18808.html) it is hard to tell from the PDB when refinement with hydrogens became hip."

Is there any foolproof way to know if a recently deposited file was refined with riding hydrogens in REFMAC, especially since some such structures do not yet have publications associated with them? How about the value of NON-BONDED CONTACTS OTHERS? If that is other than NULL, does that mean riding hydrogens were present?

Also, how about refmac version numbers?  Is there a clear demarcation when riding hydrogens a) became available, b) became default, or c) became default in the CCP4i GUI?

Thanks again,

--Paul


From: Robbie Joosten <[log in to unmask]>
To: [log in to unmask]
Sent: Tuesday, March 6, 2012 4:26 AM
Subject: Re: [ccp4bb] REFMAC Riding Hydrogens

Hi Everyone,
 
Pavel’s statement is likely a bit of an exaggeration, but he has a valid (yet hard to prove point). The default in CCP4i was (and is?) to use hydrogens only if present in the input file. This is IMO not a safe default.
Because there were some reporting errors in the past (http://proteincrystallography.org/ccp4bb/message18808.html) it is hard to tell from the PDB when refinement with hydrogens became hip. Discussions on this BB show that at the use of riding hydrogens is still not fully accepted, especially at low resolution (where they actually help most).
 
Cheers,
Robbie
 
From: CCP4 bulletin board [mailto:[log in to unmask]] On Behalf Of Pavel Afonine
Sent: Monday, March 05, 2012 21:53
To: [log in to unmask]
Subject: Re: [ccp4bb] REFMAC Riding Hydrogens
 
Dear Tim,
 
good catch, thanks; I could craft that phrase more carefully! Although often it may not be quite fair to take phrases out of context: this newsletter article was written in the context of macromolecular refinement. And yes, "recently" may be a broad term -:)
 
All the best,
Pavel
On Mon, Mar 5, 2012 at 12:45 PM, Tim Gruene <[log in to unmask]> wrote:
-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

Dear Pavel,

you may want to add to the structures mentioned in [1] one or two
organic structures present in the Cambridge Database.

"Until recently it was customary to ignore hydrogen atoms throughout the
process of crystallographic X-­ray structure determination." [1]

'recently' as in 1997 [2]? Even though 1997 is probably a poor
estimation of the corresponding year...

Cheers,
Tim


[1] "On contribution of hydrogen atoms to X-ray scattering"
http://www.phenix-online.org/newsletter/
[2] http://shelx.uni-ac.gwdg.de/SHELX/shelx.pdf

On 03/05/2012 09:14 PM, Pavel Afonine wrote:
> Hi,
>
> On Mon, Mar 5, 2012 at 11:52 AM, Matthew Franklin <[log in to unmask]>wrote:
>
>> Adding the riding hydrogens generally gives you some improvement in R
>> factors even with a good quality (i.e. stereochemically correct) model.
>>
>
> and here are the results of more or less systematic test that prove this:
>
> see "On contribution of hydrogen atoms to X-ray scattering"
> here:
> http://www.phenix-online.org/newsletter/
>
> Pavel
>
- --
Dr Tim Gruene
Institut fuer anorganische Chemie
Tammannstr. 4
D-37077 Goettingen

GPG Key ID = A46BEE1A
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.11 (GNU/Linux)
Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org/

iD8DBQFPVSXkUxlJ7aRr7hoRAm1TAJ9Hyfhkl3yhD5QSKw9I4RSK58m0fACgmlxk
YGILzeMam/3gQVmCeh0vQ8k=
=3m2J
-----END PGP SIGNATURE-----
 



Garib N Murshudov 
Structural Studies Division
MRC Laboratory of Molecular Biology
Hills Road 
Cambridge 
CB2 0QH UK
Email: garib@mrc-lmb.cam.ac.uk 
Web http://www.mrc-lmb.cam.ac.uk