- modeller_usage - salilab.org

Problem when reconstruct the protein with know structure!
by liu xin 27 Aug '07

27 Aug '07

Hi all I divided the protein with know structure (bovine rhodopsin receptor) into three part: N-terminus region, transmembrane region and C-terminus region without change their initial coordinates, and I try to use the three parts as multiple templates to reconstruct the original structure of the protein. First, I aligned the three parts as the following way: N-terminus> aaaaaaaaa-------------------------------------* transmembrane> ----------------bbbbbbbbbbbbb---------------* C-terminus> --------------------------------------ccccccccc* original> aaaaaaaaabbbbbbbbbbbbbccccccccc* then I used the following script to generate the original structure of my protein: from modeller import * from modeller.automodel import * env = environ() a = automodel(env, alnfile='bhro_mult.ali', knowns=('bhro_n','bhro_tm','bhro_c'), sequence='bhro') a.starting_model = 1 a.ending_model =1 a.make() However, the generated model was different from the original one, looking closely, I found the difference mainly occurred in the N-terminus and C-terminus regions, their coordinates have changed greatly when compared with the original structure. I was missing something important? Or, did it in a wrong way? I am new to modeller, any suggestions are appreciated! Thanks in advance yours Xin Liu

1 0

Modeling a sequence based on a low identity to a template
by ESTER FUSTE DOMINGUEZ 26 Aug '07

26 Aug '07

Dear Sir/Madam, I'm trying to model a sequence based on a low identity template. QUESTION 1: I follow the example "difficult modeling" from Modeller tutorial. However, in the last step, when I run "mod9v2 model.py" I obtain this misake described below: <<'import site' failed; use -v for traceback Traceback (most recent call last): File "model.py", line 9, in ? a.make() File "C:\Archivos de programa\Modeller9v2\modlib\modeller\automodel\automodel py", line 99, in make self.homcsr(exit_stage) File "C:\Archivos de programa\Modeller9v2\modlib\modeller\automodel\automodel py", line 412, in homcsr self.check_alignment(aln) File "C:\Archivos de programa\Modeller9v2\modlib\modeller\automodel\automodel py", line 394, in check_alignment aln.check() File "C:\Archivos de programa\Modeller9v2\modlib\modeller\alignment.py", line 182, in check libs=self.env.libs.modpt) _modeller.error: read_te_291E> Sequence difference between alignment and pdb :>> I tried to solve it, but I could'nt. Could you help me, please? mGenTHREADER PREDICTION RESULTS: HTML formatted results can be found here: http://bioinf2.cs.ucl.ac.uk/psiout/d83ffff39d5b753b.mgen.html QUESTION 2: When my sequence was submitted to the mGenThreader server to fold assignment, the server didn't returned me any significat hit, all confidence levels were GUESS (p-value >= 0.1) Then, which hit should I take? The hit which have the smallest p-value or the hit with the highest percentage identity? Or maybe If I don't have any significat hit, all models built by MODELLER will be wrong.Is that true? QUESTION 3: In "difficult example" in MODELLER's tutorial, in this alignment, I don't know which is the criterion to choose 30 :A: 209 :A:::: in structureX. How do you count it? I don't undestand the way you count it. <<Alignment between the nsp16 sequence and the 1ej0A from mGenThreader results. C; mGenThreader alignment of 30133975 and 1ej0A C; CERT significance with an e-value of 1e-04 C; Percentage Identity = 14.4% >P1;1ej0A structureX:1ej0: 30 :A: 209 :A:::: -------GLRSRAWFKL----------------------------------DEIQQSDKLFKPGMTVVDL GA------APGGWSQYVVTQIGGKGRIIACDLLPMDPIVGVDFLQGDFRDELVMKALLERVGDSKVQVVM SDMAPNMSGTPAVDIPRAMYLVELALEMCRDVLAPGGSFVVKVFQGEGFDEYLREIRSLFTKVKVRKPDS SRARSREVYIVATGRKP* >P1;30133975 sequence::::::::: ASQAWQPGVAMPNLYKMQRMLLEKCDLQNYGENAVIPKGIMMNVAKYTQLCQYLNTLTLAVPYNMRVIHF GAGSDKGVAPG--TAVLRQWLPTGTLLVDSDLNDFVSDADSTLIGDCATVH----------TANKWDLII SDMYDPRTKHVTKENDSKEGFFTYLCGFIKQKLALGGSIAVKITEHS-WNADLYKLMGHFSWWTAFVTNV NA-SSSEAFLIGANYLG* File 30133975_1ej0A_mGenThreader.ali. Red residues were manually removed from the alignment>> Thank you very much. Ester ------------------------------------------------- This mail sent through IMP: http://horde.org/imp/

2 1

ERROR MESSAGE
by ESTER FUSTE DOMINGUEZ 26 Aug '07

26 Aug '07

------------------------------------------------- This mail sent through IMP: http://horde.org/imp/

1 0

ligand modeling in binding site of protein
by mahendra awale 24 Aug '07

24 Aug '07

hi i am doing ligand modelling in binding site of protein i am using script ligand/multiple-hetero.py as u specified in tutorial but getting error pdb files read when i run script but it shows errors when doing add_restraints that "one or more atoms absent from model" below is the error file after running this script Traceback (most recent call last): File "ligand.txt", line 17, in ? a.make() File "/usr/local/modeller8v2/modlib/modeller/automodel/automodel.py", line 100, in make self.homcsr(exit_stage) File "/usr/local/modeller8v2/modlib/modeller/automodel/automodel.py", line 342, in homcsr self.mkhomcsr(aln) File "/usr/local/modeller8v2/modlib/modeller/automodel/automodel.py", line 399, in mkhomcsr self.special_restraints(aln) File "ligand.txt", line 8, in special_restraints rsr.add(atom_ids=ids, restraint_parameters=(2,1,1,22,2,2,0,3.5,0.1)) File "/usr/local/modeller8v2/modlib/modeller/model.py", line 409, in add return self.top.add_restraint('restraints.add', mdl=self.__mdl.modpt, File "/usr/local/modeller8v2/modlib/modeller/util/top.py", line 329, in add_restraint return _modeller.add_restraint(mdl, *args) _modeller.error: add_res_442E> One or more atoms absent from MODEL: N1:101 N3:13 also shown below is log file check_ali___> Checking pairwise structural superpositions. Equivalent CA pairs with distance difference larger than 6.0 angstroms: ALN_POS TMPL1 TMPL2 RID1 RID2 NAM1 NAM2 DIST ---------------------------------------------------- 95 1 2 113 94 Q S 6.894 204 1 2 222 203 R F 6.290 230 1 2 248 227 P W 7.467 233 1 2 251 230 S V 7.787 END OF TABLE check_ali___> Checking the sequence-structure alignment. Implied target CA(i)-CA(i+1) distances longer than 8.0 angstroms: ALN_POS TMPL RID1 RID2 NAM1 NAM2 DIST ---------------------------------------------- END OF TABLE read_to_681_> topology.submodel read from topology file: 3 getf_______W> RTF restraint not found in the atoms list: residue type, indices: 19 293 atom names : C +N atom indices : 2385 0 getf_______W> RTF restraint not found in the atoms list: residue type, indices: 19 293 atom names : C CA +N O atom indices : 2385 2374 0 2386 omgdel__425W> Unselected all O C +N +CA dihedrals: 306 mdtrsr__446W> A potential that relies on one protein is used, yet you have at least one known structure available. MDT, not library, potential is used. iup2crm_279W> IUPAC atom not found in topology library; residue type index: OD1 2 Possible reasons are a non-standard PDB file, or a new residue/atom types. Compare the offending residue in the PDB file with its definition in the topology library. iup2crm_279W> IUPAC atom not found in topology library; residue type index: OD2 2 Possible reasons are a non-standard PDB file, or a new residue/atom types. Compare the offending residue in the PDB file with its definition in the topology library. add_res_442E> One or more atoms absent from MODEL: N1:101 N3:13 please tell me what to do? thank u Mahendra Eknath Awale Institute Scholar M. S. (Pharm.) Dept.of Pharmacoinformatics National Institute of Pharmaceutical Education and Research NIPER S.A.S.Nagar, Mohali, Punjab-160062 INDIA. Institute Web http://www.niper.gov.in/ --------------------------------- Download prohibited? No problem. CHAT from any browser, without download.

2 1

ligand modeling in binding site of protein
by mahendra awale 24 Aug '07

24 Aug '07

hi i am doing ligand modeling in protein active site i run ligand/multiple-hetero.py script but geting errors it reads the pdb file but at the restraints line it shows one or more atoms not found in model Mahendra Eknath Awale Institute Scholar M. S. (Pharm.) Dept.of Pharmacoinformatics National Institute of Pharmaceutical Education and Research NIPER S.A.S.Nagar, Mohali, Punjab-160062 INDIA. Institute Web http://www.niper.gov.in/ --------------------------------- Get the freedom to save as many mails as you wish. Click here to know how.

2 1

ligand modeling in binding site
by mahendra awale 23 Aug '07

23 Aug '07

Hi i am doing ligand modeling in binding site of protein i am using ligand/model-multiple-hetero.py script file but geting errors it is say that alignment sequence not found in pdb file the outpu after running script shown below File "ligand.txt", line 17, in ? a.make() File "/usr/local/modeller8v2/modlib/modeller/automodel/automodel.py", line 100, in make self.homcsr(exit_stage) File "/usr/local/modeller8v2/modlib/modeller/automodel/automodel.py", line 318, in homcsr aln = self.read_alignment() File "/usr/local/modeller8v2/modlib/modeller/automodel/automodel.py", line 312, in read_alignment aln.append(file=self.alnfile, align_codes=self.knowns+[self.sequence]) File "/usr/local/modeller8v2/modlib/modeller/alignment.py", line 36, in append return self.__int_append('alignment.append', io, libs, vars) File "/usr/local/modeller8v2/modlib/modeller/alignment.py", line 68, in __int_append libs=libs.modpt, **vars) File "/usr/local/modeller8v2/modlib/modeller/util/top.py", line 33, in read_alignment return _modeller.read_alignment(aln, io, libs, *args) _modeller.error: rdpir___648E> Alignment sequence not found in PDB file: 2 1xxx.pdb the log.file out put shown below rdpir___377W> Too many fields ( 11 - should be 10) in the second line of the sequence entry (fields are separated by a colon, :) structure:1GZ8.pdb::::::::0.00:00 rdpir___648E> Alignment sequence not found in PDB file: 2 1GZ8.pdb thank u Mahendra Eknath Awale Institute Scholar M. S. (Pharm.) Dept.of Pharmacoinformatics National Institute of Pharmaceutical Education and Research NIPER S.A.S.Nagar, Mohali, Punjab-160062 INDIA. Institute Web http://www.niper.gov.in/ --------------------------------- Once upon a time there was 1 GB storage in your inbox. Click here for happy ending.

2 1

applying symmetry
by Christian Meesters 23 Aug '07

23 Aug '07

Hoi, First of all, thanks for your last explanation, Ben (was: 'code and logic question'). I just forgot to reply. My problem now is, how to create and apply arbitrary symmetry: The manual provides the nodes 28 and 103- http://www.salilab.org/modeller/manual/node28.html http://www.salilab.org/modeller/manual/node103.html - which already is a good starting point, but even after looking into the python code it is not clear to me how to distinguish between different symmetries. How, for example would I have to write: "Add restraints for a model with D6 symmetry (Schönflies notation), where all subunits have identical sequences."? Or "add restraints for a hexameric model with C3 symmetry."? Are there examples beyond the two nodes which address this question? TIA Christian

2 1

help--model_ini.py
by jinlian05＠lzu.cn 23 Aug '07

23 Aug '07

Hi: My initial model is one chain, I use model_ini.py scripts to build model, and edit my alignment file to add a starting residue number, such as: sequence:bhro:1::348:::::0.00:0.00 it also appears errors: rdpdb___303E> No atoms were read from the specified input PDB file, since the starting residue number and/or chain id in MODEL_SEGMENT (or the alignment file header) was not found; requested starting position: residue number " ", chain " " now, how does MODELLER solve it? thank you!

2 1

Reg. align_codes error in profile_build.py in Mod9v2
by saravanan 23 Aug '07

23 Aug '07

Please help me to rectify the error in mod9v2 My script from modeller import * log.verbose() env = environ() #-- Prepare the input files #-- Read in the sequence database sdb = sequence_db(env) #sdb.read(seq_database_file='C:\Program Files\Modeller9v2\modlib\pdball.pir', seq_database_format='PIR',chains_list='ALL', minmax_db_seq_len=(30, 4000), clean_sequences=True) #-- Write the sequence database in binary form #sdb.write(seq_database_file='C:\Program Files\Modeller9v2\modlib\pdball.bin', seq_database_format='BINARY', chains_list='ALL') #-- Now, read in the binary database sdb.read(seq_database_file='C:\Program Files\Modeller9v2\modlib\pdball.bin', seq_database_format='BINARY', chains_list='ALL') #-- Read in the target sequence/alignment aln = alignment(env) aln.append(file='p17787.fas', alignment_format='FASTA', align_codes='p17787') #-- Convert the input sequence/alignment into # profile format prf = aln.to_profile() #-- Scan sequence database to pick up homologous sequences prf.build(sdb, matrix_offset=-450, rr_file='${LIB}/blosum62.sim.mat', gap_penalties_1d=(-500, -50), n_prof_iterations=1, check_profile=False, max_aln_evalue=0.01) #-- Write out the profile in text format prf.write(file='build_profile.prf', profile_format='TEXT') #-- Convert the profile back to alignment format aln = prf.to_alignment() #-- Write out the alignment file aln.write(file='build_profile.ali', alignment_format='PIR') aln.check() error what i have got 'import site' failed; use -v for traceback Traceback (most recent call last): File "build_profile.py", line 21, in ? aln.append(file='p17787.fas', alignment_format='FASTA', align_codes='p17787' ) File "C:\Program Files\Modeller9v2\modlib\modeller\alignment.py", line 74, in append rewind_file, close_file) _modeller.error: read_al_373E> Protein specified in ALIGN_CODES(i) was not found in the alignment file; ALIGN_CODES( 1) = p17787 ____________________________________________________________________________________ Be a better Globetrotter. Get better travel answers from someone who knows. Yahoo! Answers - Check it out. http://answers.yahoo.com/dir/?link=list&sid=396545469

2 1

Re: [modeller_usage] Reg. comparing dope profile or both template and model
by saravanan 22 Aug '07

22 Aug '07

Now, we generated dope plot for template and model individually. How to visualize dope plot of both template and model in a single plot as mentioned in the basic tutorial.. Looking forward for reply... ---P.Saravanan ____________________________________________________________________________________ Moody friends. Drama queens. Your life? Nope! - their life, your story. Play Sims Stories at Yahoo! Games. http://sims.yahoo.com/

1 0