- modeller_usage - salilab.org

How define restrains between model and template?
by Mario A Bianchet 02 Apr '18

02 Apr '18

Hi, The bend of the modeled three strands long sheet differs too much from the template. Independently that perhaps that have meaning. I want to make closer to the template. It is not clear (or I am missing it ) if is possible restrains between model and template. Thanks in advance Mario -- ========================================================================= Mario A. Bianchet, Ph.D. | bianchet(a)jhmi.edu Assistant Professor of Neurology, | http://abaco.med.jhmi.edu/~bianchet Biophysics & Biophysical Chemistry | 725 N.Wolfe st, WBSB 608B | Baltimore, Maryland, 21205, USA Structural Enzymology | Phones: +1 (443) 591-7583 (cell) & Thermodynamics Group | +1 (410) 614-8221(work) Johns Hopkins School of Medicine | fax: +1 (410) 955-0637 =========================================================================

2 3

Re: [modeller_usage] How define restrains between model and template?
by Mario A Bianchet 28 Mar '18

28 Mar '18

Hi David Thanks. However, my problem is not a different 2ry structure. The well-modeled 3-strand sheet is bend too much compared to the template. I am looking for something like this: # Use a harmonic potential and X-Y distance group. rsr.add(forms.gaussian(group=physical.xy_distance,feature=features.distance(at['CA:77'],at['CA:418']),mean=1.0, stdev=0.5)) Is there a way to select atoms in the template something like at['{structureX/known...}:CA:47'] to restrain the model? Thanks Mario On 3/28/18 12:34 PM, David Gae wrote: > actually I have minor correction. > > modeller also has restraints for sheets as well. > > # for anti-parallel > rsr.add(secondary_structure.sheet(at['N:1'], at['O:14'], > sheet_h_bonds=-5)) > or > > # for parallel > rsr.add(secondary_structure.sheet(at['N:1'], at['O:9'], > # sheet_h_bonds=5)) > Sincerely, > David > >> On Mar 28, 2018, at 9:30 AM, David Gae <ddgae(a)ucdavis.edu >> <mailto:ddgae@ucdavis.edu>> wrote: >> >> Dear Dr. Bianchet, >> >> This may be a possible solution to keep the model similar to the >> template. You can maintain the secondary characteristics using the >> command below. >> For example, >> # beta-strand residue 1-6 should be an beta-strand >> rsr.add(secondary_structure.strand(self.residue_range('1:', '6:'))) >> >> Regards, >> David >> >> >>> On Mar 28, 2018, at 8:35 AM, Mario A Bianchet <bianchet(a)jhmi.edu >>> <mailto:bianchet@jhmi.edu>> wrote: >>> >>> Hi, >>> >>> The bend of the modeled three strands long sheet differs too much >>> from the template. Independently that perhaps that have meaning. I >>> want to make closer to the template. >>> It is not clear (or I am missing it ) if is possible restrains >>> between model and template. >>> >>> Thanks in advance >>> >>> Mario >>> -- >>> ========================================================================= >>> Mario A. Bianchet, Ph.D. |bianchet(a)jhmi.edu >>> Assistant Professor of Neurology, |http://abaco.med.jhmi.edu/~bianchet >>> Biophysics & Biophysical Chemistry | 725 N.Wolfe st, WBSB 608B >>> | Baltimore, Maryland, 21205, USA >>> Structural Enzymology | Phones: +1 (443) 591-7583 (cell) >>> & Thermodynamics Group | +1 (410) 614-8221(work) >>> Johns Hopkins School of Medicine | fax: +1 (410) 955-0637 >>> ========================================================================= >>> _______________________________________________ >>> modeller_usage mailing list >>> modeller_usage(a)salilab.org <mailto:modeller_usage@salilab.org> >>> https://salilab.org/mailman/listinfo/modeller_usage >> >

1 0

How can I change the outputted structure names?
by Ahmad Abdelzaher 21 Mar '18

21 Mar '18

Hello, I'm following this tutorial here: https://salilab.org/modeller/tutorial/basic.html. It's working great, except for not being able to control the names of the outputted structures. Any attempt to rename the alignment codes would of course give errors like "Protein specified in ALIGN_CODES(i) was not found in the alignment file". Is there a way to control the name of the outputted structures? Regards.

2 1

Evaluate protein-peptide complexes using modeller
by SIMON KØNIG 20 Mar '18

20 Mar '18

Dear all It has been pointed out to me that the evaluation method in modeller (molpdf, GA341 and DOPE) are not applicable for models of protein-peptide complexes. Is this the case? If it is, does anybody then have other recommendation for evaluating/assessing models of complexes ? Best regards Simon

2 1

residue selection by sequence
by Shyam Saladi 09 Mar '18

09 Mar '18

Hello, It is possible to specify a selection by the sequence of their residues? Usually, I use model.residue_range(), but this can be tedious. Ideally, I'd like to do something like model.residue_sequence('NNPLGTPGTS') to select all occurrences of that sequence. Thanks, Shyam

2 1

Help with rigid body minimization of protein-protein complex
by mromero＠bsc.es 13 Feb '18

13 Feb '18

Hello, I have been testing whether I could use Modeller to minimize a protein-protein complex. As a first test, I've tried to perform a rigid body minimization of the bound structures of receptor and ligand. Even though the starting structure has a Lig-RMSD of 3 A from the native complex structure (i.e., after superimposing the receptor subunits and measuring the CA RMSD between the ligand subunits), the final minimized structure does not get closer to the native complex but goes farther away, with Lig-RMSD of 12 A. I've tried with several near-native conformations and different protein complexes with similar results. The script I am running is the following: from modeller import * from modeller.scripts import complete_pdb from modeller.optimizers import conjugate_gradients, molecular_dynamics, actions env = environ() env.io.atom_files_directory = ['../atom_files'] env.edat.dynamic_sphere = True env.libs.topology.read(file='$(LIB)/top_heav.lib') env.libs.parameters.read(file='$(LIB)/par.lib') code = "test.pdb" mdl = complete_pdb(env, code) mdl.write(file=code+'.ini') # Select rigid bodies # Receptor rigid body atmsel_a = mdl.chains['A'].atoms r_a = rigid_body(atmsel_a) # Ligand rigid body atmsel_b = mdl.chains['B'].atoms r_b = rigid_body(atmsel_b) # Select all atoms: atmsel = selection(mdl) # Generate the restraints: mdl.restraints.rigid_bodies.append(r_a) mdl.restraints.rigid_bodies.append(r_b) mdl.restraints.make(atmsel, restraint_type='STEREO', spline_on_site=False) mdl.restraints.write(file=code+'.rsr') mpdf = atmsel.energy() # Create optimizer objects and set defaults for all further optimizations cg = conjugate_gradients(output='REPORT') md = molecular_dynamics(output='REPORT') # Open a file to get basic stats on each optimization trcfil = open(code+'.D00000001', 'w') # Run CG on the all-atom selection; write stats every 5 steps cg.optimize(atmsel, max_iterations=20, actions=actions.trace(5, trcfil)) # Run MD; write out a PDB structure (called '1fas.D9999xxxx.pdb') every # 10 steps during the run, and write stats every 10 steps md.optimize(atmsel, temperature=300, max_iterations=20, actions=[actions.write_structure(10, code+'.D9999%04d.pdb'), actions.trace(10, trcfil)]) # Finish off with some more CG, and write stats every 5 steps cg.optimize(atmsel, max_iterations=20, actions=[actions.trace(5, trcfil)]) mpdf = atmsel.energy() mdl.write(file=code+'.B') I have also played with several minimization schemes, with and without MD, but the final minimized model does not get closer to the native complex. Thanks in advance, Miguel http://bsc.es/disclaimer

2 1

How to make better conformational sampling of models using MODELLER?
by Mahesh VELUSAMY 30 Jan '18

30 Jan '18

Hi All, As per the suggestion given in the modeller that changing random seed would helps in producing different models. On this note I have tried to generate 100 models from single template and using different random seeds. It was mainly performed to know how changing random seed influencing on conformational sampling ?. Following the experiment I had observed that there was notable fluctuations on the sidechain orientations. But when I aligned all the models against top model (DOPE SCORE) there was no variations on the backbone conformation and also in terms of RMSD. Following these observations I have a couple of questions for the modeller experts as follows: When modeller says changing random seed helps in finding different models, does it meant for sidechain orientations or backbone? Is it possible to produce models in modeller that will have more variations on the backbone and in terms of RMSD ? Thanking you in advance

2 3

modeller error; " ValueError: too many arguments to unpack"
by Sourav Pal 19 Jan '18

19 Jan '18

I am going to generate a monomer from one single Template (tseq1) by modeller. But there one problem is arised. I have a only one template. According to the basic tutorial of modeller.. I made alignment file like this, (where tseq1 is single template file)from modeller import *env = environ()aln = alignment(env)for (pdb, chain) in (('tseq1', 'A')):m = model(env, file=pdb, model_segment=('FIRST:'+chain, 'LAST:'+chain))aln.append_model(m, atom_files=pdb, align_codes=pdb+chain)aln.malign()aln.malign3d()aln.compare_structures()aln.id_table(matrix_file='family.mat')env.dendrogram(matrix_file='family.mat', cluster_cut=-1.0)That time it was failed to run by showing the following error"ValueError:too many values to unpack" .can anyone tell me where is the problem? Thank you.Best Regards. Sourav Pal JRF

2 1

Modelling of Homotrimeric protein without biological assembly
by Sourav Pal 17 Jan '18

17 Jan '18

Actually I want to build a homology model of a human homotrimeric protein from a FASTA sequences. After BLAST I got a template which is 91% identical but it is originated from Mouse organism. and it is cryo EM structure.From there, I build a single monomeric chain of that human protein by Homology modelling. But it doesn't contain any biological assembly or asymmetric unit. So that, i am unable to build its trimeric structure. so, can you help me in two sections;1) how to build its homotrimeric structure?2) is there any way to generate the biological assembly coordinates from homology protein?Thanks in advance.Best regards.

1 0

Optimization of an existing loop model
by Claire Marks 17 Jan '18

17 Jan '18

Hi! I would like to use MODELLER to optimize some models that I have generated as part of a loop prediction protocol. The loop conformations were built onto crystal structures downloaded onto the PDB; ideally I only want the loop region to be optimized. I thought the code given on the following web page would be what I should use, however it doesn't seem to do what I want: https://salilab.org/modeller/manual/node36.html This still generates the initial 'straight' conformation that the MODELLER loop modelling algorithm starts from; it doesn't start from the loop conformation in the file I provided. Is anyone able to shed some light on this for me? Claire

2 1