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
Hello Albert,
On Tue, 06 Jul 2010 22:56:37 +0800 (CST) albert leuven@yeah.net wrote:
> 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?
<snip>
> Protein specified in ALIGN_CODES(i) was not found in the alignment file; > ALIGN_CODES( 1) = pt6
I think if you reread the error message you will see that it does not say that it cannot find the pdb file, but says that the template that you have called "pt6" is not present in the alignment file a8-pt6.ali. It appears that your alignment file does not match your script. Perhaps you labelled the sequence PT6 instead of pt6 (or something else)?
Cheers, Rob
participants (2)
-
albert
-
Robert Campbell