29 Jul
2013
29 Jul
'13
6:42 p.m.
Hi,
Do you have suggestions on how to build missing residues for a multisubunit protein while restricting non-missing residues from moving?
---
from modeller import *from modeller.automodel import * # Load the automodel classlog.verbose()env = environ()# directories for input atom filesenv.io.atom_files_directory = ['.', '../atom_files']
class MyModel(automodel): def select_atoms(self): return selection(self.residue_range('1', '1'), self.residue_range('1', '2'), self.residue_range('181','181'))
a = MyModel(env, alnfile = 'alignment.ali', knowns = 'nqo2', sequence = 'nqo2_fill')a.starting_model= 1a.ending_model = 1a.make()
---
Thank you,
Kelly