Dear All,
Please, I cannot model a beta strand (predicted as such by PSI PRED) in a protein model, even adding further restraints. I think this
depends on the fact that this region is not structured in the templates.
Please, any suggestions (see .py below)?
Many thanks.
Claudia
-------
# 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
# directories for input atom files
env.io.atom_files_directory = ['.', '../atom_files']
class MyModel(automodel):
def special_restraints(self, aln):
rsr = self.restraints
at = self.atoms
# unpick('1:', '73:')
# condense('1:', '73:')
# Add some restraints from a file:
# rsr.append(file='my_rsrs1.rsr')
# Residues 49 through 55 should be an alpha helix:
rsr.add(secondary_structure.alpha(self.residue_range('39:', '45:')))
# Two beta-strands:
# rsr.add(secondary_structure.strand(self.residue_range('1:', '5:')))
# rsr.add(secondary_structure.strand(self.residue_range('19:', '26:')))
# rsr.add(secondary_structure.strand(self.residue_range('33:', '35:')))
# rsr.add(secondary_structure.strand(self.residue_range('55:', '58:')))
rsr.add(secondary_structure.strand(self.residue_range('68:', '72:')))
# An anti-parallel sheet composed of the two strands:
# rsr.add(secondary_structure.sheet(at['N:1'], at['O:14'],
# sheet_h_bonds=-5))
# Use the following instead for a *parallel* sheet:
# rsr.add(secondary_structure.sheet(at['N:1'], at['O:9'],
# sheet_h_bonds=5))
# Restrain the specified CA-CA distance to 10 angstroms (st. dev.=0.1)
# Use a harmonic potential and X-Y distance group.
# rsr.add(forms.gaussian(group=physical.xy_distance,
# feature=features.distance(at['CA:35'],
# at['CA:40']),
# mean=10.0, stdev=0.1))
a = MyModel(env,
alnfile = 'Alignment2.ali', # alignment filename
knowns = ('pco','Tas','gen','Rap'), # codes of the templates
sequence = 'sALK1') # code of the target
a.starting_model= 1 # index of the first model
a.ending_model = 10 # index of the last model
# (determines how many models to calculate)
a.make() # do the actual homology modeling
--
Claudia Scotti, MD PhD
Department of Experimental Medicine
University of Pavia
Via Ferrata, 1
27100 Pavia
Italy
Tel. 0382 986335
Facs 0382 986893
Claudia Scotti
Dipartimento di Medicina Sperimentale
Sezione di Patologia Generale
Universita' di Pavia
Via Ferrata, 1
27100 Pavia
Italia
Tel. 0039 0382 986335/8/1
Facs 0039 0382 303673
