Dear Modellers,

I am tring to set up my script to build a multimeric model with loop refinement. The script looks like:

log.verbose()

class dope_loopmodel(dope_loopmodel): def special_patches(self, aln): # Rename both chains and renumber the residues in each self.rename_segments(segment_ids=['A', 'B'],

renumber_residues=[1, 1]) env = environ() env.io.hetatm = True a = dope_loopmodel(env, alnfile='TvLDH-1bdm.ali', knowns='1bdm', sequence='TvLDH', assess_methods=(assess.DOPE, assess.GA341)) a.starting_model = 1 a.ending_model = 2

# 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 = 2 a.max_molpdf = 1e6

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()

automodel.final_malign3d = True

It seems that the script works but I would like to check if all optimization steps have been done. How could I check it? Please, if you note errors in the script, give me suggestions

Moreover, I built a model with an HETATM residue. If I want to treat it as flexible residue: how could I do that?

Thank you in advance.

Bests

Simone