December 2004 - modeller_usage

Superposition problem
by Shamshad Zarina 29 Dec '04

29 Dec '04

Dear Modeller Community, I am trying to superpose two models using superpose command. The command runs okay and it creates a .fit file. Although the fit file looks fine but it has an unusual line connecting the last and first residues of both models. Any ideas why this is happening? thanks, Shamshad __________________________________ Do you Yahoo!? Yahoo! Mail - Easier than ever with enhanced search. Learn more. http://info.mail.yahoo.com/mail_250

2 1

Re: [modeller_usage] (no subject)
by Modeller Caretaker 28 Dec '04

28 Dec '04

Anagha Joshi wrote: > when I try to run the basic model tutorial at modeller site > I run align2d.top > I get warning > mkapsa__293W> No TOPOLOGY_LIB is in memory. Use READ_TOPOLOGY to read one. > iup2crm_280W> No TOPOLOGY_LIB in memory or assigning a BLK residue. > Default CHARMM atom type assigned: C --> CT2 > This message is written only for the first such atom. > fndatmi_285W> Number of residues <> number of atoms; atom code: 668 > 66 > 6 CA > fndatmi_285W> Number of residues <> number of atoms; atom code: 668 > 66 > 6 CA > I > don,t > understand what this warning is. It is only a warning, and can be ignored in this case. If you are going to build a comparative model, then Modeller needs to know that a 'CT2' atom is a carbon, so that it can match it against similar carbon atoms in the templates, and use molecular mechanics parameters defined for a CT2 atom. It gets this information from the topology and parameter files. However, for the case of alignment, this is not needed. 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

Problems with stereo-chemical fixing of Hydrogen atoms
by Ravi Rathore 28 Dec '04

28 Dec '04

Dear modellers, Q: Is it possible to built short-contact-free models (hydrogens included) with Modeller? I am building an all-atom model using following scripts, In modeller 7v7: SET TOPOLOGY_MODEL = 1, HYDROGEN_IO = on SET TOPLIB = $(LIB)/top_allh.lib SET PARLIB = $(LIB)/par.lib CALL ROUTINE = 'model' ----------------------------------------- In modeller6: SET TOPOLOGY_MODEL = 1, HYDROGEN_IO = on SET TOPLIB = $(LIB)/top.lib (also tried top_allh.lib) SET PARLIB = $(LIB)/par.lib CALL ROUTINE = 'model' ------------------------------------------- The models do have hydrogen atoms, however, I noticed that most of the hydrogens have short-contacts within the same atom or residue OR with neighboring residues and also have long-range clashes. Assuming that there is everything fine with topology or parameter files, I wonder, the energy minimization/MD should have corrected such clashes. Then, I tried to use other, thorough optimization protocols, but of no use: SET LIBRARY_SCHEDULE = 1 SET MAX_VAR_ITERATIONS = 300 SET MD_LEVEL = 'refine_5' SET REPEAT_OPTIMIZATION = 3 SET MAX_MOLPDF = 1E6 I would be grateful for any suggestion on proper stereo-chemical fixing of hydrogen atoms with Modeller. Sincerely, Rathore, R.S. --------------------------------- Do you Yahoo!? Send holiday email and support a worthy cause. Do good.

2 1

(no subject)
by anagha＠mbu.iisc.ernet.in 28 Dec '04

28 Dec '04

dear sir, I am new to modeller user and also to linux system. i have downloaded and installed modeller But as it doent have GUI. How will I start using it. I need help. Thank you Anagha

3 2

is it a bug?
by WHONG Tienwynn 21 Dec '04

21 Dec '04

Dear Ben Webb, the files are atteched. the configuration of my PC are as follows: Windows2000, sp4 celeron4, 1.7Ghz(but i run it at 2.04G) Physical Memery is 256M, virtual Memery is 512M. When virtual memery is 372-744M, it gave the error message as" stock overflow". i increased the virtual memery, it disapeared, but the message shown in the picture appeared again. what is wrong? tienwynn whong ===== ============================== Tienwynn(a)Yahoo.com Genejacob(a)Yahoo.com.cn ============================== _________________________________________________________ Do You Yahoo!? 注册世界一流品质的雅虎免费电邮 http://cn.rd.yahoo.com/mail_cn/tag/1g/*http://cn.mail.yahoo.com/

1 0

Regarding the running of modeller on SGI machine
by Vohra S 21 Dec '04

21 Dec '04

Hi There is a problem regarding the running of the modeller.While running the script file on the prompt an error is shown and the log file is not created. I typed the following command on the prompt > mod7 seq_search.top >seq.log01 where mod7 is the alias directed to the executable file The error at the UNIX prompt is as follows /usr28/users/gpsaro/modeller7v7/bin/mod7v7[35]: ulimit: exceeds allowable limit The log file created is empty Please can you help me with this Regards Shabana

3 2

"Precision problem"
by Oliver Hucke 21 Dec '04

21 Dec '04

Hi all, I just started to use the newest version of modeller. I noticed unfamiliar warnings in the log file: conmin__610W> Precision problem; starting from scratch. conmin__611W> Precision problem; starting from scratch. Last such mesage in this OPTIMIZE What does this mean? Should I worry about this or can it be ignored? Thanks in advance, Oliver -- ________________________________________________________________ Oliver Hucke, Dr. Biomolecular Structure Center Health Sciences Building - K418C Dept. of Biochemistry 1959 NE Pacific St. University of Washington phone: (206) 685 7046 Box 357742 fax : (206) 685 7002 Seattle, WA 98195-7742 email: ohucke(a)u.washington.edu ________________________________________________________________

2 1

Script for superimposing generated models
by Ivan Bosanac 17 Dec '04

17 Dec '04

Hello All, I have been having some problems superimposing generated models on top of my template coordinance. I have used SET FINAL_MALIGN3D = 1 command in my top file but not *_fit files were given. This is puzzling to me since using the same top script for a different modeling job worked well and *_fit files were generated. Let me say that I am very new to Modeller, and script driven software in general, and for that reason would greatly appreciate if somebody could explain to me why *_fit files were not created. Moreover, I would be very grateful if somebody could provide me with the script that could be used for the purpose of superimposing .B999000* onto the template file after top script has been executed. Thanks in advance, Ivan ps: This is my model.top file: # Homology modelling by the MODELLER TOP routine 'model'. INCLUDE # Include the predefined TOP routines SET OUTPUT_CONTROL = 1 1 1 1 1 # uncomment to produce a large log file SET ALNFILE = 'align.ali' # alignment filename SET KNOWNS = 't1b1' # codes of the templates SET SEQUENCE = 't3b1' # code of the target SET ATOM_FILES_DIRECTORY = '' # directories for input atom files SET STARTING_MODEL= 1 # index of the first model SET ENDING_MODEL = 30 # index of the last model # (determines how many models to calculate) SET MD_LEVEL = 'refine_3' # defalt is 1; 5 is very slow with MD anne. SET FINAL_MALIGN3D = 1 # superimposis template and resulting pdb SET LIBRARY_SCHEDULE = 1 SET REPEAT_OPTIMIZATION = 10 SET RSTRS_REFINED = 2 SET MD_RETURN = 'FINAL' CALL ROUTINE = 'model' # do homology modelling

2 1

array problem
by vivek sharma 17 Dec '04

17 Dec '04

Dear sir, I have included few residues in top_heav.lib and I am using MODELLER 6v2. After running it gave errors regarding array, it says array small for MAPRES, so I deleted some atoms from my defined residue and it worked a step better but now it gives another error, similar to that. Please explain what can be the reason. Error is: read_to_223E> Array too small. Increase MCNSRES current maximum, current need: 100 104 recover____E> ERROR_STATUS >= STOP_ON_ERROR: 1 1 Thanks in advance. Regards, Vivek Sharma -- VIVEK SHARMA HBG BI "You must be the change you wish to see in the world." MKG

2 1

Bad Model
by A R 16 Dec '04

16 Dec '04

Hi, Im trying to model a 58% identity protein with the template, but after alignment, the model.top generates a very different model from the template. Can anybody suggests me what can I change in my model.top to make this better? Thanks. The alignment is using MALIGN with gap penalties of -600 and -400: >P1;1xxx structureX:1xxx: 5 :A: 689 :B:undefined:undefined:-1.00:-1.00 MAVKVEYDLKRLRNIGIAAHIDAGKTTTTERILYYTGRIHKIGEVHEGAATMDFMEQERERGITITAAVTTCFWK DHRINIIDTPGHVDFTIEVERSMRVLDGAIVVFDSSQGVEPQSETVWRQAEKYKVPRIAFANKMDKTGADLWLVI RTMQERLGARPVVMQLPIGREDTFSGIIDVLRMKAYTYGNDLGTDIREIPIPEEYLDQAREYHEKLVEVAADFDE NIMLKYLEGEEPTEEELVAAIRKGTIDLKITPVFLGSALKNKGVQLLLDAVVDYLPSPLDIPPIKGTTPEGEVVE IHPDPNGPLAALAFKIMADPYVGRLTFIRVYSGTLTSGSYVYNTTKGRKERVARLLRMHANHREEVEELKAGDLG AVVGLKETITGDTLVGEDAPRVILESIEVPEPVIDVAIEPKTKADQEKLSQALARLAEEDPTFRVSTHPETGQTI ISGMGELHLEIIVDRLKREFKVDANVGKPQVAYRETITKPVDVEGKFIRQTGGRGQYGHVKIKVEPLPRGSGFEF VNAIVGGVIPKEYIPAVQKGIEEAMQSGPLIGFPVVDIKVTLYDGSYHEVDSSEMAFKIAGSMAIKEAVQKGDPV ILEPIMRVEVTTPEEYMGDVIGDLNARRGQILGMEPRGNAQVIRAFVPLAEMFGYATDLRSKTQGRGSFVMFFDH YQEVPKQVQEKLIKGQ* >P1;03883 sequence:03883: : : 696 : :::-1.00:-1.00 -MMARKFELKDYRNIGIMAHIDAGKTTTTERILFHTGKIHKIGETHDGVSQMDWMEQEKERGITITSAATTAFWK GKRINIIDTPGHVDFTVEVERSLRVLDGAVAVLDAQSGVEPQTETVWRQATNYSVPRIVYVNKMDKAGANFEASI ESVRTKLNGNAVAIQLNIGAEADFSGLIDLVEMKAYNYDGQKEEIEYEIPIPEDLFEKASQMRLALAEAVADYDE EIFNNLLEEKEILPEQLKAAIRAATITGNFFPVVCGSSFKNKGVKKMIDAVIDYLPSPVDVPPIKAFRDEEEIT- IEASDDQEFSALAFKIMNDPFVGSLTFFRVYSGVLKKGTYIINSTKGKKERVGRILAMHANSREEIDEVRTGDIG AFVGLKDTTTGDSLISEKAKTFVLERMNFPEPVISQSLEPFSKAEIEKLATALQKLANEDPTFKTWTDIETGQTI IAGMGELHLDIIVDRLKREFNVQARVGKPQVSYRETITKSAEVEGKYIKQSGGRGQYGHVWIKFEPNPE-EGFDF IDKIVGGKIPKEYIKSIQKGLEEKMQAGILAGYPLINLRATLFDGSFHEVDSSEMAFKIAASKALSRARDAVGTV LLEPIMDVSVFAPSEYAGDVMGDLSRRRGLVREQETRSDGANVIRGHVPLAEMFGYSTQLRSMTSGRGTYQMQFN HYEIVPKNISDVIVKQRAIKEDD* the model.top is: INCLUDE SET ALNFILE = 'malign.pir' SET KNOWNS = '1ktv' SET SEQUENCE = '03883' SET ATOM_FILES_DIRECTORY = './:../atom_files' SET STARTING_MODEL = 1 SET ENDING_MODEL = 5 CALL ROUTINE = 'model' ___________________________________________________________ Moving house? Beach bar in Thailand? New Wardrobe? Win £10k with Yahoo! Mail to make your dream a reality. Get Yahoo! Mail www.yahoo.co.uk/10k

2 1