May 2020 - modeller_usage

Command not found, no module named 'modeller' errors
by Xiao Ran Luo 29 May '20

29 May '20

Hello, I just downloaded version 9.24 of Modeller onto by macOS version 10.14.1 using Homebrew, but when I try to run the first few lines of the basic tutorial, I get a lot of "command not found": Luos-Air:~ luoxiaoran$ >P1;TvLDH -bash: TvLDH: command not found Luos-Air:~ luoxiaoran$ sequence:TvLDH:::::::0.00: 0.00 -bash: sequence:TvLDH:::::::0.00:: command not found Luos-Air:~ luoxiaoran$ MSEAAHVLITGAAGQIGYILSHWIASGELYGDRQVYLHLLDIPPAMNRLTALTMELEDCAFPHLAGFVATTDPKA -bash: MSEAAHVLITGAAGQIGYILSHWIASGELYGDRQVYLHLLDIPPAMNRLTALTMELEDCAFPHLAGFVATTDPKA: command not found Luos-Air:~ luoxiaoran$ AFKDIDCAFLVASMPLKPGQVRADLISSNSVIFKNTGEYLSKWAKPSVKVLVIGNPDNTNCEIAMLHAKNLKPEN -bash: AFKDIDCAFLVASMPLKPGQVRADLISSNSVIFKNTGEYLSKWAKPSVKVLVIGNPDNTNCEIAMLHAKNLKPEN: command not found Luos-Air:~ luoxiaoran$ FSSLSMLDQNRAYYEVASKLGVDVKDVHDIIVWGNHGESMVADLTQATFTKEGKTQKVVDVLDHDYVFDTFFKKI -bash: FSSLSMLDQNRAYYEVASKLGVDVKDVHDIIVWGNHGESMVADLTQATFTKEGKTQKVVDVLDHDYVFDTFFKKI: command not found Luos-Air:~ luoxiaoran$ GHRAWDILEHRGFTSAASPTKAAIQHMKAWLFGTAPGEVLSMGIPVPEGNPYGIKPGVVFSFPCNVDKEGKIHVV -bash: GHRAWDILEHRGFTSAASPTKAAIQHMKAWLFGTAPGEVLSMGIPVPEGNPYGIKPGVVFSFPCNVDKEGKIHVV: command not found Luos-Air:~ luoxiaoran$ EGFKVNDWLREKLDFTEKDLFHEKEIALNHLAQGG* -bash: EGFKVNDWLREKLDFTEKDLFHEKEIALNHLAQGG*: command not found Luos-Air:~ luoxiaoran$ Then when I try to do this (I am not sure what it is supposed to do, but have been guaranteed that it works) I get: Luos-Air:~ luoxiaoran$ cd /users/luoxiaoran/Documents/modAb101 Luos-Air:modAb101 luoxiaoran$ python3 modelAb101.py Traceback (most recent call last): File "modelAb101.py", line 2, in <module> from modeller import * # Load standard Modeller classes ModuleNotFoundError: No module named 'modeller' Luos-Air:modAb101 luoxiaoran$ I have tried troubleshooting it by downloading it via the Mac Installer instead, but that just seemed to interfere with everything so I uninstalled both versions of Modeller using brew uninstall modeller and sudo rm -rf /Library/modeller-9.24/ /usr/local/bin/mod9.24 /Library/Python/2.?/site-packages/modeller.pth /var/db/receipts/org.salilab.modeller.*, and reinstalled it using only brew install modeller. During uninstalling, I also uninstalled and reinstalled Homebrew, and manually deleted all the folders that Terminal informed me could not be deleted automatically (including the Python framework, which I later found out should never be deleted, so I reinstalled python and python(a)3.8 using both Homebrew and the official Python website's Mac Installer for Python 3.8.3) before the reinstallation of Homebrew and modeller. I also entered the modeller license key by modifying the config.py document, just as the program told me to. Do you know what's wrong, how to find out what's wrong, and/or how to solve it? Thank you so much! Best, Luna

3 9

Re: [modeller_usage] Command not found, no module named 'modeller' errors
by Xiao Ran Luo 28 May '20

28 May '20

Thanks for your reply! Unfortunately, my Mandarin is not good enough to reply in Mandarin, sorry -do you mean I should uninstall Python 3.8.3 and install Python 2.7 or below? I did that, moved Python from Applications to Trash, and did brew uninstall python and brew uninstall python(a)3.8, then I installed Python 2.7.18 from the official website (because I can't seem to install versions below 2.7, the Installer keeps telling me that "there was no software found to install"). However, when I try to check what version of Python I am using on Terminal, I get: Luos-Air:~ luoxiaoran$ python --version Python 2.7.10 Luos-Air:~ luoxiaoran$ python3 --version Python 3.8.3 Luos-Air:~ luoxiaoran$ And the above code from my first email still does not work. On Thu, May 28, 2020 at 6:14 PM Min-yi Shen <minyishen(a)gmail.com> wrote: > Modeler 其實就是 python 裡的一模塊。讓事情複雜的是：一般的機子上還會有 python 和 python3 > 這兩個，有時這破模塊會被裝到另一個 python 下，我看你用的是 python3 版，找不到。試試 python （2.7 版）下面有沒有。 > > （註明：我現在不是 modeler 開發者了，12 年前是） > > On Thu, May 28, 2020 at 18:04 Xiao Ran Luo via modeller_usage < > modeller_usage(a)salilab.org> wrote: > >> Hello, >> >> I just downloaded version 9.24 of Modeller onto by macOS version 10.14.1 >> using Homebrew, but when I try to run the first few lines of the basic >> tutorial, I get a lot of "command not found": >> >> Luos-Air:~ luoxiaoran$ >P1;TvLDH >> >> -bash: TvLDH: command not found >> >> Luos-Air:~ luoxiaoran$ sequence:TvLDH:::::::0.00: 0.00 >> >> -bash: sequence:TvLDH:::::::0.00:: command not found >> >> Luos-Air:~ luoxiaoran$ >> MSEAAHVLITGAAGQIGYILSHWIASGELYGDRQVYLHLLDIPPAMNRLTALTMELEDCAFPHLAGFVATTDPKA >> >> -bash: >> MSEAAHVLITGAAGQIGYILSHWIASGELYGDRQVYLHLLDIPPAMNRLTALTMELEDCAFPHLAGFVATTDPKA: >> command not found >> >> Luos-Air:~ luoxiaoran$ >> AFKDIDCAFLVASMPLKPGQVRADLISSNSVIFKNTGEYLSKWAKPSVKVLVIGNPDNTNCEIAMLHAKNLKPEN >> >> -bash: >> AFKDIDCAFLVASMPLKPGQVRADLISSNSVIFKNTGEYLSKWAKPSVKVLVIGNPDNTNCEIAMLHAKNLKPEN: >> command not found >> >> Luos-Air:~ luoxiaoran$ >> FSSLSMLDQNRAYYEVASKLGVDVKDVHDIIVWGNHGESMVADLTQATFTKEGKTQKVVDVLDHDYVFDTFFKKI >> >> -bash: >> FSSLSMLDQNRAYYEVASKLGVDVKDVHDIIVWGNHGESMVADLTQATFTKEGKTQKVVDVLDHDYVFDTFFKKI: >> command not found >> >> Luos-Air:~ luoxiaoran$ >> GHRAWDILEHRGFTSAASPTKAAIQHMKAWLFGTAPGEVLSMGIPVPEGNPYGIKPGVVFSFPCNVDKEGKIHVV >> >> -bash: >> GHRAWDILEHRGFTSAASPTKAAIQHMKAWLFGTAPGEVLSMGIPVPEGNPYGIKPGVVFSFPCNVDKEGKIHVV: >> command not found >> >> Luos-Air:~ luoxiaoran$ EGFKVNDWLREKLDFTEKDLFHEKEIALNHLAQGG* >> >> -bash: EGFKVNDWLREKLDFTEKDLFHEKEIALNHLAQGG*: command not found >> >> Luos-Air:~ luoxiaoran$ >> >> >> >> Then when I try to do this (I am not sure what it is supposed to do, but >> have been guaranteed that it works) I get: >> >> Luos-Air:~ luoxiaoran$ cd /users/luoxiaoran/Documents/modAb101 >> >> Luos-Air:modAb101 luoxiaoran$ python3 modelAb101.py >> >> Traceback (most recent call last): >> >> File "modelAb101.py", line 2, in <module> >> >> from modeller import * # Load standard Modeller classes >> >> ModuleNotFoundError: No module named 'modeller' >> >> Luos-Air:modAb101 luoxiaoran$ >> >> >> I have tried troubleshooting it by downloading it via the Mac Installer >> instead, but that just seemed to interfere with everything so I uninstalled >> both versions of Modeller using brew uninstall modeller and sudo rm -rf >> /Library/modeller-9.24/ /usr/local/bin/mod9.24 >> /Library/Python/2.?/site-packages/modeller.pth >> /var/db/receipts/org.salilab.modeller.*, and reinstalled it using only brew >> install modeller. During uninstalling, I also uninstalled and >> reinstalled Homebrew, and manually deleted all the folders that Terminal >> informed me could not be deleted automatically (including the Python >> framework, which I later found out should never be deleted, so I >> reinstalled python and python(a)3.8 using both Homebrew and the official >> Python website's Mac Installer for Python 3.8.3) before the reinstallation >> of Homebrew and modeller. I also entered the modeller license key by >> modifying the config.py document, just as the program told me to. >> >> >> Do you know what's wrong, how to find out what's wrong, and/or how to >> solve it? Thank you so much! >> >> >> Best, >> >> Luna >> _______________________________________________ >> modeller_usage mailing list >> modeller_usage(a)salilab.org >> https://salilab.org/mailman/listinfo/modeller_usage >> >

1 0

How to properly set the fields 5 of the second line in the .aln file?
by ZHANG Cheng 21 May '20

21 May '20

I have a two-chain protein to model, Chain A and Chain B. In the PDB file, the Chain B residues are numbered from 3-33, so only 31 residues appear in the PDB. When "1 :A:+33 :B" is used for the fields 3-6 in the second line of the .aln file, I got error: "Number of residues in the alignment and  pdb files are different". When it is changed to "1 :A:+31 :B", the same error occurred. However, when only Chain B is modelled, "3 :B:+33 :B" can be run successfully. So, it seems that, when Field 5 = n, it assumes that the last chain has in total n residues, not the last residue is numberred as n. So how should I set the Field 5 number? Do I have to renumber the Chain B residues in the PDB from 1 to 31? Thank you!

2 2

Regarding difficult modelling in modeller
by Cn rahul 20 May '20

20 May '20

I used the mGenThreader server to get a .align file as output. How to edit/ convert to respective .ali file as required by modeller python script. Do we have to do it manually ? -- Dr. C N RAHUL Ph: 09666011197 Email: cnrahul(a)gmail.com

2 1

Can the Modeller license be available as a site license?
by ZHANG Cheng 20 May '20

20 May '20

Dear Modeller, I was trying to ask our HPC staff to install the Modeller in HPC, so that we can access it using multiple cores. But I was told: The Software Manager has discovered that the licence is per user and not available as a site license. So can I ask, can a site license be available? Thank you! Yours sincerely Cheng

1 0

How to constrain the computational resources used for Modeller?
by ZHANG Cheng 18 May '20

18 May '20

I need to build a model with over 16,000 residues on Win10 which requires high demand for memory (maybe CPU as well). Can I set the maximum memory MODELLER uses, so that my PC will not freeze, and I can terminate the program without forcing the PC to restart? Thank you!

2 2

missing loop
by Qinghua Liao 16 May '20

16 May '20

Hello, In my protein, there is one missing loop, 15 residues. I used the script from the tutorials to add the missing residues while keeping other coordinates: # from modeller import * from modeller.automodel import * # Load the automodel class log.verbose() env = environ() # directories for input atom files env.io.atom_files_directory = ['.', '../atom_files'] class MyModel(automodel): def select_atoms(self): return selection(self.residue_range('133', '135'), self.residue_range('217', '230')) a = MyModel(env, alnfile = 'alignment.ali', knowns = '1qg8', sequence = '1qg8_fill') a.starting_model= 1 a.ending_model = 1 a.make() # However, the loops added are too extended, even have bad contacts with the other parts of the protein. As the figure in the Dropbox: https://www.dropbox.com/s/3obtrrvc9c9w8df/image.png?dl=0 Could someone guide me how to add the loop more properly? Thanks a lot! All the best, Qinghua

2 4

Use modeller inside class method - called from an ThreadPoolExecutor
by Pedro Guillem 07 May '20

07 May '20

Hi All, I have a class called Mutant(), its constructor generates a random FASTA sequence. A class method *writeToPDB()* uses modeller to push the FASTA into a template PDB. This is already working very well. If I do 1000 mutants this takes an eternity (one after the other)... so I thought of crafting a ThreadPoolExecutor. Modeller runs inside the class method. Concurrent calls should not share memory with other instances of Mutant(). mutants = [] #Holds an Array of mutant objects executor = ThreadPoolExecutor(max_workers=20) #Thread executor for i in range(1000): mutants.append(Mutant()) executor.submit(mutants[i].writeToPDB, templatePDB) This simply tells the executor to queue a new thread with the writeToPDB() method. When I try to do this, it seems that modeller is being executed in the first call and sharing its state with all subsequent calls... so I get like 90 models (not 1000) written in whatever state modeller was when running the other instances. How can I force a class method to use its own modeller instance?. This is the method: env = environ() env.libs.topology.read(file='$(LIB)/top_heav.lib') env.libs.parameters.read(file='$(LIB)/par.lib') aln = alignment(env) mdl = model(env, file=code) aln.append_model(mdl, atom_files=code, align_codes=code) residue = self.fasta[resid-Mutant.compensate] sel.mutate(residue_type=Mutant.res1to3[residue]) #mutate aln.append_model(mdl, align_codes='mut') mdl.clear_topology() mdl.generate_topology(aln['mut']) mdl.transfer_xyz(aln) mdl.build(initialize_xyz=False, build_method='INTERNAL_COORDINATES') name = str(self.pdbpath)+str(self.id)+".pdb" mdl.write(file=name)

2 1

Retain residue IDS when mutating and writing PDB
by Pedro Guillem 05 May '20

05 May '20

I'm using Modeller 9.24 to induce a couple of mutations into a PDB. I noticed that when I execute model.mutate() and export the PDB, the residue indexes are completely different from the original file... ebery residue looses its original position in the sequence and gets renumbered to therever it lands. Is there a way I can force modeller to stick to the original numbering? Best Pedro

2 1

Threading sequence onto template with missing residues
by Ahmad Khalifa 04 May '20

04 May '20

Hello, I have a template with lots of missing residues. I want to thread my target sequence onto the aligned parts of the template that are present without moving my target structure in these areas. In my target sequence, I also wish to know how to skip modeling the missing parts altogether, thus making the job a mere threading of aligned unknown regions onto the known ones. Is there a simple way to accomplish my task without doing too much manual alignment and residue listing? Please do post a code snippet if you can. Thanks a lot.

2 1