Hi everybody, I am really hoping that I am not asking a dumb question but I cannot seem to get the restraints to work on my protein. To try and get to grips with this I am using a target-template sequence identity of 100%. I have set the secondary structure as is in the known structure but cannot seem to 'force' the restraints through onto the target. I have posted the script below. Any help on this would be very much appreciated! Many thanks in advance, Tory from modeller import * from modeller.automodel import * class MyModel(automodel): def special_restraints(self, aln): rsr = self.restraints # Addition of 3 disulphides: self.patch(residue_type='DISU', residues=(self.residues['43'],self.residues['62'])) #Residues 20 through 30 should be an alpha helix: rsr.add(secondary_structure.alpha(self.residue_range('20:', '37:'))) rsr.add(secondary_structure.alpha(self.residue_range('54:', '65:'))) rsr.add(secondary_structure.alpha(self.residue_range('75:', '89:'))) rsr.add(secondary_structure.alpha(self.residue_range('111:', '115:'))) rsr.add(secondary_structure.alpha(self.residue_range('119:', '128:'))) rsr.add(secondary_structure.alpha(self.residue_range('137:', '154:'))) #Two beta-strands: rsr.add(secondary_structure.strand(self.residue_range('13:', '18:'))) rsr.add(secondary_structure.strand(self.residue_range('42:', '43:'))) rsr.add(secondary_structure.strand(self.residue_range('66:', '73:'))) rsr.add(secondary_structure.strand(self.residue_range('102:', '106:'))) rsr.add(secondary_structure.strand(self.residue_range('132:', '133:'))) log.verbose() env = environ() # directories for input atom files env.io.atom_files_directory = ['.', '../../atom_files'] a = MyModel(env, alnfile='../tlr1_tir_template.ali', knowns='1fyv', sequence='hTLR1TIR', assess_methods=(assess.DOPE, assess.normalized_dope, assess.DOPEHR, assess.GA341)) a.starting_model = 1 a.ending_model = 100 # Very thorough VTFM optimization: a.library_schedule = autosched.slow a.max_var_iterations = 300 # Thorough MD optimization: a.md_level = refine.very_slow # Repeat the whole cycle 2 times and do not stop unless obj.func. > 1E6 a.repeat_optimization = 5 a.max_molpdf = 1e6 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])