- modeller_usage - salilab.org

error in using mutate_model
by Anna Marabotti 24 Mar '10

24 Mar '10

Hallo, I tried to model a mutant of a protein using the mutate_model.py script available in the Wiki section of MODELLER site at: http://www.salilab.org/modeller/wiki/Mutate%20model. Therefore, after copying the script on a executable file, I opened the black window to launch Modeller commands and typed, as suggested: python mutate_model.py my_modelname my_respos my_resname my_chain > my_log.log but obtained the error: "python is not recognized as internal and/or external command, or as an executable, or as a batch file" despite I have Python installed on my PC. I also tried to launch the command without "python" and typed: mutate_model.py my_modelname my_respos my_resname my_chain > my_log.log and this time the error was: File "C:\Programmi\Modeller9v7\mutate_model.py", line 4, in <module> from modeller import * File "C:\Programmi\Modeller9v7\modlib\modeller\__init__.py", line 69, in <module> import _modeller ImportError: No module named _modeller The program ends without writing the .log file. I used this script with the previous versions of Modeller, on the same PC, without any problem. It's the first time I use this script since I installed Mod9v7 on my PC, but I've already used Modeller 9v7 on the same PC without problems. Could anybody suggest me how to do? Many thanks in advance and best regards Anna ____________________________________________________ Anna Marabotti, Ph.D. Laboratory of Bioinformatics and Computational Biology Institute of Food Science, CNR Via Roma, 64 83100 Avellino (Italy) Phone: +39 0825 299651 Fax: +39 0825 781585 Email: anna.marabotti(a)isa.cnr.it Skype account: annam1972 Web page: http://bioinformatica.isa.cnr.it/anna/anna.htm "If you think you are too small to make a difference, try sleeping with a mosquito"

2 3

Fw: mutation error
by leuven 23 Mar '10

23 Mar '10

Dear: I try to make mutation model by this scrips: http://salilab.org/modeller/wiki/Mutate%20model When I typed command in terminal: py2 1-nh-active.pdb 97 VAL A>mute-native.log It is said: ATOM 1 CB THR A 1 57.352 108.506 46.727 1.00 21.29 A ^ SyntaxError: invalid syntax I don't know where is the mistake. Thank you very much.

2 1

mutation error
by leuven 23 Mar '10

23 Mar '10

Dear: I try to make mutation model by this scrips: http://salilab.org/modeller/wiki/Mutate%20model When I typed command in terminal: py2 1-nh-active.pdb 97 VAL A>mute-native.log It is said: ATOM 1 CB THR A 1 57.352 108.506 46.727 1.00 21.29 A ^

1 0

Restraints in table 6.1
by abdullah ahmed 22 Mar '10

22 Mar '10

Dear Modeller users, I have a few very long (and perhaps naive) questions regarding "Table 6.1: List of physical Restraints types" If the Hydrogen bonding restraint is applied, does it automatically incorporate all possible Hydrogen bonds in the file, or do the specific hydrogen bonding atoms have to be specified? (Also could you maybe provide a short piece of code showing how this restraint can be applied?)Are the Phi, Psi, omega, and Chi restraints used to lock the angles of the backbone and the the side chain into specific values or are they used for something else all together?What do the accessibility, absposition, dihedral_diff and SAX restraints do? Thank you in advance for your replies, Abdullah Ahmed _________________________________________________________________ Hotmail: Trusted email with powerful SPAM protection. https://signup.live.com/signup.aspx?id=60969

2 1

Lennard Jones Restraints
by abdullah ahmed 19 Mar '10

19 Mar '10

Dear Modeller users, If Lennard Jones restraints are activated during optimization using "env.edat.dynamic_sphere = True" What does this do exactly? Thank you in advance, Abdullah _________________________________________________________________ Hotmail: Free, trusted and rich email service. https://signup.live.com/signup.aspx?id=60969

2 1

set topology for autobuilder.
by Knut J Bjuland 19 Mar '10

19 Mar '10

I used thesse commando to set toplogy information when doing an alignment. BUt it did work when I wanted to build it. What kind of command should I invoke to add thesse information in the build phase of modeller. env.libs.topology.read(file='$(LIB)/top_allh.lib') env.libs.parameters.read(file='$(LIB)/par.lib') _________________________________________________________________ Hotmail: Trusted email with powerful SPAM protection. https://signup.live.com/signup.aspx?id=60969

2 1

If I disable the default N and C-term patched, must I manually apply them?
by wayj86 wayj86 19 Mar '10

19 Mar '10

Hi all, I hope to model a protein which is only a part of a big one (i.e. not up to its real N or C-terminal). So I turn off the default N and C patches with "patch_default=false". In this case, must I manually apply neutral patches, like ACE and CT3 to my model? I mean, can Modeller work wothout any N and C-term patches? I asked this because my next step is docking, which did not require any preexisting terminal patches. Thank you. Best wishes, Stanley -- Be the change you want to see in the world.

2 1

log files and processing multiple files
by abdullah ahmed 18 Mar '10

18 Mar '10

Dear Modeller users, I would like to ask for help on two problems: 1) I do not get log files after running optimization. I am certain that it is not a matter of them taking really long to be made because I have been using Modeller for 3 months now and I have not gotten a log file. I do however, get "pdb.ini" "pdb.rsr" and "pdb.D000001" (contains info on maximal shift, energy changes and gradients). 2) I have a series of molecules that I would like to apply the same constraints to, and then optimize. Is there a script I can use to optimize all the files in one go? Thanks in advance for the help, Abdullah _________________________________________________________________ Your E-mail and More On-the-Go. Get Windows Live Hotmail Free. https://signup.live.com/signup.aspx?id=60969

2 1

modeller and virtualenv
by Will Dampier 17 Mar '10

17 Mar '10

I've been trying to modeller to play nice with the virtualenv package management system. Has anyone had any success getting this to work? thanks Will

2 1

Model refinement
by Ehab AlMoubarak 16 Mar '10

16 Mar '10

Dear All, I am modelling a potassium channel protein using a template with fairly low sequence similarity to my target, therefore, I am trying to put several restraints. I am struggling to deal with a long loop supposed to be extracellular: so its location to the rest of the protein is roughly known. As I am refining the model using MD, the loop location completely became wrong, I suppose that the absence of alignment in the refinement process is the reason (which covers only small area of the loop anyway). I don’t know if there is anyway to exclude the loop while refining the rest of the model, or to restraint the loop in an approximate location of the structure (like putting distance restraint) while refining it. I must say that the gap in alignment corresponds to the loop is well over 15 amino acids. Any thoughts are highly appreciated. My script is here: from modeller import * from modeller.automodel import * # Load the automodel class log.verbose() # Override the 'special_restraints' and 'user_after_single_model' methods: class MyModel(automodel): def special_restraints(self, aln): # Constrain the A and B chains to be identical (but only restrain # the C-alpha atoms, to reduce the number of interatomic distances # that need to be calculated): s1 = selection(self.chains['A'], self.residue_range('7:A', '30:A'),self.residue_range('107:A', '130:A')).only_atom_types('CA') #s2 = selection(self.chains['A'], self.residue_range('107:A', '130:A')).only_atom_types('CA') s2 = selection(self.chains['A'], self.residue_range('157:A', '180:A'),self.residue_range('217:A', '240:A')).only_atom_types('CA') #s4 = selection(self.chains['A'], self.residue_range('217:A', '240:A')).only_atom_types('CA') s5 = selection(self.chains['A']).only_atom_types('CA') s6 = selection(self.chains['B']).only_atom_types('CA') s7 = selection(self.chains['A'], self.residue_range('80:A', '135:A')).only_atom_types('CA') s8 = selection(self.chains['A'], self.residue_range('185:A', '250:A')).only_atom_types('CA') self.restraints.add(secondary_structure.alpha(self.residue_range('7:A', '30:A'))) #self.restraints.add(secondary_structure.alpha(self.residue_range('107:A', '130:A'))) self.restraints.add(secondary_structure.alpha(self.residue_range('157:A', '180:A'))) #self.restraints.add(secondary_structure.alpha(self.residue_range('217:A', '240:A'))) #self.restraints.add(secondary_structure.alpha(self.residue_range('73:A', '92:A'))) #self.restraints.symmetry.append(symmetry(s1, s2, 1.0)) self.restraints.symmetry.append(symmetry(s7, s8, 1.0)) #self.restraints.symmetry.append(symmetry(s2, s3, 1.0)) #self.restraints.symmetry.append(symmetry(s3, s4, 1.0)) self.restraints.symmetry.append(symmetry(s5, s6, 1.0)) def user_after_single_model(self): # Report on symmetry violations greater than 1A after building # each model: self.restraints.symmetry.report(1.0) env = environ() # directories for input atom files env.schedule_scale = physical.values(default=1.0, soft_sphere=0.7) env.io.atom_files_directory = ['.', '../atom_files'] # Be sure to use 'MyModel' rather than 'automodel' here! a = MyModel(env, alnfile = 'task3Modified8.ali' , # alignment filename knowns = '1orq_9', # codes of the templates sequence = 'TASK3h27', # code of the target assess_methods=assess.DOPE) a.starting_model= 1 # index of the first model a.ending_model = 15 # index of the last model # Very thorough VTFM optimization: a.library_schedule = autosched.slow a.max_var_iterations = 300 # Thorough MD optimization: a.md_level = refine.slow # Repeat the whole cycle 2 times and do not stop unless obj.func. > 1E6 a.repeat_optimization = 2 a.max_molpdf = 1ttt a.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: %s (DOPE score %.3f)" % (m['name'], m[key]) My best regards Ehab, Research student

1 0