- phenixbb - phenix-online.org

coot from phenix
by Andreas Förster 24 Oct '11

24 Oct '11

Dear all, I'm running phenix 1.7.2-869 on RHEL 6.1 64-bit. When I click on "Open in Coot" after refinement, coot doesn't open and I get the following error in the shell from which I started phenix: /csb/soft/Linux64/share/phenix-1.7.2-869/phenix-1.7.2-869/phenix/phenix/command_line/start_coot.sh: line 33: /csb/soft/Linux/src/coot-0.6/bin/coot: No such file or directory In the script that launches coot, $PHENIX_COOT is defined as /csb/soft/Linux/src/coot-0.6/bin/coot, which doesn't exist (and is a path that we used in the old days of 32-bit). Typing phenix.find_coot_command in the shell from which phenix was started returns /csb/soft/Linux64/src/coot-0.7-rev3689/bin/coot, which is the right coot. To me it seems the phenix script to launch coot is falling back on an internally defined $PHENIX_COOT variable (that's not defined from the shell). That's shouldn't be, should it? How do I fix this short of hardwiring $PHENIX_COOT as /csb/soft/Linux64/src/coot-0.7-rev3689/bin/coot? Thanks. Andreas -- Andreas Förster, Research Associate Paul Freemont & Xiaodong Zhang Labs Department of Biochemistry, Imperial College London http://www.msf.bio.ic.ac.uk

2 2

Unusual CCP4 to Phenix data conversion result
by R.M. Garavito 24 Oct '11

24 Oct '11

In upgrading to CCP4 6.2.0 and Phenix-1.7.2-869 (from CCP4 6.1.13 and Phenix-1.7.1-743, respectively) a bug has crept into reading of the MTZ files. We tend to scale and merge XDS processed data with Scala in CCP4, then use the resulting data in CCP4 and Phenix (with great success I might add). However, I just noticed that in testing the upgraded CCP4 and Phenix on OSX 7.2 (Lion) that the MTZ file is no longer read correctly. An example is below where after the output from scaling and merging of nonanomalous data from the CCP4i GUI gives, at the end of FREERFLAG, some 47,000 reflections. Putting this into xtriage or the reflection file editor in the Phenix GUI, one finds that the data for the F and the DANO columns are used; just selecting Fs alone is impossible. The end result is a doubling of the reflection number and the switching of the anomalous flag on. Processing with earlier versions did not do this except when actually using anomalous data. Why this has happened is not clear to me, but the column order in the MTZ file changed slightly between CCP4 6.1.13 and CCP4 6.2.0. Any ideas, am I missing something, or is this a push to only refine with intensities? Regards, Michael ***************** From Scala (3.3.20) Total number unique : 46813 From FREERFLAG Number of Reflections = 47044 (including non-observed NaN reflections) MTZ column labels H K L FreeR_flag F_S13E_X2 SIGF_S13E_X2 DANO_S13E_X2 SIGDANO_S13E_X2 F_S13E_X2(+) SIGF_S13E_X2(+) F_S13E_X2(-) SIGF_S13E_X2(-) ISYM_S13E_X2 IMEAN_S13E_X2 SIGIMEAN_S13E_X2 I_S13E_X2(+) SIGI_S13E_X2(+) I_S13E_X2(-) SIGI_S13E_X2(-) ****************** In Phenix: Xtriage output choosing "Imean" Miller array info: /Users/garavito/ccp4_projects/sus1_S13/reflections_tmp.mtz:IMEAN_S13E_X2,SIGIMEAN_S13E_X2 Observation type: xray.amplitude Type of data: double, size=46772 Type of sigmas: double, size=46772 Number of Miller indices: 46772 Anomalous flag: False Xtriage output choosing "Fs" Miller array info: /Users/garavito/ccp4_projects/sus1_S13/reflections_tmp.mtz:F_S13E_X2,SIGF_S13E_X2,DANO_S13E_X2,SIGDANO_S13E_X2 Observation type: xray.amplitude Type of data: double, size=88486 Type of sigmas: double, size=88486 Number of Miller indices: 88486 Anomalous flag: True **************************************************************** R. Michael Garavito, Ph.D. Professor of Biochemistry & Molecular Biology 513 Biochemistry Bldg. Michigan State University East Lansing, MI 48824-1319 Office: (517) 355-9724 Lab: (517) 353-9125 FAX: (517) 353-9334 Email: garavito(a)msu.edu ****************************************************************

3 6

phenix.pdbtools
by Abhinav Kumar 21 Oct '11

21 Oct '11

Hi, I had been using the pdbtools to obtain the rmsd on various geometry parameters using the command: phenix.pdbtools my.pdb --show-geometry-statistics Now pdbtools does not recognize the option "--show-geometry-statistics" What is the option now? I am using phenix-1.7.2-869. -- Thanks, Abhinav JCSG@SSRL, SLAC Phone: (650) 926-2992 Fax: (650) 926-3292

3 3

ligands at a special position
by Subhani Bandara 19 Oct '11

19 Oct '11

Hi, I am trying to solve a structure of a protein that is cocrystallized with a chelator complex using phenix (GUI). The complex is Tb(DPA)3 and DPA stands for dipicolinic acid. In the structure (C222) there are few Tb complexes on the two fold axis with symmetry coordinates 1/4X, 1/4Y, Z. In those, Tb and the center C (C4) and N (N1) of one DPA are sitting directly on the 2 fold axis. So I left only half of that DPA ligand structure on PDB file and changed the occupancy of N1 and C4 to 1/2 while keeping the occupancy of all other atoms of the half ligand structure at 1.00. Then I refined it while keeping above occupancies and xyz coordinated of the ligand fixed. After doing this the other half of the ligand could be built by symmetry, but it left a space between two halves(so it was not as one ligand). Can you please let me know what I can do to fix this structure. Also is there a way to keep only x and y coordinates fixed while varying z during refinement. Thanks Subhani

3 4

Re: [phenixbb] Phenix building residues not found in input sequence file
by Thomas C. Terwilliger 18 Oct '11

18 Oct '11

Hi Allan, I assume that there is density showing protein in the region where autobuild is building these extra residues...is that right? If so...then is the remainder of the model matched very well to sequence (i.e., is more or less what you expected to see found, plus there is extra density?) It could be that there is really something in the crystal that you didn't expect. It could also be that part of the structure is not (or not correctly) assigned to sequence, so that this part of the structure really does go with part of your sequence. Does that help? It would be interesting to hear more details of this. All the best, Tom T >> Hi everyone, >> >> Apparently my Phenix (1.6.3-473) is building in extra residues in the >> C-terminal region, which is not found in my input sequence file. >> >> It's trying to build a model which is longer than I expected. Phenix >> building 87 extra amino acids. Any clue what's wrong? >> >> I made sure that these amino acids are not residues found in the >> N-terminal region or anywhere else in the protein. >> >> Allan >> >> -- >> Allan Pang >> >> PhD Student >> >> G35 Joseph Priestley Building >> Queen Mary University of London >> London >> E1 4NS >> >> Phone number: 02078828480 >>

1 0

post-doc position at SSRL
by van den Bedem, Henry 17 Oct '11

17 Oct '11

Modeling Conformational Ensembles of Proteins and Complexes A post-doctoral position with a focus on modeling conformational ensembles of proteins and complexes is available immediately in the Structural Genomics Division of Stanford Synchrotron Radiation Lightsource at the SLAC National Accelerator Laboratory (http://ssrl.slac.stanford.edu). Macromolecules and assemblies are studied by a variety of experimental techniques, including protein crystallography and Small Angle X-ray Scattering. Nano-crystallography and single particle imaging data from the Linac Coherent Light Source (LCLS), a first-of-its-kind free electron laser, are also being targeted at these systems. Accurately modeling the conformational variability of macromolecules is essential to understanding their cellular function. This position at the forefront of computational structural biology aims to develop advanced algorithms to fit domains or proteins to experimental data from multiple sources while allowing for significant conformational changes. The goal is to automatically compute 3-D models, and to provide structural insights into functionally relevant dynamics, especially for time-resolved experiments or a series of static structures. The position offers the opportunity to closely collaborate with the Stanford Computer Science Department and work with X-ray, SAXS and LCLS nanocrystallography data. SLAC houses ample computational resources. Funding is available for up to 3 years. The starting date is as soon as possible. Candidates should have a recent PhD in computer science, applied mathematics, computational biology/crystallography or a related field with strong programming skills. Please send a CV, a short summary of your research experience and two letters of recommendation to Henry van den Bedem (vdbedem(a)slac.stanford.edu). Further details are available upon request. Relevant publications: Yao et al. Sampling-based exploration of folded state of a protein under kinematic and geometric constraints. Proteins: Structure, Function, and Bioinformatics (in press). Fraser et al. (2011) Accessing protein conformational ensembles using room-temperature X-ray crystallography. PNAS 108:16247-16252 van den Bedem et al. (2009) Modeling discrete heterogeneity in X-ray diffraction data by fitting multi-conformers. Acta Cryst. D65:1107-1117

1 0

Re: [phenixbb] STOPWIZARD
by Thomas C. Terwilliger 17 Oct '11

17 Oct '11

Hi Morten, Yes, you can use STOPWIZARD to stop a phenix wizard, including mr_rosetta. However the "restarting" part is no longer enabled for most of the wizards. For mr_rosetta, you can restart. Here is how: You can re-start mr_rosetta only at the beginning of major stages (like "place_model", "rosetta_rebuild" etc)...but not in between. Normally at the end of a major stage a .pkl file is written out with text like "type this to see all the results". You can almost always give your original command, the command "start_point=xxx" and "mr_rosetta_solutions=my_pickle_file.pkl" and it should then continue on from there. I will put this in the docs... All the best, Tom T >> Dear phenixbb, >> >> The Phenix manual mentions several times that making a file named >> STOPWIZARD >> in the wizard root directory will stop the wizard in an orderly fashion >> the >> next time it checks for STOPWIZARD. There is however no mention of how to >> properly restart the wizard after having stopped it... How? I seem to >> remember something about it being >>>phenix.appropriate_wizard params_used_in_wizard.eff >> but I'd like be sure that I don't simply start over by doing that. >> >> Also can I then change the parameters before restarting the wizard? Say, >> hypothetically of course, that I've started mr.rosetta and asked for 1000 >> prerefined models and realised that it will take ridiculously long and be >> rather redundant since I specified both the .hhr file from global and >> local >> alignment and the top hit has two chains (resulting in mr.rosetta wanting >> to >> build 4000 prerefined models) but I don't want to throw away the 700 >> models >> already built. Could I do something about that? >> >> Sincerely, >> Morten Grøftehauge, PhD >> Durham University >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://phenix-online.org/mailman/listinfo/phenixbb >>

1 0

STOPWIZARD
by Morten Grøftehauge 17 Oct '11

17 Oct '11

Dear phenixbb, The Phenix manual mentions several times that making a file named STOPWIZARD in the wizard root directory will stop the wizard in an orderly fashion the next time it checks for STOPWIZARD. There is however no mention of how to properly restart the wizard after having stopped it... How? I seem to remember something about it being >phenix.appropriate_wizard params_used_in_wizard.eff but I'd like be sure that I don't simply start over by doing that. Also can I then change the parameters before restarting the wizard? Say, hypothetically of course, that I've started mr.rosetta and asked for 1000 prerefined models and realised that it will take ridiculously long and be rather redundant since I specified both the .hhr file from global and local alignment and the top hit has two chains (resulting in mr.rosetta wanting to build 4000 prerefined models) but I don't want to throw away the 700 models already built. Could I do something about that? Sincerely, Morten Grøftehauge, PhD Durham University

1 0

mr_rosetta memory usage
by Tjaard Pijning 17 Oct '11

17 Oct '11

Hi, I am currently running mr_rosetta and found that after a while the computer (running OpenSuse/linux) slows down dramatically because its memory (8 GB) is running full. I’m not sure if the program still runs or not. Is there a way to limit the amount of memory that phenix/mr_rosetta uses ? Or should I e.g. lower the amount of models being used (currently 5) ? (Before starting mr_rosetta I ran the regression test which was fine). Kind regards, Tjaard Pijning University of Groningen This is the parameter file : mr_rosetta { input_files { seq_file = gEP_mature.txt" hhr_files = "051011_3011840.hhr" alignment_files = None model_info_file = "" data = "XDS_ASCII_C2221.mtz" labin = None search_models = None copies_in_search_models = None mr_rosetta_solutions = "" ids_to_load = None map_coeffs = "" labin_map_coeffs = None map = "" display_solutions = False fragment_files = "aat000_03_05.200_v1_3" fragment_files = "aat000_09_05.200_v1_3" use_dummy_fragment_files = False sort_fragment_files = True } output_files { log = "mr_rosetta.log" params_out = "mr_rosetta_params.eff" } directories { temp_dir = "mr_rosetta/MR_ROSETTA_7/WORK_1" workdir = "mr_rosetta/MR_ROSETTA_7" output_dir = "mr_rosetta" top_output_dir = "mr_rosetta/MR_ROSETTA_7" rosetta_path = "" 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 = 5 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 = 300 selection_criteria_rot_value = 75 fast_search_mode = True search_down_percent = 25 mr_resolution = 3 refine_after_mr = True denmod_after_refine = True find_ncs_after_mr = True min_length_ncs = 10 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 = None chain_type = *PROTEIN DNA RNA ncs_copies = 2 } control { verbose = False debug = False raise_sorry = False dry_run = False nproc = 5 group_run_command = "sh " condor = None single_run_command = "sh " background = None ignore_errors_in_subprocess = True check_run_command = False max_wait_time = 100 wait_between_submit_time = 1 n_dir_max = 100000 number_to_print = 5 write_run_directory_to_file = 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 } 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 skip_clash_guard = True correct_special_position_tolerance = None comparison_mtz = "" labin_comparison_mtz = None write_local_files = False rosetta_fixed_seed = None } }

3 2

re-building chains
by Leonid Sazanov 14 Oct '11

14 Oct '11

Hi, how exactly can one tell Autobuild to rebuild-in-place chains A and B from the input model, but build chain C de-novo? If I give it list of chains to re-build (rebuild_near_chain="A B" in GUI, dev-874 version), what happens is that either all chains are rebuild-in-place (if rebuild-in-place=True) or all chains are built de-novo (if rebuild-in-place=False). Thanks! -- Dr. Leonid A. Sazanov Research group leader Medical Research Council Mitochondrial Biology Unit Wellcome Trust / MRC Building Hills Road Cambridge CB2 0XY WEB: www.mrc-mbu.cam.ac.uk Tel: +44-1223-252910 Fax: +44-1223-252915

2 1

1.7.2 slow to open
by Jinsong Liu 12 Oct '11

12 Oct '11

Hi there, I wonder if there is anyone else experiencing a VERY VERY slow opening of the phenix GUI for the latest version of 1.7.2? We have RHEL 5.3 and CentOS 5.6, both machines show this. Most of the time, it will take minutes, sometime more than 10 minutes. There is no any warning or error message, after the long waiting the GUI will suddenly show up. In the older version, it usually takes seconds to load the GUI. Thanks. Jinsong -- Jinsong Liu, Ph.D. Principle Investigator & Head of Laboratory of Structural Biology Guangzhou Institute of Biomedicine and Health (GIBH) Chinese Academy of Sciences Guangzhou ( Canton ) 510630, China Phone: 86-20-32015317 (Office) E-mail: liu_jinsong(a)gibh.ac.cn

2 2

Re: [phenixbb] mr_rosetta
by Thomas C. Terwilliger 11 Oct '11

11 Oct '11

Hi Abhinav, This format is very unforgiving. It has to have exactly 5 lines: ## TARGET SP_3n91 # hhsearch scores_from_program: 0 1.00 0 --DDDEVTMNSK.... 0 GSDNEFPDFDYQ.... -- also note I changed SP 3n91 to SP_3n91. Your PDB file name must also start with SP_3n91 in this case. You need a pdb file name with at least 5 chars before the ".pdb" and these 5 characters have to match the 5 characters after TARGET on the first line. Note that the other format read by mr_rosetta is a lot easier: > title text for sequence of target to follow VDFNGYWKMLSNENFEEYLRALDVNVALRKIANLLKPDKEIVQDGDHMIIRTLSTFRNYIMDFQVGKEFEEDLTGIDD > title text for sequence of template (supplied PDB) to follow -AFSGTWQVYAQENYEEFLRAISLPEEVIKLAKDVKPVTEIQQNGSDFTITSKTPGKTVTNSFTIGKEAEIT--TMDG All the best, Tom T >> Somebody helped me out with the format. Thanks. >> My alignment file looks like: >> ## SP 3n91 >> >> >> 0 --DDDEVTMNSK.... >> 0 GSDNEFPDFDYQ.... >> >> >> >> But phenix says >> >> ******************* ERROR ENDING *************** >> >> Sorry, the alignment file ['xx.ali'] >> does not seem to have the correct format? >> It can either start with '##' or with '> title of target seq here' >> >> ******************* ERROR ENDING *************** >> >> I do have the file starting with '##'! >> >> Thanks, >> Abhinav >> >> JCSG@SSRL, SLAC >> Phone: (650) 926-2992 >> Fax: (650) 926-3292 >> >> >> On 10/11/2011 12:23 PM, Abhinav Kumar wrote: >>> Hi, >>> >>> I am looking for file format for alignment file (*.ali) used in >>> mr_rosetta. >>> The description on the Phenix web site is not clear. >>> Can someone please describe the file format with an example? >>> -- >>> >>> Thanks, >>> Abhinav >>> >>> JCSG@SSRL, SLAC >>> Phone: (650) 926-2992 >>> Fax: (650) 926-3292 >>> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://phenix-online.org/mailman/listinfo/phenixbb >>

2 2

Re: [phenixbb] mr_rosetta
by Thomas C. Terwilliger 11 Oct '11

11 Oct '11

Hi Abhinav, Here is an example. > cat WORK_1/alignment.ali ## TARGET coord # hhsearch scores_from_program: 0 1.00 0 AQLMDMRD 0 AQLMDMRD -- (the -- at the end is included in the file) You can also use one with just 4 lines (sculptor format): > title text for sequence of target to follow VDFNGYWKMLSNENFEEYLRALDVNVALRKIANLLKPDKEIVQDGDHMIIRTLSTFRNYIMDFQVGKEFEEDLTGIDD > title text for sequence of template (supplied PDB) to follow -AFSGTWQVYAQENYEEFLRAISLPEEVIKLAKDVKPVTEIQQNGSDFTITSKTPGKTVTNSFTIGKEAEIT--TMDG All the best, Tom T >> Hi, >> >> I am looking for file format for alignment file (*.ali) used in >> mr_rosetta. >> The description on the Phenix web site is not clear. >> Can someone please describe the file format with an example? >> -- >> >> Thanks, >> Abhinav >> >> JCSG@SSRL, SLAC >> Phone: (650) 926-2992 >> Fax: (650) 926-3292 >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://phenix-online.org/mailman/listinfo/phenixbb >>

1 0

mr_rosetta
by Abhinav Kumar 11 Oct '11

11 Oct '11

Hi, I am looking for file format for alignment file (*.ali) used in mr_rosetta. The description on the Phenix web site is not clear. Can someone please describe the file format with an example? -- Thanks, Abhinav JCSG@SSRL, SLAC Phone: (650) 926-2992 Fax: (650) 926-3292

1 1

Re: [phenixbb] phenix.mr_rosetta not running "prerefine"
by Thomas C. Terwilliger 11 Oct '11

11 Oct '11

Hi Laurie, The problem is that mr_rosetta needed a pdb file name with more than 4 characters before the ".pdb" for this shortcut of generating pre-refined models to work (the Rosetta routine requires a longer file name). So if you rename your PDB file this should run (it does on my machine.) I will fix mr_rosetta so that it can take your PDB file. All the best, Tom T >> On Tue, Oct 11, 2011 at 12:33 PM, Thomas C. Terwilliger < >> terwilliger(a)lanl.gov> wrote: >> >>> Hi Laurie, >>> >>> Yes, it is not doing anything. There is something that didn't work in >>> the >>> inputs. Can you send me your fasta and PDB files and I'll just run it >>> here >>> to see what is wrong? >>> >>> Thanks! >>> -Tom T >>> >>> >> OK - I am also sending: >>> >> >>> >> >>> /netscr/bettsl/rosetta_test/MR_ROSETTA_2/GROUP_OF_PLACE_MODEL_1/RUN_1/REBUILD_IN_SETS_1/RUN_1/WORK_1/rebuild.log >>> >> >>> >> On Tue, Oct 11, 2011 at 12:01 PM, Thomas C. Terwilliger < >>> >> terwilliger(a)lanl.gov> wrote: >>> >> >>> >>> Hi Laurie, >>> >>> >>> >>> Yes, something didn't work it all right. Can you send me, off-list, >>> this >>> >>> log file: >>> >>> >>> >>> >>> >>> >>> /netscr/bettsl/rosetta_test/MR_ROSETTA_1/GROUP_OF_PLACE_MODEL_1/RUN_1/REBUILD_IN_SETS_1/RUN_FILE_1.log >>> >>> >>> >>> which I am guessing has some kind of error message in it. >>> >>> >>> >>> Thanks! >>> >>> -Tom T >>> >>> >>> >>> >> Hi - I tried the little prerefine without any data script in your >>> >>> >> documentation - I don't think it's working right. I attach the >>> >>> script >>> >>> >> file >>> >>> >> that I made (prerefine.dat) and the RUN_FILE_1.log. >>> >>> >> >>> >>> >> I can send more info offline if need be. >>> >>> >> >>> >>> >> Thanks >>> >>> >> >>> >>> >> LAurie Betts >>> >>> >> _______________________________________________ >>> >>> >> phenixbb mailing list >>> >>> >> phenixbb(a)phenix-online.org >>> >>> >> http://phenix-online.org/mailman/listinfo/phenixbb >>> >>> >> >>> >>> >>> >>> _______________________________________________ >>> >>> phenixbb mailing list >>> >>> phenixbb(a)phenix-online.org >>> >>> http://phenix-online.org/mailman/listinfo/phenixbb >>> >>> >>> >> >>> >>> >>

1 0

refinements settings .eff
by Christian Roth 11 Oct '11

11 Oct '11

Hi Phenix gurus, I have done a refinement with an sodium ion. I had generated a phil file for the distance restraints to the surrounding water. Now i am not sure if it is really there. I removed it and removed the phil from the list in the gui. Now phenix complains about no Atom selected and stops. In the .eff file all restraints are incorporated and in I deleted therefore the geometric restraints manually from the . eff file. As a consequence an unused parameter warning appears about unused parameters or spelling errors in the input file. I guess the .eff file. My default parameter folder is empty. Does phenix load the .eff file from the previous run and would it help to remove this? I dont wanne start a new project just to be in the default configuration mode without additional information in the .eff file. I did this using the gui so I am not sure how to incorporate the --unused_ok command to which may help to convince phenix to process the job. I hope it is clear what I mean and thanks for any suggestion. Christian

2 1

phenix.mr_rosetta not running "prerefine"
by Laurie Betts 11 Oct '11

11 Oct '11

Hi - I tried the little prerefine without any data script in your documentation - I don't think it's working right. I attach the script file that I made (prerefine.dat) and the RUN_FILE_1.log. I can send more info offline if need be. Thanks LAurie Betts

2 1

Simulated annealing composite omit map conversion
by Yuri 11 Oct '11

11 Oct '11

The figures are good illustration of what seems to be the consensus: Finely sampled maps can represent the best mathematical approximation of the continuous electron density function. The weirdness I referred to probably has more to do with what Nathaniel alluded to previously. The fact that our eyes get so used to staring at coarse pointy density - caused by the default grid of building programs (in the case of "on the fly" map generation from coefficients). That of course can be changed if one changes the sampling rate and get their eyes used to the new look. that I suppose would become your new norm. -- Yuri Pompeu

1 0

MR_Rosetta and dimer of heterodimer
by Morten Grøftehauge 11 Oct '11

11 Oct '11

Hi everybody, I'm trying to phase a dimer of a heterodimer of a large, ~700aa, and small, ~150aa, component. The only good search model according to HHpred covers the first 500aa of the large component. So I am wondering whether I should add the 150aa to the end of the sequence file and add them to the fragment files in order for phenix.mr_rosetta to do density modification correctly? Or does this not matter at all? According to HHpred and others the last 200aa are beta sheets, probably a a beta sandwich or two, and the 150aa component consists of helical elements but that is as close as it really gets. Cheers, Morten

2 1

Re: [phenixbb] Simulated annealing composite omit map conversion
by Yuri 11 Oct '11

11 Oct '11

hi Nathaniel, Thats simple enough... Do I have any control over grid ( make a tighter looking mesh for publication) -- Yuri Pompeu

8 9

Phenix not reading my sca files
by Scott Classen 10 Oct '11

10 Oct '11

Why? Error is: Unknown file format for file /home/me/phenix/autosol/s4_peak.sca Unknown file format for file /home/me/phenix/autosol/s4_infl.sca Unknown file format for file /home/me/phenix/autosol/s4_rem.sca this also does not work: phenix.reflection_file_converter s4_peak.sca --mtz . error is: file_name: s4_peak.sca file_type: None No reflection data found in input file. Please advise. Thanks, Scott Classen

2 5

Phenix file organization
by R.M. Garavito 10 Oct '11

10 Oct '11

To the Phenix gurus, I am in the midst of reorganizing our little Mac computer cluster, and will have a mixed collection of machines with OS 10.6.8 or 10.7.1. While the reorganization is going fine for the most part, I am having a little trouble getting Phenix to work on one machine. The problem is integrating a machine (a MacPro) into the cluster with users having different UIDs from their NIS accounts on the rest of the cluster. Thus, the NIS UID leads to not only a different UID for the user, but a different account (and therefore a different "home" directory) name. Simply changing ownership of the user's directories and files globally in OS 10.6.8 leads to Phenix choking at some stage when used under the NIS UID, but without an informative error message. Things go on fine when used under the old local UID and user name. So my question is does Phenix store the absolute paths in its various project indices? Any suggestions would be helpful. Thanks, Michael **************************************************************** R. Michael Garavito, Ph.D. Professor of Biochemistry & Molecular Biology 513 Biochemistry Bldg. Michigan State University East Lansing, MI 48824-1319 Office: (517) 355-9724 Lab: (517) 353-9125 FAX: (517) 353-9334 Email: garavito(a)msu.edu ****************************************************************

2 1

Simulated annealing composite omit map conversion
by Yuri 10 Oct '11

10 Oct '11

"Fabricated" was probably a poorly chosen word. That implies made up data. I agree, I dont think cahnging the grid spacing would "introduce" unreal data to the maps... I just meant that it looks a little weird to have overly smooth mesh for a low res map... -- Yuri Pompeu

1 0

Computing Robs vs Rcalc for twinned data with xtriage
by Pietro Roversi 10 Oct '11

10 Oct '11

Dear all, I am trying to run phenix.xtriage to compute Robs vs Rcalc for twinned data (from A. A. Lebedev, A. A. Vagin and G. N. Murshudov Intensity statistics in twinned crystals with examples from the PDB, Acta Crystallographica, D62, 83-95, (2006) My input: phenix.xtriage xray_data.file_name=$1.mtz xray_data.obs_labels="IMEAN, SIGIMEAN" xray_data.calc_labels="FCalc, PHICalc" xray_data.high_resolution=3.1 > phenix.xtriage.$$.log The FCalc, PHICalc columns (from CCP4-sfall) are in my mtzfile: 1 ASC -27 27 0 100.00 0.7 10.4 52.72 2.60 H H 2 NONE -90 90 0 100.00 0.7 34.2 52.72 2.60 H K 3 NONE 0 15 0 100.00 5.7 5.7 52.72 2.60 H L 4 NONE 0.0 19.0 0 100.00 9.50 9.50 52.72 2.60 I FreeR_flag 5 NONE -244.3 13664.9 0 100.00 467.59 469.74 52.72 2.60 J IMEAN 6 NONE 1.8 990.9 0 100.00 43.05 43.05 52.72 2.60 Q SIGIMEAN 7 NONE 16.2 1157.9 0 100.00 173.14 173.14 52.72 2.60 F F 8 NONE 1.7 51.3 0 100.00 12.91 12.91 52.72 2.60 Q SIGF 9 NONE 0.7 7827.1 0 100.00 257.63 257.63 52.72 2.60 F FCalc 10 NONE 0.0 360.0 0 100.00 180.24 180.24 52.72 2.60 P PHICalc and yet I get the error message: creamint:~/Projects/FactorI/Tetartohedral/P1/DataFiles> ~/Scripts/phenix_xtriage.csh AUTOMATIC_DEFAULT_free.klh_sfall1 Running phenix.xtriage No matching array: scaling.input.xray_data.calc_labels=FCalc, PHICalc Possible choices: AUTOMATIC_DEFAULT_free.klh_sfall1.mtz:IMEAN,SIGIMEAN AUTOMATIC_DEFAULT_free.klh_sfall1.mtz:F,SIGF Please use scaling.input.xray_data.calc_labels to specify an unambiguous substring of the target label. Sorry: No matching array: scaling.input.xray_data.calc_labels=FCalc, PHICalc

2 1

Automated Mg & K ions picking
by Kristen Eisenberg 10 Oct '11

10 Oct '11

Hi Hasan, yes, we plan to extend water picking procedure so it can identify single-atom ions, but no specific tool is available to use right now. You need to inspect waters (for example, a water with "too small" refined B-factor is likely to be a metal), and density peaks (phenix.refine may not place a water into a density peak if it is too different from what is expected for water). Kristen Eisenberg Billige Flüge Marketing GmbH Emanuelstr. 3, 10317 Berlin Deutschland Telefon: +49 (33) 5310967 Email: utebachmeier at gmail.com Site: http://flug.airego.de - Billige Flüge vergleichen

1 0

Re: [phenixbb] PHENIX atom selection language
by Kay Diederichs 09 Oct '11

09 Oct '11

Am 09.10.2011 04:20, schrieb phenixbb-request(a)phenix-online.org: > Date: Thu, 6 Oct 2011 18:19:34 -0400 > From: Jianghai Zhu<jzhu(a)idi.harvard.edu> > To: PHENIX user mailing list<phenixbb(a)phenix-online.org> > Subject: Re: [phenixbb] PHENIX atom selection language > Message-ID:<E21B0319-CC1D-47A6-9F23-4493700B319C(a)idi.harvard.edu> > Content-Type: text/plain; charset=us-ascii > > Cool. Thanks, Pavel. It would be nice to have a list of all selection keywords and corresponding definitions somewhere. > > -- Jianghai Hi Jianghai, anybody can contribute to the CCP4 wiki, where a Phenix article already exists - check out <http://strucbio.biologie.uni-konstanz.de/ccp4wiki/index.php/Phenix#Atom_sel…> Unfortunately the wiki is hit by spammers, and we regularly delete whole batches of these pseudo-"users" to keep the userlist clean. The spammers use @gmail.com and @hotmail.com email addresses which we now do not accept any more as new users' email - such addresses are just too easy to obtain. So if you would like to contribute make sure a) to use a email adress of your institution (e.g. university) b) and to contribute to an existing article or start a new one since we won't delete "real" users HTH, Kay

1 0

Re: [phenixbb] Simulated annealing composite omit map conversion
by Yuri 09 Oct '11

09 Oct '11

I agree I like 0.25 (I believe thats what CNS spits out). Anything tighter looks fabricated -- Yuri Pompeu

1 0

Simulated annealing composite omit map conversion
by Yuri 09 Oct '11

09 Oct '11

Hello Everyone, I am calculating an SA. composite omit map to show active site density. It looks like the output format is mtz coefficients. How can I get it to ccp4 or .xplor format for PyMol? Thanks -- Yuri Pompeu

2 1

Educative Tool Idea
by Yuri 07 Oct '11

07 Oct '11

Hello Developers, This may be a silly idea after all, nonetheless: Would it be possible (worhtwhile) to have a tool that lets the user input a molecule of a certain type (probably rather small size, i.e, an amino acid or cofactor), and then input a space group. Phenix could then generate a relection file that would be displayed using either data2d or 3d viewer. Obviously there are many other factors to be considered but just an idea... -- Yuri Pompeu

6 10

Occupancy refinement issue
by AFL 07 Oct '11

07 Oct '11

Dear All, I would like to refine the occupancy of _all_ components together for the following system: - free LYS (chain A and resname LYS) - PLP-LYS complex (chain A and resname LLP) - PLP-GLU complex (chian D and resname PPE) The LYS and LLP are defined as alternative conformations. What I've tried so far is: - just phenix.refine without any additional constraints and the result is occupancy refinement of LYS vs. LLP and separately refined PPE. - constrained_group { selection = chain A and resseq 90 # position of LYS, LLP selection = chain C and resseq 1 # PPE } Result is occupancy refinement of LYS+LLP vs. PPE and total mess in the active site (=ligands moved to new positions far away from where they should be). - constrained_group { selection = chain A and resname LYS # ALYS, BLLP doesn't work selection = chain A and resname LLP selection = chain C and resname PPE } Results in occupancy refinement of LYS vs. PPE and the LLP occupancy is set to 0.00. So, the big question is: - is what I'm trying to do doable at all? - what would be the right way to do it? I will appreciate all and any suggestions on how to proceed. Regards, Andrzej Lyskowski -- alnyds(a)gmail.com e-mail home

2 1

PHENIX atom selection language
by Jianghai Zhu 06 Oct '11

06 Oct '11

Hi, Where can I find a complete description of PHENIX atom selection language? The examples in the PHENIX.REFINE manual are not enough. Thanks. -- Jianghai

2 6

Fwd: Really good stats... Unreal!??
by Yuri 05 Oct '11

05 Oct '11

Hi everyone, So I am done refining a structure with around 760 residues in the model to a resolution of 2.29A. Rwork 0.18 Rfree 0.24 RMSD angle 1.1 Bond 0.008, clashscore 28. As a I last round I ran optmize_weights and got some interesting numbers. Rwork 0.17 Rfree 0.23 RMSD Angle 0.68! Bond 0.003! clashscore 10.2 I dont know how real these are. That would make this a new pdb record for tight geometry!!! Any ideas? Best, -- Yuri Pompeu -- Yuri Pompeu

2 1

local refinement of target weights?
by Yarrow Madrona 04 Oct '11

04 Oct '11

Hello, I am wondering if there is a way to restrain geometry target weights for a ligand seperate from the rest of the molecule. I would like to impose higher restraints for the ligand than that of the protein in the same refinement run. Thanks for your help. -Yarrow -- Yarrow Madrona Graduate Student Molecular Biology and Biochemistry Dept. University of California, Irvine Natural Sciences I, Rm 2403 Irvine, CA 92697

2 2

combinding two CIF's and creating links in elbow
by Yarrow Madrona 04 Oct '11

04 Oct '11

Hello, 1. I have a question on creating links using elbow. I have a CIF of NO and I want to define a chemical link between NO and Fe of a HEME group. I need to input the geometry as NO is acting as a metal ligand. I really don't know where to begin. Any help would be greatly appreciated. 2. Is there a quick way of combining CIF files into one large file without a lot of copy and pasting? Thank you. -- Yarrow Madrona Graduate Student Molecular Biology and Biochemistry Dept. University of California, Irvine Natural Sciences I, Rm 2403 Irvine, CA 92697

2 1

Phenix version 1.7.2 released
by Paul Adams 30 Sep '11

30 Sep '11

The Phenix developers are pleased to announce that version 1.7.2 of Phenix is now available. Binary installers for Linux, and Mac OSX platforms are available at the download site: http://phenix-online.org/download/ Some of the new features in this version are: Graphical interface: ==================== - phenix.data_viewer added under Reflection Tools - New GUIs: HySS, phenix.cut_out_density, phenix.reciprocal_space_arrays General: ======== - New commands: phenix.cut_out_density, phenix.data_viewer, phenix.morph_model - Default map output format changed to CCP4 - Native support for CIF-format reflection files (e.g. from PDB) in many commands and GUI Refinement: =========== - New Xray/restraints weight optimization. Keywords optimize_wxc and optimize_wxu replaced with optimize_xyz_weight and optimize_adp_weight, respectively. - Grid searches parallelized: use 'nproc' parameter to set number of CPUs to use - this affects bulk solvent mask optimization, and XYZ and ADP weight optimization (in addition to TLS identification) - can reduce runtime by 75% or more depending on settings (see manual for details) - Modified non-bonded parameters resulting in improved geometry, especially at low resolution - Rigid-body refinement mode defaults to grouping by chain - Printout of validation statistics during coordinate refinement - Bug fix in phase-combined map calculation results in significantly improved maps - Bug fix in X-ray/ADP restraints weight calculation results in smaller B-factors deviations for bonded atoms - Experimental feature (still in testing): - Torsion-space NCS (ncs.type=torsion), automatic parameterization, similar to reference model restraints - Uses top-out potential to allow deviations between NCS related atoms New commands: ============= - phenix.cut_out_density: - carve out a region from map coefficients (MTZ file) - density may be specified as a sphere, a box, or masked around a PDB file - phenix.data_viewer: - visualization of reciprocal-space reflection data - 3D OpenGL view of all data, or 2D view of planes (pseudo-precession photograph) - phenix.morph_model: - Use after molecular replacement to improve model for rebuilding - Morphs model to match a prime-and-switch map Molecular Replacement: ====================== - Phaser (and AutoMR wizard): - significant speedup of fast rotation function Miscellaneous: ============== - phenix.reel: - better integration with main GUI, especially eLBOW - phenix.fobs_minus_fobs_map: - option to disable strict unit cell isomorphism check - phenix.ramalyze: - optional output of PNG images of Ramachandran plots; add --plot to command line arguments - phenix.fetch_pdb: - new flags: --all fetches PDB, CIF, and FASTA files simultaneously --mtz runs phenix.cif_as_mtz on downloaded CIF file For a full list of changes see: http://www.phenix-online.org/documentation/CHANGES Please note that this publication should be used to cite use of Phenix: PHENIX: a comprehensive Python-based system for macromolecular structure solution. P. D. Adams, P. V. Afonine, G. Bunkóczi, V. B. Chen, I. W. Davis, N. Echols, J. J. Headd, L.-W. Hung, G. J. Kapral, R. W. Grosse-Kunstleve, A. J. McCoy, N. W. Moriarty, R. Oeffner, R. J. Read, D. C. Richardson, J. S. Richardson, T. C. Terwilliger and P. H. Zwart. Acta Cryst. D66, 213-221 (2010). Full documentation is available here: http://www.phenix-online.org/documentation/ There is a Phenix bulletin board: http://www.phenix-online.org/mailman/listinfo/phenixbb/ Please consult the installer README file or online documentation for installation instructions. Direct questions and problem reports to the bulletin board or: help(a)phenix-online.org and bugs(a)phenix-online.org Commercial users interested in obtaining access to Phenix should visit the Phenix website for information about the Phenix Industrial Consortium. The development of Phenix is principally funded by the National Institute of General Medical Sciences (NIH) under grant P01-GM063210. We also acknowledge the generous support of the members of the Phenix Industrial Consortium. -- Paul Adams Deputy Division Director, Physical Biosciences Division, Lawrence Berkeley Lab Division Deputy for Biosciences, Advanced Light Source, Lawrence Berkeley Lab Adjunct Professor, Department of Bioengineering, U.C. Berkeley Vice President for Technology, the Joint BioEnergy Institute Head, Berkeley Center for Structural Biology Laboratory Research Manager, ENIGMA Science Focus Area Building 64, Room 248 Tel: 1-510-486-4225, Fax: 1-510-486-5909 http://cci.lbl.gov/paul Lawrence Berkeley Laboratory 1 Cyclotron Road BLDG 64R0121 Berkeley, CA 94720, USA. Executive Assistant: Louise Benvenue [ LBenvenue(a)lbl.gov ][ 1-510-495-2506 ] --

1 0

Twinned Refinement in P 1 21 1
by Yuri 29 Sep '11

29 Sep '11

I dont see spcgrp number 20 on here? Should I be concerned? >>> Unit cell: (44.767, 130.64, 88.33, 90, 104.6, 90) Space group: P 1 (No. 1) Angular tolerance: 5.000 degrees Similar symmetries ================== Symmetry in minimum-lengths cell: C m m m (x+y,2*y,z) (No. 65) Input minimum-lengths cell: (44.767, 88.33, 130.64, 90, 90, 104.6) Symmetry-adapted cell: (44.767, 88.36, 130.64, 90, 90, 104.674) Conventional setting: C m m m (No. 65) Unit cell: (44.767, 170.956, 130.64, 90, 90, 90) Change of basis: x-1/2*z,1/2*z,-y Inverse: x+y,-z,2*y Maximal angular difference: 0.079 degrees Symmetry in minimum-lengths cell: C 1 2/m 1 (x-y,2*x,z) (No. 12) Input minimum-lengths cell: (44.767, 88.33, 130.64, 90, 90, 104.6) Symmetry-adapted cell: (44.767, 88.36, 130.64, 90, 90, 104.674) Conventional setting: C 1 2/m 1 (No. 12) Unit cell: (170.956, 44.767, 130.64, 90, 90, 90) Change of basis: 1/2*z,-x+1/2*z,-y Inverse: x-y,-z,2*x Maximal angular difference: 0.079 degrees Symmetry in minimum-lengths cell: C 1 2/m 1 (x+y,2*y,z) (No. 12) Input minimum-lengths cell: (44.767, 88.33, 130.64, 90, 90, 104.6) Symmetry-adapted cell: (44.767, 88.36, 130.64, 90, 90, 104.674) Conventional setting: C 1 2/m 1 (No. 12) Unit cell: (44.767, 170.956, 130.64, 90, 90, 90) Change of basis: x-1/2*z,1/2*z,-y Inverse: x+y,-z,2*y Maximal angular difference: 0.079 degrees Symmetry in minimum-lengths cell: P 1 1 2/m (No. 10) Input minimum-lengths cell: (44.767, 88.33, 130.64, 90, 90, 104.6) Symmetry-adapted cell: (44.767, 88.33, 130.64, 90, 90, 104.6) Conventional setting: P 1 2/m 1 (No. 10) Unit cell: (44.767, 130.64, 88.33, 90, 104.6, 90) Change of basis: x,y,z Inverse: x,y,z Maximal angular difference: 0.000 degrees Symmetry in minimum-lengths cell: P -1 (No. 2) Input minimum-lengths cell: (44.767, 88.33, 130.64, 90, 90, 104.6) Symmetry-adapted cell: (44.767, 88.33, 130.64, 90, 90, 104.6) Conventional setting: P -1 (No. 2) Unit cell: (44.767, 88.33, 130.64, 90, 90, 104.6) Change of basis: x,z,-y Inverse: x,-z,y Maximal angular difference: 0.000 degrees -- Yuri Pompeu

3 2

Twin Refinement in P 1 21 1
by Yuri 27 Sep '11

27 Sep '11

So in other words, if I input the reflection and the space group as P 1 21 1 (lower symmetry vs the initially attempted C 2 2 21), and I specified the operator in the refinement parameters I should be ok, right? Thank you -- Yuri Pompeu

2 1

Twin Refinement in P 1 21 1
by Yuri 27 Sep '11

27 Sep '11

Hello everyone, I have a twinned data set with PM twin operator h,-k,-h-l I was able to input the twin law into the refinement parameters, however I do not see a place that asks me for the law when generating the Rfree flags. Could someone point me in the right direction? Best, -- Yuri Pompeu

2 1

renaming an holo PDB after the alo one
by Francois Berenger 27 Sep '11

27 Sep '11

Hello, I have one bound complexe (a ligand + a protein in holo conformation). I also have the apo structure for a very similar protein. Is there a tool to create a new PDB, whose coordinates are taken from the holo structure but residue names and numbers are taken from the alo structure (by looking at the corresponding residues in the output of a sequence alignment program)? I thought about doing a script using clustalw for the alignment part, but the task seems not trivial. Thanks a lot, Francois.

1 0

real space residuals
by Bradley Hintze 26 Sep '11

26 Sep '11

Hey phenix, I see a util for real space correlation. Real space residuals are slightly different. Is the a tool to calculate RSR in phenix? -- Bradley J. Hintze Duke University Graduate Student Department of Biochemistry

3 3

map CC given a model and map coefficients
by Francis E Reyes 26 Sep '11

26 Sep '11

Hi all I'd like to compute a resid by resid and overall map CC between a set of coefficients I specify and a model. Any utils to get this to me quickly? phenix.model_vs_data doesn't seem to take in custom coefficients. Thanks! F --------------------------------------------- Francis E. Reyes M.Sc. 215 UCB University of Colorado at Boulder

2 1

Pseudo-merohedral twinning refinement in C2, when beta = 94 degrees
by John Pak 22 Sep '11

22 Sep '11

Hi all, I currently have a dataset that scales well in C2. (a = 224,b = 129,c = 392. beta= 94.5 degrees). Previously, I had thought that these crystals were in P6322 (posted this a few months ago), but I've since grown better crystals, where the beta angle is now 94.5, and not close to 90.0. I'm reasonably confident with the scaling - however, it has always been really difficult to index these crystals in HKL2000. Like the previous crystals, this new C2 dataset looks like it might be perfectly twinned, based on Xtriage's L test. I've solved the structure by MR, and refined it to 3.0 angstroms to R=30%, Rfree=33%, but I am running out of things to fix. I'd like to perform twinning refinement, based on this Xtriage L test result. The problem is that Xtriage cannot find any pseudomerohedral twinning laws, likely because the beta angle is too far from 90 degrees?? For kicks and giggles, I started to enter arbitrary beta angles in Xtriage, and at a beta angle below 93.4, Xtriage can find two twinning laws (closer to 90 degrees, Xtriage can find 5 twin laws). Using "-1/2*h+3/2*k, 1/2*h+1/2*k, -l" in Phenix refine (the beta angle is still 94.5 degrees, I haven't changed this) results in the R and Rfree dropping a lot to 25.5 and 28.6, respectively. Used the same random seed as previous runs to generate the test set from the scalepack file, with the "use lattice symmetry to generate test set" flag enabled. The maps look a lot noisier. So my questions are: 1. Is this twinning refinement valid? I'm really thinking no - since the maps didn't improve. Suppose the maps didn't get worse though - would this twinning refinement then be valid? 2. How do I go about reindexing this dataset in I2? i.e. to see if the beta angle with this unit cell selection is closer to 90 degrees. Would this even help, with respect to pseudomerohedral twinning refinement? Thanks!

2 2

Maps Bias
by Yuri 21 Sep '11

21 Sep '11

I was having a discussion with a friend about maps and model bias. Is the Fobs/PHImodel more biased than 2FOFCWT/PHI2FOFCWt? What about an omit region of these two? Thanks for any opinion on this -- Yuri Pompeu

2 1

Ligand / Loop occupancy refinement
by Antony Oliver 19 Sep '11

19 Sep '11

Dear all, Apologies if this has been covered already in a previous post — but I can't find a suitable response in the archive. We have the following situation... We have soaked a ligand into an apo-crystal, collected diffraction data, and solved the structure. What is apparent from the electron density, is that the ligand isn't in the crystal at full occupancy, it's roughly about 70-80%. However, as the ligand binds, it causes a small conformation change in the active site of the protein, altering the position of around 4 amino acids. What I want to do, and can't quite get phenix.refine to do is the following... Refine the occupancy of the ligand (chain C resname LIG) with the "A" conformation of the active site residues (which should be the same, as they are mutually dependent) — and then refine the occupancy of the "B" conformation of the active site residues — which should all theoretically add up to a total of 1. Could anyone help me with how the occupancies / constrained_group parameters should be set up in this case? With thanks, Tony. --- Dr Antony W Oliver Senior Research Fellow CR-UK DNA Repair Enzymes Group Genome Damage and Stability Centre Science Park Road University of Sussex Falmer, Brighton, BN1 9RQ email: antony.oliver(a)sussex.ac.uk<mailto:[email protected]> tel (office): +44 (0)1273 678349 tel (lab): +44 (0)1273 677512

2 1

bulk solvent mask question
by Laurie Betts 16 Sep '11

16 Sep '11

I notice that during refinement of a MR solution for a protein (search model is 100% homologous to target), a very large ligand-binding cavity (as much as 20 A wide in some parts) in my protein gives weak positive difference density for a ligand soaked into the crystals, unless I turn off bulk solvent modeling. Then the fo-fc difference density is quite strong (2.6 Angstrom data). I think this is due to part of the cavity being included in bulk solvent? Can I edit the mask easily, or flag parts of the protein as NOT exposed to solvent? Thanks Laurie Betts

2 1

Refinement with optimize weights
by Ankit Gupta 16 Sep '11

16 Sep '11

Hello Forum, I am using Phenix 1.7.1-743 to refine my structure. My R/Rfree values for a 2.5 A structure were 0.22/0.305. I used optimize weights along with TLS and the R/Rfree values improved to 0.23/0.29. I wanted to know what parameters should I use for subsequent refinement. I read on phenix bb that one needs to fix the wxc_scale and wxu_scale. Is this correct? What values should be used for these parameters after running optimize wxc/wxu weights. Is there any place where I can get starting values for my future refinement? Also, I would like to know what Fix X-ray/stereochemistry weight and Fix X-ray/ADP weights are for? Thanks for the help.

2 1

Phenix.refine distorts geometry during rigid body refinement
by Ian Slaymaker 14 Sep '11

14 Sep '11

I am bringing a PDB refined in Refmac back to phenix.refine. After an initial rigid body refinement to compare numbers (one chain per group), the geometry (specifically the bond length rmsd) had surprisingly increased. Phenix.validate confirmed the Refmac (initial) geometry numbers and then afterwards confirmed the distorted geometry. What is strange is that when opened in COOT, the initial and final models are completely super imposed. REMARK Start: r_work = 0.3375 r_free = 0.3624 bonds = 0.006 angles = 0.768 REMARK Final: r_work = 0.3379 r_free = 0.3626 bonds = 0.145 angles = 1.019 can you help me figure out what is going on here? These are the final few refinements of a 4.1 A structure. Thanks, Ian

4 4

Linux kernel version 3.0 detection
by Ben Eisenbraun 14 Sep '11

14 Sep '11

Hi, Robert Campbell at Queens pointed out to me that the "machine_type" script bundled with PHENIX doesn't understand the new version 3.0 linux kernel strings. He tweaked his local copy for testing purposes on a Debian testing machine, and just forcing the 3.0 version to use the 2.6 binaries worked: case "$KERNEL" in 2.4.*) kernel_version=2.4 ;; 2.6.*) kernel_version=2.6 ;; 3.0.*) kernel_version=2.6 ;; *) kernel_version="$KERNEL" ;; esac I haven't tested it myself. It looks like Fedora 16 and Ubuntu 11.10, both targeted for November releases, will be using the new 3.0 kernel. -ben -- | Ben Eisenbraun | SBGrid Consortium | http://sbgrid.org | | Harvard Medical School | http://hms.harvard.edu |

2 1

Multiple conformations of disulfide
by Matthew Merski 13 Sep '11

13 Sep '11

Hello, I am trying to refine a surface cys that is alkylated by beta-mercaptoethanol (BME). The BME appears to have two poses and I am having trouble with getting phenix to understand the restraints. When I give the BME two conformations I get the error: Sorry: No atom selected: "chain A and resname ABME and resid 165 and name S2" My link_params file looks like this: refinement.geometry_restraints.edits { drgA_selection1 = chain A and resname ABME and resid 165 and name S2 bond { action = *add atom_selection_1 = $drgA_selection1 atom_selection_2 = $cys215_selection distance_ideal = 2.05 sigma = 0.02 slack = None } Not calling it ABME of course results in a non-unique atom error. Thanks for your help. Matthew Merski UCSF Shoichet Group

2 1

phenix.refine restraint
by Simon Kolstoe 12 Sep '11

12 Sep '11

Dear phenixbb, I'm trying to constrain a ligand-metal bond that is continuously getting altered by phenix.refine (my data's only 3.2A). I've added the following restraint but it just doesn't seem to be holding the bond length during refinement. Is there a different way of telling phenix.refine not to alter this distance? (I'm doing it all on the command line & by editing the .def file). refinement.geometry_restraints.edits { bond { action = *add atom_selection_1 = name CA and chain X and resname CA and resseq 1 atom_selection_2 = name O1 and chain L and resname ROC and resseq 1 distance_ideal = 2.380000 sigma = 0.100 } } Thanks, Simon --------------------------------------------------------------- Dr Simon Kolstoe Laboratory for Protein Crystallography Centre for Amyloidosis & Acute Phase Proteins UCL Medical School Rowland Hill Street, London NW3 2PF Tel: 020 7433 2765 http://www.ucl.ac.uk/~rmhasek ---------------------------------------------------------------

2 1