[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [modeller_usage] SequenceMismatchError when adding missing residues in LoopModel()



Dear Ben,

Okay, I think I see where I’m making the mistake then!  Currently, I am following the tutorial (), and my interpretation was that “1qg8_fill” would be the template.  It seems like that is NOT the case.  Instead, we can ignore 1qg8_fill for our current purposes.
>P1;1qg8
structureX:1qg8:   2 :A: 256 :A:undefined:undefined:-1.00:-1.00
PKVSVIMTSYNKSDYVAKSISSILSQTFSDFELFIMDDNSNEETLNVIRPFLNDNRVRFYQSDISGVKERTEKTR
YAALINQAIEMAEGEYITYATDDNIYMPDRLLKMVRELDTHPEKAVIYSASKTYHLN---DIVKETVRPAAQVTW
NAPCAIDHCSVMHRYSVLEKVKEKFGSYWDESPAFYRIGDARFFWRVNHFYPFYPLDEELDLNYIT---------
-----EFVRNLPPQRNCRELRESLKKLGMG*
>P1;1qg8_fill
sequence:::::::::
PKVSVIMTSYNKSDYVAKSISSILSQTFSDFELFIMDDNSNEETLNVIRPFLNDNRVRFYQSDISGVKERTEKTR
YAALINQAIEMAEGEYITYATDDNIYMPDRLLKMVRELDTHPEKAVIYSASKTYHLNENRDIVKETVRPAAQVTW
NAPCAIDHCSVMHRYSVLEKVKEKFGSYWDESPAFYRIGDARFFWRVNHFYPFYPLDEELDLNYITDQSIHFQLF
ELEKNEFVRNLPPQRNCRELRESLKKLGMG*
I mistook the above and created a “6z43” template myself which contained the dashes like in 1qg8 in the tutorial.  I’ll replace the dashes in that case with actual residues then.  Thank you for the help!

Many thanks,
Steven Truong


On 4/5/21 5:42 PM, Steven Truong wrote:
I think I have the template sequence alignment matching exactly to the PDB file, in the same way the second entry “1qg8_fill” has a full sequence.  However, I am indicating missing residues with “1qg8.”   Should I fill in those dashes with missing residues?  I’m likely misunderstanding how AutoModel/LoopModel works then.  When you refer to “template sequence,” is that “1qg8” or “1qg8_fill” in the align.pir file on the tutorial?

The template is the known structure, i.e. the PDB file, 1qg8 in the example.

>P1;6z43
structureM:6z43:27:A:+3600:C:::-1.00:-1.00
AYTNSFTR*-*VYYPDKVFRSSVLHSTQDLFLPFFSNVTWFH--------------NPVLPF
NDGVYFAST

Well, this simply isn't correct. The "missing" glycine here isn't missing in the 6z43 PDB file - it is there (residue 35 in chain A). What are you trying to accomplish here?

Ben Webb, Modeller Caretaker
--
" class="">             https://salilab.org/modeller/
Modeller mail list: https://salilab.org/mailman/listinfo/modeller_usage