- modeller_usage - salilab.org

Loop modelling
by Bishwajit Das 23 Nov '11

23 Nov '11

Dear modellers when i use automodel class in modelling a protein then is it necessary to run another loop model program?? please help me about this . -- * BISHWAJIT DAS* *+91-9433237920* KOLKATA

2 1

loop modelling
by Bishwajit Das 23 Nov '11

23 Nov '11

Dear modellers when i use automodel class in modelling a protein then is it necessary to run another loop model program?? please help me about this . -- *BISHWAJIT DAS* *+91-9433237920* KOLKATA

1 0

Loop Model ranking
by flavio seixas 23 Nov '11

23 Nov '11

Hi Ben First of all I would like to apologise for a so basic question. However, I am not a pyton expert and do not have another one to ask here. I run the loop.py scritp to generate 200 models. However, my attempt to rank the models by DOPE score fails. I am using the script below: ############################################ from modeller import * from modeller.automodel import * log.verbose() env = environ() env.io.atom_files_directory = ['.', '../atom_files'] # Read in HETATM records from template PDBs env.io.hetatm = True # Create a new class based on 'loopmodel' so that we can redefine # select_loop_atoms class MyLoop(loopmodel): # This routine picks the residues to be refined by loop modeling def select_loop_atoms(self): # Two residue ranges (both will be refined simultaneously) return selection(self.residue_range('23:A', '25:A'), self.residue_range('73:A', '75:A')) m = MyLoop(env, inimodel='Modeloutput23.pdb', sequence='Target1', loop_assess_methods=(assess.DOPE, assess.GA341)) m.loop.starting_model= 1 # index of the first loop model m.loop.ending_model = 200 # index of the last loop model m.loop.md_level = refine.very_fast # loop refinement method m.make() # Get a list of all successfully built models from a.outputs ok_models = filter(lambda x: x['failure'] is None, a.outputs) # Rank the models by DOPE score key = 'DOPE score' ok_models.sort(lambda a,b: cmp(a[key], b[key])) # Get top model m = ok_models[0] print "Top model 1: %s (DOPE score %.3f)" % (m['name'], m[key]) m = ok_models[1] print "Top model 2: %s (DOPE score %.3f)" % (m['name'], m[key]) m = ok_models[2] print "Top model 3: %s (DOPE score %.3f)" % (m['name'], m[key]) m = ok_models[3] print "Top model 4: %s (DOPE score %.3f)" % (m['name'], m[key]) m = ok_models[4] print "Top model 5: %s (DOPE score %.3f)" % (m['name'], m[key]) ############################## I realise that "loop assess methods" could not be the same sintax for "model assess methods", but I do not know how to correct it in the script Could you give me any help? Thanks in advance, Flavio

2 3

error in align codes
by si yang 17 Nov '11

17 Nov '11

Hi, i am Yang here.. I try to build a model by using basic tutorial, that begin with only protein sequence, to look for potential templates. However, i have encountered a problem during build_profile.py loading, it generate problem as below: Protein specified in ALIGN_CODES<i>, was not found in the alignment file;ALIGN_CODES< 1>= ALL The alignment file i saved in plain text format, with suffix .ali (bleg1_2439.ali), is that affected the reading by modeller? I don't specified the templates are used, so just follow the python file from basic example, the only thing i change, is only the align file that attached above.. Your feedback is important and appreciated.. Thank you

2 1

Protein-DNA loop model
by Rubben Torella 17 Nov '11

17 Nov '11

Hi all, I used Modeller in the past to create small loops from scratch. Now I have a protein-DNA complex, and I have to build and model 4 loops of my protein. My first problem is that MODELLER does not support DNA, so I have to strip the DNA from the PDB and use the protein alone. Considering the problem that I have to "paste" the DNA to the modelled structure, I use the options *a.md_level= None* to keep the rest of the protein fixed and the option * a.loop.md_level* to model the loops. Now, these two commands work on three of the 4 loops, but I have another problem now: I have a very long loop to build (35 aa) and using these two commands, MODELLER create a very strange hairpin loop to connect these two positions. I know which will be the correct structure of this loop because I have a similar protein that present this loop, and the modelling is wrong. Finally, if I don't usage a.md_level=None, MODELLER builds a very good loop, but all the protein is in a complete different position with respect of the starting one. Any suggestion? Best regards Rubben -- PhD Student, Prof.Glen Group University of Cambridge, Unilever Center, Lensfield Road CB2 1EW Cambridge, Cambridgeshire, United Kingdom University Phone : +44 (0)1223 336432 Cell Phone: +44 (0)774 3358973

2 1

Error in align codes
by si yang 17 Nov '11

17 Nov '11

Hi, i am Yang here.. I am trying to build my model by using basic tutorial, however during build_profile, i have encountered a problem,Modeller error: read_al_373E> Protein specified in ALIGN_CODES<i> was not found in the alignment files; ALIGN_CODES< 1>= ALL, because i am looking for potential templates to my target. I have saved my alignment file (bleg1_2439.ali) in plain text format with suffix .ali, is that affect the reading? Your feedback is important and appreciated.. Thank you

1 0

Problem in running script file
by Anushree Tripathi 16 Nov '11

16 Nov '11

I have installed Modeller and python2.3 version.After making all three files (i.e.,atom,alignment and script files in a folder placed on desktop) and changing directory, when I run command i.e., mod9v8 model-default.py.<http://model-default.py.it/>It is showing result i.e.,'mod9v8' is not recognised as an internal and external command,operable file or batch file.Please guide me to troubleshoot this problem as soon as possible.

2 1

error with loops in transmembrane receptor model
by Vasillis Myrianthopoulos 16 Nov '11

16 Nov '11

Hi everybody, I am trying to build a model of a transmembrane receptor and I get serious problems and obvious steric errors on the loops which connect the TM helices. The alignment is correct on a biological and phylogenetic basis however there are large insertions on the loops. As a result, there are two main errors that appear in the model. First, the loops are modelled as packed between the TM helices and not as free moving segments that can float in the solvent environment and second, in some cases the loops form kind of knots which are obviously wrong. I guess that there must be some kind of parameter within MODELLER that directs the construction towards the maximum intramolecular contact, this is the only way I can explain this wrong output. Does anybody have an idea on the matter? many thanks, Vassilis

2 1

error with loops in transmembrane receptor model
by Vasillis Myrianthopoulos 16 Nov '11

16 Nov '11

Hi everybody, I am trying to build a model of a transmembrane receptor and I get serious problems and obvious steric errors on the loops which connect the TM helices. The alignment is correct on a biological and phylogenetic basis however there are large insertions on the loops. As a result, there are two main errors that appear in the model. First, the loops are modelled as packed between the TM helices and not as free moving segments that can float in the solvent environment and second, in some cases the loops form kind of knots which are obviously wrong. I guess that there must be some kind of parameter within MODELLER that directs the construction towards the maximum intramolecular contact, this is the only way I can explain this wrong output. Does anybody have an idea on the matter? many thanks, Vassilis

1 0

Lucky16: Ordered search across multiple models?
by R K Belew 10 Nov '11

10 Nov '11

hi Ben, i'm finding an odd pattern across what i think of as random trials, and i wonder why. Suppose i run N independently seeded runs from the same alignment and generate K models each time; ie, i have a distribution of N * K models. Now consider a small fraction EPSILON of these with good GA341 scores. i would expect the fraction of models which happened to be model "target.B ..._k" ( ie, the model that gets generated the k-th time by Modeller) to be uniform over choice of k, wouldn't you? instead, i find that for some values of k, ~50% have very good GA341 scores (across various random seed instances), while some have none. here's a sample distribution with K=20: > ModNum Freq > 1 0.01 > 3 0.03 > 4 0.02 > 5 0.01 > 6 0.09 > 7 0.01 > 8 0.01 > 9 0.05 > 11 0.06 > 12 0.08 > 15 0.01 > 16 0.50 > 19 0.05 > 20 0.01 what makes model#16 so consistently good?1 rik

2 1