Hi Modeller users,
I want to add some additional restraints during homology modelling, but
it seems to be not working. And I did not found any error complaining
the restraints.
Here is my script:
from modeller import *
from modeller.automodel import *
log.verbose()
env = environ()
env.io.atom_files_directory = '/home/qinghua/Projects/Tycko/HM'
class MyModel(automodel):
def special_restraints(self, aln):
rsr = self.restraints
at = self.atoms
# CONSTRAINTS ON HISTIDINE 6
rsr.add(forms.gaussian(group=physical.xy_distance,
feature=features.distance(at['CG:13:B'],at['CG:14:B']),
mean=4.00, stdev=0.05))
# rsr.add(forms.gaussian(group=physical.xy_distance,
# feature=features.distance(at['NE2:14:B'],at['O:13:B']),
# mean=3.03, stdev=0.05))
env = environ()
a = automodel(env, alnfile='2LMO-aln.ali',
knowns='abcdB', sequence='2LMO',
assess_methods=(assess.DOPE, assess.GA341))
a.starting_model = 1
a.ending_model = 5
a.make()
This is just a test, I want to make sure that the restraints work here.
However, I checked the output pdb files, and I found that the distance
between Atom CG and CG is larger
than 4 angstrom.
Could someone give me some suggestion for it? Thanks very much!
All the best,
Qinghua