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
jo hanna wrote: > 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!
In order to be able to help you, you need to post your full inputs. See http://salilab.org/modeller/manual/node11.html It is essentially impossible to reproduce your problem if you only post part of your alignment.
If your sequences are confidential or just rather large, you can send them directly to me rather than the list.
> 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 also be interested in seeing the input files that result in this error, as I'm not sure what you're doing wrong here (but Modeller should certainly be fixed to give a more understandable error message in this case).
Ben Webb, Modeller Caretaker
participants (2)
-
jo hanna
-
Modeller Caretaker