- modeller_usage - salilab.org

Modeller Model Building Inquiry using existing PDB cooridnates
by Joel Subach 13 Aug '20

13 Aug '20

Hi I want to build a C and N tail and one cytoplasmic loop that are missing from the PDB Cryo-EM 3-D coordinate file of a transmembrane protein. These regions were too dynamic during their Cryo-EM process and therefore are left out of the structural file that the investigators input into the PDB. Is there a way to model just these missing tails and one loop WITH the rest of the protein meaning not modeling the regions separately and then stitching them into the existing PDB coordinate file subsequently? I did not see this methodology amongst the Modeller tutorials but think it is best to maintain the fidelity of the protein as a whole i.e. non-stitching. Thanks:) Joel

2 1

Change ROOT_NAME for automodel.make()
by Brian Wiley 07 Aug '20

07 Aug '20

Hi, I thought this would be extremely easy to find if you can do, but I guess not. I would like to be able to change the ROOT_NAME of the model from automodel.make(). For instance: 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 Is this possible? Brian Wiley JHU

3 5

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