Script is wrong,
Use knowns='5fd1A' (with ChainID) because this is a protein.
For ligand bound residue put a . in the alignment file, so that MODELLER can include the attached ligand too in the model.
Then it will correctly call
assess_methods=(assess.DOPE, assess.GA341)) (Last line of your initial set of lines).
If you are using 5FD1 as the PDB file to model your sequence 1fdx, then you should definitely have the DOPE Score calculated.
Ashish Runthala,
Lecturer, Structural Biology Cell,
Biological Sciences Group,
BITS, Pilani
Rajasthan, INDIA
----- Original Message -----
From: "modeller usage-request" modeller_usage-request@salilab.org
To: "modeller usage" modeller_usage@salilab.org
Sent: Tuesday, April 10, 2012 11:12:30 AM GMT +05:30 Chennai, Kolkata, Mumbai, New Delhi
Subject: modeller_usage Digest, Vol 11, Issue 44
Send modeller_usage mailing list submissions to
modeller_usage@salilab.org
To subscribe or unsubscribe via the World Wide Web, visit
https://salilab.org/mailman/listinfo/modeller_usage
or, via email, send a message with subject or body 'help' to
modeller_usage-request@salilab.org
You can reach the person managing the list at
modeller_usage-owner@salilab.org
When replying, please edit your Subject line so it is more specific
than "Re: Contents of modeller_usage digest..."
Today's Topics:
1. Re: modeller_usage Digest, Vol 11, Issue 42 (Modeller Caretaker)
2. calculating DOPE score (NUTAN CHAUHAN)
----------------------------------------------------------------------
Message: 1
Date: Mon, 09 Apr 2012 20:35:37 -0700
From: Modeller Caretaker modeller-care@salilab.org
To: Ashish Runthala ashishr@bits-pilani.ac.in
Cc: modeller usage modeller_usage@salilab.org
Subject: Re: [modeller_usage] modeller_usage Digest, Vol 11, Issue 42
Message-ID: 4F83AA89.9080408@salilab.org
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
On 4/9/12 8:00 PM, Ashish Runthala wrote:
> MODELLER based on your constructed alignment file
> parsing the reliable template residues mapped against the target
> ones, constructs and fixes the average alpha carbon backbone
Actually all atoms that have the same name in target and template will
be simply copied into the model in this case. Usually this means all
backbone atoms, not just alpha carbons.
> On increased sampling,
> MODELLER tries best to satisfy the Lennard Jones and all the
> energetic equations composing its DOPE energy function.
In fact the optimization tries to minimize the molecular PDF, which
includes contributions from soft sphere (not Lennard Jones). DOPE is a
separate, atomistic pairwise statistical potential. It is not used in
comparative modeling.
Ben Webb, Modeller Caretaker
--
modeller-care@salilab.org
http://www.salilab.org/modeller/
Modeller mail list:
http://salilab.org/mailman/listinfo/modeller_usage
------------------------------
Message: 2
Date: Tue, 10 Apr 2012 13:29:27 +0800 (SGT)
From: NUTAN CHAUHAN
nutanbt@yahoo.co.in
To: "modeller_usage@salilab.org"
modeller_usage@salilab.org
Subject: [modeller_usage] calculating DOPE score
Message-ID:
1334035767.17687.YahooMailNeo@web192502.mail.sg3.yahoo.com
Content-Type: text/plain; charset="iso-8859-1"
Hello all,
I am trying to build a protein-ligand model using align-ligand.ali file. . Below is ?the python script to generate 5 model and to ?calculate the Dope and GA score.?but I am unable to calculate the DOPE score and GA score. Plz help
from modeller import *
from modeller.automodel import *
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']
env = environ()
a = automodel(env, alnfile='alignment.ali',
? ? ? ? ? ? ? knowns='5fd1', sequence='1fdx',
? ? ? ? ? ? ? assess_methods=(assess.DOPE, assess.GA341))
a.starting_model = 1
a.ending_model = 5
a.make()
?
Regards,
Nutan Chauhan
Research Scholar
Department of Biotechnology
BIT, Mesra
Ranchi-835215, Jharkhand
________________________________
From: "modeller_usage-request@salilab.org"
modeller_usage-request@salilab.org
To: modeller_usage@salilab.org
Sent: Tuesday, 10 April 2012 8:58 AM
Subject: modeller_usage Digest, Vol 11, Issue 43
Send modeller_usage mailing list submissions to
??? modeller_usage@salilab.org
To subscribe or unsubscribe via the World Wide Web, visit
???
https://salilab.org/mailman/listinfo/modeller_usage
or, via email, send a message with subject or body 'help' to
??? modeller_usage-request@salilab.org
You can reach the person managing the list at
??? modeller_usage-owner@salilab.org
When replying, please edit your Subject line so it is more specific
than "Re: Contents of modeller_usage digest..."
Today's Topics:
? 1. Re: In what are the initial structures different?
? ? ? (Modeller Caretaker)
? 2. Re: modeller_usage Digest, Vol 11, Issue 42 (Ashish Runthala)
? 3. Re: modeller_usage Digest, Vol 11, Issue 42 (Ashish Runthala)
----------------------------------------------------------------------
Message: 1
Date: Mon, 09 Apr 2012 11:46:02 -0700
From: Modeller Caretaker
modeller-care@salilab.org
To: BIN ZHANG
zhngbn@gmail.com
Cc: modeller_usage@salilab.org
Subject: Re: [modeller_usage] In what are the initial structures
??? different?
Message-ID:
4F832E6A.5020001@salilab.org
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
On 04/08/2012 11:58 PM, BIN ZHANG wrote:
> I have a question about the differences of the initial structure built
> for homology modeling. Can anyone tell me how exactly are these initial
> structures different? Do they differ in the secondary structures in
> uncertain regions? I cannot find this information in the manual.
See
http://salilab.org/modeller/9.10/manual/node470.html, in particular
step 3 and substep 3. The initial structures are simply randomized by
calling selection.randomize_xyz().
??? Ben Webb, Modeller Caretaker
--
modeller-care@salilab.org? ? ? ? ? ?
http://www.salilab.org/modeller/
Modeller mail list:
http://salilab.org/mailman/listinfo/modeller_usage
------------------------------
Message: 2
Date: Tue, 10 Apr 2012 08:33:02 +0530 (IST)
From: Ashish Runthala??? ?
ashishr@bits-pilani.ac.in
To: modeller usage
modeller_usage@salilab.org
Subject: Re: [modeller_usage] modeller_usage Digest, Vol 11, Issue 42
Message-ID:
???
15240798.16751334026982399.JavaMail.root@mailserver.bits-pilani.ac.in
???
Content-Type: text/plain; charset=utf-8
Focus on the statement with ALIGN_CODES, as it is inclusive of the PDB_ID and the PDB_CHAIN used in the alignment file.
So it could probably be '2AHNA' or '2AHNB'.
Ashish Runthala,
Lecturer, Structural Biology Cell,
Biological Sciences Group,
BITS, Pilani
Rajasthan, INDIA
------------------------------
Message: 2
Date: Mon, 9 Apr 2012 17:21:12 +0530
From: Sumedha Roy
sumedharoy@gmail.com
To: modeller usage
modeller_usage@salilab.org
Subject: [modeller_usage] Error in ALIGN_CODES(i)
Message-ID:
???
CABFmE0dNMawUwZgsh07bsRd43+-O3Yrzq6VhfkSQJOTjp=NECA@mail.gmail.com
Content-Type: text/plain; charset="iso-8859-1"
Hello all,
When I try to model my protein-DNA complex, I get the error:
_modeller.ModellerError: read_al_373E> Protein specified in ALIGN_CODES(i)
was not found in the alignment file; ALIGN_CODES(? ? ? 2) =? 2AHN
But I see no mistake either in the alignment or the format. I've attached
both the ALI file and model-ligand.py. Please help!!!
--
Sumedha Roy
Ashish Runthala
Modeller Caretaker