Sorry for the delay, it's hard for me to do non-trivial mails while I'm 
travelling.

On 23/11/2017 15:52, Toon Van Thillo wrote:


For my master thesis I am busy with neutron crystallography. For refinement, the software used is /phenix.refine /and COOT for manual corrections.

However, it seems that deuterium atoms will not correctly refine because there are no restraints available in the COOT dictionary.

Yes. It's using the standard Refmac dictionary.


Thus we set out to introduce restraints to COOT by creating cif files for each amino acid with restraints for deuteriums.

This seems to work for perdeuterated proteins, but not for those that are 
H/D-exchanged.

OK, makes sense.


The cause of this problem seems to be in the way H/D positions are defined: /phenix.refine/ assigns H and D to alternative conformations A and B,

while the other part of the residue remains normal (eg. has no alternative 
conformation assigned).


ATOM   1080  N          ALA A  67      16.629  -3.415   9.989  1.00  7.92       
    N
ATOM   1081  CA        ALA A  67      15.560  -3.601   9.012  1.00 10.82        
   C
ATOM   1082  C          ALA A  67      14.894  -4.957   9.211  1.00  8.09       
    C
ATOM   1083  O          ALA A  67      14.918  -5.812   8.329  1.00  6.06       
    O
ATOM   1084  CB        ALA A  67      16.119  -3.459   7.599  1.00  7.64        
   C
ATOM   1085  HA       ALA A  67      14.801  -2.842   9.167  1.00  9.05         
  H
ATOM   1086  HB1     ALA A  67      15.405  -3.583   6.954  1.00  9.12          
 H
ATOM   1087  HB2     ALA A  67      16.502  -2.575   7.490  1.00  9.97          
 H
ATOM   1088  HB3     ALA A  67      16.806  -4.129   7.454  1.00  9.88          
 H
ATOM   1089  H       AALA A  67      17.253  -4.214  10.089  0.50  9.73         
  H
ATOM   1090  D       BALA A  67      17.253  -4.214  10.089  0.50  9.73         
  D

Yes.


The restraints we introduced in the cif files do work when deuterium and the corresponding hydrogen have not been assigned to an alternative conformation.

Also makes sense.


When they are, they (both hydrogen and deuterium) refuse to refine. Does 
someone know a way around this problem?



Phenix neutron cif files are not well handled by coot. The expectation that there is only one dictionary atom for any given atom name is deeply entrenched in the way Coot deals with cif files, and stores and uses dictionary information. If it had been straightforward to fix this I would have done (Pavel had told me about this a while ago) - but it's not :-(. I did make Deuterium atoms pink though.

Paul.

Reply via email to