August 2014 - modeller_usage

malign3d error with different structures but the same sequence
by Felipe Villanelo 03 Sep '14

03 Sep '14

Hello Modellers users: I want to analyze a set of PDB models, all with same sequence (84 Aa peptide) but different conformations (obteined by different methods). So I used the compare.py script from tutorial to do a first classification from structure but I get the error: _modeller.ModellerError: fit2xyz_296E> Number of equivalent positions < 3: 0 I search this error, and most of the time is related to too different proteins that do not align well in sequence. But in my case all the protein are the same sequence!. Even I add a write step between malign() and malign3d() in the script so I see a perfect sequence alignment. Some idea? could be the reason that some structures are too different? Thanks! Felipe Villanelo Lizana Bioquímico Laboratorio de Biología Estructural y Molecular Universidad de Chile "La curiosidad no mató al gato, sino que la falta de ella"

2 1

Models based on the multiple templates
by James Starlight 03 Sep '14

03 Sep '14

Dear Modeler users! Could you provide me with the ideas bases on which the modeler build the final models having several templates as well as its multiple-structural alignment as the inputs? In my case I'd like to build model of the GPCR having 5 templates of other receptors (with known X-ray structures) having low sequence identity against my target. What will be better i) to include all structures as the templates or ii) choose only one template with the highest sequence identify to my target? In the first case we will obtain trivial result: model will looks like the template but in the second case ( model built based on the multiple templates) - what type of averaging between templates will be used by modeler to produce final model? Please provide me with some ideas and reading about this topic. Thanks for suggestions, James

2 1

secondary structure restrain error
by Albert 18 Aug '14

18 Aug '14

Hello: I am trying to use the following script for restrain my secondary structure: from modeller import * from modeller.automodel import * # Load the automodel class log.verbose() env = environ() # directories for input atom files env.io.atom_files_directory = [’.’, ’../atom_files’] class MyModel(automodel): def special_restraints(self, aln): rsr = self.restraints at = self.atoms # Add some restraints from a file: # rsr.append(file=’my_rsrs1.rsr’) # Residues 20 through 30 should be an alpha helix: rsr.add(secondary_structure.alpha(self.residue_range(’20:’, ’30:’))) # Two beta-strands: rsr.add(secondary_structure.strand(self.residue_range(’1:’, ’6:’))) rsr.add(secondary_structure.strand(self.residue_range(’9:’, ’14:’))) a = MyModel(env, alnfile='CBM-mult.ali', knowns=('2z73A','3uonA'), sequence='CBM', assess_methods=(assess.DOPE, assess.GA341)) a.starting_model= 1 a.ending_model = 1 a.make() But it always failed with messages: File "model_mult.py", line 36 SyntaxError: Non-ASCII character '\xe2' in file model_mult.py on line 36, but no encoding declared; see http://www.python.org/peps/pep-0263.html for details 02-modeller/work> python test.py File "test.py", line 7 SyntaxError: Non-ASCII character '\xe2' in file test.py on line 7, but no encoding declared; see http://www.python.org/peps/pep-0263.html for details Does anybody have any idea where is the problem? THX ALbert

2 1

distance restraint
by flavio seixas 13 Aug '14

13 Aug '14

Hi Ben, I try to make a distance restraint between two atoms (protein and cofator) in a homodimer protein. However, I face another error: [flavios@localhost]$ mod9.13 model-mult-simmetry-bond.py 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 Traceback (most recent call last): File "model-mult-simmetry-bond.py", line 40, in ? a.make() File "/usr/lib/modeller9.13/modlib/modeller/automodel/automodel.py", line 108, in make self.multiple_models(atmsel) File "/usr/lib/modeller9.13/modlib/modeller/automodel/automodel.py", line 212, in multiple_models self.outputs.append(self.single_model(atmsel, num)) File "/usr/lib/modeller9.13/modlib/modeller/automodel/automodel.py", line 308, in single_model self.model_analysis(atmsel, filename, out, num) File "/usr/lib/modeller9.13/modlib/modeller/automodel/automodel.py", line 350, in model_analysis self.user_after_single_model() File "model-mult-simmetry-bond.py", line 34, in user_after_single_model rsr.add(forms.gaussian(group=physical.xy_distance, feature=features.distance(at['NZ:38'], at['C4:382']), mean=1.5, stdev=0.1)) NameError: global name 'rsr' is not defined I am certain that it must be an script error, but I do not know how to fix (I am not familiar with python script). My script is: # Homology modeling with single template from modeller import * # Load standard Modeller classes from modeller.automodel import * # Load the automodel class 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 = ['.', '../atom_files'] # Read in HETATM records from template PDBs env.io.hetatm = True env.io.water = True class MyModel(automodel): def special_restraints(self, aln): rsr = self.restraints at = self.atoms # Constrain the A, B, C and D 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']).only_atom_types('CA') s2 = selection(self.chains['B']).only_atom_types('CA') self.restraints.symmetry.append(symmetry(s1, s2, 1)) def user_after_single_model(self): # Report on symmetry violations greater than 1A after building # each model: self.restraints.symmetry.report(1.0) #Restrain between LYS and PLP (st. dev.=0.1) #Use a harmonic potential and X-Y distance group. rsr.add(forms.gaussian(group=physical.xy_distance, feature=features.distance(at['NZ:38'], at['C4:382']), mean=1.5, stdev=0.1)) rsr.add(forms.gaussian(group=physical.xy_distance, feature=features.distance(at['NZ:424'], at['C4:768']), mean=1.5, stdev=0.1)) a=MyModel(env, alnfile='AlRacem-mult.ali', knowns=('4A3Q', '1BD0'), sequence='AlRacem', assess_methods=(assess.DOPE, assess.GA341)) a.starting_model = 1 a.ending_model = 10 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 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]) ------------------------------------- Flavio Augusto Vicente Seixas Laboratory of Structural Biochemistry Department of Biochemistry Universidade Estadual de Maringá, PR, Brazil http://www.uem.br

2 1

User restriction
by flavio seixas 11 Aug '14

11 Aug '14

Hi, I was trying to make a restriction to generate a covalent bond between NZ atom from Lisine with C4 atom of PLP cofator. I search Manual and modeler-List to find some clues about how to do, but USER restriction appears to be not supported by 9.13 version: [flavios@localhost AlaninaRacemase-Saureus-T9]$ mod9.13 model-mult-simmetry-restr.py & 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 Traceback (most recent call last): File "model-mult-simmetry-restr.py", line 35, in ? a.make() File "/usr/lib/modeller9.13/modlib/modeller/automodel/automodel.py", line 103, in make self.rd_restraints() File "/usr/lib/modeller9.13/modlib/modeller/automodel/automodel.py", line 191, in rd_restraints self.restraints.append(file=self.csrfile) File "/usr/lib/modeller9.13/modlib/modeller/restraints.py", line 77, in append return _modeller.mod_restraints_read(self.__mdl.modpt, fh.file_pointer) _modeller.ModellerError: rdcsrs__703E> USER restraints format is no longer supported. My script contains model-mult-simmetry-restr.py contains: # Homology modeling from modeller import * # Load standard Modeller classes from modeller.automodel import * # Load the automodel class 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 = ['.', '../atom_files'] # Read in HETATM records from template PDBs env.io.hetatm = True env.io.water = True class MyModel(automodel): def special_restraints(self, aln): rsr = self.restraints at = self.atoms # Constrain the A, 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']).only_atom_types('CA') s2 = selection(self.chains['B']).only_atom_types('CA') self.restraints.symmetry.append(symmetry(s1, s2, 1)) def user_after_single_model(self): # Report on symmetry violations greater than 1A after building # each model: self.restraints.symmetry.report(1.0) a=MyModel(env, alnfile='AlRacem-mult.ali', knowns=('4A3Q', '1BD0'), sequence='AlRacem', csrfile='my.rsr', assess_methods=(assess.DOPE, assess.GA341)) a.starting_model = 1 a.ending_model = 1 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 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]) my.rsr file contains: MODELLER5 VERSION: USER FORMAT R 3 1 1 1 2 2 0 1.5000 0.2000 NZ:38:A C4:382:A R 3 1 1 1 2 2 0 1.5000 0.2000 NZ:424:B C4:768:B I appreciate any help. Best regards, Flavio ------------------------------------- Flavio Augusto Vicente Seixas Laboratory of Structural Biochemistry Department of Biochemistry Universidade Estadual de Maringá, PR, Brazil http://www.uem.br

2 1