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

[modeller_usage] Build missing residues for multisubunit protein



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