- modeller_usage - salilab.org

Re: [modeller_usage] modeller_usage Digest, Vol 11, Issue 54
by Ashish Runthala 04 May '12

04 May '12

No, Mary Varughese. Sir just meant that you should parse your sequence mentioned as M in his alignment (with / being the chain delimiter), (M: your target Model residues and not Missing residue) He simply pointed that You can easily cover up the missing residues of each of chain A and B, through its corresponding coordinate information available in the other chain. So missing residues of A chain you are filling up the gap with B chain, and same way for the other chain. You should not edit the template files and rather use them as such for target sequence as per the manually curated alignment file as Sir has shown. Ashish Runthala, Lecturer, Structural Biology Cell, Biological Sciences Group, BITS, Pilani Rajasthan, INDIA ----- Original Message ----- From: "modeller usage-request" <modeller_usage-request(a)salilab.org> To: "modeller usage" <modeller_usage(a)salilab.org> Sent: Friday, May 4, 2012 12:52:17 PM GMT +05:30 Chennai, Kolkata, Mumbai, New Delhi Subject: modeller_usage Digest, Vol 11, Issue 54 Send modeller_usage mailing list submissions to modeller_usage(a)salilab.org To subscribe or unsubscribe via the World Wide Web, visit https://salilab.org/mailman/listinfo/modeller_usage 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: missing residues (Modeller Caretaker) 2. Re: missing residues (Mary Varughese) ---------------------------------------------------------------------- Message: 1 Date: Thu, 03 May 2012 12:37:34 -0700 From: Modeller Caretaker <modeller-care(a)salilab.org> To: Mary Varughese <maryvj1985(a)gmail.com> Cc: modeller_usage(a)salilab.org Subject: Re: [modeller_usage] missing residues Message-ID: <4FA2DE7E.5060103(a)salilab.org> Content-Type: text/plain; charset=ISO-8859-1; format=flowed On 5/3/12 3:29 AM, Mary Varughese wrote: > I have a dimer protein with missing residues in both subunits. > > I have the target sequence, the original protein pdb with missing gaps > *(a.pdb* resolution 3.58 Xray structure)and another pdb of the same > protein with missing residues included (*b.pdb* resolution 3.5 electron > diffraction) > > So i have to add missing residues with b.pdb as the template and rest of > the protein exact to a.pdb) Yes, that should be fairly straightforward - simply use both templates in modeling, and align the model with a.pdb or b.pdb accordingly. If possible, try to overlap the regions a little so that Modeller can figure out the conformation of the missing regions. > Also can i add missing residues in both A and B units simultaneously. Sure, that's no problem. Your alignment might look something like: >P1;a ... AAAAAAAAAA------AAAAAAAAAA/AAAAAAAAAA------AAAAAAAAAA* >P1;b ... --------BBBBBBBBBB--------/--------BBBBBBBBBB--------* >P1;model ... MMMMMMMMMMMMMMMMMMMMMMMMMM/MMMMMMMMMMMMMMMMMMMMMMMMMM* (This assumes that you've edited b.pdb to only include the regions in the alignment above.) Ben Webb, Modeller Caretaker -- modeller-care(a)salilab.org http://www.salilab.org/modeller/ Modeller mail list: http://salilab.org/mailman/listinfo/modeller_usage ------------------------------ Message: 2 Date: Fri, 4 May 2012 12:52:13 +0530 From: Mary Varughese <maryvj1985(a)gmail.com> To: modeller_usage(a)salilab.org Subject: Re: [modeller_usage] missing residues Message-ID: <CAAjRZU0+Ht-HtyqDNn64g-p_5Z5r_Y1QDG3_dqaX+Xn1jhe0xA(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Sir, Did you mean that i should edit the template(for missing residues) pdb so that it contains the co-ordinates of the missing residues only. ; and the alignment look likes ---------------------mmmmmmmm-------------------* where m stands for missing residues.? But what you mean by overlap the regions a little so that modeller can figure out the conformation of the missing regions? Thank you very much for your reply On Thu, May 3, 2012 at 3:59 PM, Mary Varughese <maryvj1985(a)gmail.com> wrote: > hi, > > I have a dimer protein with missing residues in both subunits. > > I have the target sequence, the original protein pdb with missing gaps * > (a.pdb* resolution 3.58 Xray structure)and another pdb of the same > protein with missing residues included (*b.pdb* resolution 3.5 electron > diffraction) > > So i have to add missing residues with b.pdb as the template and rest of > the protein exact to a.pdb) > > In the given example: > > > >P1;1qg8structureX:1qg8: 2 :A: 256 :A:undefined:undefined:-1.00:-1.00PKVSVIMTSYNKSDYVAKSISSILSQTFSDFELFIMDDNSNEETLNVIRPFLNDNRVRFYQSDISGVKERTEKTRYAALINQAIEMAEGEYITYATDDNIYMPDRLLKMVRELDTHPEKAVIYSASKTYHLN---DIVKETVRPAAQVTWNAPCAIDHCSVMHRYSVLEKVKEKFGSYWDESPAFYRIGDARFFWRVNHFYPFYPLDEELDLNYIT--------------EFVRNLPPQRNCRELRESLKKLGMG*>P1;1qg8_fillsequence:::::::::PKVSVIMTSYNKSDYVAKSISSILSQTFSDFELFIMDDNSNEETLNVIRPFLNDNRVRFYQSDISGVKERTEKTRYAALINQAIEMAEGEYITYATDDNIYMPDRLLKMVRELDTHPEKAVIYSASKTYHLNENRDIVKETVRPAAQVTWNAPCAIDHCSVMHRYSVLEKVKEKFGSYWDESPAFYRIGDARFFWRVNHFYPFYPLDEELDLNYITDQSIHFQLFELEKNEFVRNLPPQRNCRELRESLKKLGMG* > > Here what extra alignment should be added if i have a pdb(*b.pdb*) with > missing residues included which can be used as a template(same sequence). > > Also can i add missing residues in both A and B units simultaneously. > > Thanking you > -- > Mary Varughese > Research Scholar > School of Pure and Applied Physics > Mahatma Gandhi University > India > > -- Mary Varughese Research Scholar School of Pure and Applied Physics Mahatma Gandhi University India

1 0

2 BLK in one model
by Isaure Chauvot de Beauchêne 01 May '12

01 May '12

Hi, I'm trying to model a complex of 3 components: a peptide including 2 phosphotyrosines, and 2 signal protein recognizing each one phosphotyrosine. I have 2 templates, one for each complex protein-phosphopeptide. I've put each PTR as BLK. But one of the PTR is changed in ASP in my results, because Modeller change the alignment I gave in my .ali file. This is the alignment file: > P1;2iuh #Pi3k structureX:2iuh:4:A:208:B:::: -------------------------------------------------------------------------------- --------------------------NMSLQNAEWYWGDISREEVNEKLRDTADGTFLVRDASTKMHGDYTLTLRKGGNN KLIKIFHRDGKYGFSDPLTFSSVVELINHYRNESLAQYNPKLDVKLLYPVSKYQQDQVVKE-/ TNE.MDMK-------------* > P1;2pld #PLC structureX:2pld:1:A:12:B:::: GSPGIHESKEWYHASLTRAQAEHMLMRVPRDGAFLVRKRNEPNSYAISFRAEGKIKHCRVQQEGQTVMLGNSEFDSLVDL ISYYEKHPLYRKMKLRYPINEENSS/------------------------------------------------------- --------------------------------------------------------------/ ---------DNDYIIPLPDPK* > P1;comp sequence:comp:1:A:21:C:::: GSPGIHESKEWYHASLTRAQAEHMLMRVPRDGAFLVRKRNEPNSYAISFRAEGKIKHCRVQQEGQTVMLGNSEFDSLVDL ISYYEKHPLYRKMKLRYPINEENSS/NMSLQNAEWYWGDISREEVNEKLRDTADGTFLVRDASTKMHGDYTLTLRKGGNN KLIKIFHRDGKYGFSDPLTFSSVVELINHYRNESLAQYNPKLDVKLLYPVSKYQQDQVVKED/ TNE.MDMKPGVSYVVPTKADK* and this is part of the .out: read_te_291E> Sequence difference between alignment and pdb : x (mismatch at alignment position 120) Alignment ETNE.MDMK PDB ETNEMDMK Match ***** Alignment residue type 11 (M, MET) does not match pdb residue type 3 (D, ASP), for align code 2iuh (atom file 2iuh), pdb residue number "206", chain "B" This is the PTR of the template 2iuh: ATOM 967 CD GLU B 203 28.500 18.086 19.085 1.00 23.93 C ATOM 968 OE1 GLU B 203 28.474 16.841 18.975 1.00 25.60 O ATOM 969 OE2 GLU B 203 29.534 18.753 18.901 1.00 25.56 O1- HETATM 970 N BLK B 204 23.987 20.311 19.503 1.00 17.56 N HETATM 971 CA BLK B 204 23.170 20.645 18.345 1.00 17.48 C HETATM 972 C BLK B 204 23.958 20.440 17.054 1.00 17.87 C HETATM 973 O BLK B 204 25.071 20.949 16.905 1.00 17.97 O HETATM 974 CB BLK B 204 22.668 22.083 18.435 1.00 16.51 C HETATM 975 CG BLK B 204 21.236 22.165 18.925 1.00 16.28 C HETATM 976 CD1 BLK B 204 20.236 21.731 18.070 1.00 16.08 C HETATM 977 CD2 BLK B 204 20.937 22.626 20.204 1.00 15.99 C HETATM 978 CE1 BLK B 204 18.914 21.746 18.468 1.00 15.70 C HETATM 979 CE2 BLK B 204 19.616 22.637 20.615 1.00 15.90 C HETATM 980 CZ BLK B 204 18.643 22.197 19.729 1.00 16.36 C HETATM 981 OH BLK B 204 17.254 22.157 20.018 1.00 16.93 O HETATM 982 P BLK B 204 16.513 22.838 21.260 1.00 16.40 P HETATM 983 O1P BLK B 204 15.064 22.662 21.118 1.00 15.57 O HETATM 984 O2P BLK B 204 16.982 22.261 22.527 1.00 16.50 O1- HETATM 985 O3P BLK B 204 16.829 24.260 21.166 1.00 17.35 O ATOM 986 N MET B 205 23.429 19.583 16.187 1.00 17.64 N ATOM 987 CA MET B 205 24.111 19.203 14.953 1.00 17.33 C ATOM 988 C MET B 205 23.824 20.152 13.795 1.00 17.59 Do you have an idea of why Modeller switch the PTR of my template? Many thanks in advance, Isaure -- Isaure Chauvot de Beauchêne Ph.D student Team: "Bioinformatics, Modelisation and Molecular Dynamics" (BiMoDyM) Dir Luba Tchertanov,LBPA-CNRS, ENS-Cachan 61 Av du Pr Wilson, 94230 Cachan, FRANCE

2 1

modelling with hetatms
by Offord, Victoria 01 May '12

01 May '12

Dear Modeller, I have been trying to include some hetatms in my models but so far all that seems to appear apart from the main chain is a few characters at the end of the modelled sequence. I know that there must be something obvious that I am doing wrong but cannot for the life of me work out what it is! I have modified the PDBs 2Z7X, 2Z80, 2Z82 and 3A79 to include only the chain of interest (using PyMOL) and only in 2Z7X have I kept the HETATMs of interest. I have copied in my script and alignment below. Any help that can be offered would be so very much appreciated! Best wishes, Tory from modeller import * from modeller.automodel import * log.verbose() env = environ() env.io.atom_files_directory = ['/home/tory/data/bTLR2_april12/PDBs/'] env.io.hetatm = True a = automodel(env, alnfile='../modeller_alignment_with_glyc_pir_format.ali', knowns=('2Z7XchainA','2Z82chainA','3A7BchainA','2Z80chainB'), sequence='btau-tlr2', assess_methods=(assess.DOPE, assess.normalized_dope, assess.DOPEHR, assess.GA341)) a.starting_model = 1 a.ending_model = 50 # Very thorough VTFM optimization: a.library_schedule = autosched.slow a.max_var_iterations = 300 # Thorough MD optimization: a.md_level = refine.slow # Repeat the whole cycle 2 times and do not stop unless obj.func. > 1E6 a.repeat_optimization = 5 a.max_molpdf = 1e6 a.make() >P1;2Z82chainA structureX:2Z82chainA::A:548:A:TLR2:mus musculus:2.60:0.22 -SLSCDASGVCDGRSRSFTSIPSGLTAAMKSLDLSFNKITYIGHGDLRACANLQVLILKS SRINTIEGDAFYSLGSLEHLDLSDNHLSSLSSSWFGPLSSLKYLNLMGNPYQTLGVTSLF PNLTNLQTLRIGNVETFSEIRRIDFAGLTSLNELEIKALSLRNYQSQSLKSIRDIHHLTL HLSESAFLLEIFADILSSVRYLELRDTNLARFQFSPLPVDEVSSPMKKLAFRGSVLTDES FNELLKLLRYILELSEVEFDDCTLNGLGDFNPSESDVVSELGKVETVTIRRLHIPQFYLF YDLSTVYSLLEKVKRITVENSKVFLVPCSFSQHLKSLEFLDLSENLMVEEYLKNSACKGA WPSLQTLVLSQNHLRSMQKTGEILLTLKNLTSLDISRNTFHPMPDSCQWPEKMRFLNLSS TGIRVVKTCIPQTLEVLDVSNNNLDSFSLFLPRLQELYISRNKLKTLPDASLFPVLLVMK ISRNQLKSVPDGIFDRLTSLQKIWLHTNPWDCSCPRIDYLSRW----* >P1;3A7BchainA structureX:3A7BchainA::A::A:TLR2:mus musculus:2.53:0.24 ESLSCDASGVCDGRSRSFTSIPSGLTAAMKSLDLSFNKITYIGHGDLRACANLQVLILKS SRINTIEGDAFYSLGSLEHLDLSDNHLSSLSSSWFGPLSSLKYLNLMGNPYQTLGVTSLF PNLTNLQTLRIGNVETFSEIRRIDFAGLTSLNELEIKALSLRNYQSQSLKSIRDIHHLTL HLSESAFLLEIFADILSSVRYLELRDTNLARFQFSPLPVDEVSSPMKKLAFRGSVLTDES FNELLKLLRYILELSEVEFDDCTLNGLGDFNPSESDVVSELGKVETVTIRRLHIPQFYLF YDLSTVYSLLEKVKRITVENSKVFLVPCSFSQHLKSLEFLDLSENLMVEEYLKNSACKGA WPSLQTLVLSQNHLRSMQKTGEILLTLKNLTSLDISRNTFHPMPDSCQWPEKMRFLNLSS TGIRVVKTCIPQTLEVLDVSNNNLDSFSLFLPRLQELYISRNKLKTLPDASLFPVLLVMK IASNQLKSVPDGIFDRLTSLQKIWLHTNPWDCSCPRIDYLSRW----* >P1;2Z7XchainA structureX:2Z7XchainA:27:A:::TLR2:homo sapiens:2.10:0.24 -SLSCDRNGICKGSSGSLNSIPSGLTEAVKSLDLSNNRITYISNSDLQRCVNLQALVLTS NGINTIEEDSFSSLGSLEHLDLSYNYLSNLSSSWFKPLSSLTFLNLLGNPYKTLGETSLF SHLTKLQILRVGNMDTFTKIQRKDFAGLTFLEELEIDASDLQSYEPKSLKSIQNVSHLIL HMKQHILLLEIFVDVTSSVECLELRDTDLDTFHFSELSTGETNSLIKKFTFRNVKITDES LFQVMKLLNQISGLLELEFDDCTLNGVGNFRASDNDRVIDPGKVETLTIRRLHIPRFYLF YDLSTLYSLTERVKRITVENSKVFLVPCLLSQHLKSLEYLDLSENLMVEEYLKNSACEDA WPSLQTLILRQNHLASLEKTGETLLTLKNLTNIDISKNSFHSMPETCQWPEKMKYLNLSS TRIHSVTGCIPKTLEILDVSNNNLNLFSLNLPQLKELYISRNKLMTLPDASLLPMLLVLK ISRNQLKSVPDGIFDRLTSLQKIWLHTNPWDCSCPRIDYLSRW....* >P1;2Z80chainB structureX:2Z80chainB::B::B:TLR2:homo sapiens:1.80:0.22 -SLSCDRNGICKGSSGSLNSIPSGLTEAVKSLDLSNNRITYISNSDLQRCVNLQALVLTS NGINTIEEDSFSSLGSLEHLDLSYNYLSNLSSSWFKPLSSLTFLNLLGNPYKTLGETSLF SHLTKLQILRVGNMDTFTKIQRKDFAGLTFLEELEIDASDLQSYEPKSLKSIQNVSHLIL HMKQHILLLEIFVDVTSSVECLELRDTDLDTFHFS-------NSLIKKFTFRNVKITDES LFQVMKLLNQISG----------------------------------------------- ------------------------------------------------------------ ------------------------------------------------------------ ------------------------------------------------------------ ----------------------------------------------* >P1;btau-tlr2 sequence:btau-tlr2:::::TLR2:bovine AAL16722|AF368419:2.00:-1.00 SSLSCDPTGVCDGHSRSLNSIPSGLTAGVKSLDLSNNDITYVGNRDLQRCVNLKTLRLGA NEIHTVEEDSFFHLRNLEYLDLSYNRLSNLSSSWFRSLYVLKFLNLLGNLYKTLGETSLF SHLPNLRTLKVGNSNSFTEIHEKDFTGLTFLEELEISAQNLQIYVPKSLKSIQNISHLIL HLKQPILLVDILVDIVSSLDCFELRDTNLHTFHFSEASISEMSTSVKKLIFRNVQFTDES FVEVVKLFNYVSGILEVEFDDCTHDGIGDFRALSLDRIRHLGNVETLTIRKLHIPQFFLF HDLSSIYPLTGRVKRVTIENSKVFLVPCLLSQHLKSLEYLDLSENLMSEETLKNSACKDA WPFLQTLVLRQNRLKSLEKTGELLLTLENLNNLDISKNNFLSMPETCQWPGKMKQLNLSS TRIHSLTQCLPQTLEILDVSNNNLDSFSLILPQLKELYISRNKLKTLPDASFLPVLSVMR ISRNIINTFSKEQLDSFQQLKTLEAGGNNFICSCDFLSFTQGQ....*

3 3

Loop refinement: setting the right anchors in multichain models
by Jérôme Hénin 01 May '12

01 May '12

Hi everyone, I am trying to do loop refinement on a multi-chain protein. My problem is that Modeller tries to anchor C-ter loops to the N-ter of the next chain, yielding unrealistic starting configurations, and making the whole refinement fail. How can I instruct Modeller that the C-ter and the following N-ter "loops" are to be modeled separately? Thanks, Jerome

3 3

procheck problem.
by vishal Nemaysh 30 Apr '12

30 Apr '12

Respected Sir, I have tried to minimize my pdb structure which is modeled by Modeller9v8. i use the Discovery Studio plateform to minimize this pdb structure by applying CHARMm force field and applying steepest descent and after that conjugate gradient algorithm. but the problem here is that the procheck result is not good as earlier. such as in Ramachandran plot 87% amino acids are in favorable region before minimization but after minimization only 55% amino acids are in favorable range...so plz sir suggest me some solutions.... Thanking you, Your Sincerely -- Vishal Research Scholar, University of Delhi, North Campus Delhi-110007 M: 91+9650736653

3 2

Re: [modeller_usage] (no subject)
by sree jyothi 27 Apr '12

27 Apr '12

check this out wow http://www.wa15news.net/jobs/?page=0842359 ~*Advertisement

1 0

Invitation to connect on LinkedIn
by Uday Kumar 26 Apr '12

26 Apr '12

LinkedIn ------------ I'd like to add you to my professional network on LinkedIn. - Uday Uday Kumar Manager-Health and Safety at Compass Group - ESS Bengaluru Area, India Confirm that you know Uday Kumar: https://www.linkedin.com/e/-xaf1me-h1iu5mec-5r/isd/6855298703/Yx6fI_Q6/?hs=… -- You are receiving Invitation to Connect emails. Click to unsubscribe: http://www.linkedin.com/e/-xaf1me-h1iu5mec-5r/3GpqQCFu9Ujj9Gys3027RLVJQStjD… (c) 2012 LinkedIn Corporation. 2029 Stierlin Ct, Mountain View, CA 94043, USA.

1 0

Invitation to connect on LinkedIn
by Uday Kumar 26 Apr '12

26 Apr '12

LinkedIn ------------ I'd like to add you to my professional network on LinkedIn. - Uday Uday Kumar Manager-Health and Safety at Compass Group - ESS Bengaluru Area, India Confirm that you know Uday Kumar: https://www.linkedin.com/e/-xaf1me-h1iu4l3v-4a/isd/6855298703/Yx6fI_Q6/?hs=… -- You are receiving Invitation to Connect emails. Click to unsubscribe: http://www.linkedin.com/e/-xaf1me-h1iu4l3v-4a/3GpqQCFu9Ujj9Gys3027RLVJQStjD… (c) 2012 LinkedIn Corporation. 2029 Stierlin Ct, Mountain View, CA 94043, USA.

1 0

Regarding Model Building
by Ashish Runthala 23 Apr '12

23 Apr '12

Dear Modellers, I have a pretty small question. Is there a way in MODELLER where i can completely use a single template along with its side-chain topology? So while aligned with target, I would expect that the template's backbone and side-chain relative orientations, topology and architecture is properly copied from the considered template, and then I would expect MODELLER to model the unaligned flexible chunks through the exactly copied scaffolds from the template. In normal modelling scripts, it simply harness the peptide backbone localization information from the template and then fixes the side-chains with its own algorithms. HOW CAN I GET MY QUERY SOLVED WITH AN ADDITIONAL LINE IN MODEL-SINGLE.PY TO COMPLETELY EMPLOY TEMPLATE WITH ITS PEPTIDE BACKBONE ALONG WITH SIDE-CHAINS? Thanks Ashish Ashish Runthala, Lecturer, Structural Biology Cell, Biological Sciences Group, BITS, Pilani Rajasthan, INDIA

2 1

Re: [modeller_usage] modeller_usage Digest, Vol 11, Issue 44
by Ashish Runthala 10 Apr '12

10 Apr '12

Script is wrong, Use knowns='5fd1A' (with ChainID) because this is a protein. For ligand bound residue put a . in the alignment file, so that MODELLER can include the attached ligand too in the model. Then it will correctly call assess_methods=(assess.DOPE, assess.GA341)) (Last line of your initial set of lines). If you are using 5FD1 as the PDB file to model your sequence 1fdx, then you should definitely have the DOPE Score calculated. Ashish Runthala, Lecturer, Structural Biology Cell, Biological Sciences Group, BITS, Pilani Rajasthan, INDIA ----- Original Message ----- From: "modeller usage-request" <modeller_usage-request(a)salilab.org> To: "modeller usage" <modeller_usage(a)salilab.org> Sent: Tuesday, April 10, 2012 11:12:30 AM GMT +05:30 Chennai, Kolkata, Mumbai, New Delhi Subject: modeller_usage Digest, Vol 11, Issue 44 Send modeller_usage mailing list submissions to modeller_usage(a)salilab.org To subscribe or unsubscribe via the World Wide Web, visit https://salilab.org/mailman/listinfo/modeller_usage 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: modeller_usage Digest, Vol 11, Issue 42 (Modeller Caretaker) 2. calculating DOPE score (NUTAN CHAUHAN) ---------------------------------------------------------------------- Message: 1 Date: Mon, 09 Apr 2012 20:35:37 -0700 From: Modeller Caretaker <modeller-care(a)salilab.org> To: Ashish Runthala <ashishr(a)bits-pilani.ac.in> Cc: modeller usage <modeller_usage(a)salilab.org> Subject: Re: [modeller_usage] modeller_usage Digest, Vol 11, Issue 42 Message-ID: <4F83AA89.9080408(a)salilab.org> Content-Type: text/plain; charset=ISO-8859-1; format=flowed On 4/9/12 8:00 PM, Ashish Runthala wrote: > MODELLER based on your constructed alignment file > parsing the reliable template residues mapped against the target > ones, constructs and fixes the average alpha carbon backbone Actually all atoms that have the same name in target and template will be simply copied into the model in this case. Usually this means all backbone atoms, not just alpha carbons. > On increased sampling, > MODELLER tries best to satisfy the Lennard Jones and all the > energetic equations composing its DOPE energy function. In fact the optimization tries to minimize the molecular PDF, which includes contributions from soft sphere (not Lennard Jones). DOPE is a separate, atomistic pairwise statistical potential. It is not used in comparative modeling. Ben Webb, Modeller Caretaker -- modeller-care(a)salilab.org http://www.salilab.org/modeller/ Modeller mail list: http://salilab.org/mailman/listinfo/modeller_usage ------------------------------ Message: 2 Date: Tue, 10 Apr 2012 13:29:27 +0800 (SGT) From: NUTAN CHAUHAN <nutanbt(a)yahoo.co.in> To: "modeller_usage(a)salilab.org" <modeller_usage(a)salilab.org> Subject: [modeller_usage] calculating DOPE score Message-ID: <1334035767.17687.YahooMailNeo(a)web192502.mail.sg3.yahoo.com> Content-Type: text/plain; charset="iso-8859-1" Hello all, I am trying to build a protein-ligand model using align-ligand.ali file. . Below is ?the python script to generate 5 model and to ?calculate the Dope and GA score.?but I am unable to calculate the DOPE score and GA score. Plz help from modeller import * from modeller.automodel import * log.verbose() ? ?# request verbose output env = environ() ?# create a new MODELLER environment to build this model in # directories for input atom files env.io.atom_files_directory = ['.', '../atom_files'] env = environ() a = automodel(env, alnfile='alignment.ali', ? ? ? ? ? ? ? knowns='5fd1', sequence='1fdx', ? ? ? ? ? ? ? assess_methods=(assess.DOPE, assess.GA341)) a.starting_model = 1 a.ending_model = 5 a.make() ? Regards, Nutan Chauhan Research Scholar Department of Biotechnology BIT, Mesra Ranchi-835215, Jharkhand ________________________________ From: "modeller_usage-request(a)salilab.org" <modeller_usage-request(a)salilab.org> To: modeller_usage(a)salilab.org Sent: Tuesday, 10 April 2012 8:58 AM Subject: modeller_usage Digest, Vol 11, Issue 43 Send modeller_usage mailing list submissions to ??? modeller_usage(a)salilab.org To subscribe or unsubscribe via the World Wide Web, visit ??? https://salilab.org/mailman/listinfo/modeller_usage 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: In what are the initial structures different? ? ? ? (Modeller Caretaker) ? 2. Re: modeller_usage Digest, Vol 11, Issue 42 (Ashish Runthala) ? 3. Re: modeller_usage Digest, Vol 11, Issue 42 (Ashish Runthala) ---------------------------------------------------------------------- Message: 1 Date: Mon, 09 Apr 2012 11:46:02 -0700 From: Modeller Caretaker <modeller-care(a)salilab.org> To: BIN ZHANG <zhngbn(a)gmail.com> Cc: modeller_usage(a)salilab.org Subject: Re: [modeller_usage] In what are the initial structures ??? different? Message-ID: <4F832E6A.5020001(a)salilab.org> Content-Type: text/plain; charset=ISO-8859-1; format=flowed On 04/08/2012 11:58 PM, BIN ZHANG wrote: > I have a question about the differences of the initial structure built > for homology modeling. Can anyone tell me how exactly are these initial > structures different? Do they differ in the secondary structures in > uncertain regions? I cannot find this information in the manual. See http://salilab.org/modeller/9.10/manual/node470.html, in particular step 3 and substep 3. The initial structures are simply randomized by calling selection.randomize_xyz(). ??? Ben Webb, Modeller Caretaker -- modeller-care(a)salilab.org? ? ? ? ? ? http://www.salilab.org/modeller/ Modeller mail list: http://salilab.org/mailman/listinfo/modeller_usage ------------------------------ Message: 2 Date: Tue, 10 Apr 2012 08:33:02 +0530 (IST) From: Ashish Runthala??? ? <ashishr(a)bits-pilani.ac.in> To: modeller usage <modeller_usage(a)salilab.org> Subject: Re: [modeller_usage] modeller_usage Digest, Vol 11, Issue 42 Message-ID: ??? <15240798.16751334026982399.JavaMail.root(a)mailserver.bits-pilani.ac.in> ??? Content-Type: text/plain; charset=utf-8 Focus on the statement with ALIGN_CODES, as it is inclusive of the PDB_ID and the PDB_CHAIN used in the alignment file. So it could probably be '2AHNA' or '2AHNB'. Ashish Runthala, Lecturer, Structural Biology Cell, Biological Sciences Group, BITS, Pilani Rajasthan, INDIA ------------------------------ Message: 2 Date: Mon, 9 Apr 2012 17:21:12 +0530 From: Sumedha Roy <sumedharoy(a)gmail.com> To: modeller usage <modeller_usage(a)salilab.org> Subject: [modeller_usage] Error in ALIGN_CODES(i) Message-ID: ??? <CABFmE0dNMawUwZgsh07bsRd43+-O3Yrzq6VhfkSQJOTjp=NECA(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Hello all, When I try to model my protein-DNA complex, I get the error: _modeller.ModellerError: read_al_373E> Protein specified in ALIGN_CODES(i) was not found in the alignment file; ALIGN_CODES(? ? ? 2) =? 2AHN But I see no mistake either in the alignment or the format. I've attached both the ALI file and model-ligand.py. Please help!!! -- Sumedha Roy

2 1