Something was wrong with that script, I have taken model-outputs.py modified it and it worked fine.
Best regards, Andrew
-------- Пересылаемое сообщение -------- 08.11.2011, 16:41, "Andrew Voronkov" drugdesign@yandex.ru:
Does it have name molpdf? I can't find DOPE or GA341 words in the log file. I have also got an error while ranking the solutions. Like this? But I have both DOPE and GA341, which of them is this?
>> Summary of successfully produced models:
Filename molpdf ---------------------------------------- chimera.B99990001.pdb 1580.34619 chimera.B99990002.pdb 1593.52393 chimera.B99990003.pdb 1693.15637 chimera.B99990004.pdb 1743.10718 chimera.B99990005.pdb 1586.12122 chimera.B99990006.pdb 1652.49170 chimera.B99990007.pdb 1477.14709 chimera.B99990008.pdb 1681.59668 chimera.B99990009.pdb 1625.87500 chimera.B99990010.pdb 1632.77686 chimera.B99990011.pdb 1689.23584 chimera.B99990012.pdb 1557.99890 chimera.B99990013.pdb 1622.84631 chimera.B99990014.pdb 1805.83228 chimera.B99990015.pdb 1525.45911
Also I have got an error while trying to select best solutions based on DOPE score (script is below):
ok_models.sort(lambda a,b:cmp(a[key]), b[key])) KeyError: 'DOPE score'
Here is the script, which was used:
from modeller import * from modeller.automodel import * import sys
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 = ['.'] # Read in HETATM records from template PDBs env.io.hetatm = True
a = automodel(env, alnfile = '2rf53mhj.pir', # alignment filename knowns = ('3mhj', '2rf5'), # codes of the templates sequence = 'chimera') # code of the target assess_methods=(assess.DOPE, assess.GA341) # assess each loop with DOPE a.starting_model= 1 # index of the first model a.ending_model = 15 # index of the last model # (determines how many models to calculate) # Thorough MD optimization: a.md_level = refine.slow
a.make() # do the actual homology modeling
# Get a list of all successfully built models from a.outputs ok_models = [x for x in a.outputs if x['failure'] is None]
# Rank the models by DOPE score key = 'DOPE score' if sys.version_info[:2] == (2,3): # Python 2.3's sort doesn't have a 'key' argument ok_models.sort(lambda a,b: cmp(a[key], b[key])) else: ok_models.sort(key=lambda a: a[key])
# Get top model m = ok_models[0] print("Top model: %s (DOPE score %.3f)" % (m['name'], m[key]))
Best regards, Andrew
07.11.2011, 19:51, "Modeller Caretaker" modeller-care@salilab.org:
> On 11/7/11 5:21 AM, Andrew Voronkov wrote: >> Dear Ben, very simple question . I have used DOPE assessment, but in >> which file can I find the assessment results for the models? It's not >> in resulted PDB, there is only objective function.Where can I find >> those DOPE and GA341 values? > In the log file. They are not put in the PDB file. > > Ben Webb, Modeller Caretaker > > -- > modeller-care@salilab.org http://www.salilab.org/modeller/ > Modeller mail list: http://salilab.org/mailman/listinfo/modeller_usage -------- Завершение пересылаемого сообщения --------