December 2009 - phenixbb - phenix-online.org

16 Dec '09

Hi, has anyone had any experience using phenix with refmac5? I seem to be having issues with the Rfree flag. Phenix sets aside 8.1% of the reflection for Rfree. However when i use refmac for refinement, it reports that its using 91.4% of the reflections for Rfree........? thanks

3 2

polygon output (dev-256)
by Bryan W. Lepore 16 Dec '09

16 Dec '09

[ phenix dev 256 ] where - or should we - be looking for the polygon output .png (the one with the polygon)? its obvious in 1.5-2, so i think its dev-256 only... -bryan

2 3

.cif use/restraints
by Bryan W. Lepore 16 Dec '09

16 Dec '09

a ligand is distorting - as if it is ignoring the .cif and i already have the angles tighted up in the .def. what else should i check? what i did : generated some elbow output. the .cif was provided to phenix.refine on cmd line. the .pdb had the H's and CONNECT's removed, but the .cif retains them. the .pdb was incorporated into a single .pdb input file. -bryan

2 1

Firefox problem
by Phil Evans 15 Dec '09

15 Dec '09

Dear Gurus If I click the Help button on the Phenix GUI I get the following fail messages. What am I doing wrong? I'm running 64-bit Phenix on a 64-bit machine, but our Firefox (as I understand it) is 32-bit, since there is no 64-bit Java Best wishes Phil hex_pre> phenix Using browser /usr/bin/firefox (Gecko:3603): Gtk-WARNING **: /lmb/home/public/xtal/Phenix/Phenix64/phenix-dev-2 49/build/intel-linux-2.6-x86_64/base/lib/gtk-2.0/2.10.0/engines/libclearlooks.so : wrong ELF class: ELFCLASS64 (Gecko:3603): Pango-WARNING **: /lmb/home/public/xtal/Phenix/Phenix64/phenix-dev -249/build/intel-linux-2.6-x86_64/base/lib/pango/1.6.0/modules/pango-basic-fc.so : wrong ELF class: ELFCLASS64 Failed to load Pango module for id: 'BasicScriptEngineFc' (Gecko:3603): Pango-WARNING **: /lmb/home/public/xtal/Phenix/Phenix64/phenix-dev -249/build/intel-linux-2.6-x86_64/base/lib/pango/1.6.0/modules/pango-basic-fc.so : wrong ELF class: ELFCLASS64 Failed to load Pango module for id: 'BasicScriptEngineFc' (Gecko:3603): Pango-WARNING **: /lmb/home/public/xtal/Phenix/Phenix64/phenix-dev -249/build/intel-linux-2.6-x86_64/base/lib/pango/1.6.0/modules/pango-basic-fc.so : wrong ELF class: ELFCLASS64 Failed to load Pango module for id: 'BasicScriptEngineFc' (Gecko:3603): Pango-WARNING **: /lmb/home/public/xtal/Phenix/Phenix64/phenix-dev -249/build/intel-linux-2.6-x86_64/base/lib/pango/1.6.0/modules/pango-basic-fc.so : wrong ELF class: ELFCLASS64 .. snip ... Failed to load Pango module for id: 'BasicScriptEngineFc' (Gecko:3603): Pango-WARNING **: /lmb/home/public/xtal/Phenix/Phenix64/phenix-dev -249/build/intel-linux-2.6-x86_64/base/lib/pango/1.6.0/modules/pango-basic-fc.so : wrong ELF class: ELFCLASS64 Failed to load Pango module for id: 'BasicScriptEngineFc' (Gecko:3603): GLib-GObject-CRITICAL **: g_object_get_qdata: assertion `G_IS_OBJE CT (object)' failed (Gecko:3603): Pango-CRITICAL **: pango_font_describe: assertion `font != NULL' f ailed /usr/lib/firefox-3.0.15/run-mozilla.sh: line 131: 3603 Segmentation fault "$prog" ${1+"$@"}

2 1

phenix and refmac5
by anon 15 Dec '09

15 Dec '09

Hi, recently i used phenix (AutoSol and AutoBuild) on a SAD dataset containing Iodine. I was wanting to do my refinement using Refmac instead of phenix.refine. I was wondering where would i be able to locate the correct .mtz file to input into Refmac5 for rounds of refinement? I was under the impression the .mtz file used for refinement was the one produced after scaling, but im not sure where phenix puts this. Also I was wondering if anyone had any issues with phenix.xtriage carrying out its Matthews coefficient calculation? When using Phenix, it gave me a solution with 3 mols/AU, despite there being 4 mols/AU according to the Matthews coeffiecient calculation in CCP4. Looking back at the xtriage log in autosol it has n_residues= 200 (despite me not specifying anywhere the number of residues, my protein has 147 residues). Realising this I re-ran Autosol, inputting my .mtz file along with my sequence this time, and specifying the number of NCS copies as 4. Looking at the phenix.xtriage log this time results in it recognising the correct number of residues, however the matthews coefficient calculation is still incorrect, claiming 1 mol/Au. Any ideas what may be the problem? ##-------------------------------------------## ## Unit cell defined manually, will ignore ## specification in reflection file: ## From file : None ## From input: (140.124, 77.954, 54.519, 90, 98.549, 90) ##-------------------------------------------## ##-------------------------------------------## ## Space group defined manually, will ignore ## specification in reflection file: ## From file : C 1 2 1 ## From input: C 1 2 1 ##-------------------------------------------## Effective parameters: #phil __ON__ scaling { input { asu_contents { n_residues = 588 n_bases = None n_copies_per_asu = 4 } Matthews coefficient and Solvent content statistics ---------------------------------------------------------------- | Copies | Solvent content | Matthews Coef. | P(solvent cont.) | |--------|-----------------|----------------|------------------| | 1 | 0.448 | 2.226 | 1.000 | ---------------------------------------------------------------- | Best guess : 1 copies in the asu | ----------------------------------------------------------------

2 1

input/output error
by Geoffrey Feld 14 Dec '09

14 Dec '09

Greetings fellow crystallographers, This past weekend, I collected some new data on my new project, indexed, scaled it, and found a phaser solution. The resolution is 3.6 A and I need to find one of the protein components, so I'm slowly refining piece-by-piece so I don't overdo it and flatten out the missing protein. To do this, I simply do rigid body and tls refinements for decreasingly smaller pieces of the protein I've already found. This method worked for me in the past for a similar structure. However, this time, as the refinement is finishing up, I get an input/output error and the program shuts down. Sometimes I get similar errors, and I can fix them by simply running the job again, overwriting the failed job, but this time that doesn't seem to be working. Here's the output from my .log file: Traceback (most recent call last): File "/programs/phenix-1.5-2/phenix/phenix/command_line/refine.py", line 11, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/programs/phenix-1.5-2/phenix/phenix/refinement/command_line.py", line 89, in run call_back_handler=call_back_handler) File "/programs/phenix-1.5-2/phenix/phenix/refinement/driver.py", line 1121, in run call_back_handler = call_back_handler) File "/programs/phenix-1.5-2/phenix/phenix/refinement/strategies.py", line 428, in refinement_machine monitors = monitors) File "/programs/phenix-1.5-2/cctbx_project/mmtbx/refinement/rigid_body.py", line 559, in __init__ out = log) File "/programs/phenix-1.5-2/cctbx_project/mmtbx/refinement/rigid_body.py", line 619, in show translations=t_vec) File "/programs/phenix-1.5-2/cctbx_project/mmtbx/refinement/print_statistics.py", line 161, in show_rigid_body_rotations_and_translations % euler_angle_convention, suffix=frame) File "/programs/phenix-1.5-2/cctbx_project/libtbx/utils.py", line 652, in write file_object.write(str) IOError: [Errno 5] Input/output error Again, the only thing I've done is made a few rigid body and tls groups. Please advise! Thanks, Geoff -- Geoffrey K. Feld Department of Chemistry 492 Stanley Hall University of California, Berkeley "Vigilia pretium libertatis"

2 1

TLS-derived anisoU restraint
by Bryan W. Lepore 14 Dec '09

14 Dec '09

is wxu the easiest - perhaps indirect - way to restrain anisoU's derived from TLS? under the target_weights scope i see : wxu_scale fix_wxu optimize_wxu <-already using this .... b/c i want to pull back some anisoU deviations across C-N. -bryan

2 1

Re: [phenixbb] How to extend R-free set to higher resolution by new phenix-1.5.2
by Thomas C. Terwilliger 14 Dec '09

14 Dec '09

Hi PC, >> >> Hi, >> >> I know this topic was discussed before, but I am still confused and still >> unable to refine my map merged from hires data because phenix.refine >> complains "No array of R-free flags found" >> >> The command was used to generate a hires map. >> phenix.autobuild data=resolve_1.mtz hires_file=hires.sca >> \seq_file=sequence.dat maps_only=True >> >> So, now I have "exptl_fobs_phases_freeR_flags.mtz", "cycle_best_1.mtz", >> "overall_best_denmod_map_coeffs.mtz"...several files in my autobuild >> folder, What this command has done is to take the phases in resolve_1.mtz and combined them with amplitudes from hires.sca and added free R flags to generate exptl_fobs_phases_freeR_flags.mtz . Additionally the file overall_best_denmod_map_coeffs.mtz now contains the above amplitudes along with some kind of density-modified phases. I would caution that this is probably not what you really wanted to do. Instead: Almost always you should use for "data=" a file that contains experimental phases only (not density-modified phases). So it would be better to use the command: phenix.autobuild data=phaser_1.mtz hires_file=hires.sca \ eq_file=sequence.dat maps_only=True which would instead density modify using hl coeffs from modify phaser_1.mtz, extending phase information to the data in hires.sca . Now you will have exptl_fobs_phases_freeR_flags.mtz with (a) experimental amplitudes from hires.sca (b) experimental phases and hl coefficients (where available) from phaser_1.mtz (c) free R flags to the full resolution of hires.sca and overall_best_denmod_map_coeffs.mtz with coefficients for a density-modified map at the resolution of hires.sca. >> >> How do I merge low and high resolution data by >> "phenix.reflection_file_editor"? What files (input and output arrays) >> should I choose? I think reflection_file_editor will only copy arrays (or columns) of data, and not merge several columns of data into one. Note that merging high and low resolution data requires (1) a scale factor between the datasets and (2) possibly a b-factor or anisotropic b-factor relating the two datasets. It cannot be properly done simply by merging reflection files. The reflection_flle_editor will add free R flags if they are missing or extend them for you if they are present only to low resolution, however. This allows you to take a low-res dataset with free R flags and a high-res dataset with no free R and to transfer the existing low-res free R flags to the high-res dataset and to extend them to full resolution. If you want, you can use phenix.solve to scale low and high resolution data by pretending that they are part of a mad dataset: make a file "scale.csh" like this, after editing to match your case: #!/bin/csh -f cp $PHENIX/solve_resolve/ext_ref_files/c2.sym . phenix.solve<<EOD resolution 20 3.0 cell 76 28 42 90 103 90 symfile c2.sym readdenzo ! readformatted/readdenzo/readtrek/readccp4_unmerged premerged ! premerged/ unmerged read_intensities ! read_intensities/read_amplitudes lambda 1 mad_atom c fprprv_mad 1. fprimv_mad 1. rawmadfile low.sca rawmadfile high.sca skip_solve sad EOD Then run with chmod +x scale.csh ./scale.csh Now solve_1.sca should have your scaled, merged data. I hope that helps! -Tom T

3 2

How to extend R-free set to higher resolution by new phenix-1.5.2
by Pei-Chun Lin 12 Dec '09

12 Dec '09

Hi, I know this topic was discussed before, but I am still confused and still unable to refine my map merged from hires data because phenix.refine complains "No array of R-free flags found" The command was used to generate a hires map. phenix.autobuild data=resolve_1.mtz hires_file=hires.sca \seq_file=sequence.dat maps_only=True So, now I have "exptl_fobs_phases_freeR_flags.mtz", "cycle_best_1.mtz", "overall_best_denmod_map_coeffs.mtz"...several files in my autobuild folder, How do I merge low and high resolution data by "phenix.reflection_file_editor"? What files (input and output arrays) should I choose? Detail description will be very appreciate. Thanks for answer my question. PC _________________________________________________________________ Windows Live: Make it easier for your friends to see what you’re up to on Facebook. http://www.microsoft.com/middleeast/windows/windowslive/see-it-in-action/so…

1 0

Side chain at special position
by Nate Nicely 11 Dec '09

11 Dec '09

I have a structure wherein the data show alternate conformations for a particular Trp side chain. The problem is that each of these conformations also occupies a special position: rotamer A clashes with rotamer A of the symmetry object, and similar for rotamer B. I'm attaching a screenshot wherein I tried to capture what's going on. The prerefined structure has a yellow carbon backbone, and you can see how refinement pushes rotamer A out of its rightful place in the map. I readjusted it and added the alternate conformation as seen in the green carbon backbone, but refinement after that just pushes apart both rotamers (not shown). I'd appreciate advice how to handle refining this structure in phenix. Thanks, Nate N

4 3