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