[modeller_usage] Problems with modelling with ligand
To:
Subject: [modeller_usage] Problems with modelling with ligand
From: jo hanna <>
Date: Tue, 3 Jan 2006 09:19:51 +0000
Hello
I apologise if this message is posted more than once, I am having difficulties with my email system.
My questions regards the process for modelling with ligands.
>From the tutorial i understand that you have the ligands in the pdb
file template and then you manually alter the alignment to take account
of these ligands, for me that means adding /..w* to the end of the
structure and seqeunce alignments (i'm unsure if you add it to the
seqeunce, if you want it included). This also means changing the number
of residues in the structure file in the alignment to include the 3
BLK/water items, as even if they are included in the pdb prior to
alignment they are not included in the residue number line in the
alignment (i.e. in alignment below changing original 514 to 517).
I then do the modelling using the alignment, as appended below and the
model.top script (7v7 user!) also below and get the following error
which i do not understand.
check_ali___> Checking the sequence-structure alignment.
Implied target CA(i)-CA(i+1) distances longer than 8.0 angstroms:
ALN_POS TMPL RID1 RID2 NAM1 NAM2 DIST
----------------------------------------------
134 1 134
138 G
E 8.407
END OF TABLE
delete__442E> One or more atoms absent from MODEL: O: 522: C: 522: N: 523: CA: 523:
delete__442E> One or more atoms absent from MODEL: O: 523: C: 523: N: 524: CA: 524:
delete__442E> One or more atoms absent from MODEL: O: 524: C: 524: N: 525: CA: 525:
delete__443_> Restraints marked for deletion were removed.
Total number of restraints before, now:
51247 47550
nselat__596E> Atom index is out of range (iatm,natm): 4033 4032
recover____E> ERROR_STATUS >=
STOP_ON_ERROR:
1 1
Dynamically allocated memory
at finish
[B,kB,MB]: 26700639
26074.844 25.464
Starting
time
: 2006/01/02 21:50:15.029
Closing
time
: 2006/01/02 21:50:35.218
Total CPU time
[seconds]
: 20.18
I unerstand distance table and know which region of my alignment this
points to but I do not understand the missing atoms or the out of range
error, as I have 4142 atoms and no missing atoms. Could anyone help
with where I am going wrong?! And also tell me if my procedure is
correct, should the alignment automatically include the 3 BLK/water
residues in the resdiue number if it is in the pdb file? and if i want
the ligands in the output model is it correct to state /..w* at the end
of this sequence in the alignment too? Also in the alignment is it
correct for the seqeunce line to read:
sequence:CcHh1: : : :
:::-1.00:-1.00 ?
INCLUDE
SET ALNFILE = 'CcHh1-Bamu.ali'
SET KNOWNS = 'Bamu'
SET SEQUENCE = 'CcHh1'
SET STARTING_MODEL = 1
SET ENDING_MODEL = 1
SET HETATM_IO = ON, WATER_IO = ON
CALL ROUTINE = 'model'