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 class
log.verbose()
env = environ()
# directories for input atom files
env.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= 1
a.ending_model = 1
a.make()
---
Thank you,
Kelly