[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [ccp4bb]: warp/refmac5 bug?



***  For details on how to be removed from this list visit the  ***
***    CCP4 home page http://www.dl.ac.uk/CCP/CCP4/main.html    ***

> Dear all,
>         I'm getting an error message when auto-building hydrogens into
> my structure with refmac5 and then cycling through to warp to find new
> waters. The message,
>
> Illegal atom name: 1HG1 ILE A
>
> arrives while reading in the pdb generated by refmac5 to warp. Any
> ideas?
>

That's Victor's code so I am not sure, but after a quick look at source
I can see that the valid atom names are:

      DATA AT/' H','HE','LI','BE',' B',' C',' N',' O',' F','NE',
     *        'NA','MG','AL','SI',' P',' S','CL','AR',' K','CA',
     *        'SC','TI',' V','CR','MN','FE','CO','NI','CU','ZN',
     *        'GA','GE','AS','SE','BR','KR','RB','SR',' Y','ZR',
     *        'NB','MO','TC','RU','RH','PD','AG','CD','IN','SN',
     *        'SB','TE',' I','XE','CS','BA','LA','CE','PR','ND',
     *        'PM','SM','EU','GD','TB','DY','HO','ER','TM','YB',
     *        'LU','HF','TA',' W','RE','OS','IR','PT','AU','HG',
     *        'TL','PB','BI','PO','AT','RN','FR','RA','AC','TH',
     *        'PA',' U','NP','PU','AM','CM','BK','CF','ES','FM'/

ARP the checks:

IF(ATNAM(1:2).EQ.AT(I)) THEN
    .....

ie '1HG1' is an illegal name.
' HG1' is OK.

I guess we have to fix that to get compatible with refmac5.

Again, thanx to Garib, refmac5 is/will be handling arp/warp issues
much better. We also have to make changes in arp_warp,
which have started but are not finished or ready for release.
So, for a couple of months more I would restrain myself to refmac4,
when used together with arp_warp, at least for 'unusual'
options like that or at the first sight of arp/warp - refmac5
compatibility trouble.

        Tassos