- phenixbb - phenix-online.org

extend R-free set to higher resolution
by NATHAN THOMSEN 25 Sep '09

25 Sep '09

Is there a phenix utility to extend a phenix generated R-free set to higher resolution (while still using lattice symmetry etc.)? Thanks! Nathan Thomsen UC Berkeley

2 1

Re: [phenixbb] extend R-free set to higher resolution
by Thomas C. Terwilliger 25 Sep '09

25 Sep '09

Hi Nathan, Yes, you can use the GUI to do this (you need a recent nightly build to do it). Just click on "Reflection file editor" and you get an option to extend a free R set, maintaining the free R set you already have. All the best, Tom T >> Is there a phenix utility to extend a phenix generated R-free set to >> higher resolution (while still using lattice symmetry etc.)? Thanks! >> >> Nathan Thomsen >> UC Berkeley >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://www.phenix-online.org/mailman/listinfo/phenixbb >>

1 0

wxu_scale is ignored with TLS
by Joe Krahn 25 Sep '09

25 Sep '09

When using TLS, the wxu_scale is changed, even if it is given in the input file. IMHO, TLS should only modify the default value, and not a user-defined wxu_scale. Thanks, Joe Krahn

3 3

Various problems with phenix.refine
by Joe Krahn 25 Sep '09

25 Sep '09

Here are a few problems I have found using phenix.refine 1.4: If NCS restraint selection in the .def file have newline characters, the REMARK 3 formatting is messed up. It should probably replace all white-space regions with a single space. When NCS restraints contain a residue with alternate conformations, that residue gets excluded from all equivalences, even if the selections contain ALTID to avoid the redundant atoms. If the .def file has a syntax error, phenix.refine dies with an error that the file format is not recognized. Why not print syntax errors? The optimize_wxc and optimize_wxu features are useful, but it is difficult to use the result to define an improved fixed scale on subsequent rounds, without having to optimize every time. It is not obvious whether the scale written in the log represents the new wxc_scale, or if it is a multiplier versus the current wxc_scale. The log also truncates to 2 decimal places. In my case, I only get one significant digit from scale=0.05. I propose that the optimizations write something like "Optimal wxc_scale = 0.533". The input data is written as "refine_data.mtz", not using the serial number, and there seems to be no flag to disable writing it. That means I have to force overwrite to run iterations in the same directory. If each run is in its own directory, the serial number is not very useful for filenames. So, the output MTZ should have the serial number, and you might also consider an option to disable the serial number when separate directories are used. Thanks, Joe Krahn

4 5

specifying column labels for input MTZ file
by Tirumala Kumar Chowdary 24 Sep '09

24 Sep '09

Hi, I am trying to run phenix.refine with a sharp output MTZ file (column labels H K L NUMORB FB PHIB FOM HX HY FH PHIH FP SIGFP HLA HLB HLC HLD HL0 Fcent PHIcent FOMcent Fmap PHImap FHref PHIHref Fav SIGFav). Phenix.refine stopped at 'Multiple equally suitable arrays of observed xray data found. Possible choices: sharp_sad_se_eden.mtz:FP,SIGFP sharp_sad_se_eden.mtz:Fav,SIGFav Please use refinement.input.xray_data.labels' Can someone tell me on how to use this 'refinement.input.xray_data.labels' and how to specify which column labels to use. There is nothing written in the documentation about it. Thanks Tirumal

15 26

question about rfree flag transfer
by crystallogrphy 24 Sep '09

24 Sep '09

Hi, I have built a model against a 3.7A data set by AutoBuild. The resulting model is overall_best.pdb. The data with the rfree test flags is exptl_fobs_phases_freeR_flags.mtz Later on, I got a 3.4A data, which is isomorphous with the 3.7A data. Now I want to transfer the rfree flag from the 3.7A to 3.4A data for refinement. I know one strategy if the data file is in CNS format. That is, use DATAMAN to transfer the flags and generate new test set flags in the high resolution range: 3.7-3.4A. However, when I convert the format by phenix.reflection_file_converter, it fails: phenix.reflection_file_converter exptl_fobs_phases_freeR_flags.mtz --cns=exptl_fobs_phases_freeR_flags.cv --label="PHIM" --label="FP,SIGFP" --label="FOMM" --label="FreeR_flag" --label="FC" However, the resulting data does not have any FP, SIGFP, and FreeR flags. { file: exptl_fobs_phases_freeR_flags.cv.cns } { sg=C2 a=175.385 b=102.472 c=92.057 alpha=90 beta=111.312 gamma=90 } { source of data: exptl_fobs_phases_freeR_flags.mtz:FC } NREFlections=15307 ANOMalous=FALSe DECLare NAME=DATA DOMAin=RECIprocal TYPE=REAL END INDEx 0 0 5 DATA= 0 INDEx 0 0 7 DATA= 0 INDEx 0 0 8 DATA= 0 INDEx 0 0 9 DATA= 0 INDEx 0 0 10 DATA= 0 INDEx 0 0 11 DATA= 0 INDEx 0 0 12 DATA= 0 INDEx 0 0 13 DATA= 0 I also tried mtz2various in CCP4i, the resulting hkl file does not have rfree flags at all. Does anyone know who to convert the exptl_fobs_phases_freeR_flags.mtz into cns format or simple transfer the flags by phenix?

2 1

supported compilers for source builds
by Ben Eisenbraun 23 Sep '09

23 Sep '09

For compiling PHENIX from source, what are the supported compilers? Are the Intel or Portland compilers supported? There are notes in the list archives about using OpenMP with gcc, and hints that the Intel compiler is supported, but I didn't see any docs on building with icc/ifort. In one email Ralf mentions that you might actually lose performance compiling your own version of Solve/Resolve because gfortran optimizes poorly, which leads me to think that the distributed binaries are using ifort for at least the Solve/Resolve bits. Did I overlook some notes on how to do this? -ben -- | Ben Eisenbraun | Software Sysadmin | | Structural Biology Grid | http://sbgrid.org | | Harvard Medical School | http://hms.harvard.edu |

2 3

question about composite SA_omit map in phenix
by crystallogrphy 22 Sep '09

22 Sep '09

Hi, I am using molecular replacement to build a model against a 2.2A resolution data set. However, the template model only has 23% sequence identity. First, I use chainsaw in CCP4i to prune all side chain to be poly-Ala except the identical residues. Then I use AutoMR and rigid body refinement, I got Rfree=50%. The continuous electron density is not bad showing second structures. However, the model has some out-of-registration problem, for example, two proline are in a alpha-helix. In other words, it has model bias. I want to use composite SA-omit to get a less model bias map, not a simple composite omit map or a SA_omit map ommiting a part of the region. Does someone know who to do this? Thanks! Qun Wan postdoc fellow Case Western Reserve University

10 19

heavy atom in pdb for refinement
by crystallogrphy 20 Sep '09

20 Sep '09

Hi, I have Zn atoms in my pdb for refinement. However, it failed with the following scripts. Do I have to create a libary of this heavy atom? Number of atoms with unknown nonbonded energy type symbols: 6 "ATOM 1 Zn Zn Z 1 .*. Zn " "ATOM 2 Zn Zn Z 2 .*. Zn " "ATOM 3 Zn Zn Z 3 .*. Zn " "ATOM 4 Zn Zn Z 4 .*. Zn " "ATOM 5 Zn Zn Z 5 .*. Zn " "ATOM 6 Zn Zn Z 6 .*. Zn " Time building chain proxies: 1.78, per 1000 atoms: 0.45 Here is the format of my pdb: ..... ATOM 1 Zn Zn Z 1 32.028 47.603 13.955 1.00 30.00 Zn ATOM 2 Zn Zn Z 2 29.017 71.434 8.164 1.00 30.00 Zn ATOM 3 Zn Zn Z 3 39.223 49.330 -6.810 1.00 30.00 Zn ATOM 4 Zn Zn Z 4 24.811 89.144 -4.406 1.00 30.00 Zn ATOM 5 Zn Zn Z 5 23.652 72.429 6.125 1.00 30.00 Zn ATOM 6 Zn Zn Z 6 30.132 63.054 -31.918 1.00 30.00 Zn TER 7 Zn 6 ATOM 10 CL CL L 1 20.125 63.405 10.719 1.00 30.00 CL TER 11 CL 1

5 7

RuntimeError : min()
by RAVINDRA D MAKDE 19 Sep '09

19 Sep '09

Hi all, I will highly appreciate your help to resolve the following bugs encountered during running phenix.refine on GUI as well as on command lines. My version of Phenix and platform are: PHENIX: Python-based Hierarchical ENvironment for Integrated Xtallography Version: 1.4 Release tag: 3 cctbx tag: 2008_12_07_1353 Platform: mac-intel-osx osx-10.5.8 The error encountered using phenix refinement GUI (phenix.refine_gui) is as follows: A Python error was detected. This is almost certainly a bug; please send an email containing this message to phenixbb(a)phenix-online.org so we can fix it. Thanks!RuntimeError : min() argument is an empty arrayTraceback:Traceback (most recent call last): File "/usr/local/phenix-1.4-3/phenix/wxGUI2/Base.py", line 274, in OnBrowse self.file_text_field.SetValue(new_path) File "/usr/local/phenix-1.4-3/phenix/wxGUI2/Controls.py", line 974, in SetValue self.phil_handler.propagate_update(self.phil_def_name, new_value) File "/usr/local/phenix-1.4-3/phenix/wxGUI2/PhilInterfaceHandler.py", line 103, in propagate_update self._update_handlers[def_name](self, def_name, value) File "/usr/local/phenix-1.4-3/phenix/wxGUI2/Extensions.py", line 28, in extract_xray_neutron_params phil_def_name = phil_object_name) File "/usr/local/phenix-1.4-3/phenix/wxGUI2/Reflections.py", line 155, in __init__ self.load_reflection_file(file_name, phil_def_name) File "/usr/local/phenix-1.4-3/phenix/wxGUI2/Reflections.py", line 159, in load_reflection_file self.hkl_server = self.hkl_handler.add_hkl_file(file_name, phil_def_name) File "/usr/local/phenix-1.4-3/phenix/wxGUI2/Reflections.py", line 66, in add_hkl_file err=sys.stderr) File "/usr/local/phenix-1.4-3/iotbx/reflection_file_utils.py", line 393, in __init__ force_symmetry=self.force_symmetry)) File "/usr/local/phenix-1.4-3/iotbx/reflection_file_reader.py", line 197, in as_miller_arrays base_array_info=base_array_info) File "/usr/local/phenix-1.4-3/iotbx/mtz/__init__.py", line 434, in as_miller_arrays and flex.min(column_group.sigmas()) > 0): The error encountered using phenix refinement on command line (phenix.refine coordinate.pdb reflection.mtz) is as follows: Traceback (most recent call last): File "/usr/local/phenix-1.4-3/phenix/phenix/command_line/refine.py", line 9, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/usr/local/phenix-1.4-3/phenix/phenix/refinement/command_line.py", line 68, in run log=log) File "/usr/local/phenix-1.4-3/phenix/phenix/refinement/command_line.py", line 231, in __init__ reflection_file_server = self.reflection_file_server() File "/usr/local/phenix-1.4-3/phenix/phenix/refinement/command_line.py", line 566, in reflection_file_server err=self.log) File "/usr/local/phenix-1.4-3/iotbx/reflection_file_utils.py", line 393, in __init__ force_symmetry=self.force_symmetry)) File "/usr/local/phenix-1.4-3/iotbx/reflection_file_reader.py", line 197, in as_miller_arrays base_array_info=base_array_info) File "/usr/local/phenix-1.4-3/iotbx/mtz/__init__.py", line 434, in as_miller_arrays and flex.min(column_group.sigmas()) > 0):RuntimeError: min() argument is an empty array N.B. Earlier, I had encountered the similar error (using command line) using PHENIX version 1.3 on the same platform. Thanks in advance. Regards, Ravi Ravindra D. Makde, PhD Postdoctoral fellow, Tan lab, The Department of Biochemistry and Molecular Biology, The Pennsylvania State University, University Park, PA 16802 USA

2 1