October 2007 - modeller_usage

ImportError: No module named socket (or math)
by Greipel.Joachim＠MH-Hannover.DE 21 Oct '07

21 Oct '07

Hi all, I just installed Modeller 9v2 on my linux box and ran the examples. Occasionally I got the Error Messages ImportError: No module named socket resp. ImportError: No module named math. Does anybody know how to cope with this? Regards, Joachim (Here is an example of the output) mod9v2 center_of_mass.py Traceback (most recent call last): File "center_of_mass.py", line 3, in ? import math ImportError: No module named math mod9v2 selection.py mod9v2 parallel-msg.py Traceback (most recent call last): File "parallel-msg.py", line 2, in ? from modeller.parallel import * File "/network/apps/share/modeller-9v2/modlib/modeller/parallel/__init__.py", line 1, in ? from job import job File "/network/apps/share/modeller-9v2/modlib/modeller/parallel/job.py", line 1, in ? import socket, random, os, select ImportError: No module named socket mod9v2 parallel-task.py Traceback (most recent call last): File "parallel-task.py", line 2, in ? from modeller.parallel import * File "/network/apps/share/modeller-9v2/modlib/modeller/parallel/__init__.py", line 1, in ? from job import job File "/network/apps/share/modeller-9v2/modlib/modeller/parallel/job.py", line 1, in ? import socket, random, os, select ImportError: No module named socket

2 1

regarding N terminus region
by simran jeet 17 Oct '07

17 Oct '07

I have a protein sequence whose N terminus region does not matches with the template. Is it neccessary that N and C terminus regions of a protein should be fixed for homology modeling? The n terminus region contains 12 residues.  There are long gaps also within the alignment. Can these regions be modelled by loop modelling? I am also attaching the alignment file.

2 1

Warning messages
by csreddy＠ncbs.res.in 17 Oct '07

17 Oct '07

Hi All, I recently installed modller9v2 on the linux (UBUNTU) machine. I could modeller successfully, but I am getting few warning messages at the beginning of the execution. Below I am pasting the warning messages. Could not find platform independent libraries <prefix> Could not find platform dependent libraries <exec_prefix> Consider setting $PYTHONHOME to <prefix>[:<exec_prefix>] 'import site' failed; use -v for traceback I am thankful if some one can help me to correct it. Thnks in advance Regards Chandu -------------------------------------------------- C.CHANDRA SEKHAR REDDY(Research Scholar) Dr.SOWDHAMINI'S LAB NATIONAL CENTRE FOR BIOLOGICAL SCIENCES TATA INSTITUTE OF FUNDAMENTAL RESEARCH UAS-GKVK CAMPUS BELLARY ROAD BANGALORE-560065 INDIA PHONE NO.#+91-(0)80- 2366 6251 E.MAIL# csreddy(a)ncbs.res.in chandubio(a)yahoo.co.in URL# http://www.ncbs.res.in/ --------------------------------------------------

2 1

regarding N terminus region
by simran jeet 17 Oct '07

17 Oct '07

I have a protein sequence whose N terminus region does not matches with the template. Is it neccessary that N and C terminus regions of a protein should be fixed for homology modeling? The n terminus region contains 12 residues.  There are long gaps also within the alignment. Can these regions be modelled by loop modelling? I am also attaching the alignment file.

1 0

Re: [modeller_usage] RE : Re: Class mymodel: How to randomize coordinates
by Modeller Caretaker 11 Oct '07

11 Oct '07

Baptiste Legrand wrote: > Thanks for your answer, I didn't know that the default deviation is 4.0. > Is it interessant to set it to 8.0 or more or the models will be too far > from the template? Is there a theorical (estimation) limit to explore > the conformationnal space without deviate totally from the template? The deviation parameter controls the spread of initial models, not the final optimized models - if the optimization is functioning correctly, the optimized models should look like the template anyway. (However, optimization is not a 'solved' problem, which is why the initial model is the sequence threaded onto the template structure rather than, say, the extended chain or randomly placed points.) So in practice changing deviation will not have a large effect on your model generation (unless you set it to 0, in which case no randomization is done). Ben Webb, Modeller Caretaker -- modeller-care(a)salilab.org http://www.salilab.org/modeller/ Modeller mail list: http://salilab.org/mailman/listinfo/modeller_usage

1 0

Class mymodel: How to randomize coordinates
by Baptiste Legrand 10 Oct '07

10 Oct '07

Hi all, I try to make models from a template, with constraints (Helix) and a randomization of the coordinates for each models. I set deviation to 4.0 but MODELLER dont use it. A part of my script: # Add Helix class mymodel(automodel): def special_restraints(self, aln): rsr = self.restraints at = self.atoms rsr.add(secondary_structure.alpha(self.residue_range('2:', '13:'))) # Build models, and assess with both DOPE and GA341 a = mymodel(env, alnfile = 'f125-1s35.ali', # alignment filename (use default: 'a = automodel', 'loopmodel' for loop, 'allhmodel' to add protons) knowns = '1j18', # codes of the templates (for multiple template, just add codes: knowns = ('1s35', '1quu')) sequence = 'f3s5', # code of the target deviation = 4.0, assess_methods =(assess.DOPE, assess.GA341)) # assessment scripts #inifile = 'my-initial.pdb') # use my initial structure What should I do to randomize in this case? Thanks in advance, Regards Baptiste --------------------------------- Ne gardez plus qu'une seule adresse mail ! Copiez vos mails vers Yahoo! Mail

2 1

Re: [modeller_usage] exploding GTP
by Modeller Caretaker 09 Oct '07

09 Oct '07

Jonathan Winger wrote: > I'm using Modeller 9v2. My protein is a homodimer, and my template is > also a homodimer. There is one molecule of AMPCPP and 2 metals bound to > each chain of the template. I adjusted the PDB file (2BW7) so that the > metals and nucleotides followed directly after the appropriate chain > TER's. I should add that the modeling runs fine in the absence of > nucleotide. You could try using the BLK residue type (.) instead of GTP (y) in your model sequence - it is entirely possible that the CHARMM forcefield parameters for GTP are not doing a great job. As for your metal, your model contains MG but your template contains CA, so since you have no equivalent template atom, the MG gets no homology-derived restraints on it and it can wander (Modeller is not smart enough to recognize automatically that one metal could be considered equivalent to another here.) To fix that, you could either create manual distance restraints on your metal ion to hold it in place, or edit your template PDB and change the metal ion there to MG. Ben Webb, Modeller Caretaker -- modeller-care(a)salilab.org http://www.salilab.org/modeller/ Modeller mail list: http://salilab.org/mailman/listinfo/modeller_usage

1 0

exploding GTP
by Jonathan Winger 09 Oct '07

09 Oct '07

When modeling GTP and metal ions into my protein, using a template with 2 metal ions and an ATP analog, the GTP atoms end up disconnected from each other and scattered around the active site, and one of the two metals is launched out of the protein. What can I do to fix this? Thanks much, Jon

2 1

qsub and cis restraints
by Gordon Wells 08 Oct '07

08 Oct '07

Do classes inheriting from automodel get propagated when using the sge_qsub parallel job? I can get a parallel run to work with automodel, but not with a new class including a cis-peptide restraint defined in special_restraints (which I can get to run in single proc mode). I get the following error: # Traceback (most recent call last): File "psamdc.py", line 77, in ? a.make() File "/home/namd/lib/modeller9v2/modlib/modeller/automodel/automodel.py", line 110, in make self.multiple_models(atmsel) File "/home/namd/lib/modeller9v2/modlib/modeller/automodel/automodel.py", line 214, in multiple_models self.parallel_multiple_models(atmsel) File "/home/namd/lib/modeller9v2/modlib/modeller/automodel/automodel.py", line 234, in parallel_multiple_models self.outputs.extend(job.run_all_tasks()) File "/home/namd/lib/modeller9v2/modlib/modeller/parallel/job.py", line 110, in run_all_tasks self.finish_next_task() File "/home/namd/lib/modeller9v2/modlib/modeller/parallel/job.py", line 127, in finish_next_task self.process_event(obj, s) File "/home/namd/lib/modeller9v2/modlib/modeller/parallel/job.py", line 141, in process_event task = obj.task_results() File "/home/namd/lib/modeller9v2/modlib/modeller/parallel/slave.py", line 54, in task_results r = self.get_data(allow_heartbeat=True) File "/home/namd/lib/modeller9v2/modlib/modeller/parallel/communicator.py", line 78, in get_data (cmdtype, obj) = self._recv() File "/home/namd/lib/modeller9v2/modlib/modeller/parallel/communicator.py", line 116, in _recv raise RemoteError(obj, self) modeller.parallel.communicator.RemoteError: AttributeError: 'module' object has no attribute 'cismodel' from <SGE qsub slave, ID 107.1> Is there a convenient way around this? Thanks in advance Regards -- I think perhaps the most important problem is that we are trying to understand the fundamental workings of the universe via a language devised for telling one another when the best fruit is. -- Terry Pratchett, alt.fan.pratchett -- Gordon Wells Bioinformatics and Computational Biology Unit Department of Biochemistry University of Pretoria

2 1

Re: [modeller_usage] Mutations
by laura.buttazzo＠email.it 04 Oct '07

04 Oct '07

Hi, I tried the example in http://salilab.org/modeller/wiki/Mutate_model, but my model hasn't chains. It has been calculated with Modeller and the residue that I would like to change is the 27: ATOM 203 N ASP 27 16.536 9.555 15.986 1.00 91.26 1SG 204 ATOM 204 CA ASP 27 16.891 10.947 16.025 1.00 91.26 1SG 205 ATOM 205 CB ASP 27 18.072 11.218 16.974 1.00 91.26 1SG 206 ATOM 206 CG ASP 27 18.307 12.718 17.052 1.00 91.26 1SG 207 ATOM 207 OD1 ASP 27 17.470 13.485 16.504 1.00 91.26 1SG 208 ATOM 208 OD2 ASP 27 19.335 13.117 17.662 1.00 91.26 1SG 209 ATOM 209 C ASP 27 17.319 11.354 14.653 1.00 91.26 1SG 210 ATOM 210 O ASP 27 17.069 12.473 14.212 1.00 91.26 1SG 211 However I tried to change the following parts in the script: >#first argument >modelname, respos, restyp, chain, = sys.argv[1:] I put: modelname, respos, restyp, = sys.argv[1:] >s = selection(mdl1.chains[chain].residues[respos]) I put s = selection(mdl1.residues[respos]) >s = selection(mdl1.chains[chain].residues[respos]) s = selection(mdl1.residues[respos]) But I have obtained a model very similar as that in input. Should I try to change some parameter in the script? Thank you in advance. Laura --------- Original Message -------- Da: Modeller Caretaker <modeller-care(a)salilab.org> To: Cc: modeller_usage(a)salilab.org Oggetto: Re: [modeller_usage] Mutations Data: 25/09/07 19:46 > > > > laura.buttazzo(a)email.it wrote: > > I'm a new user of Modeller and I have a question about the model computing. > > If I take one/or more models with a given alignment and I change just one or > > more amino acids maintaning the same alignment, should I obtain a different > > model? > > In other words, Can Modeller be used for calculating mutations? > > Certainly - there has been some work done for SNPs using Modeller; see > http://salilab.org/modeller/wiki/Mutate_model for an example script. > > Ben Webb, Modeller Caretaker > -- > modeller-care(a)salilab.org http://www.salilab.org/modeller/ > Modeller mail list: http://salilab.org/mailman/listinfo/modeller_usage > > > -- Email.it, the professional e-mail, gratis per te: http://www.email.it/f Sponsor: In REGALO 'Meravigliosa Creatura' la super hit di GIANNA NANNINI Clicca qui: http://adv.email.it/cgi-bin/foclick.cgi?mid=6615&d=20071004

2 1