- modeller_usage - salilab.org

mutate_model on protein with ligand
by Anna Marabotti 18 Nov '13

18 Nov '13

Dear Modeller caretaker, I used mutate_model.py to introduce a single point mutation into the structure of my protein with a ligand bound into the active site. The mutation I introduced was near the ligand, and when I had a look at the resulting structure, I saw that the mutated residue has been located in too strict proximity with the ligand (i.e. the mutated residue and the ligand collide with each other). My suspicion is that this script does not take into account the presence of ligands or other elements marked as HETATM in the original protein. Is this true? In order to encompass this problem, if this is true, can I simply replace HETATM with ATOM in order to make mutate_model.py correctly functioning? Many thanks in advance and best regards Anna Marabotti -- ______________________________________________ Anna Marabotti, Ph.D. Assistant Professor Department of Chemistry and Biology University of Salerno Via Giovanni Paolo II, 132 84084 Fisciano (SA) Italy Phone: +39 089 969583 Fax: +39 089 969603 E-mail: amarabotti(a)unisa.it Skype: annam1972 Web page: http://www.unisa.it/docenti/annamarabotti/index "Indifference is the eighth deadly sin" (don Andrea Gallo)

2 1

Multiple models for a single template
by vigneshwari subramanian 15 Nov '13

15 Nov '13

Dear Modeller users, I am interested in creating models for 100 different targets based on a single template. Is there an easy way to do this in modeller? Thanks, vigneshwari

2 1

How to improve sequence alignment?!
by Amjad Farooq 13 Nov '13

13 Nov '13

Hi everyone, I am noticing that Modeller does not do a good job when it comes to sequence alignment of short peptides or DNA sequences. For example, please see the the comparison below: (1) Modeller alignment: Template GTPPPP - - YTVG Peptide TVLPPPPYRHRN (1) ClustalW alignment: Template -GTPPPPYTVG- Peptide TVLPPPPYRHRN As you can see, Modeller is unable to align up properly the central PPPPY motif in my peptide with the template!! In contrast, ClustalW is spot on in its ability to recognize a key pattern in the sequence. Is there a way to improve sequence alignment in Modeller? How can I make changes to the following Modeler file that I am using for my sequence alignment: # Make a pairwise alignment of two sequences aln = alignment(env, file='template.ali') aln.salign(overhang=30, gap_penalties_1d=(-450, -50), alignment_type='tree', output='ALIGNMENT') aln.write(file='template.pir', alignment_format='PIR') aln.write(file='template.pap', alignment_format='PAP') What parameters or settings, if any, could be changed in the above setup to improve sequence alignment in Modeller? Thank you so much for your advice and any comments that you may have. *AMJAD FAROOQ PhD DIC | Associate ProfessorDept of Biochemistry & Molecular Biology | Miller School of Medicine | University of Miami | Miami | FL 33136Located @ Gautier Building #217 | Mail @ 1011 NW 15th Street #217, Miami , FL 33136 amjad(a)farooqlab.net <amjad(a)farooqlab.net> | off 305-243-2429 | lab 305-243-9799 | fax 305-243-3955 | www.farooqlab.net <http://www.farooqlab.net>*

2 2

help regarding some problems related to modelling
by DHWANI DHOLAKIA 12 Nov '13

12 Nov '13

Dear sir/mam I am Dhwani Dholakia Doing M.Tech from Delhi Technological University. I am Using advanced homology modelling for my project work but while doing pdb search i could not find much good E value or % identity or query coverage.. I am confused that should i use a invivo modelling software like I-TASSER or advanced modelling . And also if can help me with the threshold values of various parameters to be used for advanced homology modelling.. Thanking You -- Dhwani Dholakia M.Tech bioinformatics Delhi Technological University

2 1

questions about multichains modeller
by Yi Kong 12 Nov '13

12 Nov '13

Dear colleagues, I have a problem in multichain modeller. according the template file in automodel folder, I created the input file and script file. however when I　run this script, I recieved warnings listed below. Who knows how to reslove the problem? Thanks a lot! You can find many useful example scripts in the examples\automodel directory. It is recommended that you use Python to run Modeller scripts. However, if you don't have Python installed, you can type 'mod9.12' to run them instead. d:\Program Files (x86)\Modeller9.12> d:\Program Files (x86)\Modeller9.12>cd examples d:\Program Files (x86)\Modeller9.12\examples>cd automodel d:\Program Files (x86)\Modeller9.12\examples\automodel>python model-multichain-d ec1.py > model-multichain-dec1.log Traceback (most recent call last): File "model-multichain-dec1.py", line 37, in <module> a.make() # do homology modeling File "d:\Program Files (x86)\Modeller9.12\modlib\modeller\automodel\automodel. py", line 96, in make self.homcsr(exit_stage) File "d:\Program Files (x86)\Modeller9.12\modlib\modeller\automodel\automodel. py", line 410, in homcsr aln = self.read_alignment() File "d:\Program Files (x86)\Modeller9.12\modlib\modeller\automodel\automodel. py", line 400, in read_alignment aln.append(file=self.alnfile, align_codes=self.knowns+[self.sequence]) File "d:\Program Files (x86)\Modeller9.12\modlib\modeller\alignment.py", line 78, in append allow_alternates) _modeller.ModellerError: read_al_373E> Protein specified in ALIGN_CODES(i) was n ot found in the alignment file; ALIGN_CODES( 1) = 4h10 d:\Program Files (x86)\Modeller9.12\examples\automodel>

2 1

multiple templates - loop modelling
by Simone Fulle 08 Nov '13

08 Nov '13

Dear all, I would like to model a protein using two templates and (!) define which parts of the structure are modeled by which template. I.e. the core of the structure by template_1 and a specific loop by template_2. The template structures are already aligned. Any suggestions and comments will be highly appreciated. I am aware that I could generate a "Fused_template_structure" manually but it would be really nice if I could do it in an automatic fashion. Kind regards, Simone

2 1

homology modelling
by pr＠tibh＠ 07 Nov '13

07 Nov '13

How can i build a homology model using template,my template is a crystal structure of a chimeric protein. though it has a 94% query coverage and 48% identity with my target sequence. sequence alignment of non chimeric protein as a whole is very much identical.Can i use the chimeric crystal structure for modelling?

2 1

Hardware?
by Istvan Jonyer 06 Nov '13

06 Nov '13

What hardware is everyone using for Modeller? What do you think would be the ideal configuration? Anyone using any cloud providers? ᐧ

1 0

Protein specified in ALIGN_CODES<i> was not found in the alignment file; ALIGN_CODES< 1> = 1N11
by listek_wfrp＠interia.pl 30 Oct '13

30 Oct '13

Hello everybody ! I am new to modeller, so please excuse me if my question seems trivial. I have a problem building a model from multiple templates. When i try to run build model script I get an error: "Protein specified in ALIGN_CODES<i> was not found in the alignment file; ALIGN_CODES< 1> = 1N11" I followed all the previous steps from advanced modeling tutorial and all of them worked fine. Now I am stuck, and all solutions I`ve tried dont work. I attatch my script and alignment file. Thank you for your help in advance !

2 1

Fwd: Extract Chain ID
by Amjad Farooq 28 Oct '13

28 Oct '13

I am not scanning against PDB95 but rather working on multi-chain PDBs locally. For instance, if a PDB has 4 chains (say A, D, F, P) and I am interested in extracting the chain ID of the second segment (D) that best aligns with my query sequence, I do not believe you can extract that from PIR alignment file, which only provides chain IDs of the first (A) and the last chain (P). Also, how would I know if my query sequence aligned up against chain D unless I look at the alignment file manually. I want to automate this procedure. In other words, have a python script extract the chain ID of the aligned segment in a multi-chain PDB and automatically proceed to next step for homology modeling--that is I only want to use segment D as a template during homology modeling rather than the whole multi-chain PDB, which compromises the quality of the alignment and hence the model. Thanks! On Mon, Oct 28, 2013 at 2:36 PM, Modeller Caretaker < modeller-care(a)salilab.org> wrote: > On 10/27/2013 11:08 AM, Amjad Farooq wrote: > >> Is it possible to perform an automatic search of a multi-chain PDB >> template with a known protein sequence that is to be modeled (say >> 100-200 residues in length) and then extract the chain ID of the segment >> of the PDB template that best matches with the query sequence. >> > > If I understand you correctly, this is straightforward. If you're scanning > against our PDB95 database, for example, all the templates are single > chain, and the chain ID is simply the last character in the align code. > > If you're trying to align against a multi-chain template (presumably a > local alignment) you could extract the chain ID trivially from the > alignment object. > > Ben Webb, Modeller Caretaker > -- > modeller-care(a)salilab.org http://www.salilab.org/**modeller/<http://www.salilab.org/modeller/> > Modeller mail list: https://salilab.org/mailman/**listinfo/modeller_usage<https://salilab.org/mailman/listinfo/modeller_usage> >

2 3