Thijs Beuming wrote: > I have some problems making a disulphide bridge in MODELLER. Any help > would be greatly appreciated. > > Thanks, > > Thijs Beuming > > My script: > > # Homology modeling by the automodel class > from modeller import * # Load standard Modeller classes > from modeller.automodel import * # Load the automodel class > > # 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 residues 176 and 185: > self.patch(residue_type='DISU', residues=(self.residues['176'], > self.residues > ['185']))
According to your output, you are running Modeller 8v2. But this script you have written is intended for Modeller 9v1. Either upgrade to 9v1, or look in the 8v2 manual (http://salilab.org/modeller/8v2/manual/node124.html) for the old way to patch residues.
> env.io.atom_files_directory = '~/NTProject2005/NEW_TCA/MODELLER'
Neither Python nor Modeller will expand '~' in path names. You can use Python's os.path module to do this (http://www.python.org/doc/2.3.5/lib/module-os.path.html) or Modeller will expand environment variables such as ${HOME}.
Ben Webb, Modeller Caretaker