Take the 2-minute tour ×
Biology Stack Exchange is a question and answer site for biology researchers, academics, and students. It's 100% free, no registration required.

I am doing a basic modelling following the instructions here. But different target / template sequences.


Do I need to convert the FASTA to PIR with a script something like this


or in the 3rd step in basic tutorial

from modeller import *
env = environ()
aln = alignment(env)
mdl = model(env, file='1bdm', model_segment=('FIRST:A','LAST:A'))
aln.append_model(mdl, align_codes='1bdmA', atom_files='1bdm.pdb')
aln.append(file='TvLDH.ali', align_codes='TvLDH')
aln.write(file='TvLDH-1bdmA.ali', alignment_format='PIR')
aln.write(file='TvLDH-1bdmA.pap', alignment_format='PAP')

Changing alignment_format='PIR' to alignment_format='FASTA' is sufficient enough?

share|improve this question
From the second link you posted, it looks like simply changing the line of code will not work -- in the caretaker's message, it sounds like it needs to have the empty header information. However, why don't you try it both ways and report back if both give the same result? –  dd3 May 13 '13 at 22:02
add comment

Know someone who can answer? Share a link to this question via email, Google+, Twitter, or Facebook.

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.