- phenixbb - phenix-online.org

Using a parameter file in GUI
by Subhani Bandara 24 Feb '12

24 Feb '12

Hi everyone, I need to input a parameter file describing the bond lengths of a metal chelator complex in my protein. I used a text file containing the commands with .params extension. when I add this file in GUI inputs in refinement, it gives a message "Phenix did not recognize the file type". What is the file type I should use? The file contained commands similar to the following with my metal and ligand atom names. refinement.geometry_restraints.edits { zn_selection = chain X and resname ZN and resid 200 and name ZN his117_selection = chain X and resname HIS and resid 117 and name NE2 asp130_selection = chain X and resname ASP and resid 130 and name OD1 bond { action = *add atom_selection_1 = $zn_selection atom_selection_2 = $his117_selection symmetry_operation = None distance_ideal = 2.1 sigma = 0.02 slack = None } bond { action = *add atom_selection_1 = $zn_selection atom_selection_2 = $asp130_selection symmetry_operation = None distance_ideal = 2.1 sigma = 0.02 slack = None } } Then I used command line to run "phenix.refine model.pdb data.mtz elbow.cif restraints_edits.params" it again gave error saying .params is an unknown file type. Am I doing anything wrong and is there any other way to fix the distances of the metal chelator complex. Thanks Subhani

7 35

importing projects problem
by Haytham Wahba 24 Feb '12

24 Feb '12

i moved my project from computer to another i did import project under Files. i got my projects with all jop history but when i try to restor jop i get this report (see below) after that i see only .eff file under the input file, No .Sca or the pdb i used before. even the result of refinment is not there and refinment setting which i did not use. i tried to change the path in .eff file to the new location for .sca and .pdb but i did not help could you help me in that. thank you Haytham (((ImportError : No module named f_model_info Traceback (most recent call last): File "/usr/local/phenix-1.7-650/phenix/wxGUI2/Home/Jobs.py", line 327, in OnDoubleClick self.list_frame.restore_job(job) File "/usr/local/phenix-1.7-650/phenix/wxGUI2/Home/Jobs.py", line 241, in restore_job app.restore_job(self.project.get_id(), current_job) File "/usr/local/phenix-1.7-650/phenix/wxGUI2/App.py", line 544, in restore_job use_template_gui=use_template_gui) File "/usr/local/phenix-1.7-650/phenix/wxGUI2/App.py", line 501, in restore_session frame.restore_session(output_dir, config_file, result) File "/usr/local/phenix-1.7-650/phenix/wxGUI2/Programs/Refine/Main.py", line 473, in restore_session (result, last_status_file) = self._load_result_files(output_dir) File "/usr/local/phenix-1.7-650/phenix/wxGUI2/Programs/Refine/Main.py", line 467, in _load_result_files refine_result = easy_pickle.load(result_file) File "/usr/local/phenix-1.7-650/cctbx_project/libtbx/easy_pickle.py", line 21, in load return cPickle.loads(_open(file_name, "rb").read()) ImportError: No module named f_model_info)))

2 1

question about phenix rosetta with two chains
by Elena Seiradake 23 Feb '12

23 Feb '12

Dear Phenixians, I would like to use mr_rosetta to re-build a model with two copies per AU. Although the two copies have the same sequence, they adopt different conformations and one is partially disordered. So I don't want to impose NCS. Could someone please send me an example input file for running mr_rosetta with phenix for multiple chains of the same protein? I think this probably would have to contain "fragment_files_chain_list" and "fragment_files_3_mer_by_chain" etc. (?) I have tried including the above keys in my script, in various ways, but I always get the error: "Sorry, the file path/test_1.pdb representing group 1 has no atoms? This could indicate a problem with sequence alignment in rebuilding" When I run it with just one of the two chains in my pdb (using simple "fragment_file ="), it rebuilds happily. Also, when use both chains and impose NCS=2 it works, but then one of the two chains is forced into the wrong conformation. Any advice on how to change my script would be welcome! Below is part of my simple script, the one that runs ok, using "fragment_file =". seq_file= Seq.dat \ data=test.mtz \ search_models=test.pdb \ already_placed=True \ fragment_files= frag_A_3.gz \ fragment_files= frag_A_9.gz \ rescore_mr.relax=False \ rosetta_models=20 \ ncs_copies= 1\ Cheers Elena

2 1

composite omit map
by Ursula Schulze-Gahmen 22 Feb '12

22 Feb '12

I tried to calculate a composite omit map in autobuild for 3 complexes in the au. The calculation finished without any problems, but the maps look not right. There is density for the molecule, but there is a lot of density in the solvent areas, which I never saw when I calculated regular 2fo-fc maps or fo-fc maps. I also have calculated omit maps in autobuild before for single problem region and never saw anything like this. Any idea what could have gone wrong? Where should I start looking? Ursula -- Ursula Schulze-Gahmen, Ph.D. Assistant Researcher UC Berkeley, QB3 356 Stanley Hall #3220 Berkeley, CA 94720-3220

2 1

autoMR fixed ensemble error
by Alexander Batyuk 22 Feb '12

22 Feb '12

Dear all, I have the following error while trying to put fixed ensemble into AutoMR: Adding Ensemble c1 as fixed solution with: Euler angles: [359.98, 0.0, 0.0] Coords (orthogonal A): [-1.0, 0.476, -1.001] *************ERROR ENDING ******************* **************************************** AutoMR Input failed Python argument types in SOLU.addSOLU_6DIM_ENSE(InputMR_AUTO, str, list, bool, list, float, bool, bool, bool) did not match C++ signature: addSOLU_6DIM_ENSE(phaser::SOLU {lvalue}, std::string, scitbx::vec3<double>, bool, scitbx::vec3<double>, double, bool, bool, bool, bool, scitbx::vec3<double>, double) **************************************** *************ERROR ENDING ******************* What could be the reason for that? I'm running dev-986 on OS X 10.6.8. Stable 1.7.3 gives the same error. Best wishes, Alex -- Alex Batyuk The Plueckthun Lab www.bioc.uzh.ch/plueckthun

4 5

Re: [phenixbb] CCN Jan 2012 article : find_alt_origin_sym_mate
by R.D. Oeffner 22 Feb '12

22 Feb '12

Dear Bryan, The factor 0.9 is a purely empirical term. It serves two functions. One is to soften the minimas when atoms coalesce and to also prevent taking the logarithm of a negative number. The closer it gets to zero the more accurately it pinpoints the global minimum when the MLAD score is used for identifying the best superposition betweent two molecules. But if it becomes two small chances are that the landscape of MLAD values will be dotted with many very deep local minima which might lead to wrong results. The square boxes under the min, max symbols are typos to be ignored, sorry. There are currently no plans to publish further details about find_alt_origin_sym_mate. Regards, Robert Oeffner On Wed, Feb 15, 2012 at 4:33 PM, Bryan Lepore <bryanlepore at gmail.com> wrote: > will the find_alt_origin_sym_mate article in CCN Jan. 2012 be > described in any other publications? ... e.g. Acta Cryst., because I am interested in particular : 1. what is the significance of the factor of 0.9 2. is the MLAD a novel formula also, the 'min' and 'max' on p. 6 have small boxes under them [*] as if missing a character - are these integrals? -Bryan [*] tried evince and Adobe Reader on Linux, Adobe Acrobat on OSX 10.5.8 -- Robert Oeffner, Ph.D. Research Associate Department of Haematology, University of Cambridge Cambridge Institute of Medical Research Wellcome Trust / MRC Building, Hills Road, Cambridge, CB2 0XY www-structmed.cimr.cam.ac.uk, tel:01223763234, mobile:07712 887162

5 5

error in running autobuild
by intekhab alam 21 Feb '12

21 Feb '12

hi i am working with phenix version 1.7.1-743. after running auto mr i tried to run autobuild but it failed and generated the following message. i am unable to run autobuild now. kindly help. OError : [Errno 5] Input/output error Traceback (most recent call last): File "/usr/local/phenix-1.7.1-743/phenix/wxGUI2/Runtime.py", line 172, in OnUpdate self.propagate_event(event, "OnUpdate") File "/usr/local/phenix-1.7.1-743/phenix/wxGUI2/Runtime.py", line 228, in propagate_event fn(event) File "/usr/local/phenix-1.7.1-743/phenix/wxGUI2/Programs/Phaser/Output.py", line 194, in OnUpdate self.process_callback(callback.message, callback.data) File "/usr/local/phenix-1.7.1-743/phenix/wxGUI2/Programs/Phaser/Output.py", line 209, in process_callback print "CALLBACK: %s %s" % (str(message), str(data)) IOError: [Errno 5] Input/output erro -- INTEKHAB ALAM LABORATORY OF STRUCTURAL BIOINFORMATICS KOREA UNIVERSITY, SEOUL

3 2

FFT map coefficients only for certain chains
by Zhang yu 20 Feb '12

20 Feb '12

Dear Phenixers, Is that possible to generate map coefficients only for certain chains? For example, I have two chains, A and B, and I would like to output a map file only contains coefficients for chain A. The "isomesh" command in Pymol could generate similar images. But my purpose is not for presentation, I need a map file only contains coefficient for certain chains. In the interface of "FFT" tool in Phenix or CCP4, there is an option to include a PDB file and define atom selections. It describe that "If a PDB is supplied, the output map will cover the model plus a buffer on all sides. The atom selection parameters can be used to specify a smaller region" . If I define the selection as chain A when I run the FFT, the output map still covers a rectangular block containing chain A, instead of regions only surrounding chain A. Thanks. Yu Zhang

3 4

CCN Jan 2012 article : find_alt_origin_sym_mate
by Bryan Lepore 17 Feb '12

17 Feb '12

will the find_alt_origin_sym_mate article in CCN Jan. 2012 be described in any other publications? -Bryan

2 2

MRage
by Graeme Winter 17 Feb '12

17 Feb '12

Hi Folks, Following from the newsletter article I am trying to test out phaser.MRage. I have a well solved case which I know is in the pdb and a sequence file: phaser.MRage hklin=data.mtz composition.component.sequence=sequence.pir symmetry_exploration=pointgroup mode=full search.services=ncbi composition.count=3 however it generates a second component then complains it knows nothing about this: Sorry: Molecular weight needed for component 2. Any suggestions how to get around this? Many thanks, Graeme [gw56@ws050 rage]$ phaser.MRage hklin=data.mtz composition.component.sequence=sequence.pir symmetry_exploration=pointgroup mode=full search.services=ncbi composition.count=3 ############################################################################### # MRage version 0.1.0 # ############################################################################### All configuration options: -------------------------- hklin = data.mtz labin = None resolution_cutoff = 2.5 crystal_symmetry { unit_cell = None space_group = None } mode = quick *full symmetry_exploration = *pointgroup enantiomorph dataset output { root = "mrage" max_solutions_to_write = 0 gui_base_dir = None save_aniso_data = False job_title = None } composition { count = 3 component { sequence = sequence.pir stoichiometry = 1 ensemble { coordinates { pdb = None identity = None } } model_collection { trim = False coordinates { pdb = None identity = None } } template { pdb = None alignment = None } homology { file_name = None max_hits = 3 } search { services = local ncbi max_hits = 3 } } component { sequence = None stoichiometry = 1 ensemble { coordinates { pdb = None identity = None } } model_collection { trim = False coordinates { pdb = None identity = None } } template { pdb = None alignment = None } homology { file_name = None max_hits = 3 } search { services = local *ncbi max_hits = 3 } } } queue { technology = lsf threading sge *multiprocessing pbs cpus = 1 submission_command = None qslot_cpus = 1 } packing_pool = 20 rotation_peaks_cutoff = 0.75 post_refinement_cutoff = 0.75 final_selection_cutoff = 0.75 b_factor_refinement = True significant_peak_threshold = 7.0 sculptor_protocols = 11 10 13 12 1 3 2 5 4 7 6 9 8 *all minimal template_equivalence = False exclude_pdb_ids = None simple_run { enable = False component_sequence = None template_model = None blast_services = local ncbi max_blast_hits = 3 } Sorry: Molecular weight needed for component 2.

2 2