Hello
I have been trying for some time now to model with 2 ligands and 1 water molecule and I feel that I have tried everything from the list, faq and tutorials and not got very far. I have posted a number of questions to the list but can still not get this to work, properly! I have managed to include the ligands and the water O into my output pdb by using the script and alignment below, however there is not recognition that this oxygen is a water and in the output model it is listed as BLK along with the other ligands. From reading posts in the user list I was expecting this water oxygen to be recognised as such, since i have used the w character to represent it and also HOH in the pdb file (the file has been edited to include only the protain, ligand and water that I am wishing to model around). I cannot type w at the end of the sequence alignment however as this causes modeller to fail with an out of range error. I would really appreciate some help as I have completelty exhaused all of my options now.
Thanks Jo
top INCLUDE SET ALNFILE = 'CcHh1-Bamu.pap', ALIGNMENT_FORMAT = 'PAP' SET KNOWNS = 'Bamu' SET SEQUENCE = 'CcHh1' SET STARTING_MODEL = 1 SET ENDING_MODEL = 1 SET HETATM_IO = on, WATER_IO = on CALL ROUTINE = 'model'
end of alignment (pap format) alignment _aln.pos 480 490 500 510 520 Bamu RKQLPEPDLTF/..w CcHh1 RNALPEPVIPA/...
produce model with water included as:
HETATM 4032 O BLK 525 37.988 89.184 33.629 1.00 25.09 2SG4034 END