- modeller_usage - salilab.org

loop model error
by albert 18 Feb '11

18 Feb '11

Hi: I am using model loop based on initial PDB files, but there is a error: Traceback (most recent call last): File "model.py", line 41, in ? a.make() File "/usr/lib64/python2.4/site-packages/modeller/automodel/loopmodel.py", line 36, in make automodel.make(self, exit_stage) File "/usr/lib64/python2.4/site-packages/modeller/automodel/automodel.py", line 98, in make self.homcsr(exit_stage) File "/usr/lib64/python2.4/site-packages/modeller/automodel/automodel.py", line 411, in homcsr aln = self.read_alignment() File "/usr/lib64/python2.4/site-packages/modeller/automodel/automodel.py", line 401, in read_alignment aln.append(file=self.alnfile, align_codes=self.knowns+[self.sequence]) File "/usr/lib64/python2.4/site-packages/modeller/alignment.py", line 81, in append allow_alternates) IOError: openf______E> Cannot open file : No such file or directory

2 1

how seed influence the Modeller results?
by albert 18 Feb '11

18 Feb '11

Hi: I've got a question about how seed influence the Modeller results? If I start from the same seed in the same computer, does it means it should get the same results? If two job both start from random seed or different seeds, they will get different results, and those results is comparable with running model building jobs with doubled model number building? THX

2 1

plot_profiles.py errors
by hari krishna 18 Feb '11

18 Feb '11

I installed MatPlotLib 2.6v. (python2.6) Still it is showing me the error: File "plot_profiles.py", line 1, in ? import pylab Importerror: No module named pylab Please help me to over come this problem. *Thanks in advance..*

2 1

modeller and win64
by Igor Dvortsov 17 Feb '11

17 Feb '11

so the subj. could authors of modeller build x64 version of _modeller.pyd? kind regards, Igor

2 1

Protein - DNA Modelling
by Koushik A C 17 Feb '11

17 Feb '11

Hi, I wanna perform a protein - DNA modelling. But I don't know how to proceed. Can anyone help me with any tutorial or steps? --Kaushik

2 1

Re: [modeller_usage] modeller_usage Digest, Vol 10, Issue 33
by Ashish Runthala 16 Feb '11

16 Feb '11

You remove the chain identifier from both align2d.py and model_multiple.py (Or, whatever you are using) This is how i did it in one of the run. Thus, Align2d.py becomes mdl = model(env, file='Run2BestModel', model_segment=('1','292')) #No ChainID is here with model residue positions aln.append_model(mdl, align_codes='Run2BestModel', atom_files='Run2BestModel') So, Model_multiple.py will be a = automodel(env, alnfile = 'alignment.ali', knowns = ('Run2BestModel'), sequence = 'target', assess_methods=(assess.DOPE, assess.GA341,assess.normalized_dope)) a.starting_model= 1 a.ending_model = 1000 a.make() Ashish Ashish Runthala, Lecturer, Structural Biology Cell, Birla Institute of Technology & Science, Pilani Rajasthan, INDIA

1 0

mutiple CPU error
by albert 16 Feb '11

16 Feb '11

Hi I try to build model with my own constrain(beta sheet) with multiple CPU, but it said: ----------------log-------------------- raise RemoteError(obj.exc, self) modeller.parallel.communicator.RemoteError: AttributeError: 'module' object has no attribute 'MyModel' from <Slave on localhost> I ues the following script for multiple CPU: -------------script------------------------ from modeller import * from modeller.automodel import * from modeller.parallel import * j = job() j.append(local_slave()) j.append(local_slave()) j.append(local_slave()) j.append(local_slave()) env = environ() log.verbose() ---------------------------------------------- THX

2 1

initial model error
by albert 16 Feb '11

16 Feb '11

Hi: I would like to build a model staring from my own model. but it didn't work, it is said: -------log-------------- _modeller.ModellerError: rdpdb___303E> No atoms were read from the specified input PDB file, since the starting residue number and/or chain id in MODEL_SEGMENT (or the alignment file header) was not found; requested starting position: residue number " ", chain " "; atom file name: init.pdb -------------------------- In fact in the alignment file I already add something like this: -------------------align---------------------------- >P1;seq sequence:seq: : : : ::: 0.00: 0.00 aaaaaaaaaaaaaaaaaaaaaaaaaaaaaa >P1;init structureX:init.pdb:1 :A: :A:::-1.00:-1.00 aaaaaaaaaaaaaaaaaaaaaaaaaaaaaa >P1;temp structureX:temp.pdb:1 :A: :A:::-1.00:-1.00 aaaaaaaaaaaaaaaaaabbbbbbbbbbbb ----------script------------------------------------- class MyModel(automodel): def special_restraints(self, aln): rsr = self.restraints at = self.atoms # Two beta-strands: rsr.add(secondary_structure.strand(self.residue_range('146:', '150:'))) rsr.add(secondary_structure.strand(self.residue_range('156:', '161:'))) a = MyModel(env, alnfile='mult.ali', knowns=('tempA'), sequence='my', inifile = 'init.pdb') a.starting_model = 1 a.ending_model = 2 ---------------------------------------------------------------------------- the seq and init share the same sequences. The starting model init.pdb contains a china named A. could you give me some advices? THX

2 1

Modeling circular permutation
by bharat lal 16 Feb '11

16 Feb '11

Hi, I am trying to model circular permutation in my protein which is a beta barrel and I modeled it the following way-- AAAA-----BBB--------CCC------- A and B comprise the two halves of beta barrel and ---- refers to a cut in the loop region of beta barrel where the other domain has to be attached ... I generated some 20 models but none of them was satisfactory as the beta-barrel gets distorted in each of those models.. Can any body tell how can I achieve good models for modeling circular permutation ??? ---------- Bharat Ph.D. Candidate Room No. : 7202A, 2nd Floor Biomolecular Engineering Laboratory Division of Chemical Engineering and Polymer Science Pusan National University Busan -609735 South Korea Lab phone no. - +82-51-510-3680, +82-51-583-8343 E-mail : monu46010(a)yahoo.com

2 1

Differences using Xeon vs Opteron systems
by Robert Healey 15 Feb '11

15 Feb '11

Hello. I'm the systems administrator for a research lab using Modeller 9v3. I've installed the x86_64 RPM file on a pair of Xeon X5355 based systems and some Opteron 6136 based systems. The researchers are reporting getting radically different results between these two systems. I've done some testing and Opteron 252 based systems give the same results as the Opteron 6136 systems and Xeon X5570 systems give the same results as the Xeon X5355. I'm working on getting a set of results I can post, but has anyone seen anything like this before? I was being shown summaries where the final results were different in value by 1.5 Thank you. -- Bob Healey Systems Administrator Biocomputation and Bioinformatics Constellation and Molecularium healer(a)rpi.edu (518) 276-4407

2 3