- phenixbb - phenix-online.org

problem with refinement in AutoBuild
by Michael Thompson 15 Aug '13

15 Aug '13

Hello Everyone, I am having a problem running autobuild, which appears to be happening during the refinement step. I start autobuild using : phenix.autobuild data=data.mtz model=model.mtz seq_file=sequence.dat rebuild_from_fragments=True i_ran_seed=124889 find_ncs=False refine_with_ncs=False nproc=4 The runs are failing during the first cycle with the following error message: AutoBuild_build_cycle AutoBuild Run 2 Thu Aug 15 13:37:34 2013 Build cycle 1 of 16 method:rebuild Initial rebuild cycle: just combining previous models or copying over initial PDB file Using NO_LIG_NO_WAT_CloMP_AutoBuild_input_ed.pdb for rebuild, copied to unrefined.pdb_1 Refining model: unrefined.pdb_1 Refinement strategy: ['individual_sites', 'individual_adp', 'occupancies'] No refinement carried out (may need free R flags) Error message from phenix.refine: Refinement failed... Sorry: ******************************************************************************** Failed to carry out AutoBuild_build_cycle: Initial phenix refinement of AutoBuild_run_2_/TEMP0/NO_LIG_NO_WAT_CloMP_AutoBuild_input_ed.pdb with aniso_refinement_data_PHX.mtz failed...quitting... /home/miket/projects/CloMP/Oxidized/APSoct2012/CloMP-TMAOcryo_B2_3/MRwithChainsAB/RefineCChalf/AutobuildPrep/AutoBuild_run_2_/LAST.LOG may have additional information... ******************************************************************************** I also checked the LAST.LOG file that is referenced in the error message. That appears to be a phenix.refine logfile, however it just cuts off after the "refinement start" header, with no information about the cause of the failure (no error message). My input MTZ file does indeed contain R-free flags. I have used this MTZ for refinement with phenix.refine (outside of autobuild) with no problems. Furthermore, earlier in the autobuild logfile, I see: Copying refinement_PHX.mtz to exptl_fobs_phases_freeR_flags.mtz Columns used: FP=FP SIGFP=SIGFP FreeR_flag=FreeRflag Which seems to me like it is recognizing my R-free flags column in the input MTZ. Does anyone have any idea what might be going wrong here?? Thanks, Mike -- Michael C. Thompson Graduate Student Biochemistry & Molecular Biology Division Department of Chemistry & Biochemistry University of California, Los Angeles miket(a)chem.ucla.edu

2 1

Phenix.refine suddenly stopped to refine B-factors
by krimmer＠staff.uni-marburg.de 12 Aug '13

12 Aug '13

Dear all, Phenix.refine suddenly stopped to refine B-factors, even though according to the settings it should refine them (checked box “individual B-factors”). This problem only occurs for one job. The weird thing is that this was always working and I did not do anything different, at least as far as I know. Also, this just happened to a colleague, totally independent from me. We never had this problem before. Here is a part of my log file: |-----overall-----|---macromolecule----|------solvent-------| stage b_max b_min b_ave b_max b_min b_ave b_max b_min b_ave 0 : 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 1_bss: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 1_ohs: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 1_xyz: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 1_adp: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 1_occ: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 2_bss: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 2_ohs: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 2_xyz: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 2_adp: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 2_occ: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 3_bss: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 3_ohs: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 3_xyz: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 3_adp: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 3_occ: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 4_bss: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 4_ohs: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 4_xyz: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 4_adp: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 4_occ: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 5_bss: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 5_ohs: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 5_xyz: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 5_adp: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 5_occ: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 5_bss: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 5_ohs: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 The values are not changing at all. The structure is almost done, but it still should change a bit during the refinement cycles. I noticed the problem because I added another water molecule to the structure, and the B-factor stayed at 30.00 which resulted in a negative density blob after refinement. Thank you for your help! Stefan

3 3

How to correct translational psuedosymmetry
by Appu kumar 12 Aug '13

12 Aug '13

Dear all, I have processed a protein-DNA complex data with imosflm. Phenix.xtriage suggest that the is a off origin peak of hight 41.702 hight at fractional coordiantes of 0.500 0.021 0.500 with distance of 40.019 from origin. I am new to such kind of problem, so would you please suggest me how to solve the psuedosymmetry problem encountered in data process. The present unit cells are 55.41 109.72 66.61 90 98.74 90 in space group P1211 with p_value is 3.116e-04. Your help would be much appreciated Thank you Appu

2 1

is phenix trying to polymerize b-octyl glucoside?
by Edward A. Berry 09 Aug '13

09 Aug '13

I'm refining a model with some molecules of the detergent b-octyl glucoside (BOG). I got BOG.cif from a recent ccp4 library, so not sure if this is a suggestion for phenix or ccp4, but: If two BOG molecules are consecutive in the pdb file, phenix stops in setup because it can't find atom HO1 to delete. Even if I had H's in my model there would be no HO1 because O1 is in glycosidic linkage with the octyl chain. I'm guessing phenix somehow realizes there is a sugar and when two are consecutive it autogenerates a link to make a glycosidic bond between them - not appropriate in the case of BOG. easy work around is to intersperse with non-sugar molecules, but if there is a better way, I'll take it. ed RuntimeError: delete_atom_in_place: unknown atom_id: HO1 apply_mod failure: pdbres="BOG R 18 " segid="R " comp id: BOG mod id: DEL-O1 Version: dev Release tag: 1420 Platform: intel-linux-2.6-x86_64 linux

2 2

Ensemble refinement?
by Joseph Noel 08 Aug '13

08 Aug '13

Is anyone having issues executing ensemble refinement in Phenix? I've tried several times now to run things and it always crashes without returning any maps and just returning a set of large PDB files such as: "TEAS-W273E-FPP-A_ensemble_0.6.pdb, TEAS-W273E-FPP-A_ensemble_0.8.pdb, TEAS-W273E-FPP-A_ensemble_0.9.pdb and TEAS-W273E-FPP-A_ensemble_0.1.pdb". The log file always ends with: " Time building chain proxies: 3.01, per 1000 atoms: 0.35 Number of scatterers: 8638 At special positions: 0 Unit cell: (125.75, 125.75, 120.07, 90, 90, 90) Space group: P 41 21 2 (No. 92) Number of sites at special positions: 0 Number of scattering types: 7 Type Number sf(0) S 19 16.00 P 2 15.00 Mg2+ 3 10.00 O 835 8.00 N 717 7.00 C 2765 6.00 H 4297 1.00 sf(0) = scattering factor at diffraction angle 0. AttributeError" I only entering a single chain with solvent and ions and hydrogens. All of the other parameters from the GUI are kept the same. Any suggestions or comments would be much appreciated. Joe ______________________________________________________________________________________ Joseph P. Noel, Ph.D. Arthur and Julie Woodrow Chair Investigator, Howard Hughes Medical Institute Professor, The Jack H. Skirball Center for Chemical Biology and Proteomics The Salk Institute for Biological Studies 10010 North Torrey Pines Road La Jolla, CA 92037 USA Phone: (858) 453-4100 extension 1442 Cell: (858) 349-4700 Fax: (858) 597-0855 E-mail: noel(a)salk.edu Publications & Citations: http://scholar.google.com/citations?user=xiL1lscAAAAJ Homepage Salk: http://www.salk.edu/faculty/noel.html Homepage HHMI: http://hhmi.org/research/investigators/noel.html ______________________________________________________________________________________

1 0

Restraints file for a modified peptide
by Michael Feldkamp 05 Aug '13

05 Aug '13

I have a peptide that has a fluorescent tag attached to its N-terminus along with a unnatural amino acid incorporated into the sequence. I know how to generate .cif files for the individual ligands using elbow, but I'm at a loss as to get Phenix to recognize these modifications to the peptide so that it refines the structure as one molecule with covalent linkages between the fluorescent tag, peptide and modified amino acid. Any help would be much appreciated. Thanks, Michael

3 2

MR Rosetta Error
by Masaki UNNO 05 Aug '13

05 Aug '13

Dear all I am trying to determine a protein structure using Rosetta MR. It is the first time to use Rosetta (Phenix). When I performed Rosetta MR via Phenix, I received the Error message as follows: What is wrong? What should I do to solve the problem? I use Ubuntu 12.10 as a platform. Starting mr_rosetta Date: Mon Aug 5 11:06:04 2013 Directory: /home/pad1/Work/20130222_PF/PAD1034/PHENIX/REF01/MR_ROSETTA_1 mr_rosetta { input_files { seq_file = "" hhr_files = None alignment_files = None model_info_file = "" data = "/home/pad1/Work/20130222_PF/PAD1034/process/PAD1034_NW12A_P61_scale2.mtz" data_labels = "F_PAD1034,SIGF_PAD1034" free_r_data = "" free_r_labels = "FreeR_flag" labin = None search_models = "/home/pad1/Work/20130222_PF/PAD1034/MR/PAD4_Arg_cut.pdb" copies_in_search_models = None mr_rosetta_solutions = "" ids_to_load = None map_coeffs = "" labin_map_coeffs = None map = "" refinement_params = "" display_solutions = False fragment_files = None fragment_files_chain_list = None fragment_files_9_mer_by_chain = None fragment_files_3_mer_by_chain = None use_dummy_fragment_files = False sort_fragment_files = True } output_files { log = "mr_rosetta.log" params_out = "mr_rosetta_params.eff" } directories { temp_dir = "/home/pad1/Work/20130222_PF/PAD1034/PHENIX/REF01/MR_ROSETTA_1/WORK_1" workdir = "/home/pad1/Work/20130222_PF/PAD1034/PHENIX/REF01/MR_ROSETTA_1" output_dir = "/home/pad1/Work/20130222_PF/PAD1034/PHENIX/REF01" gui_output_dir = "/home/pad1/Work/20130222_PF/PAD1034/PHENIX/REF01" top_output_dir = "/home/pad1/Work/20130222_PF/PAD1034/PHENIX/REF01/MR_ROSETTA_1" rosetta_path = "/home/unno19-frt/programs/rosetta-3.5" rosetta_binary_dir = "rosetta_source/bin" rosetta_binary_name = "mr_protocols.default" rosetta_script_dir = "rosetta_source/src/apps/public/electron_density" rosetta_pilot_script_dir = "rosetta_source/src/apps/pilot/frank/" rosetta_database_dir = "rosetta_database" } read_hhpred { number_of_models = 1 number_of_models_to_skip = 0 copies_to_extract = None only_extract_proper_symmetry = False } place_model { run_place_model = True prerefine { run_prerefine = False number_of_prerefine_models = 1000 number_of_models_in_ensemble = 1 } model_already_placed = False model_already_aligned = False number_of_output_models = 5 align_with_sculptor = True identity = None identity_for_scoring_only = 25 use_all_plausible_sg = True overlap_allowed = 10 selection_criteria_rot_value = 75 fast_search_mode = True search_down_percent = 25 mr_resolution = 4.2 refine_after_mr = True denmod_after_refine = True ps_in_rebuild = False find_ncs_after_mr = True min_length_ncs = 10 fixed_model = "" fixed_model_identity = None fixed_ensembles { fixed_ensembleID_list = None fixed_euler_list = 0 0 0 fixed_frac_list = 0 0 0 fixed_frac_list_is_fractional = True } copies_of_search_model_to_place = None } rescore_mr { run_rescore_mr = True nstruct = 5 relax = False include_unrelaxed_in_scoring = False align = True edit_model = False stage_to_rescore = "mr_solution" } rosetta_rebuild { run_rosetta_rebuild = True stage_to_rebuild = "rescored_mr_solution" max_solutions_to_rebuild = 5 min_solutions_to_rebuild = 1 llg_percent_of_max_to_keep = 50 rosetta_models = 100 chunk_size = 1 edit_model = True superpose_model = False } rosetta_rescore { run_rosetta_rescore = True percentage_to_rescore = 20 min_solutions_to_rescore = 2 } similarity { run_similarity = False required_cc = 0.2 number_of_required_cc = 5 } refine_top_models { run_refine_top_models = True stage_to_refine = None sort_score_type = None percent_to_refine = 20 denmod_after_refine = True } average_density_top_models { run_average_density_top_models = True percent_to_average = 100 } relax_top_models { run_relax_top_models = True stage_to_relax = None number_to_relax = 2 nstruct = 5 } autobuild_top_models { run_autobuild_top_models = True number_to_autobuild = 2 quick = False phase_and_build = False macro_cycles = None morph = False edit_model = True use_map_coeffs = True } setup_repeat_mr_rosetta { run_setup_repeat_mr_rosetta = True repeats = 1 template_repeats = 0 morph_repeats = 0 number_to_repeat = 1 acceptable_r = 0.25 minimum_delta_r = None } repeat_mr_rosetta { run_repeat_mr_rosetta = True copies_in_new_search_group = 1 update_map_coeffs_with_autobuild = True } rosetta_modeling { map_resolution = 3 map_grid_spacing = 1.5 map_weight = 1 map_window = 5 include_solvation_energy = True weights_file = "" } crystal_info { resolution = 0 space_group = "P 61" chain_type = *PROTEIN DNA RNA ncs_copies = Auto } control { verbose = False debug = False raise_sorry = False dry_run = False nproc = 1 group_run_command = "sh " queue_commands = None condor_universe = "vanilla" add_double_quotes_in_condor = True condor = None one_subprocess_level = None single_run_command = "sh " last_process_is_local = True background = None ignore_errors_in_subprocess = True check_run_command = False max_wait_time = 100 wait_between_submit_time = 1 wizard_directory_number = 1 n_dir_max = 100000 number_to_print = 5 write_run_directory_to_file = None rosetta_command = None resolve_command_list = None start_point = *place_model rescore_mr rosetta_rebuild rosetta_rescore \ similarity refine_top_models average_density_top_models \ relax_top_models autobuild_top_models \ setup_repeat_mr_rosetta repeat_mr_rosetta stop_point = place_model rescore_mr rosetta_rebuild rosetta_rescore \ similarity refine_top_models average_density_top_models \ relax_top_models autobuild_top_models \ setup_repeat_mr_rosetta repeat_mr_rosetta clean_up = False add_id = True } non_user_params { file_base = None print_citations = True highest_id = 0 is_sub_process = False dummy_autobuild = False dummy_refinement = False dummy_rosetta = False prerefine_only = False skip_clash_guard = True correct_special_position_tolerance = None ncs_in_refinement = *torsion cartesian None comparison_mtz = "" labin_comparison_mtz = None write_local_files = False rosetta_fixed_seed = None } } Best regards ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Masaki UNNO, Ph.D. Frotier Research Center for Applied Atomic Sciences, Ibaraki University Ibaraki Quantum Beam Research Center 162-1 Shirakata, Tokai, Naka, Ibaraki 319-1106, Japan Tel: 029-352-3239, Fax: 029-287-7872 E-mail: unno19(a)mx.ibaraki.ac.jp HP: http://www.fas.ibaraki.ac.jp/?page_id=961 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

2 1

Molecular replacement pipeline MRage
by Morten Grøftehauge 05 Aug '13

05 Aug '13

Hi guys, So I've been trying out MRage. I used the GUI since I didn't care too much and just wanted it started. I noticed two things. It doesn't show you the best results so far. It doesn't even save the best results so far. Everything is hidden in a giant log file in the .phenix/project_data folder (.phenix is hidden so I'm guessing that means you don't want people snooping around in there?). I can grep the TFZ= lines with grep 'TFZ=' mr_pipeline.log > grep-TFZ.txtand sort them with sort -r -k 13 grep-TFZ.txt but I think this is a bit excessive for the average user. ... Oh, and now I tried to re-run the highest scoring model in Phaser and it crashed and took MRage with it. Are you sure MRage is out of beta? Trying to be facetious, not mean. Cheers, Morten -- Morten K Grøftehauge, PhD Pohl Group Durham University

3 2

How does phenix handle deuterium atoms in neutron diffraction data?
by Javier Gonzalez 02 Aug '13

02 Aug '13

Hello, I'm performing a joint neutron-xray refinement with phenix, for a dataset collected from a crystal exchanged with deuterated water, ie, all exchangeable H atoms are replaced by D. I was wondering whether phenix has libraries of .cif dictionaries (or equivalent) in order to handle the name "D" of some "H" atoms in the .pdb files? I would like to use those dictionaries to input into Coot, so that it can handle deuterated molecules during geometry regularization, real space refinement, etc. Any help will be greatly appreciated. Thanks, Javier -- Javier M. Gonzalez Protein Crystallography Station Bioscience Division Los Alamos National Laboratory TA-3, Building 4200, Room 202 Mailstop T001 work: (505) 667-9376 (temporary)

2 1

feature enhanced maps
by Osman Asghar Mirza 02 Aug '13

02 Aug '13

Dear all i just installed phenix-1.8.2-1309 (OSX) and was looking for the "phenix.development.fem" but haven't been able to find it. any ideas? best osman Osman Mirza, Associate Professor PhD Department of Drug Design and Pharmacology Faculty of Health and Medical Sciences University of Copenhagen Universitetsparken 2, DK-2100 Denmark Phone: +45 3533 6175 Fax: +45 3533 6040 E-mail: om(a)farma.ku.dk<mailto:[email protected]>

2 1