I have a doubt in plotting the graph for calculating the DOPE Score using GNUPLOT. In GNUPLOT using the command 'plot "TvLDH.profile" using 1:42 with lines'. (From the Tutorial ) Only a single graph is shown but in the tutorial 2 lines of 2 profiles are shown.(TVLDH.profile showing red & green )
My question is how dow we compare the DOPE Score of 2 or more profiles using GNU plot. Please give me a posetive reply
Regards Salam Pradeep Bioinformatics Center North-East Institute of Science & Technology
Jorhat, Assam, India
When replying, please edit your Subject line so it is more specific than "Re: Contents of modeller_usage digest..."
Today's Topics:
1. Rino Ragno ? fuori ufficio/ is out of office
(">) 2. Rino Ragno ? fuori ufficio/ is out of office (">) 3. Re: problem with ligand modelling (Modeller Caretaker)
4. Re: modeller_usage Digest, Vol 7, Issue 54 (sanjay singh)
Message: 3 Date: Mon, 18 Aug 2008 09:21:51 -0700 From: Modeller Caretaker <">> Subject: Re: [modeller_usage] problem with ligand modelling
To: sanjay singh <">> Cc: "> Message-ID: <">>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
sanjay singh wrote: ... > my alignment looks like this > > >P1;1iax > structureX:1iax:11 :A :438 :B :ferredoxin:Azotobacter vinelandii:
> 1.90: 0.19 > ILSKLATNE-----SPYFDGWKAYDSDPFHPLKNPNGVIQMGLAENQLCLDLIEDWIKRN // > NAGLFCWMDLRPLLRESTFDSEMSLWRVIINDVKLNVSPGSSFECQEPGWFRVCFANMDD > GTVDIALARIRRFVGVEK/..* ... > read_te_290E> Number of residues in the alignment and pdb files are
> different: 836 834 > For alignment entry: 1 1iax > x (mismatch at alignment position 835) > Alignment VDIALARIRRFVGVEK.. > PDB VDIALARIRRFVGVEK
> Match **************** > > Please check your alignment file header to be sure you correctly specified > the starting and ending residue numbers and chains. The alignment sequence > must match that from the atom file exactly.
You should follow the advice in the error message, and check your alignment file header. Your header tells Modeller to read residues 11:A through 438:B from your PDB file, so it's not going to read the HETATMs,
which come after 438:B. Fix your alignment header.
BTW, your PDB file contains 4 HETAMs (two SO4 and two PLP) but your alignment sequence contains only two.
? hi all, i got success to build a 2chain model with two ligand.I am sincerely thankful to Ben Webb,modeller care taker for all this thing .but now I have another question :P...last portion of my final PDB look like this
ATOM 6730 C ASN B 852 5.358 41.203 4.902 1.00 44.07 2SG6732 ATOM 6731 O ASN B 852 5.201 40.108 4.300 1.00 44.07 2SG6733 ATOM 6732 OXT ASN B 852 4.554 41.656 5.760 1.00 44.07 2SG6734
TER 6732 ASN 852 2SG6735 HETATM 6733 S SO4 C 853 2.345 45.702 73.000 0.40 17.18 3SG6736 HETATM 6734 O1 SO4 C 853 1.549 46.638 73.796 1.00 17.18 3SG6737
HETATM 6735 O2 SO4 C 853 3.775 45.879 73.250 1.00 17.18 3SG6738 HETATM 6736 O3 SO4 C 853 2.020 44.340 73.405 1.00 17.18 3SG6739 HETATM 6737 O4 SO4 C 853 2.051 45.892 71.587 1.00 17.18 3SG6740
HETATM 6738 S SO4 C 854 1.214 44.486 35.947 0.40 11.34 3SG6741
here you can see that modeller take HEATATMs as chain C.then my question is whether ir hamper the model or not.will it acceptable???
regards sanjay
On Fri, 15 Aug 2008 "> wrote : >Send modeller_usage mailing list submissions to > ">
> >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
> "> > >You can reach the person managing the list at > ">
> >When replying, please edit your Subject line so it is more specific >than "Re: Contents of modeller_usage digest..." > > >Today's Topics: > > 1. Re: using modeller to refine a docking conformation
> (Modeller Caretaker) > 2. problem with model-multiple-hetero.py (sanjay singh) > > >---------------------------------------------------------------------- > >Message: 1
>Date: Thu, 14 Aug 2008 11:27:49 -0700 > From: Modeller Caretaker <">> >Subject: Re: [modeller_usage] using modeller to refine a docking
> conformation >To: Jake Gunn-Glanville <">> >Cc: "> >Message-ID: <">>
>Content-Type: text/plain; charset=ISO-8859-1; format=flowed > >Jake Gunn-Glanville wrote: > > Does modeller provide a means of generating models that differ only > > slightly from an input template?
> >Sure - you could build some models using a simple 1:1 alignment to the >template. The final models will look very similar to the input (by >construction) but you should expect to see some variability,
>particularly in the regions that are not strongly restrained (e.g. >sidechains). You could add in some extra restraints around the docked >region, if you have some extra data or intuition about the interactions.
>But Modeller is designed to generate good solutions to your scoring >function, not to sample from an ensemble - if you want to do that you >might have more luck with an MD package. > > Ben Webb, Modeller Caretaker
>-- >">http://www.salilab.org/modeller/ >Modeller mail list: http://salilab.org/mailman/listinfo/modeller_usage
> > >------------------------------ > >Message: 2 >Date: 15 Aug 2008 08:03:14 -0000 > From: "sanjay singh" <">>
>Subject: [modeller_usage] problem with model-multiple-hetero.py >To: "> >Message-ID: <">>
>Content-Type: text/plain; charset="iso-8859-1" > >My aim is to model single chain protein from a 2 chain template with 2 ligands ( PLP).following the advance tutorial I have edited template PDB ( i.e.1IAX) as per my requirement.I have deleted chain B in the ligand which will interact with it.now template is look like this-
> >ATOM 2020 N SER A 277 6.836 50.703 64.530 1.00 50.02 N >ATOM 2021 CA SER A 277 5.736 50.683 65.488 1.00 52.57 C >ATOM 2022 C SER A 277 6.158 50.953 66.920 1.00 53.76 C
>ATOM 2023 O SER A 277 5.656 51.877 67.558 1.00 53.78 O >ATOM 2024 CB SER A 277 5.023 49.340 65.432 1.00 53.52 C >ATOM 2025 OG SER A 277 4.399 49.181 64.168 1.00 60.54 O
>ATOM 2026 N LYS A 278 7.083 50.145 67.425 1.00 55.18 N >ATOM 2027 CA LYS A 278 7.540 50.304 68.795 1.00 56.49 C >ATOM 2028 C LYS A 278 8.430 51.529 69.009 1.00 56.08 C
>ATOM 2029 O LYS A 278 8.170 52.321 69.915 1.00 56.54 O >ATOM 2030 CB LYS A 278 8.218 49.013 69.274 1.00 59.46 C >ATOM 2031 CG LYS A 278 7.268 47.807 69.271 1.00 63.38 C
>ATOM 2032 CD LYS A 278 6.048 48.015 70.191 1.00 69.44 C >ATOM 2033 CE LYS A 278 4.984 46.952 69.962 1.00 73.69 C >ATOM 2034 NZ LYS A 278 5.598 45.765 69.377 1.00 82.22 N
>TER 2035 LYS A 278 >HETATM 6623 S SO4 600 2.001 43.793 71.949 0.40 47.42 S >HETATM 6624 O1 SO4 600 1.250 44.661 72.856 1.00 49.68 O >HETATM 6625 O2 SO4 600 3.441 43.960 72.132 1.00 48.63 O
>HETATM 6626 O3 SO4 600 1.702 42.399 72.263 1.00 51.95 O >HETATM 6627 O4 SO4 600 1.623 44.087 70.571 1.00 47.72 O >HETATM 6633 N1 PLP A 500 6.942 41.456 67.024 1.00 91.96 N
>HETATM 6634 C2 PLP A 500 7.237 41.844 68.301 1.00 92.40 C >HETATM 6635 C2A PLP A 500 8.277 41.064 69.104 1.00 91.84 C >HETATM 6636 C3 PLP A 500 6.579 42.967 68.880 1.00 91.52 C
>HETATM 6637 O3 PLP A 500 6.876 43.351 70.176 1.00 90.45 O >HETATM 6638 C4 PLP A 500 5.613 43.677 68.118 1.00 90.22 C >HETATM 6639 C4A PLP A 500 4.901 44.884 68.713 1.00 86.77 C
>HETATM 6640 C5 PLP A 500 5.327 43.248 66.794 1.00 91.55 C >HETATM 6641 C6 PLP A 500 6.014 42.124 66.269 1.00 92.40 C >HETATM 6642 C5A PLP A 500 4.287 43.977 65.936 1.00 92.05 C
>HETATM 6643 O4P PLP A 500 4.919 44.983 65.137 1.00 91.48 O >HETATM 6644 P PLP A 500 3.850 45.746 64.212 0.40 89.98 P >HETATM 6645 O1P PLP A 500 2.942 46.559 65.059 1.00 90.28 O
>HETATM 6646 O2P PLP A 500 3.052 44.763 63.428 1.00 89.69 O >HETATM 6647 O3P PLP A 500 4.581 46.636 63.284 1.00 90.10 O >HETATM 6663 O HOH 1 -2.199 65.598 27.605 1.00 37.63 O
>HETATM 6664 O HOH 2 13.146 43.997 69.980 1.00 28.90 O >HETATM 6665 O HOH 3 3.514 34.147 84.461 1.00 40.65 O >HETATM 6666 O HOH 4 14.245 36.668 40.667 1.00 33.18 O
>HETATM 6667 O HOH 5 6.748 58.168 71.500 1.00 37.09 O >HETATM 6668 O HOH 6 -3.541 20.765 51.638 1.00 34.79 O >HETATM 6669 O HOH 7 -4.102 61.235 21.691 1.00 38.94 O
>HETATM 6670 O HOH 8 10.326 70.845 65.572 1.00 35.43 O >HETATM 6671 O HOH 9 -4.995 51.193 45.379 1.00 38.86 O >HETATM 6672 O HOH 10 -6.557 33.928 65.811 1.00 39.72 O
>HETATM 6673 O HOH 11 -1.108 49.863 77.220 1.00 41.22 O >HETATM 6674 O HOH 12 -26.403 30.845 53.499 1.00 43.93 O >HETATM 6675 O HOH 13 -6.662 49.811 60.632 1.00 20.93 O
>HETATM 6676 O HOH 14 0.524 46.344 51.157 1.00 30.77 O >HETATM 6677 O HOH 15 4.444 27.533 53.723 1.00 45.47 O >HETATM 6678 O HOH 16 -9.377 50.142 40.921 1.00 38.67 O
>HETATM 6679 O HOH 17 16.732 42.270 53.224 1.00 43.17 O >HETATM 6680 O HOH 18 36.203 41.041 67.666 1.00 41.06 O >HETATM 6681 O HOH 19 -23.753 42.588 51.308 1.00 39.07 O
>CONECT 2034 6639 >CONECT 5345 6654 >CONECT 6623 6624 6625 6626 6627 >CONECT 6624 6623 >CONECT 6625 6623 >CONECT 6626 6623 >CONECT 6627 6623 >CONECT 6628 6629 6630 6631 6632
>CONECT 6629 6628 >CONECT 6630 6628 >CONECT 6631 6628 >CONECT 6632 6628 >CONECT 6633 6634 6641 >CONECT 6634 6633 6635 6636 >CONECT 6635 6634 >CONECT 6636 6634 6637 6638 >CONECT 6637 6636
>CONECT 6638 6636 6639 6640 >CONECT 6639 2034 6638 >CONECT 6640 6638 6641 6642 >CONECT 6641 6633 6640 >CONECT 6642 6640 6643 >CONECT 6643 6642 6644 >CONECT 6644 6643 6645 6646 6647
>CONECT 6645 6644 >CONECT 6646 6644 >CONECT 6647 6644 >CONECT 6648 6649 6656 >CONECT 6649 6648 6650 6651 >CONECT 6650 6649 >CONECT 6651 6649 6652 6653 >CONECT 6652 6651 >CONECT 6653 6651 6654 6655
>CONECT 6654 5345 6653 >CONECT 6655 6653 6656 6657 >CONECT 6656 6648 6655 >CONECT 6657 6655 6658 >CONECT 6658 6657 6659 >CONECT 6659 6658 6660 6661 6662 >CONECT 6660 6659 >CONECT 6661 6659
>CONECT 6662 6659 >MASTER 342 0 4 47 25 0 0 6 6679 2 42 66 >END >II have made alignment in PIR using chain breaks "/" and block residues "." for SO4,PLP AND HOH.My alignment looks like this:
> > > >P1;TvLDH >sequence:TvLDH: : : : ::: 0.00: 0.00 >MRIYGEEHPNQQILSRIATNDGHGENSSYFDGWKAYEKDPFHLTDNPTGVIQMGLAENQL >SLDLIRDWMKKNPQASICTEEGVSEFKAIANFQDYHGLPTFRKAIAQFMEKVRGGRTRFD
>PDRIVMSGGATGAQETIAFCLADPGEAFLIPTPYYPGFDR-FRWRTGVQLLPIHCHSSNK >FKITQAALETAYRKARNSHIRVKGILVTNPSNPLGTTMDRETLRTLVSFVNEKRMHLVCD >EIFSGTAFDKPSYVSVSEVIEDE--PYCDRDLIHIAYSLSKDLGVPGFRVGVIYSYNDAV >VSCARKMSSFGLVSSQTQHLLASMLGDEELTTSFLATSRTGLCGRRRVFTDGLKRVGIHC
>LDGNAGLFCWMDLRPLLKEATVEAELRLWRVIINDVKLNISPGSSFHCSEPGWFRVCFAN >MDDTAMKIALRRIESFVYRENDAAVQAKNKRRWDEALRLSLPRRRFEDPTIMTPHLMSPH >SPLVQAAT/..* > >P1;TvLDH_model >structureX:TvLDH-loop:1 : :485 : ::: 0.00: 0.00
>MRIYGEEHPNQQILSRIATNDGHGENSSYFDGWKAYEKDPFHLTDNPTGVIQMGLAENQL >SLDLIRDWMKKNPQASICTEEGVSEFKAIANFQDYHGLPTFRKAIAQFMEKVRGGRTRFD >PDRIVMSGGATGAQETIAFCLADPGEAFLIPTPYYPGFDR-FRWRTGVQLLPIHCHSSNK >FKITQAALETAYRKARNSHIRVKGILVTNPSNPLGTTMDRETLRTLVSFVNEKRMHLVCD
>EIFSGTAFDKPSYVSVSEVIEDE--PYCDRDLIHIAYSLSKDLGVPGFRVGVIYSYNDAV >VSCARKMSSFGLVSSQTQHLLASMLGDEELTTSFLATSRTGLCGRRRVFTDGLKRVGIHC >LDGNAGLFCWMDLRPLLKEATVEAELRLWRVIINDVKLNISPGSSFHCSEPGWFRVCFAN >MDDTAMKIALRRIESFVYRENDAAVQAKNKRRWDEALRLSLPRRRFEDPTIMTPHLMSPH
>SPLVQAAT/--* > >P1;1IAX >structureX:1IAX: 126 :A : 500 : :undefined:undefined:-1.00:-1.00 >------------------------------------------------------------ >------------------------------------------------------------
>--------GATGANETIIFCLADPGDAFLVPSPYYPAFNRDLRWRTGVQLIPIHCESSNN >FKITSKAVKEAYENAQKSNIKVKGLILTNPSNPLGTTLDKDTLKSVLSFTNQHNIHLVCD >EIYAATVFDTPQFVSIAEILDEQEMTYCNKDLVHIVYSLSK------------------- >------------------------------------------------------------
>------------------------------------------------------------ >------------------------------------------------------------ >--------/..* > > >Then I used advanced example as a template to write my input script:
> > from modeller import * > from modeller.automodel import * > >class mymodel(automodel): > def special_restraints(self, aln): > rsr = self.restraints > for ids in (('NZ:278:A','C4A:500:A')):
> > > atoms = [self.atoms[i] for i in ids] > rsr.add(forms.upper_bound(group=physical.upper_distance, > feature=features.distance(*atoms),
> mean=3.5, stdev=0.1)) > >env = environ() >env.io.hetatm = True >a = mymodel(env, alnfile='align-ligand.ali', > knowns=('TvLDH_model','1IAX'), sequence='TvLDH')
>a.starting_model = 1 >a.ending_model = 5 >a.make() > > >then program starts calculating, and then I always get : > > > > >check_ali___> Checking pairwise structural superpositions.
> >Equivalent CA pairs with distance difference larger than 6.0 angstroms: > >ALN_POS TMPL1 TMPL2 RID1 RID2 NAM1 NAM2 DIST >---------------------------------------------------- >END OF TABLE
> >check_ali___> Checking the sequence-structure alignment. > >Implied target CA(i)-CA(i+1) distances longer than 8.0 angstroms: > >ALN_POS TMPL RID1 RID2 NAM1 NAM2 DIST >----------------------------------------------
>END OF TABLE >read_to_681_> topology.submodel read from topology file: 3 > >getf_______W> RTF restraint not found in the atoms list: > residue type, indices: 17 485
> atom names : C +N > atom indices : 3847 0 > >getf_______W> RTF restraint not found in the atoms list: > residue type, indices: 17 485
> atom names : C CA +N O > atom indices : 3847 3843 0 3848 > > >fndatmi_285W> Only 153 residues out of 155 contain atoms of type CA
> (This is usually caused by non-standard residues, such > as ligands, or by PDB files with missing atoms.) >mdtrsr__446W> A potential that relies on one protein is used, yet you have at
> least one known structure available. MDT, not library, potential is used. >iup2crm_280W> No topology library in memory or assigning a BLK residue. > Default CHARMM atom type assigned: S --> S
> This message is written only for the first such atom. >20 atoms in HETATM residues constrained >to protein CA atoms within 10.00 angstroms >20 atoms in residues without defined topology
>constrained to be rigid bodies >indxatm_278E> No ":" in ATOM:RESID[:CHAINID] atom identifier: N > > > >but when I used model-ligand.py command it generate model including ligands successfully.
>So my questions are: >1. What am I doing wrong with model-multiple-hetero.py command? >2.how 2 solve this problem?? > >What do you think? Thank you for your suggestions.Sorry for my bad english.
>sanjay > > > > > >sanjay Kumar singh kolkata >Bose Institute >Department of Botany >Main Campus >93\1 A.P.C.Road >Kolkata 700 009 >India >-------------- next part --------------
>An HTML attachment was scrubbed... >URL: http://salilab.org/archives/modeller_usage/attachments/20080815/07a82fbb/attachment.html
> >------------------------------ > >_______________________________________________ >modeller_usage mailing list >">
>https://salilab.org/mailman/listinfo/modeller_usage > > >End of modeller_usage Digest, Vol 7, Issue 54 >*********************************************