- modeller_usage - salilab.org

Re: [modeller_usage] modeller_usage Digest, Vol 19, Issue 52
by Brian Wiley 07 Aug '20

07 Aug '20

Nevermind. I just created a new function reset_defaults() to rewrite self.set_defaults() in automodel and I can get it. Thanks for your help. -----Original Message----- From: modeller_usage <modeller_usage-bounces(a)salilab.org> On Behalf Of modeller_usage-request(a)salilab.org Sent: Friday, August 7, 2020 6:47 PM To: modeller_usage(a)salilab.org Subject: modeller_usage Digest, Vol 19, Issue 52 External Email - Use Caution Send modeller_usage mailing list submissions to modeller_usage(a)salilab.org To subscribe or unsubscribe via the World Wide Web, visit https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fsalilab.o… or, via email, send a message with subject or body 'help' to modeller_usage-request(a)salilab.org You can reach the person managing the list at modeller_usage-owner(a)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: Change ROOT_NAME for automodel.make() (Modeller Caretaker) 2. Re: Change ROOT_NAME for automodel.make() (Brian Wiley) ---------------------------------------------------------------------- Message: 1 Date: Fri, 7 Aug 2020 15:33:39 -0700 From: Modeller Caretaker <modeller-care(a)salilab.org> To: Brian Wiley <bwiley4(a)jhu.edu>, "modeller_usage(a)salilab.org" <modeller_usage(a)salilab.org> Subject: Re: [modeller_usage] Change ROOT_NAME for automodel.make() Message-ID: <9ffc2904-4103-32e6-3c62-a4f8fd11a4f7(a)salilab.org> Content-Type: text/plain; charset=windows-1252; format=flowed On 8/7/20 2:21 PM, Brian Wiley wrote: > a = automodel(env, > alnfile = 'TSC2_6ssh.ali', > knowns = ('6ssh'), > sequence = 'TSC2_HUMAN') > > a.make() # Here I want the prefix to be 'TSC2_6ssh' and not > 'TSC2_HUMAN' so all files and pdb file at end would be > TSC2_6ssh.B99990001.pdb automodel is written under the assumption that files are named for the sequence. If you want to name the output PDB files differently, that's straightforward; just override get_model_filename() as per https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fsalilab.o… If you want to name the other output files differently, there isn't currently a mechanism for that, but it would be easy enough to modify the Python code to do so. Ben Webb, Modeller Caretaker -- modeller-care(a)salilab.org https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fsalilab.o… Modeller mail list: https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fsalilab.o… ------------------------------ Message: 2 Date: Fri, 7 Aug 2020 22:47:13 +0000 From: Brian Wiley <bwiley4(a)jhu.edu> To: Modeller Caretaker <modeller-care(a)salilab.org>, "modeller_usage(a)salilab.org" <modeller_usage(a)salilab.org> Subject: Re: [modeller_usage] Change ROOT_NAME for automodel.make() Message-ID: <BYAPR01MB3622652122731C17B45282EAEF490(a)BYAPR01MB3622.prod.exchangelabs.com> Content-Type: text/plain; charset="us-ascii" Thanks, I would like to change the root_name for all the output file. It's kind of silly there is not class root_name attribute in automodel. Then you could just name that. I am having trouble overriding the __init__ in automodel because it calls self.set_defaults(). Brian -----Original Message----- From: Modeller Caretaker <modeller-care(a)salilab.org> Sent: Friday, August 7, 2020 6:34 PM To: Brian Wiley <bwiley4(a)jhu.edu>; modeller_usage(a)salilab.org Subject: Re: [modeller_usage] Change ROOT_NAME for automodel.make() External Email - Use Caution On 8/7/20 2:21 PM, Brian Wiley wrote: > a = automodel(env, > alnfile = 'TSC2_6ssh.ali', > knowns = ('6ssh'), > sequence = 'TSC2_HUMAN') > > a.make() # Here I want the prefix to be 'TSC2_6ssh' and not > 'TSC2_HUMAN' so all files and pdb file at end would be > TSC2_6ssh.B99990001.pdb automodel is written under the assumption that files are named for the sequence. If you want to name the output PDB files differently, that's straightforward; just override get_model_filename() as per https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fsalilab.o… If you want to name the other output files differently, there isn't currently a mechanism for that, but it would be easy enough to modify the Python code to do so. Ben Webb, Modeller Caretaker -- modeller-care(a)salilab.org https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fsalilab.o… Modeller mail list: https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fsalilab.o…

1 0

Re: [modeller_usage] Treading onto known PDB template
by Modeller Caretaker 20 Jul '20

20 Jul '20

On 7/6/20 2:49 PM, Geoffrey Woollard wrote: > Can I do the modelling in the context of the tetramer? Otherwise I would > have to just do it with the protomer, and then delete out parts that > were clashing. > > Are there options for applying symmetry during modelling? Sure, this is also covered in the manual: https://salilab.org/modeller/9.24/manual/node29.html Ben Webb, Modeller Caretaker -- modeller-care(a)salilab.org https://salilab.org/modeller/ Modeller mail list: https://salilab.org/mailman/listinfo/modeller_usage

2 5

Re: [modeller_usage] reg: modeling against docked protein-ligand complex
by Anwesha Mohapatra 16 Jul '20

16 Jul '20

Hello Sir, Could you tell me where am going wrong? On Wed, Jul 15, 2020, 2:10 PM Anwesha Mohapatra < anwesha.mohapatra11(a)gmail.com> wrote: > > > <https://mail.google.com/mail/u/0?ui=2&ik=2650c6341a&attid=0.1&permmsgid=msg…> > Hello Sir, > > I have a protein (6eqe.pdb) which I have docked to a ligand using autodock > vina. I want to model my query gene taking this receptor-ligand complex as > the template. The ligand is identified as UNK in the pdb file. On modeling > I am not getting all the atoms of the ligand in the modelled output (Only 1 > carbon and 1 oxygen).Attached to this mail are the template file which has > been docked to the ligand *(6eqe_bhet.pdb) , *the modelled protein with > incomplete/erroneous ligand > *(ARN19002.1_Rhizobacter_gummiphilus_NS21.B99990001.pdb)* and the > alignment file for the same* > (ARN19002.1_Rhizobacter_gummiphilus_NS21.ali) . * > Could you please look into the files? > > Kindly let me know where I am going wrong. > > Thanks and regards > Anwesha >

1 0

Multiple CPUs
by Casey Becker 06 Jul '20

06 Jul '20

Hello again, I'm trying to use modeller directly through python's IDLE interface to run a parallel job. I've attached my script and also the python core file that's created up to a point, because it stops before it actually makes the models. Three files with .slave0, .slave1, and .slave2 are created, but when edited with notepad, show an error. I've attached one of them as well. Thanks.

2 1

Treading onto known PDB template
by Geoffrey Woollard 06 Jul '20

06 Jul '20

I'm studying a transmembrane ion channel (protein 1). There are structures of protein 1 in a closed conformation, and I am interested in the open conformation. There is a very similar protein in an open conformation (protein 2). So I am trying to thread the sequence of protein 1 onto the atomic model of protein 2. I know the template. I tried robetta and swiss prot for this. How should I do things with modeller? Could someone point me to a tutorial? The sequences are quite similar, except for some insertions. I know my protein very well and am looking for a protein modelling approach where I can have low level control, and manual manipulation at the sequence alignment and structure level. I don't need something that does everything automagically. -- Geoffrey Woollard, MSc PhD student Department of Medical Biophysics, University of Toronto Princess Margaret Cancer Centre, University Health Network

2 2

Re: [modeller_usage] z-DOPE score
by Mariana GONZALEZ MEDINA 22 Jun '20

22 Jun '20

Makes sense now, thank you! Mariana ________________________________ From: Min-yi Shen <minyishen(a)gmail.com> Sent: Friday, June 19, 2020 5:27:08 PM To: Modeller Caretaker Cc: Mariana GONZALEZ MEDINA; modeller_usage(a)salilab.org Subject: Re: [modeller_usage] z-DOPE score Original DOPE author here. The Z-dope is effectively an estimated z-score of your protein model based on the protein composition and size, as the raw DOPE scores are size and protein dependent. Z-dope used a linear model to estimate the mean and sd of a protein’s DOPE score based on their sizes and compositions. It was trained by millions of models from ModBase. -1.5 means your model’s DOPE score is -1.5 standard deviation lower than the average... usually a good indicator of good native overlap. I can’t believe I can still recall all these... On Thu, Jun 18, 2020 at 11:59 Modeller Caretaker <modeller-care(a)salilab.org<mailto:modeller-care@salilab.org>> wrote: On 6/18/20 9:48 AM, Mariana GONZALEZ MEDINA wrote: > I have a question regarding the zDOPE score. I have read some paper that > state this: "The Z-DOPE scores ranged from −1.63 to −1.85, where a score > of less than −1 indicates a “reliable” model (i.e., 80% of its Cα atoms > are within 3.5 Å of their correct positions)" > > I am wondering how is it that a score of less than -1 indicates that 80% > of its Cα atoms are within 3.5 Å of their correct positions? I would say "implies" rather than "indicates". DOPE is a model evaluation function. Of course we don't know the true structure so we can only guess with some degree of confidence (the more negative the score, the more confident I would be) whether a good DOPE score really means a correct model. Benchmarking does show that this is typically the case though; see the DOPE paper. Ben Webb, Modeller Caretaker -- modeller-care(a)salilab.org<mailto:modeller-care@salilab.org> https://salilab.org/modeller/ Modeller mail list: https://salilab.org/mailman/listinfo/modeller_usage _______________________________________________ modeller_usage mailing list modeller_usage(a)salilab.org<mailto:modeller_usage@salilab.org> https://salilab.org/mailman/listinfo/modeller_usage

1 0

z-DOPE score
by Mariana GONZALEZ MEDINA 18 Jun '20

18 Jun '20

Hi, I have a question regarding the zDOPE score. I have read some paper that state this: "The Z-DOPE scores ranged from −1.63 to −1.85, where a score of less than −1 indicates a “reliable” model (i.e., 80% of its Cα atoms are within 3.5 Å of their correct positions)" I am wondering how is it that a score of less than -1 indicates that 80% of its Cα atoms are within 3.5 Å of their correct positions? How is it possible to know the percent of atoms in the "correct position" based on how small (or not) the z-DOPE is? What can I say about the magnitude of the negative values besides stating that since the models have negative values they are acceptable? Best, Mariana

2 1

Modeller Outputs TER field wrong in PDB
by Pedro Guillem 18 Jun '20

18 Jun '20

I had some issues today using some PDBs generated by Modeller, and I tracked the error to the fact the TER fields inside the PDB have more columns (which should not be there) For example ATOM 1393 C2' DT D 35 81.341 17.516 55.945 1.00 48.63 C ATOM 1394 C3' DT D 35 82.152 18.800 55.754 1.00 48.81 C ATOM 1395 O3' DT D 35 83.036 19.108 56.819 1.00 49.40 O *TER 1396 DT D 35* ATOM 1397 N SER A 20 33.719 48.929 107.140 1.00 49.89 N ATOM 1398 CA SER A 20 34.080 48.490 108.522 1.00 50.05 C ATOM 1399 CB SER A 20 33.010 48.924 109.527 1.00 50.13 C I understand that the TER separator is not supposed to have anything afterwards. Is there a particular reason why Modeller is writing my TER records like that? It should read only "TER" Best Pedro

2 3

Can I calculate the "Objective Function" based on a PDB without doing the comparative modelling?
by ZHANG Cheng 14 Jun '20

14 Jun '20

The "Objective Function" is automatically shown on the second line of remark after doing the MODELLER comparative modelling. Now I have a structure generated by RosettaCM, and I want to use the MODELLER's "Objective Function" to evaluate its reliability. Would it be possible? Thank you!

2 1

Can't really edit my script
by Casey Becker 11 Jun '20

11 Jun '20

Hello, I'm new to modeller and I'm eventually trying to make several hundred models using multiple cpus, however it seems I can't change anything in my script except number of models without causing the system to break and not work. No error pops up, it simply says the task completed instantaneously, doesn't make any models, and doesn't pop up the ok_model window. My script is below from modeller import * from modeller.automodel import * # ---------------------- namelist.dat -------------------------------- # A "namelist.dat" file contains the file names, which was output from # the Modeller dialog in Chimera based on user's selection. # The first line it the name of the target sequence, the remaining # lines are name of the template structures namelist = open( 'namelist.dat', 'r' ).read().split('\n') tarSeq = namelist[0] template = tuple( [ x.strip() for x in namelist[1:] if x != '' ] ) # ---------------------- namelist.dat -------------------------------- # This instructs Modeller to display all log output. log.verbose() # create a new Modeller environment env = environ() # Directory of atom/PDB/structure files. It is a relative path, inside # a temp folder generated by Chimera. User can modify it and add their # absolute path containing the structure files. env.io.atom_files_directory = ['.', './template_struc'] # create a subclass of automodel or loopmodel, MyModel. # user can further modify the MyModel class, to override certain routine. class MyModel(automodel): def customised_function(self): pass #code overrides the special_restraints method #def special_restraints(self, aln): #code overrides the special_patches method. # e.g. to include the addtional disulfides. #def special_patches(self, aln): a = MyModel(env, sequence = tarSeq, # alignment file with template codes and target sequence alnfile = 'alignment.ali', # PDB codes of the templates knowns = template) # index of the first model a.starting_model = 1 # index of the last model a.ending_model = 5 loopRefinement = False # Assesses the quality of models using # the DOPE (Discrete Optimized Protein Energy) method (Shen & Sali 2006) # and the GA341 method (Melo et al 2002, John & Sali 2003) a.assess_methods = (assess.GA341, assess.normalized_dope) # ------------------------- build all models -------------------------- a.make() I've tried changing the MyModel part to just be automodel, because I was told it was a custom subclass and needed special treatment, which broke the code. I also tried changing the knowns = template to knowns = my template ID. It also broke the code. Thanks.

2 1