Hi all,
I'm tring to align an amino acid sequence with a protein (solved by X-ray) which conteins 6 chains in its PDB file. for my job I use the script:
>SET ATOM_FILES_DIRECTORY = './../atom_files' >READ_MODEL FILE = '1icf.pdb' >SEQUENCE_TO_ALI ALIGN_CODES = '1icf', ATOM_FILES = '1icf' >READ_ALIGNMENT FILE = 't.ali', ALIGN_CODES = ALIGN_CODES 't', ADD_SEQUENCE = ON >ALIGN2D >WRITE_ALIGNMENT FILE='t-1fh0.ali', ALIGNMENT_FORMAT='PIR' >WRITE_ALIGNMENT FILE='t-1fh0.pap', ALIGNMENT_FORMAT='PAP
but it don't work end the output is:
>iatmcls_286W> MODEL atom not classified: THR:OT1 THR >iatmcls_286W> MODEL atom not classified: THR:OT2 THR >iatmcls_286W> MODEL atom not classified: VAL:OT1 VAL >iatmcls_286W> MODEL atom not classified: VAL:OT2 VAL >iatmcls_286W> MODEL atom not classified: THR:OT1 THR >iatmcls_286W> MODEL atom not classified: THR:OT2 THR >iatmcls_286W> MODEL atom not classified: VAL:OT1 VAL >iatmcls_286W> MODEL atom not classified: VAL:OT2 VAL >iatmcls_286W> MODEL atom not classified: SER:OT1 SER >iatmcls_286W> MODEL atom not classified: SER:OT2 SER >iatmcls_286W> MODEL atom not classified: SER:OT1 SER >iatmcls_286W> MODEL atom not classified: SER:OT2 SER >read_al_375E> Unknown residue type,position,sequence: 1 2 >recover____E> ERROR_STATUS >= STOP_ON_ERROR: 1 1
any ideas?
Thanks!
G. **********************************
Dr. Gianluca Croce Ph.D. Student DISTA - University Amedeo Avogadro Corso Borsalino 54-15100 Alessandria-Italy Tel_office: +390131287414 Fax: +390131287416