- modeller_usage - salilab.org

Best workflow to modeling a disordered loop in a complicated structure
by Zachary W Carpenter 19 Feb '14

19 Feb '14

Hello, I have used modeller to attempt to model a loop region in a tetrameric structure. To speed up computational time I've started with just the dimer ( which I think also makes sense biologically ) and ran the following script with an .ali file I created. The region of interest is a disordered loop AA~399-417. Following this I used the evaluate model script to pick the best model, which I then used as the template input for the loop.py script. I run the loop.py script at many iterations and then when it finishes I take the top 20 models (out of 1000) by DOPE score and attempt to analyze them. Any advice on how I could improve my method would be great.. I really am at the edge of my understanding here. I've edited the name of the protein. I can include it as soon as I can verify that this can stay anonymous. Thanks!!! Scripts +files : First )using a template pdb and .ali file I made myself of the Dimer (BIO_headers in chimera and deletion of one dimer of the tetramer) Input : >P1;ABCD structureX:ABCD.pdb: 3 :A:+990:B:::-1.00:-1.00 --TSWSDRLQNAADMPANMDKHALKKYRREAYHRVFVNRSLAMEKIKCFGFNMDYTLAVYKSPEYESLGFELTVE RLVSIGYPQELLSFAYDSTFPTRGLVFDTLYGNLLKVDAYGNLLVCAHGFNFIRGPETREQYPNKFIQRDDTERF YILNTLFNLPETYLLACLVDFFTNCPRYTSCETGFKDGDLFMSYRSMFQDVRDAVDWVHYKGSLKEKTVENLEKY VVKDGKLPLLLSRMKEVGKVFLATNSDYKYTDKIMTYLFDFPHGPKPGSSHRPWQSYFDLILVDARKPLFFGEGT VLRQVDTKTGKLKIGTYTGPLQHGIVYSGGSSDTICDLLGAKGKDILYIGDHIFGDILKSKKRQGWRTFLVIPEL AQELHVWTDKSSLFEELQSLDIFLAS----------------SIQRRIKKVTHDMDMCYGMMGSLFRSGSRQTLF ASQVMRYADLYAASFINLLYYPFSYLFRAAHVLMPHES/--TSWSDRLQNAADMPANMDKHALKKYRREAYHRVFVNRSLAMEKIKCFGFNMDYTLAVYKSPEYESLGFELTVE RLVSIGYPQELLSFAYDSTFPTRGLVFDTLYGNLLKVDAYGNLLVCAHGFNFIRGPETREQYPNKFIQRDDTERF YILNTLFNLPETYLLACLVDFFTNCPRYTSCETGFKDGDLFMSYRSMFQDVRDAVDWVHYKGSLKEKTVENLEKY VVKDGKLPLLLSRMKEVGKVFLATNSDYKYTDKIMTYLFDFPHGPKPGSSHRPWQSYFDLILVDARKPLFFGEGT VLRQVDTKTGKLKIGTYTGPLQHGIVYSGGSSDTICDLLGAKGKDILYIGDHIFGDILKSKKRQGWRTFLVIPEL AQELHVWTDKSSLFEELQSLDIFLAS----------------SIQRRIKKVTHDMDMCYGMMGSLFRSGSRQTLF ASQVMRYADLYAASFINLLYYPFSYLFRAAHVLMPHES* >P1;X sequence:ABCD: : : : ::: 0.00: 0.00 MSTSWSDRLQNAADMPANMDKHALKKYRREAYHRVFVNRSLAMEKIKCFGFDMDYTLAVYKSPEYESLGFELTVE RLVSIGYPQELLSFAYDSTFPTRGLVFDTLYGNLLKVDAYGNLLVCAHGFNFIRGPETREQYPNKFIQRDDTERF YILNTLFNLPETYLLACLVDFFTNCPRYTSCETGFKDGDLFMSYRSMFQDVRDAVDWVHYKGSLKEKTVENLEKY VVKDGKLPLLLSRMKEVGKVFLATNSDYKYTDKIMTYLFDFPHGPKPGSSHRPWQSYFDLILVDARKPLFFGEGT VLRQVDTKTGKLKIGTYTGPLQHGIVYSGGSSDTICDLLGAKGKDILYIGDHIFGDILKSKKRQGWRTFLVIPEF AQELHVWTDKSSLFEELQSLDIFLAELYKHLDSSSNERPDISSIQRRIKKVTHDMDMCYGMMGSLFRSGSRQTLF ASQVMRYADLYAASFINLLYYPFSYLFRAAHVLMPHES/MSTSWSDRLQNAADMPANMDKHALKKYRREAYHRVFVNRSLAMEKIKCFGFDMDYTLAVYKSPEYESLGFELTVE RLVSIGYPQELLSFAYDSTFPTRGLVFDTLYGNLLKVDAYGNLLVCAHGFNFIRGPETREQYPNKFIQRDDTERF YILNTLFNLPETYLLACLVDFFTNCPRYTSCETGFKDGDLFMSYRSMFQDVRDAVDWVHYKGSLKEKTVENLEKY VVKDGKLPLLLSRMKEVGKVFLATNSDYKYTDKIMTYLFDFPHGPKPGSSHRPWQSYFDLILVDARKPLFFGEGT VLRQVDTKTGKLKIGTYTGPLQHGIVYSGGSSDTICDLLGAKGKDILYIGDHIFGDILKSKKRQGWRTFLVIPEF AQELHVWTDKSSLFEELQSLDIFLAELYKHLDSSSNERPDISSIQRRIKKVTHDMDMCYGMMGSLFRSGSRQTLF ASQVMRYADLYAASFINLLYYPFSYLFRAAHVLMPHES* # Homology modeling by the automodel class # # Demonstrates how to build multi-chain models, and symmetry restraints # from modeller import * from modeller.automodel import * # Load the automodel class log.verbose() # Override the 'special_restraints' and 'user_after_single_model' methods: class MyModel(automodel): def special_restraints(self, aln): # Constrain the A, B, C and D chains to be identical s1 = selection(self.chains['A']).only_atom_types('CA') s2 = selection(self.chains['B']).only_atom_types('CA') self.restraints.symmetry.append(symmetry(s1, s2, 1)) def user_after_single_model(self): # Report on symmetry violations greater than 1A after building # each model: self.restraints.symmetry.report(1) env = environ() # directories for input atom files env.io.atom_files_directory = ['.', '../atom_files'] # Be sure to use 'MyModel' rather than 'automodel' here! a = MyModel(env, alnfile = 'WT_dimer_NoCterm.ali' , # alignment filename knowns = 'ABCD', # codes of the templates sequence = 'ABCD') # code of the target a.starting_model= 1 # index of the first model a.ending_model = 20 # index of the last model # (determines how many models to calculate) a.make() # do homology modeling # # class MyModel(automodel): 2nd) Run eval model and pick top model to use as input for the loop.py script below ---> # Loop refinement of an existing model from modeller import * from modeller.automodel import * log.verbose() env = environ() # directories for input atom files env.io.atom_files_directory = './:../atom_files' # Create a new class based on 'loopmodel' so that we can redefine # select_loop_atoms (necessary) class MyLoop(loopmodel): # This routine picks the residues to be refined by loop modeling def select_loop_atoms(self): # 10 residue insertion return selection(self.residue_range('399:A', '417:A')) m = MyLoop(env, inimodel='ABCD_Top_Model.pdb', # initial model of the target sequence='ABCD') # code of the target m.loop.starting_model= 1 # index of the first loop model m.loop.ending_model = 1000 # index of the last loop model m.loop.md_level = refine.slow # loop refinement method; this yields # models quickly but of low quality; # use refine.slow for better models m.make() 3rd) Run model_energies script Then I rank them using this perl one liner ( just cats and displays the dope and the file ) can run in any unix terminal. cat model*.log |perl -lane' print if /DOPE\ score/ || /ABCD.BL/'| perl -pe 's/ +/\t/g'|perl -pe 's/:\t//g' |cut -f3,4|perl -pe 's/\n/\t/'|perl -pe 's/(\d)\t(ABCD)/$1\n$2/g'|sort -k2,2n|head -20 4) Take top 20 models and analyze Any help with where I'm doing things wrong or how I can improve would be awesome. How many models in loop.py would you recommend to get a good answer? Is there away to interpret the models created as a whole ....for example say something like "60% of the time residue A is next to residue B and 20% of the time its over in this pocket , therefor".... etc My major aim is to compare different ligand bound crystal states, or mutation states, and show that the models created are different between them. Then use this to design wet lab experiments to test. Thanks again

1 0

Run Modeller on multiple cpu
by Lalith Kumar 17 Feb '14

17 Feb '14

Dear Modeller team, I need your quick suggestions and corrections on my scripts that I have been using to run on multiple cpus. I am not getting any errors, when I run the script parallel-task.py with python. But it is not utilizing multiple cpus. Can you please look in to the attached scripts and let me know where I might have went wrong. Please... Thank you lot Regards, Lalith

2 1

Modeller 9.13 release
by Modeller Caretaker 12 Feb '14

12 Feb '14

The new version of Modeller, 9.13, is now available for download! Please see the download page at http://salilab.org/modeller/ for more information. If you have a license key for Modeller 8 or 9, there is no need to reregister for Modeller 9.13 - the same license key will work. (It won't do any harm to reregister if you want to, though!) 9.13 is primarily a bugfix release relative to the last public release (9.12). Major user-visible changes include: # Modeller now includes a variety of SOAP (statistically optimized atomic potential) scores for assessing proteins, loops, and interfaces. # The Lennard-Jones interaction energy is now artificially truncated at very short distance; this makes simulations with poor starting conditions much less likely to 'blow up'. # model.get_insertions(), model.get_deletions() and model.loops() now have an include_termini option; if False, residue ranges that include chain termini are excluded from the output. See the Modeller manual for a full change log: http://salilab.org/modeller/9.13/manual/node39.html If you encounter bugs in Modeller 9.13, please see http://salilab.org/modeller/9.13/manual/node10.html for information on how to report them. Note: you are receiving this email either because you subscribed to the modeller_usage mailing list, or you provided this address when you requested a Modeller license and you ticked the "notify me of new releases" box. In the latter case, if you no longer wish to receive announcements of new Modeller releases, simply reply to this email and let us know. 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

Reading heteroatoms
by Lucas 12 Feb '14

12 Feb '14

Dear all, I'm willing to model a protein containing zinc atoms, something I've done successfully before using earlier versions of modeller. The template has protein and two zinc ions in chain A, and so I added two dots in the end of my alignment sequences, but I get the "Number of residues in the alignment and pdb files are different" message (the difference is actually the two additional residues). The FAQ (question #17) suggests running a simple script to see what Modeller is reading from the PDB and, running it on PDB 1HRA, I get this: >P1;1hra structureX:1hra: 1 :A:+80 :A:MOL_ID 1; MOLECULE RETINOIC ACID RECEPTOR; CHAIN A; ENGINEERED YES:MOL_ID 1; ORGANISM_SCIENTIFIC HOMO SAPIENS; ORGANISM_COMMON HUMAN; ORGANISM_TAXID 9606:-1.00: -1.00 MPRVYKPCFVCQDKSSGYHYGVSACEGCKGFFRRSIQKNMIYTCHRDKNCVINKVTRNRCQYCRLQKCFEVGMSK ESVRN* I.e., for some reason modeller only reads the amino acids and not the two zinc ions, and the last two dots are missing. What should I do to read the heteroatoms? Lucas

2 1

Best workflow to modeling a disordered loop in a complicated structure
by Zachary W Carpenter 11 Feb '14

11 Feb '14

Hello, I have used modeller to attempt to model a loop region in a tetrameric structure. To speed up computational time I've started with just the dimer ( which I think also makes sense biologically ) and ran the following script with an .ali file I created. The region of interest is a disordered loop AA~399-417. Following this I used the evaluate model script to pick the best model, which I then used as the template input for the loop.py script. I run the loop.py script at many iterations and then when it finishes I take the top 20 models (out of 1000) by DOPE score and attempt to analyze them. Any advice on how I could improve my method would be great.. I really am at the edge of my understanding here. I've edited the name of the protein. I can include it as soon as I can verify that this can stay anonymous. Thanks!!! Scripts +files : *First )using a template pdb and .ali file I made myself of the Dimer (BIO_headers in chimera and deletion of one dimer of the tetramer) * *Input : * >P1;ABCD structureX:ABCD.pdb: 3 :A:+990:B:::-1.00:-1.00 --TSWSDRLQNAADMPANMDKHALKKYRREAYHRVFVNRSLAMEKIKCFGFNMDYTLAVYKSPEYESLGFELTVE RLVSIGYPQELLSFAYDSTFPTRGLVFDTLYGNLLKVDAYGNLLVCAHGFNFIRGPETREQYPNKFIQRDDTERF YILNTLFNLPETYLLACLVDFFTNCPRYTSCETGFKDGDLFMSYRSMFQDVRDAVDWVHYKGSLKEKTVENLEKY VVKDGKLPLLLSRMKEVGKVFLATNSDYKYTDKIMTYLFDFPHGPKPGSSHRPWQSYFDLILVDARKPLFFGEGT VLRQVDTKTGKLKIGTYTGPLQHGIVYSGGSSDTICDLLGAKGKDILYIGDHIFGDILKSKKRQGWRTFLVIPEL AQELHVWTDKSSLFEELQSLDIFLAS----------------SIQRRIKKVTHDMDMCYGMMGSLFRSGSRQTLF ASQVMRYADLYAASFINLLYYPFSYLFRAAHVLMPHES/--TSWSDRLQNAADMPANMDKHALKKYRREAYHRVFVNRSLAMEKIKCFGFNMDYTLAVYKSPEYESLGFELTVE RLVSIGYPQELLSFAYDSTFPTRGLVFDTLYGNLLKVDAYGNLLVCAHGFNFIRGPETREQYPNKFIQRDDTERF YILNTLFNLPETYLLACLVDFFTNCPRYTSCETGFKDGDLFMSYRSMFQDVRDAVDWVHYKGSLKEKTVENLEKY VVKDGKLPLLLSRMKEVGKVFLATNSDYKYTDKIMTYLFDFPHGPKPGSSHRPWQSYFDLILVDARKPLFFGEGT VLRQVDTKTGKLKIGTYTGPLQHGIVYSGGSSDTICDLLGAKGKDILYIGDHIFGDILKSKKRQGWRTFLVIPEL AQELHVWTDKSSLFEELQSLDIFLAS----------------SIQRRIKKVTHDMDMCYGMMGSLFRSGSRQTLF ASQVMRYADLYAASFINLLYYPFSYLFRAAHVLMPHES* >P1;X sequence:ABCD: : : : ::: 0.00: 0.00 MSTSWSDRLQNAADMPANMDKHALKKYRREAYHRVFVNRSLAMEKIKCFGFDMDYTLAVYKSPEYESLGFELTVE RLVSIGYPQELLSFAYDSTFPTRGLVFDTLYGNLLKVDAYGNLLVCAHGFNFIRGPETREQYPNKFIQRDDTERF YILNTLFNLPETYLLACLVDFFTNCPRYTSCETGFKDGDLFMSYRSMFQDVRDAVDWVHYKGSLKEKTVENLEKY VVKDGKLPLLLSRMKEVGKVFLATNSDYKYTDKIMTYLFDFPHGPKPGSSHRPWQSYFDLILVDARKPLFFGEGT VLRQVDTKTGKLKIGTYTGPLQHGIVYSGGSSDTICDLLGAKGKDILYIGDHIFGDILKSKKRQGWRTFLVIPEF AQELHVWTDKSSLFEELQSLDIFLAELYKHLDSSSNERPDISSIQRRIKKVTHDMDMCYGMMGSLFRSGSRQTLF ASQVMRYADLYAASFINLLYYPFSYLFRAAHVLMPHES/MSTSWSDRLQNAADMPANMDKHALKKYRREAYHRVFVNRSLAMEKIKCFGFDMDYTLAVYKSPEYESLGFELTVE RLVSIGYPQELLSFAYDSTFPTRGLVFDTLYGNLLKVDAYGNLLVCAHGFNFIRGPETREQYPNKFIQRDDTERF YILNTLFNLPETYLLACLVDFFTNCPRYTSCETGFKDGDLFMSYRSMFQDVRDAVDWVHYKGSLKEKTVENLEKY VVKDGKLPLLLSRMKEVGKVFLATNSDYKYTDKIMTYLFDFPHGPKPGSSHRPWQSYFDLILVDARKPLFFGEGT VLRQVDTKTGKLKIGTYTGPLQHGIVYSGGSSDTICDLLGAKGKDILYIGDHIFGDILKSKKRQGWRTFLVIPEF AQELHVWTDKSSLFEELQSLDIFLAELYKHLDSSSNERPDISSIQRRIKKVTHDMDMCYGMMGSLFRSGSRQTLF ASQVMRYADLYAASFINLLYYPFSYLFRAAHVLMPHES* # Homology modeling by the automodel class # # Demonstrates how to build multi-chain models, and symmetry restraints # from modeller import * from modeller.automodel import * # Load the automodel class log.verbose() # Override the 'special_restraints' and 'user_after_single_model' methods: class MyModel(automodel): def special_restraints(self, aln): # Constrain the A, B, C and D chains to be identical s1 = selection(self.chains['A']).only_atom_types('CA') s2 = selection(self.chains['B']).only_atom_types('CA') self.restraints.symmetry.append(symmetry(s1, s2, 1)) def user_after_single_model(self): # Report on symmetry violations greater than 1A after building # each model: self.restraints.symmetry.report(1) env = environ() # directories for input atom files env.io.atom_files_directory = ['.', '../atom_files'] # Be sure to use 'MyModel' rather than 'automodel' here! a = MyModel(env, alnfile = 'WT_dimer_NoCterm.ali' , # alignment filename knowns = 'ABCD', # codes of the templates sequence = 'ABCD') # code of the target a.starting_model= 1 # index of the first model a.ending_model = 20 # index of the last model # (determines how many models to calculate) a.make() # do homology modeling # # class MyModel(automodel): *2nd) Run eval model and pick top model to use as input for the loop.py script below --->* # Loop refinement of an existing model from modeller import * from modeller.automodel import * log.verbose() env = environ() # directories for input atom files env.io.atom_files_directory = './:../atom_files' # Create a new class based on 'loopmodel' so that we can redefine # select_loop_atoms (necessary) class MyLoop(loopmodel): # This routine picks the residues to be refined by loop modeling def select_loop_atoms(self): # 10 residue insertion return selection(self.residue_range('399:A', '417:A')) m = MyLoop(env, inimodel='ABCD_Top_Model.pdb', # initial model of the target sequence='ABCD') # code of the target m.loop.starting_model= 1 # index of the first loop model m.loop.ending_model = 1000 # index of the last loop model m.loop.md_level = refine.slow # loop refinement method; this yields # models quickly but of low quality; # use refine.slow for better models m.make() *3rd) Run model_energies script* Then I rank them using this perl one liner ( just cats and displays the dope and the file ) can run in any unix terminal. cat model*.log |perl -lane' print if /DOPE\ score/ || /ABCD.BL/'| perl -pe 's/ +/\t/g'|perl -pe 's/:\t//g' |cut -f3,4|perl -pe 's/\n/\t/'|perl -pe 's/(\d)\t(ABCD)/$1\n$2/g'|sort -k2,2n|head -20 *4) Take top 20 models and analyze * Any help with where I'm doing things wrong or how I can improve would be awesome. How many models in loop.py would you recommend to get a good answer? Is there away to interpret the models created as a whole ....for example say something like "60% of the time residue A is next to residue B and 20% of the time its over in this pocket , therefor".... etc My major aim is to compare different ligand bound crystal states, or mutation states, and show that the models created are different between them. Then use this to design wet lab experiments to test. Thanks again

1 0

Using 2 templates, one to just model the protein and the other to just model the ligand
by Juan Munoz-Garcia 06 Feb '14

06 Feb '14

Dear MODELLER users, I want to do homology modelling of a protein based on just one template because it is so far the only template of a particular protein state of its family I' m interest in. However, that template does not contain the natural ligand of the target. Then, I want to use another template of the protein-ligand(natural) complex of interest but just to say Modeller to use the ligand coordinates. So, it is a multi-alignment using the HET option but I don't know how to say Modeller that I just want to model the target sequence based on just one template and considering the ligand of the other template. Is this possible? I look forward to hearing from you. Juan C. Munoz-Garcia, PhD

3 3

Reg modeller
by Manjula Mummadisetti 05 Feb '14

05 Feb '14

I would like to understand the link below. I am not quite sure what each of those numbers mean, although the link does mention, can somebody clearly explain what those numbers are in restraints panel? http://salilab.org/modeller/8v1/manual/node157.html#SECTION:atomid Also, I would like to restrain N-N distance rather than C-alpha C-alpha. How do I make those changes? Thanks a lot Manju

2 1

model.patch
by Roland Rivera 31 Jan '14

31 Jan '14

Hi there, I'm having some trouble using the model_patch() command to create a disulfide bond. If I use the following alignment file: >P1;1PRX_edited structure:1PRX_edited.pdb:2 :B:224 :B: : : : -PGGLLLGDVAPNFEANTTVGRIRFHDFLGDSWGILFSHPRDFT--.-----------PE FAKRNVKLIALSIDS---------------SEEPTEKLPFPIIDDRNRELAILLGMLDPA E-----MPVTARVVFVFGPDKKLKLSILYPATTGRNFDEILRVVISLQLTAEKRVATPVD WKDGDSVMVLPTIPEEEAKKLFPKGVFTKELPSGKKYLRYTPQP* >P1;1PRX_WT sequence:1PRX_WT:2 :A:224 :A: : : : MPGGLLLGDVAPNFEANTTVGRIRFHDFLGDSWGILFSHPRDFTPVCTTELGRAAKLAPE FAKRNVKLIALSIDSVEDHLAWSKDINAYNCEEPTEKLPFPIIDDRNRELAILLGMLDPA EKDEKGMPVTARVVFVFGPDKKLKLSILYPATTGRNFDEILRVVISLQLTAEKRVATPVD WKDGDSVMVLPTIPEEEAKKLFPKGVFTKELPSGKKYLRYTPQP* in conjunction with the following script: # Modeling a cysteine bond in a new model from modeller import * from modeller.automodel import * # Load the automodel class from modeller.scripts import complete_pdb log.verbose() # Verbose error description - helps with troubleshooting! env = environ() # Defines the environment # directories for input atom files env.io.atom_files_directory = ['.','../atom-files'] #Enables MODELLER to read HETATM entries from a PDB file env.io.hetatm = True env.libs.topology.read(file='$(LIB)/top_heav.lib') env.libs.parameters.read(file='$(LIB)/par.lib') # Make the disulfide bond: def patches (mdl): mdl.patch(residue_type='DISU', residues=(mdl.residues['47:A'], mdl.residues['91:A'])) # Read the sequence: code = '1PRX' mdl = complete_pdb(env, code, special_patches=patches) # Create the stereochemical restraints: sel = selection(mdl) mdl.restraints.make(sel, restraint_type='stereo', spline_on_site=False) # Calculate the energy to test the disulfide: sel.energy() a = automodel(env, alnfile = '1PRX_edited-1PRX_WT.ali', # alignment filename knowns = '1PRX_edited', # codes of the templates sequence = '1PRX_WT', # code of the target assess_methods=(assess.DOPE, assess.GA341)) a.starting_model= 1 # index of the first model a.ending_model = 1 # index of the last model # (together, they determine how many models to calculate) a.make() # do homology modeling It generates a PDB file that has the desired sequence, but no disulfide bond. If I use the following alignment: >P1;1PRX structure:1PRX.pdb:2 :B:224 :B: : : : -PGGLLLGDVAPNFEANTTVGRIRFHDFLGDSWGILFSHPRDFTPV.TTELGRAAKLAPE FAKRNVKLIALSIDSVEDHLAWSKDINAYNSEEPTEKLPFPIIDDRNRELAILLGMLDPA E-----MPVTARVVFVFGPDKKLKLSILYPATTGRNFDEILRVVISLQLTAEKRVATPVD WKDGDSVMVLPTIPEEEAKKLFPKGVFTKELPSGKKYLRYTPQP* >P1;1PRX_WT sequence:1PRX_WT:2 : :224 : : : : : MPGGLLLGDVAPNFEANTTVGRIRFHDFLGDSWGILFSHPRDFTPVCTTELGRAAKLAPE FAKRNVKLIALSIDSVEDHLAWSKDINAYNCEEPTEKLPFPIIDDRNRELAILLGMLDPA EKDEKGMPVTARVVFVFGPDKKLKLSILYPATTGRNFDEILRVVISLQLTAEKRVATPVD WKDGDSVMVLPTIPEEEAKKLFPKGVFTKELPSGKKYLRYTPQP* Combined with the following script: # Modeling a cysteine bond in a new model from modeller import * from modeller.automodel import * # Load the automodel class from modeller.scripts import complete_pdb log.verbose() # Verbose error description - helps with troubleshooting! env = environ() # Defines the environment # directories for input atom files env.io.atom_files_directory = ['.','../atom-files'] #Enables MODELLER to read HETATM entries from a PDB file env.io.hetatm = True env.libs.topology.read(file='$(LIB)/top_heav.lib') env.libs.parameters.read(file='$(LIB)/par.lib') # Special distance restraints class MyModel (automodel): # Make the disulfide bond: def patches (mdl): mdl.patch(residue_type='DISU', residues=(mdl.residues['47:'], mdl.residues['91:'])) # Read the sequence: code = '1PRX' mdl = complete_pdb(env, code, special_patches=patches) # Create the stereochemical restraints: sel = selection(mdl) mdl.restraints.make(sel, restraint_type='stereo', spline_on_site=False) # Calculate the energy to test the disulfide: sel.energy() a = MyModel(env, alnfile = '1PRX-1PRX_WT.ali', # alignment filename knowns = '1PRX', # codes of the templates sequence = '1PRX_WT', # code of the target assess_methods=(assess.DOPE, assess.GA341)) a.starting_model= 1 # index of the first model a.ending_model = 1 # index of the last model # (together, they determine how many models to calculate) a.make() # do homology modeling It returns a KeyError: 'No such residue'. I researched this in the answer archives and found the suggestion to create a .ini file to determine the model residue numbers, which are exactly correspondent to those found in the script. I would appreciate any help deciphering what exactly I'm doing wrong. Thanks, -- Roland F. Rivera-Santiago Ph.D. Candidate Biochemistry & Molecular Biophysics Graduate Group Speicher Lab, Wistar Institute E-mail: roland.f.rivera.santiago(a)gmail.com

2 1

Re: [modeller_usage] problem in generating models
by Nupur Bansal 29 Jan '14

29 Jan '14

Hi, Thanks for replying. Now, after making that change, I am getting the following error: Alignment sequence not found in PDB file: 1 atom_files/2hzv.pdb (You didn't specify the starting and ending residue numbers and chain IDs in the alignment, so Modeller tried to guess these from the PDB file.) But, I don’t know how to specify the starting and ending residue numbers. On Jan 29, 2014, at 12:11 PM, Nupur Bansal <nupur.b298(a)gmail.com> wrote: > Hi, > > I am new to Modeller and has been trying to use it for generating models for my system. > > I have done sequence alignment using Clustalw in .pir format. But, the error I am getting is : > check_a_282E> Unable to read structural information for alignment > sequence: 1 2hzv > This command requires structures. > > > My alignment file is : 2ca9_2hzv.ali > > >P1;2hzv > NMR:2hzv.pdb:FIRST:@ END:::::: > --MQRVTITLDDDLLETLDSLSQRRGYNNRSEAIRDILRSALAQEATQQHGTQG---FAV > LSYVYEHEKRDLASRIVSTQHHHHDLSVATLHVHINHDDCLEIAVLKGDMGDVQHFADDV > IAQRGVRHGHLQCLP--KEDMQRVTITLDDDLLETLDSLSQRRGYNNRSEAIRDILRSAL > AQEATQQHGTQG---FAVLSYVYEHEKRDLASRIVSTQHHHHDLSVATLHVHINHDDCLE > IAVLKGDMGDVQHFADDVIAQRGVRHGHLQCLP----KEDMQRVTITLDDDLLETLDSLS > QRRGYNNRSEAIRDILRSALAQEATQQHGTQG---FAVLSYVYEHEKRDLASRIVSTQHH > HHDLSVATLHVHINHDDCLEIAVLKGDMGDVQHFADDVIAQRGVRHGHLQCLP--KEDMQ > RVTITLDDDLLETLDSLSQRRGYNNRSEAIRDILRSALAQEATQQHGTQG---FAVLSYV > YEHEKRDLASRIVSTQHHHHDLSVATLHVHINHDDCLEIAVLKGDMGDVQHFADDVIAQR > GVRHGHLQCLPKED--* > > >P1;2CA9_apo > sequence:2CA9_apo:FIRST:@ END:::::: > DSIIRFSVSLQQNLLDELDNRIIKNGYSSRSELVRDMIREKLVEDNWAEDNPNDESKIAV > LVVIYDHHQRELNQRMIDIQHASGTHVLCTTHIHMDEHNCLETIILQGNSFEIQRLQLEI > GGLRGVKFAKLTKASSFEYSIIRFSVSLQQNLLDELDNRIIKNGYSSRSELVRDMIREKL > VEDNWAEDNPNDESKIAVLVVIYDHHQRELNQRMIDIQHASGTHVLCTTHIHMDEHNCLE > TIILQGNSFEIQRLQLEIGGLRGVKFAKLTKASSFEYNDSIIRFSVSLQQNLLDELDNRI > IKNGYSSRSELVRDMIREKLVEDNWAEDNPNDESKIAVLVVIYDHHQRELNQRMIDIQHA > SGTHVLCTTHIHMDEHNCLETIILQGNSFEIQRLQLEIGGLRGVKFAKLTKASSFEYSII > RFSVSLQQNLLDELDNRIIKNGYSSRSELVRDMIREKLVEDNWAEDNPNDESKIAVLVVI > YDHHQRELNQRMIDIQHASGTHVLCTTHIHMDEHNCLETIILQGNSFEIQRLQLEIGGLR > GVKFAKLTKASSFEYN* > > > I am also attaching my pdb file for reference. > > Any help would be appreciated. > > Thanks, > Nupur > <2hzv.pdb>

2 1

Refinement step questions
by gbanoczi＠mail.bme.hu 21 Jan '14

21 Jan '14

2 1