I want to model ONLY magnesium ion (from protein B) on an already crystallized structure (protein A). So the sequence of my target and template is exactly the same. My alignment file is ---------------------------------------------------------------------------------------------------------------------------- >P1;proteinA structureX:proteinA: 37 :B:+266 :B:::-1.00:-1.00 -QYQVGPLLGSGGFGSVYSGIRVSDNLPVAIKHVEKDRISDWGEGTRVPMEVVLLKKVSSGFSGVIRLLDWFERPDSFVLILERPEPVQDLFDFITERGALQEELARSFFWQVLEAVRHCHNCGVLHRDIKDENILIDLNRGELKLIDFGSGALLKDTVYTDFDGTRVYSPPEWIRYHRYHGRSAAVWSLGILLYDMVCGDIPFEHDEEIIGGQVFFRQRVSSECQHLIRWCLALRPSDRPTFEEIQNHPWMQDVLLPQETAEIHLH*
>P1;proteinB structureX:proteinB: 401 :A:+1 :A:::-1.00:-1.00
$--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------*
>P1;target sequence:target: : : : ::: 0.00: 0.00
$QYQVGPLLGSGGFGSVYSGIRVSDNLPVAIKHVEKDRISDWGEGTRVPMEVVLLKKVSSGFSGVIRLLDWFERPDSFVLILERPEPVQDLFDFITERGALQEELARSFFWQVLEAVRHCHNCGVLHRDIKDENILIDLNRGELKLIDFGSGALLKDTVYTDFDGTRVYSPPEWIRYHRYHGRSAAVWSLGILLYDMVCGDIPFEHDEEIIGGQVFFRQRVSSECQHLIRWCLALRPSDRPTFEEIQNHPWMQDVLLPQETAEIHLH*
when i run this script
from modeller import *
from modeller.automodel import *
log.verbose()
env = environ()
env.io.atom_files_directory = ['.', '../atom_files'] env.io.hetatm = True
a = automodel(env, alnfile = 'alignment.ali',
knowns=('proteinA','proteinB'), sequence='target') a.starting_model= 1
a.ending_model = 1
a.make()
I am getting segmentation fault. Could you please help me? Regards Renu