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

[modeller_usage] Problem modeling B-sheet



Hi.

I am trying to model a protein that has a C-terminal region without template match. However, the secondary structure prediction strongly suggests that part of this terminal region forms a beta-strand (residues 125 to 128).
I am trying to generate an antiparallel beta sheet from this terminal segment, bonding with a pre-existing beta strand (residues 79-82).
However, the generated model does not answer this command, being equal to the model generated without the command.

Could you point out what I am doing wrong?

Below are the script file.

Thanks in advance.


# Homology modeling with multiple templates
from modeller import *              # Load standard Modeller classes
from modeller.automodel import *    # Load the automodel class

log.verbose()    # request verbose output
env = environ()  # create a new MODELLER environment to build this model in
env.io.water = True

# Redefine the special_patches routine to include the additional disulfides
# (this routine is empty by default):
class MyModel(automodel):
    def special_patches(self, aln):
        rsr = self.restraints
        at = self.atoms

#       A disulfide between residues:
        self.patch(residue_type='DISU', residues=(self.residues['9'], self.residues['30']))
        self.patch(residue_type='DISU', residues=(self.residues['53'], self.residues['74']))

#       Two beta-strands:
        rsr.add(secondary_structure.strand(self.residue_range('79:', '82:')))
        rsr.add(secondary_structure.strand(self.residue_range('125:', '128:')))
#       An anti-parallel sheet composed of the two strands:
        rsr.add(secondary_structure.sheet(at['N:79'], at['O:128'],
                                          sheet_h_bonds=-4))


a = MyModel(env, alnfile = 'Maycon3-1template.ali', knowns = '1template', sequence = 'Maycon3', assess_methods=(assess.DOPE, assess.GA341))

a.starting_model = 1
a.ending_model = 100
a.make()