- modeller_usage - salilab.org

radii.lib
by Joaquim Rui Rodrigues 09 Oct '08

09 Oct '08

Dear Modellers, I am trying to model a non-standard residue, with non-standard lib files, by following the examples in the wiki page "Using non-standard parameter files". This works fine for toplib and parlib. However, I also need to read my extended versions of radii.lib and radii14.lib. Here is the relevant part of my script: env = environ(restyp_lib_file='/Users/jrui/modeller/LIBS/restyp.lib') class MyModel(automodel): .............. a = MyModel(env, ....) a.toplib = '/Users/jrui/modeller/LIBS/top_heav.lib' a.parlib = '/Users/jrui/modeller/LIBS/par.lib' I tried the obvious a.radlib = '/Users/jrui/modeller/LIBS/radii.lib' but I get: runcmd_____W>: creation of new member 'radlib' in <Model containing 0 chains, 0 residues, and 0 atoms>: possible typo! I searched for toplib and parlib in the manual but could not anything about these keywords. How can I tell modeller to read my radii files? Many thanks, Rui Rodrigues -- Webmail ESTG de Leiria (http://webmail.estg.ipleiria.pt)

2 1

PCA problem
by Julien Pele 06 Oct '08

06 Oct '08

Dear all, My question deals with the command which performs Principal Component Analysis (PCA). I use an input matrix of pairwise distances calculated from percentage residue identity of a multiple protein sequence alignment. I executed the following script : ################# env = environ() aln = alignment(env, file='alignment.pir') aln.id_table(matrix_file='id.mat') env.principal_components(matrix_file='id.mat', file='alignment_mod.princ') ################# The file 'alignment_mod.princ' contains projected coordinates of the first two principal components. I compared these results with other programs which hold PCA and go by the same input matrix. For instance, ade4, a package of the R software. I executed the following commands under R : ################# require (ade4) x <- read.table("id.mat", sep="") y <- dudi.pca(x, scan = F, nf = 2) write.table(y$li, file = "alignment_ade4.princ") ################# The file 'alignment_ade4.princ' also contains projected coordinates of the first two principal components. However, the coordinates acquired with these two softwares are not the same. I would like to know if the PCA command of Modeller uses internal procedures which could explain this disagreement ? How does the treatment of input matrix work with Modeller ? Any help on this will be greatly appreciated. Best regards, Julien

1 0

clockwise pentamer to anticlockwise model
by plmallip＠mail.uh.edu 29 Sep '08

29 Sep '08

Dear colleagues, I want to model a pentamer. I have two templates. One of the templates is arranged anti-clockwise in the order ABCDE. The other template runs clockwise (AEDCB) and has identical chains. Thus if I overlap the two templates, subunit B overlaps with subunit E and so on. I want to do modelling in 2 steps. a. I want to get an initial anticlockwise model (ABCDE) from clockwise template (identical chains AEDCB). b. The model obtained from step 1 will be overlapped with the anticlockwise template to build some missing residues. I tried to get through step a by 1. Changing the alignment file 2. by changing chain identification number in the template 3. by changing the order of subunits arranged in the template Each time, I end up with the same clockwise model. Thanks in advance for any kind of suggestion in this regard. Prema, Graduate student, UOH.

2 1

Cluster PDBs
by Wang Yanli 29 Sep '08

29 Sep '08

Dear modeller, I had build 500 models, and want cluster the build models. I had try the Cluster_PDBs(http://www.salilab.org/modeller/wiki/Cluster%20PDBs), but I just got the average-structure. I want analysis the clusters of the build models. How to do it? Thanks very much!

2 1

Re: [modeller_usage] Need help in not modeling the propep feature in final model
by Modeller Caretaker 29 Sep '08

29 Sep '08

Saravanan wrote: > I am modeling a protein in which two subunits are linked by a "propep" > which is not present in crystal structurei.e, its cleaved when protein > is active. I want to generate models without the "PROPEP". > > Crystal Structure available looks like this: (chain a and b) > > aaaaaaaaaabbbbbbbbbbb > > Model sequence is like this : (chain a and b and p-propep) > > aaaaaaaaaaapppppppppppppppbbbbbbbbbbb > > Model structure should be like this : (chain a and b AND propep chain p > should not present) > > aaaaaaaaaaabbbbbbbbbbb I don't think I fully understand your question, because the model structure always has the same sequence as the "model sequence", by definition. If you don't want to have the propep sequence in your model, just don't include it in the model sequence in your alignment. If, on the other hand, you want to build a model containing the propep sequence but you don't want to write it into the output PDB file, you can either edit the output PDB file by hand, or you can use selection() to select the non-propep parts and then selection.write() to write out a PDB. Ben Webb, Modeller Caretaker -- modeller-care(a)salilab.org http://www.salilab.org/modeller/ Modeller mail list: http://salilab.org/mailman/listinfo/modeller_usage

1 0

W> MODEL atom not classified: GLY:OXT GLY
by Joaquim Rui Rodrigues 24 Sep '08

24 Sep '08

Dear Modellers, I am trying to build some missing loops in a PDB structure with modeller 9v4. For this, I am using the script found in examples/automodel/model-loop-define.py for which I only changed residue_range's, and alnfile, knowns and sequence. The protein is a dimer, and both subunits end with a Gly residue. I am getting this warning at the DOPE potential assessment: >> Model assessment by DOPE potential iatmcls_286W> MODEL atom not classified: GLY:OXT GLY iatmcls_286W> MODEL atom not classified: GLY:OXT GLY A related message in the archives suggests editing modlib/atmcls-melo.lib (http://www.salilab.org/archives/ modeller_usage/2004/msg00169.html). This file seems to be used for loop modelling, but not for DOPE evaluation, as suggested in the log file: openf___224_> Open ${LIB}/atmcls-mf I have also tried to use the dope_loopmodel class, instead of loopmodel. Here, I get the same warnings for both stages (loop modelling and DOPE assessment). How can I avoid such messages? Should I somehow edit atmcls-mf.lib? Many thanks, Rui Rodrigues -- Webmail ESTG de Leiria (http://webmail.estg.ipleiria.pt)

2 1

from ylwang: loop conformation sampling
by Wang Yanli 19 Sep '08

19 Sep '08

Dear modeller, Firt, thanks for Ben Webb's help on my question, I follow your advice and found that error in MASS line. Now, I want to ask questioins about loop conformation. A loop in involved in the ligand binding and also involved in the dimmer interface. In the templet, the loop is short, but for my protein, it is longer than templet by 12 residues. The protein can be monomer or dimmer, but only be funcational in dimmer form. So, I want study: 1: the impact of dimmer on the loop conformation. 2: the impact of ligand on the loop conformation. Can I do loop-refine process times for dimmer and monomer to the study 1 issue 1? And Can I do modeller by "modeling ligands in the binding site" to study the issue 2? Or, is there other ways for my issue. I had read some paper about "ligand-supported homology model", but can't understand it well. It seems a method, but no tool-package availble? Any advices, please! -- PH.D Student Shanghai Key Lab. of Chemical Biology School of Pharmacy East China University of Science and Technology

2 1

how to constrain automodel to remain closer to the template
by Jake Gunn-Glanville 19 Sep '08

19 Sep '08

I have been modeling a protein with a template I am quite confident in, and a loop that I would like to use a separate template for. When I generate 200 models, they all seem to show a great deal of structural diversity. Curious, I developed 20 models of a solved structure against itself, and found a bit of diversity in all of these resulting models as well. Is there a good way to constrain automodel so that it diverses less from the templates when thre are no obvious residue clashes forcing it to make adjustments? -- Jacob Glanville Informatics Associate Rinat Biotherapeutics & Bioinnovation Center; Pfizer Inc. 415 497 3552

2 1

How to introduce a covalent bond in a model?
by DimitryASuplatov 19 Sep '08

19 Sep '08

Hello, I want to link N-end of the chain alpha with C-end of the chain beta in my peptide. Terminal acids are quite near to each other (4 A), thus it is possible. I have made up some workarounds like switching the chains in pdb and seq files before making the alignment for modeller etc.; I dont like it though. I suggest that there should be some universal function built-in modeller for this kind of purpose. Do you know how to manage this? Thank you. I appreciate your help. SDA

2 1

CHARMM atom type is not recognized (not in RADII_LIB): 15.999
by Wang Yanli 09 Sep '08

09 Sep '08

Hi, Modeller When I add the ligand (Allosamidin) as new residue type. I use insightII to generate the rtf file, ad add information in *.lib files. But I still get bellow error. How can I do? Thanks! env.libs.topology.read(file='$(LIB)/top_heav.lib') File "/home/ylwang/bin/modeller9v4/modlib/modeller/topology.py", line 23, in read return self.append(file) File "/home/ylwang/bin/modeller9v4/modlib/modeller/topology.py", line 27, in append return _modeller.mod_topology_read(self._modpt, self.__libs.modpt, file) _modeller.ModellerError: read_to_231E> CHARMM atom type is not recognized (not in RADII_LIB): 15.999

2 1