| Dear All, I have a problem with restraints after loop refeinement, in particular the distances between some atoms do not respect the restraint values set up. These values are respected after model building (without a slow MD) but not after loop refinements. I didn't use any refinement of imput model. I show the used scripts: This is for automodel: from modeller import * from modeller.automodel import * # Load the automodel class log.verbose() env = environ() env.io.atom_files_directory = ['.', '../atom_files'] class MyModel(automodel): def special_restraints(self, aln): rsr = self.restraints at = self.atoms rsr.add(forms.gaussian(group=physical.xy_distance, feature=features.distance(at['NZ:114'], at['OD1:288']), mean=2.6, stdev=0.1)) rsr.add(forms.gaussian(group=physical.xy_distance, feature=features.distance(at['NE:136'], at['OD1:260']), mean=2.6, stdev=0.3)) a = MyModel (env, alnfile='allinea.pir', knowns='AB, BC', sequence='CD') a.starting_model = 1 a.ending_model = 1 a.md_level = refine.very_fast a.make() This is for loop refinement: log.verbose() env = environ() a = loopmodel(env, alnfile = 'allinea.pir', # alignment filename knowns = ''AB, BC', # codes of the templates sequence = 'CD') # code of the target) a.starting_model = 1 # index of the first model a.ending_model = 1 # index of the last model # (determines how many models to calculate)a. md_level = None # No refinement of model a.loop.starting_model = 1 # First loop model a.loop.ending_model = 1 # Last loop model a.loop.md_level = refine.fast # Loop model refinement level a.make() Thank you very much in advance Mario |