[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[modeller_usage] optimization protocol and HETATM flexibility



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