You are here

Threading a sequence onto a template with ligands/metal ions included

3 posts / 0 new
Last post
Threading a sequence onto a template with ligands/metal ions included
#1

I want to include the GDP/GTP and magnesium ions in the threading process for a tubulin homology model I am trying to make (since I do believe this is the ideal way to start formulating the homology model). For some reason, however, only the GTP is included in the output pdb. I get an error message that says: 

Error: potential mismatch between sequence from alignment  and sequence from PDB!

However, as far as I can tell, the sequence for the template in the alignment is the same as the sequence in the PDB. Here are the associated files, command line input, and error message. Please let me know if there are any other files you need to help solve this problem.

Sincerely,

Andrew Kalenkiewicz

AttachmentSize
threading_command.txt318 bytes
fixaln.grishin.txt1.84 KB
threading_error_msg.txt2.08 KB
Post Situation: 
Wed, 2017-04-19 11:53
ajkal

Hi Andrew,

One potential solution to your problem is that you don't need to include the GDP/GTP/Mg in the alignment for the initial threading. Add them after you have generated your initial threaded models (renaming the atoms from HETATM to ATOM), then add a forward slash to the end of your sequence file, and then run rosettaCM.

Wed, 2017-04-19 20:01
Derek

Derek's method (adding back the ligand after threading) is a good one, but if you want to try threading with the ligand present, what you probably want to try is using the non-annotated sequence in your grishin files. Grishin alignment files take plain one-letter-code sequences, rather than annotated sequences. Remove the parts in square brackets, and I think that should work for you.  (Your fasta file will contain the full annotated sequence, but the alignment files should just be the un-annotated one letter codes.)

Thu, 2017-04-20 13:57
rmoretti