#redefine the special_patches routine to include the additional disulfides
#This routine is empty by default
class mymodel(automodel):
def special_patches(self, aln):
# A disulfide between 199-205
self.patch(residue_type='DISU', residues=(self.residues['199'],
self.residues['205']))
log.verbose() # request verbose output
env = environ() # create a new MODELLER environment to build this model in
# directories for input atom files
env.io.atom_files_directory = './:..'
a = mymodel(env,
alnfile = 'align2.ali', # alignment filename
knowns = 'e791', # codes of the templates
sequence = 'spnG') # 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.make
() # do the actual homology modeling