Hello: I am running modeller job by "python model-parall.py" based on the following scrips. But it is said in the termial which seems that modeller coundn't find the template pdb file. However, my template pdb file is already in the same working directory which named pt6A.pdb. I don't why it said couldn't find it. Would you please give me some advices?
File "/soft/python2.5/lib/python2.5/site-packages/modeller/automodel/automodel.py", line 98, in make self.homcsr(exit_stage) File "/soft/python2.5/lib/python2.5/site-packages/modeller/automodel/automodel.py", line 411, in homcsr aln = self.read_alignment() File "/soft/python2.5/lib/python2.5/site-packages/modeller/automodel/automodel.py", line 401, in read_alignment aln.append(file=self.alnfile, align_codes=self.knowns+[self.sequence]) File "/soft/python2.5/lib/python2.5/site-packages/modeller/alignment.py", line 81, in append allow_alternates) _modeller.ModellerError: read_al_373E> Protein specified in ALIGN_CODES(i) was not found in the alignment file; ALIGN_CODES( 1) = pt6
---------------------------------------------------------------------------------------- # Homology modeling by the automodel class, using multiple processors from modeller import * from modeller.automodel import * # Load the automodel class from modeller.parallel import *
# Use 2 CPUs in a parallel job on this machine j = job() j.append(local_slave()) j.append(local_slave())
log.verbose() # request verbose output env = environ() # create a new MODELLER environment to build this model in
# directories for input atom files env.io.atom_files_directory = ['.', '../atom_files']
a = automodel(env, alnfile = 'a8-pt6.ali', # alignment filename knowns = 'pt6', # codes of the templates sequence = 'a8') # code of the target a.starting_model= 1 # index of the first model a.ending_model = 5 # index of the last model # (determines how many models to calculate) a.use_parallel_job(j) # Use the job for model building a.make() # do the actual homology modeling