Hello
while performing the modellin using multitemplate I am facing problem . the model-multiple.py , align.ali and generated log files are pasted below . plz help memodel-multiple.py file
# 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']
a = automodel(env, alnfile = 'align-multiple.ali', # alignment filename knowns = ('2hyd', '3g5u'), # codes of the templates sequence = 'pgp_ref') # code of the target a.starting_model= 1 # index of the first model a.ending_model = 1 # index of the last model # (determines how many models to calculate) a.make() # do the actual homology modeling
C; A multiple alignment in the PIR format;
>P1;2hyd structureX:2hyd:1: A:578 :A:ferredoxin:Azotobacter vinelandii: 1.90: 0.19 ----------------------------------MIKRYLQFVKPYKYRIFATIIVGIIKFGIPMLIPLLIKYAIDGVIN ------------NHALTTDEKVHHLTIAIGIALFIFVIVRP-----PIEFIRQYLAQWTSNKILYDIRKKLYNHLQALSA RFYANNQVGQVISRVINDVEQTKDFILTGLMNIWLDCITIFGGFIIGFTRGWKLTLVILAISPVLGLSAGIWAKILSSFT DKELHAYAKAQGFLHERVQGISVVKSFAIEDNEAKNFDKKNTNFLTRALKHTRWNAYSFAAINTVTDIGPIIVIGVGAYL AISGSITVGTLAAFVGYLELLFGPLRRLVASFTTLTQSFASMDRVFQLIDEDYDIKNG-VGAQPIEIKQGRIDIDHVSFQ YND-NEAPILKDINLSIEKGETVAFVGMSGGGKSTLINLIPRFYDVTSGQILIDGHNIKDFLTGSLRNQIGLVQQDNILF SDTVKENILLGRPTATDEEV VEAAKMANAHDFIMNLPQGYDTEVGERGVKLSGGQKQRLSIARIFLNNPPILILDEATS ALDLESESIIQEALDVLSKDRTTLIVAHRLSTITHADKIVV IENGHIVETGTHRELIAKQGAYEHLYSIQNL------- -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- ----------*
>P1;3g5u structureX:3g5u:1:A: 1271 :A:ferredoxin:Bacillus schlegelii:-1.00:-1.00 MELEEDLKGRADK-NFSKMGKKSKKEKKEKKPAVSVLTMFRYAGWLDRLY MLVGTLAAIIHGVALPLMMLIFGDMTDSF ASVGNVS---KNSTNMSEADK RAMFAKLEEEMTTYAYYYTGIGAGVLIVAYIQVSFWCLAAGRQIHKIRQKFFHAIMNQ EIGWFDVHDVGELNTRLTDDVSKINEGIGDKIGMFFQAMATF FGGFIIGFTRGWKLTLVILAISPVLGLSAGIWAKILS SFTDKELHAYAKAGAVAEEVLAAIRTVIAFGGQKKELERYNNNLEEAKRLGIKKAITANISMGAAFLLIYASYALAFWYG TSLVISKEYSIGQVLTVFFSVLIGAFSVGQASPNIEAFANARGAAYEVFKIIDNKPSIDSFSKSGHKPDNIQGNLEFKNI HFSYPSRKEVQILKGLNLKVKSGQTVALVGNSGCGKSTTVQLMQRLYDPLDGMVSIDGQDIRTINVRYLREIIGVVSQEP VLFATTIAENIRYGREDVTMDEIEKAVKEANAYDFIMKLPHQFDTLVGERGAQLSGGQKQRIAIARALVRNPKILLLDEA TSALDTESEAVVQAALDKAREGRTTIVIAHRLSTVRNADVIAGFDGGVIVEQGNHDELMREKGIYFKLVMTQTAGNEIEL GNEACKSKDEIDNLDMSSKDSGSSLIRRRSTRKSICGPHDQDRKLSTKEALDEDVPPASFWRILKLNSTEWPYFVVGIFC AIINGGLQPAFSVIFSKVVGVFTNGGPPETQRQNSNLFSLLFLILGIISFITFFLQGFTFGKAGEILTKRLRYMVFKSML RQDVSWFDDPKNTTGALTTRLANDAAQVKGATGSRLAVIFQNIANLGTGIIISLIYGWQLTLLLLAIVPIIAIAGVVEMK MLSGQALKDKKELEGSGKIATEAIENFRTVVSLTREQKFETMYAQSLQIPYRNAMKKAHVFGITFSFTQAMMYFSYAAAF RFGAYLVTQQLMTFENVLLVFSAIVFGAMAVGQVSSFAPDYAKATVSASHIIRIIEKTPEIDSYSTQGLKPNMLEGNVQF SGVVFNYPTRPSIPVLQGLSLEVKKGQTLALVGSSGCGKSTVVQLLERFYDPMAGSVFLDGKEIKQLNVQWLRAQLGIVS QEPILFDCSIAENIAYGDNSRVVSYEEIVRAAKEANIHQFIDSLPDKYNTRVGDKGTQLSGGQKQRIAIARALVRQPHIL LLDEATSALDTESEKVVQEALDKAREGRTCIVIAHRLSTIQNADLIVVIQNGKVKEHGTHQQLLAQKGIYFSMVSVQAGA KRSYVHHHHHH*
>P1;pgp_ref sequence:pgp_ref:1 : :1280 : :ferredoxin:Peptococcus aerogenes: 2.00:-1.00 MDLEGDRNGGAKKKNFFKLNNKSEKDKKEKKPTVSVFSMFRYSNWLDKLYMVVGTLAAIIHGAGLPLMMLVFGEMTDIFA NAGNLEDLMSNITNRSDINDTGFFMNLEEDMTRYAYYYSGIGAGVLVAAYIQVSFWCLAAGRQIHKIRKQFFHAIMRQEI GWFDVHDVGELNTRLTDDVSKINEGIGDKIGMFFQSMATFFTGFIVGFTRGWKLTLVILAISPVLGLSAAVWAKILSSFT DKELLAYAKAGAVAEEVLAAIRTVIAFGGQKKELERYNKNLEEAKRIGIKKAITANISIGAAFLLIYASYALAFWYGTTL VLSGEYSIGQVLTVFFSVLIGAFSVGQASPSIEAFANARGAAYEIFKIIDNKPSIDSYSKSGHKPDNIKGNLEFRNVHFS YPSRKEVKILKGLNLKVQSGQTVALVGNSGCGKSTTVQLMQRLYDPTEGMVSVDGQDIRTINVRFLREIIGVVSQEPVLF ATTIAENIRYGRENVTMDEIEKAVKEANAYDFIMKLPHKFDTLVGERGAQLSGGQKQRIAIARALVRNPKILLLDEATSA LDTESEAVVQVALDKARKGRTTIVIAHRLSTVRNADVIAGFDDGVIVEKGNHDELMKEKGIYFKLVTMQTAGNEVELENA ADESKSEIDALEMSSNDSRSSLIRKRSTRRSVRGSQAQDRKLSTKEALDESIPPVSFWRIMKLNLTEWPYFVVGVFCAII NGGLQPAFAIIFSKIIGVFTRIDDPETKRQNSNLFSLLFLALGIISFITFFLQGFTFGKAGEILTKRLRYMVFRSMLRQD VSWFDDPKNTTGALTTRLANDAAQVKGAIGSRLAVITQNIANLGTGIIISFIYGWQLTLLLLAIVPIIAIAGVVEMKMLS GQALKDKKELEGSGKIATEAIENFRTVVSLTQEQKFEHMYAQSLQVPYRNSLRKAHIFGITFSFTQAMMYFSYAGCFRFG AYLVAHKLMSFEDVLLVFSAVVFGAMAVGQVSSFAPDYAKAKISAAHIIMIIEKTPLIDSYSTEGLMPNTLEGNVTFGEV VFNYPTRPDIPVLQGLSLEVKKGQTLALVGSSGCGKSTVVQLLERFYDPLAGKVLLDGKEIKRLNVQWLRAHLGIVSQEP ILFDCSIAENIAYGDNSRVVSQEEIVRAAKEANIHAFIESLPNKYSTKVGDKGTQLSGGQKQRIAIARALVRQPHILLLD EATSALDTESEKVVQEALDKAREGRTCIVIAHRLSTIQNADLIVVFQNGRVKEHGTHQQLLAQKGIYFSMVSVQAGTKRQ --------*
this the log file
ate and time of compilation : 2010/04/21 09:52:56 MODELLER executable type : i386-w32 Job starting time (YY/MM/DD HH:MM:SS): 2010/06/16 12:26:20
openf___224_> Open $(LIB)/restyp.lib openf___224_> Open ${MODINSTALL9v8}/modlib/resgrp.lib rdresgr_266_> Number of residue groups: 2 openf___224_> Open ${MODINSTALL9v8}/modlib/sstruc.lib
Dynamically allocated memory at amaxlibraries [B,KiB,MiB]: 3234072 3158.273 3.084
Dynamically allocated memory at amaxlibraries [B,KiB,MiB]: 3234600 3158.789 3.085 openf___224_> Open ${MODINSTALL9v8}/modlib/resdih.lib
Dynamically allocated memory at amaxlibraries [B,KiB,MiB]: 3283200 3206.250 3.131 rdrdih__263_> Number of dihedral angle types : 9 Maximal number of dihedral angle optima: 3 Dihedral angle names : Alph Phi Psi Omeg chi1 chi2 chi3 chi4 chi5 openf___224_> Open ${MODINSTALL9v8}/modlib/radii.lib
Dynamically allocated memory at amaxlibraries [B,KiB,MiB]: 3292440 3215.273 3.140 openf___224_> Open ${MODINSTALL9v8}/modlib/radii14.lib openf___224_> Open ${MODINSTALL9v8}/modlib/af_mnchdef.lib rdwilmo_274_> Mainchain residue conformation classes: APBLE openf___224_> Open ${MODINSTALL9v8}/modlib/mnch.lib rdclass_257_> Number of classes: 5 openf___224_> Open ${MODINSTALL9v8}/modlib/mnch1.lib openf___224_> Open ${MODINSTALL9v8}/modlib/mnch2.lib openf___224_> Open ${MODINSTALL9v8}/modlib/mnch3.lib openf___224_> Open ${MODINSTALL9v8}/modlib/xs4.mat rdrrwgh_268_> Number of residue types: 21 openf___224_> Open align-multiple.ali
Dynamically allocated memory at amaxalignment [B,KiB,MiB]: 3303853 3226.419 3.151
Dynamically allocated memory at amaxalignment [B,KiB,MiB]: 3305303 3227.835 3.152
Dynamically allocated memory at amaxalignment [B,KiB,MiB]: 3308203 3230.667 3.155
Dynamically allocated memory at amaxalignment [B,KiB,MiB]: 3314003 3236.331 3.160
Dynamically allocated memory at amaxalignment [B,KiB,MiB]: 3325603 3247.659 3.172
Dynamically allocated memory at amaxalignment [B,KiB,MiB]: 3348803 3270.315 3.194
Dynamically allocated memory at amaxalignment [B,KiB,MiB]: 3800163 3711.097 3.624
Dynamically allocated memory at amaxsequence [B,KiB,MiB]: 3802479 3713.358 3.626
Dynamically allocated memory at amaxsequence [B,KiB,MiB]: 3807623 3718.382 3.631
Dynamically allocated memory at amaxsequence [B,KiB,MiB]: 3812739 3723.378 3.636 read_al_375E> Unknown residue type,position,sequence: 501 1
On 6/16/11 1:05 AM, shikha agarwal wrote: > while performing the modellin using multitemplate I am facing problem ... > read_al_375E> Unknown residue type,position,sequence: 501 1
Modeller cannot identify the 501st residue in the first sequence in your alignment. If you look at your alignment, the first sequence is 2hyd, and if you look at the sequence around position 501, it is:
TATDEEV VEAAKMANAHDFIMNLPQ
A space character is not a valid residue type. Remove the space between EEV and VEA and the problem should go away.
Ben Webb, Modeller Caretaker
participants (2)
-
Modeller Caretaker
-
shikha agarwal