- phenixbb - phenix-online.org

b factors reporting
by Maia Cherney 22 Apr '11

22 Apr '11

To: PHENIX user mailing list<phenixbb(a)phenix-online.org> Date: Thu, 21 Apr 2011 15:09:19 -0600 Subject: Re: [phenixbb] B factor for DNA and Protein Pavel Afonine wrote: > > Yes, this trivial combination of two commands will do it: > > phenix.pdbtools model.pdb keep=protein output.file_name=protein_only.pdb > phenix.pdbtools model.pdb keep=rna output.file_name=rna_only.pdb > phenix.pdbtools model.pdb keep=water output.file_name=water_only.pdb > Dear Pavel, These commands did not work for me (version 1.7-650). I need to report average b factors separately for the protein, water, several ligands and ions, all of them in one chain A. Traceback (most recent call last): File "/raid/programs/linux/ccp4-6.1.0/ccp4-6.1.3/lib/cctbx/cctbx_sources/mmtbx/mmtbx/command_line/pdbtools.py", line 3, in<module> from mmtbx import pdbtools File "/raid/programs/linux/ccp4-6.1.0/ccp4-6.1.3/lib/cctbx/cctbx_sources/mmtbx/mmtbx/__init__.py", line 3, in<module> import cctbx.array_family.flex File "/raid/programs/linux/ccp4-6.1.0/ccp4-6.1.3/lib/cctbx/cctbx_sources/cctbx/cctbx/array_family/flex.py", line 1, in<module> import scitbx.array_family.flex File "/raid/programs/linux/ccp4-6.1.0/ccp4-6.1.3/lib/cctbx/cctbx_sources/scitbx/scitbx/array_family/flex.py", line 2, in<module> import boost.optional File "/raid/programs/linux/ccp4-6.1.0/ccp4-6.1.3/lib/cctbx/cctbx_sources/boost_adaptbx/boost/optional.py", line 1, in<module> import boost.python File "/raid/programs/linux/ccp4-6.1.0/ccp4-6.1.3/lib/cctbx/cctbx_sources/boost_adaptbx/boost/python.py", line 45, in<module> ext = import_ext("boost_python_meta_ext") File "/raid/programs/linux/ccp4-6.1.0/ccp4-6.1.3/lib/cctbx/cctbx_sources/boost_adaptbx/boost/python.py", line 24, in import_ext + [" "+p for p in sys.path])) ImportError: __import__("boost_python_meta_ext"): /raid/programs/linux/ccp4-6.1.0/ccp4-6.1.3/src/phaser/phaser-2.1.4/build/intel-linux/lib/libboost_python.so: undefined symbol: PyUnicodeUCS2_AsWideChar sys.path: /raid/programs/linux/ccp4-6.1.0/ccp4-6.1.3/lib/cctbx/cctbx_sources/mmtbx/mmtbx/command_line /raid/programs/linux/ccp4-6.1.0/ccp4-6.1.3/lib/cctbx/cctbx_sources/mmtbx /raid/programs/linux/ccp4-6.1.0/ccp4-6.1.3/lib/cctbx/cctbx_sources/iotbx /raid/programs/linux/ccp4-6.1.0/ccp4-6.1.3/lib/cctbx/cctbx_sources/PyCifRW /raid/programs/linux/ccp4-6.1.0/ccp4-6.1.3/src/phaser/source/phaser /raid/programs/linux/ccp4-6.1.0/ccp4-6.1.3/lib/cctbx/cctbx_sources/cctbx /raid/programs/linux/ccp4-6.1.0/ccp4-6.1.3/lib/cctbx/cctbx_sources/scitbx /raid/programs/linux/ccp4-6.1.0/ccp4-6.1.3/lib/cctbx/cctbx_sources/boost_adaptbx /raid/programs/linux/ccp4-6.1.0/ccp4-6.1.3/lib/cctbx/cctbx_sources/libtbx /raid/programs/linux/ccp4-6.1.0/ccp4-6.1.3/src/phaser/phaser-2.1.4/build/intel-linux/lib /raid/programs/linux/ccp4-6.1.0/ccp4-6.1.3/share/python /usr/lib/python2.6 /usr/lib/python2.6/plat-linux2 /usr/lib/python2.6/lib-tk /usr/lib/python2.6/lib-old /usr/lib/python2.6/lib-dynload /usr/lib/python2.6/dist-packages /usr/lib/python2.6/dist-packages/Numeric /usr/lib/python2.6/dist-packages/PIL /usr/lib/python2.6/dist-packages/gst-0.10 /var/lib/python-support/python2.6 /usr/lib/python2.6/dist-packages/gtk-2.0 /var/lib/python-support/python2.6/gtk-2.0 /usr/local/lib/python2.6/dist-packages

3 2

Creating cif file of ligands by Elbow/ReadySet
by Ngan Nguyen Bich 21 Apr '11

21 Apr '11

Dear Phenix developers, I'm quite new in crystallography and use PhenixGUI for refinement my protein. I got a problem when I had to create cif file of the new ligand to feed in phenix_refine by Elbow/ReadySet I saw that program creates some bonds shorter or longer compared to the statistic value from the CSD for example. Here are some examples when I play with Elbow (same problem with ReadySet) 1. when I feed program the chemical string, (template PDB file is the mylidand_idea.pdb) CR2 C3 CA3 single 1.569 0.038 and if I choose AM1 optimization CR2 C3 CA3 single 1.602 0.040 2. when I feed program the .pdb file, choose Use geometry for Initial geometry. GYC N CA1 single 1.409 0.020 * GYC N HN1 single 0.902 0.020 GYC N HN2 single 0.901 0.020 GYC SG1 CB1 single 1.797 0.020 GYC SG1 HG1 single 1.204 0.020 GYC CB1 CA1 single 1.500 0.020 GYC CB1 HB11 single 1.006 0.020 GYC CB1 HB12 single 0.972 0.020 GYC CA1 C1 single 1.528 0.020 GYC CA1 HA1 single 1.010 0.020 GYC C1 N2 aromatic 1.412 0.020 GYC C1 N3 aromatic 1.304 0.020 GYC N2 CA2 aromatic 1.505 0.020 ** GYC N3 C2 aromatic 1.349 0.020 GYC N3 CA3 single 1.423 0.020 GYC C2 O2 double 1.370 0.020 GYC C2 CA2 aromatic 1.509 0.020 GYC CA2 CB2 single 1.532 0.020 GYC CA3 C single 1.496 0.020 It's strange to me when comparing * and ** and if I run AM1 GYC CB1 CA1 single 1.533 0.035 GYC CA3 C single 1.564 0.035 as I know, the C-C single bond should be about 1.51. Does anyone get a problem like me? I wonder if I didn't do properly or some thing else, or any suggestion for my case. One more questions is that afterward I found cif file (seemly mmcif file) of my ligand in phenix-1.7-650/chem_data/chemical_components And I'm wondering that How the program use these cif file, or I can use them in which way? Thank you in advance. Best regards, Ngan Nguyen Bich KULeuven Belgium

2 1

unit cell discrepancies
by Michael Hothorn 21 Apr '11

21 Apr '11

Dear all, I am sure there is a simple explanation for this, but I cannot find a solution. I am reading in an .mtz file in phenix.refine. mtzdump reports the following unit cell constants on the file test.mtz itself: 175.0930 67.2510 119.0500 90.0000 121.5500 90.0000 The refined .pdb file and .maps report however a slightly different unit cell, as does the logfile >snip Miller array info: test.mtz:FP,SIGFP Observation type: xray.amplitude Type of data: double, size=39702 Type of sigmas: double, size=39702 Number of Miller indices: 39702 Anomalous flag: False Unit cell: (175.17, 67.236, 119.246, 90, 121.416, 90) What am I missing? Thanks! Michael

6 16

Installation issue
by Simon Kolstoe 21 Apr '11

21 Apr '11

Hi, I've just installed the latest stable version of phenix on a new macbook pro using phenix-installer-1.7-650-mac-intel-osx.tar However, as I like running phenix.refine on the command line, I normally source the script phenix_env.sh which on my other computer is in the directory /Applications/PHENIX-1.7-650/Contents/phenix-1.7-650/ However the latest install doesn't seem to have the final phenix-1.7-650 directory! Has phenix_env.sh moved somewhere else and if so where? Thanks, Simon

2 2

Out of Memory Errors
by Heather Condurso 21 Apr '11

21 Apr '11

Hello all, I am currently running the Phenix GUI version 1.7-650 on an iMac with an 3.2Ghz Intel Core i3 processor, 8GB of memory operating under OS 10.6. During many of my runs (mostly AutoMR runs but sometimes refine) I get out of memory errors. I have my system booted in 64 bit mode and I watch the activity monitor and find that phenix doesn't use much more than 2GB of memory. My search model is very large. My unit cell contains 4 trimers (each monomer is 48kDa). I hadn't encountered this issue until I began to increase my search model from one trimer to four. Do you have any suggestions as to how I could avoid this error without minimizing my search model? Often my AutoMR runs get the out of memory error at the start of refinement, but just now the memory was exhausted during the fast rotation function. During this run I used all default parameters with a max res of 4.0 and RMSD of 1.0. Any advice you can give would be greatly appreciated! Heather

3 2

R free test set size
by Yuri 21 Apr '11

21 Apr '11

hello everyone, I just noticed that the latest model I am working on has only 2% of reflections marked as free R after refinement! I usually give the GUI the .sca files out of HKL2000 that do not have an Rfree array pre-defined. The GUI lets you know about the lack of it and asks me if I want generate a test set of 10%. I say yes. However the pdb header shows the only 2.14% are chosen for free R. I see that some of my other models (lower resolution) have test arrays of around 7-10% as I would expect. How can I explicitly define the size I want? Thanks -- Yuri Pompeu

2 1

optimizing weights
by Alexandra Deaconescu 20 Apr '11

20 Apr '11

Dear all: I would like to optimize my weights during Phenix refinement. I am running Phenix 1.6.4-486 using the gui, and I have checked the "optimize weights", but I cannot see the results of the optimization in the output, so I do not know the relationship between Rfree/geometry and wc and wxc_scale (I might not want to pick the automatically selected wxc_scale for the lowest Rfree, but would like to find a reasonable values to give me good geometry). Some earlier posts indicated that this is not output in the log file, but I am running a more recent version... Also, in case I wanted to adjust the weights manually, do I have set wxc_scale AND check the fix wc at the same time in the gui? Or just one of them? What would be typical values to start with at 2.8A (even if the answer depends on the particular data)? Same question for the ADP refinement.... I am working at 2.8A and my B factors are high (higher I think than what one would expect at this resolution). Thanks for your help, Alex -- Alexandra M. Deaconescu, Ph.D. Research Associate, HHMI c/o Grigorieff Laboratory Rosenstiel Center MS029 Brandeis University 415 South St. Waltham, MA 02454 phone: 781-736-2498 For FEDEX Brandeis University Kalman Receiving Dock 415 South St. Waltham, MA 02454

2 1

obtaining f'/f'' from solved structures and comparison to a fluorescence scan
by Francis E Reyes 20 Apr '11

20 Apr '11

Hi phenixer's I'm in need of reliable f'/f'' for my heavy atom (that isn't Sasaki tables which I believe is for free, unbound atoms). I've used group_anomalous in phenix.refine before so I'm sure a phenix utility is capable of extracting it. I wanted to if it's possible to extract it from several solved structures (I have the DANO/SIGDANO, F(+)/F(-), or if needed the original I(+)/I(-)) I'm curious to see if there's wide variation and whether it'd be as reliable (or even better) than actually taking a fluorescence scan (and of course use the values for my current phasing challenge as a starting point). F --------------------------------------------- Francis E. Reyes M.Sc. 215 UCB University of Colorado at Boulder

2 1

AutoSol
by Nora Cronin 18 Apr '11

18 Apr '11

Is have 2 difference anomalous datasets with two different heavy atoms. I can input 2 SAD datasets in AutoSol using the Wavelength and Group = keywords. How do I input 2 different crystal_info keywords since this keyword does not accept the Group= option? cheers Nora Dr N B Cronin X-Ray Laboratory Manager Tel. (+44) 020 7153 5445 Section of Structural Biology Fax. (+44) 020 7153 5457 Chester Beatty Laboratories mobile: (+44) 07787554059 The Institute of Cancer Research 237 Fulham Road London SW3 6JB The Institute of Cancer Research: Royal Cancer Hospital, a charitable Company Limited by Guarantee, Registered in England under Company No. 534147 with its Registered Office at 123 Old Brompton Road, London SW7 3RP. This e-mail message is confidential and for use by the addressee only. If the message is received by anyone other than the addressee, please return the message to the sender by replying to it and then delete the message from your computer and network.

1 0

unusual maps and Rfree for 1.75 structure
by Yuri 18 Apr '11

18 Apr '11

>>you did not mention the resolution... From your email I guess it is 1.9A, is it correct? It is 1.75A. But like I said the maps right now do not look like those of 1.75 resolution. my 1.9 maps looked way better... The P1 strategy is valid. I will also try processing without the screw axis P 4 2 2, or maybe even P 4 21 2. It looks like the other systematic absences are ok (they are pretty weak reflections). -- Yuri Pompeu

1 0

unusual maps and Rfree for 1.75 structure
by Yuri 17 Apr '11

17 Apr '11

I am solving a structure of an enzyme that I have solved many times before, the only difference is a single mutation in the active site. The crystal grew in basically exact conditions (macroscopically speaking they look different) as the other mutants I had. I scaled the data P 43 21 2 Rsym 0.068. Like all the other data sets I had before. This time however, the initial model has Rfree of about 0.27 and the maps look more like those of a 2.5 - 2.8 structure. My experinence with this enzyme is that I should be looking at around 0.22 or lower Rfree, and pretty clear maps as I saw for my 1.9 structure. Xtriage does not detect any pathologies. I rececked my scaling log file, and it lists the systematic absences for the data. I see some that are not that "weak" for instance : h k l I/sigma 0 0 1 6.0 0 0 2 5.8 0 0 13 7.8 So I am inclined to reprocess my data in P 4 2 2. Are these "symptoms" of a possible wrong space group assignment? Is there something else obvious I am overlooking? Thank you -- Yuri Pompeu

3 2

Re: [phenixbb] Autobuild error -- merging results of nearly-finished OMIT and autobuild runs
by Thomas C. Terwilliger 15 Apr '11

15 Apr '11

Dear Nick, I'm sorry for the error ("[Errno 5] Input/output error") that you are getting with your i7-based PC's with Fedora 14. I have not seen this error before, but we can try to debug it with "debug=true" and also I have a work-around for combining the files that you already have created so you don't need to rerun your whole jobs. Can you possibly run this quick test job to try and get it to fail: phenix.autobuild $PHENIX/phenix_regression/wizards/perfect.mtz \ input_labels="FP SIGFP" \ map_file=$PHENIX/phenix_regression/wizards/perfect.mtz \ input_map_labels="FP PHIC FOM" \ seq_file=$PHENIX/phenix_regression/wizards/seq.dat \ resolution=3.0 \ solvent_fraction=.6 \ model=$PHENIX/phenix_regression/wizards/coords.pdb \ rebuild_in_place=True \ n_cycle_rebuild_min=1 \ n_cycle_rebuild_max=1 \ multiple_models=True \ multiple_models_number=1 \ multiple_models_group_number=2 \ highest_resno=200 \ debug=true This takes about 5 minutes. If it fails, can you send me (terwilliger(a)lanl.gov, not the list) the end out the output where it gives a traceback of where the error occurred? Combining files from a nearly-complete autobuild run with rebuild-in-place=true You can combine models with autobuild, providing all models have exactly the same set of atoms (aside from any solvent) by editing this script (basically you run a dummy autobuild run to create a directory and database entries, then you copy your files there, then you run autobuild and tell it to carry on and do the combine step): -------------------------------------------------------- #!/bin/csh -f #COMBINE_MODELS SCRIPT if (-d PDS || -d AutoBuild_run_1_) then echo "Please run in a directory without PDS or AutoBuild_run_1_" exit 1 endif echo "Setting up combine models with a dummy run. NOTE: multiple_models_group_number must be correct" phenix.autobuild fobs.mtz multiple_models=true seq_file=seq.dat combine_only=true multiple_models_group_number=2 multiple_models_number=1>& dummy_autobuild.log echo "Copying files to AutoBuild_run_1_/MULTIPLE_MODELS" mkdir AutoBuild_run_1_/MULTIPLE_MODELS cp coords1.pdb AutoBuild_run_1_/MULTIPLE_MODELS/initial_model.pdb_1_1 cp coords2.pdb AutoBuild_run_1_/MULTIPLE_MODELS/initial_model.pdb_1_2 cp map_coeffs_1.mtz AutoBuild_run_1_/MULTIPLE_MODELS/initial_model.mtz_1_1 cp map_coeffs_2.mtz AutoBuild_run_1_/MULTIPLE_MODELS/initial_model.mtz_1_2 ls AutoBuild_run_1_/MULTIPLE_MODELS/ echo "Running autobuild to combine files in AutoBuild_run_1_/MULTIPLE_MODELS" phenix.autobuild combine_only=true seq_file=seq.dat carry_on=true run=1 > autobuild_combine.log # END OF COMBINE_MODELS SCRIPT ------------------------------------------------------- Combining OMIT runs that have been run on separate computers or that come from a nearly-completed omit run: You can do this by copying all the individual mtz files with map coefficients for omit regions to a single directory and running a script like this one: (NOTE: you need to ensure that the OMIT regions are defined the same in the runs where you got your overall_best_denmod_map_coeffs.mtz_OMIT_REGION_1 etc files and this run. You ensure that with the n_xyz command that sets the grid. You can copy this from one of your resolve log files created when you ran your omit (i.e., AutoBuild_run_1_/TEMP0/AutoBuild_run_1_/TEMP0/resolve.log will have a line "nu nv nw: 32 32 32 " and you copy those numbers). ------------------------------------ #!/bin/csh -f # COMBINE OMIT SCRIPT phenix.resolve<<EOD hklin exptl_fobs_phases_freeR_flags.mtz labin FP=FP SIGFP=SIGFP n_xyz 32 32 32 # YOU MUST SET THIS BASED ON THE nu nv nw in a resolve log file. solvent_content 0.85 no_build ha_file NONE combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_1 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_2 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_3 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_4 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_5 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_6 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_7 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_8 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_9 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_10 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_11 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_12 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_13 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_14 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_15 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_16 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_17 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_18 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_19 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_20 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_21 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_22 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_23 combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_24 omit EOD # END OF COMBINE OMIT SCRIPT All the best, Tom T >> Dear Phenix Developers, >> >> I have been experiencing problems with autobuild whenever I try to rebuild >> an MR model or calculate a composite OMIT map. These tasks always end the >> same way: with "[Errno 5] Input/output error". The failure always happens >> when the run is basically finished and trying to combine the output from >> multiple tasks. My first thought was that my hard drive was full or fixing >> to crash, but these problems happen across several machines and have persisted through multiple official and nightly Phenix releases (currently >> running 1.7-634). The computers that cannot run autobuild successfully are >> both core i7-based PCs running Fedora 14 (though these problems also occurred running earlier Fedora versions). My older xeon-based PC runs autobuild without any trouble, which makes me think this a core i7-related >> problem. I have searched the PhenixBB archive and have not turned up any >> similar problems - have any other users experienced trouble with Phenix and core i7? Is there a work-around I ca >> n use to manually combine the various outputs? The core i7 machines are >> so much faster, I would really like to solve this problem so that I can stop relying on my aging xeon PC for these long calculations. I have included the end of two log files for reference. Also, autosol runs fine >> on the core i7 computers and successfully builds wonderful models from my >> SeMet SAD data. Thank you for Phenix and for any advice you can give. >> >> Best wishes, >> Nick Silvaggi >> >> P.S. I apologize if this message appears twice - I sent it once from an address that was not registered with the BB and decided to re-send from a >> registered address. >> >> Autobuild after MR: >> ----------------------------------------------------------------------------------------------------- Standard build in parallel starting with refine.pdb_1 >> This is the first try at building this model >> Try: 1 building 1 model >> Try: 2 building 1 model >> Try: 3 building 1 model >> Running up to 3 jobs in parallel... with total of 3 jobs >> >> Splitting work into 3 jobs and running with 3 processors using sh background=True in >> /home/nicks/projects/r51/phenix_mr/AutoBuild_run_2_/TEMP0 >> >> Starting job 1...Log will be: >> /home/nicks/projects/r51/phenix_mr/AutoBuild_run_2_/TEMP0/RUN_FILE_1.log Starting job 2...Log will be: >> /home/nicks/projects/r51/phenix_mr/AutoBuild_run_2_/TEMP0/RUN_FILE_2.log Starting job 3...Log will be: >> /home/nicks/projects/r51/phenix_mr/AutoBuild_run_2_/TEMP0/RUN_FILE_3.log Collecting all runs now >> *************ERROR ENDING ******************* >> >> ******************************************************************************** Failed to carry out AutoBuild_build_cycle: >> >> [Errno 5] Input/output error >> ******************************************************************************** >> >> *************ERROR ENDING ******************* >> >> >> Calculating SA composite omit map: >> ----------------------------------------------------------------------------------------------------- Starting job 34...Log will be: >> /home/nicks/projects/r51/phenix_mr/AutoBuild_run_3_/TEMP0/RUN_FILE_34.log Starting job 35...Log will be: >> /home/nicks/projects/r51/phenix_mr/AutoBuild_run_3_/TEMP0/RUN_FILE_35.log Starting job 36...Log will be: >> /home/nicks/projects/r51/phenix_mr/AutoBuild_run_3_/TEMP0/RUN_FILE_36.log Collecting all runs now >> *************ERROR ENDING ******************* >> >> ******************************************************************************** Failed to carry out AutoBuild_parallel_omit: >> >> [Errno 5] Input/output error >> ******************************************************************************** >> >> *************ERROR ENDING ******************* >> >> Nicholas R Silvaggi, Ph.D. >> University of Wisconsin, Milwaukee >> Department of Chemistry and Biochemistry >> 3210 North Cramer St >> Milwaukee, WI 53211 >> >> Office: Chemistry 272B >> Phone: 414-229-2647 >> silvaggi(a)uwm.edu >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://phenix-online.org/mailman/listinfo/phenixbb >>

1 0

Autobuild error
by Nicholas R Silvaggi 15 Apr '11

15 Apr '11

Dear Phenix Developers, I have been experiencing problems with autobuild whenever I try to rebuild an MR model or calculate a composite OMIT map. These tasks always end the same way: with "[Errno 5] Input/output error". The failure always happens when the run is basically finished and trying to combine the output from multiple tasks. My first thought was that my hard drive was full or fixing to crash, but these problems happen across several machines and have persisted through multiple official and nightly Phenix releases (currently running 1.7-634). The computers that cannot run autobuild successfully are both core i7-based PCs running Fedora 14 (though these problems also occurred running earlier Fedora versions). My older xeon-based PC runs autobuild without any trouble, which makes me think this a core i7-related problem. I have searched the PhenixBB archive and have not turned up any similar problems - have any other users experienced trouble with Phenix and core i7? Is there a work-around I can use to manually combine the various outputs? The core i7 machines are so much faster, I would really like to solve this problem so that I can stop relying on my aging xeon PC for these long calculations. I have included the end of two log files for reference. Also, autosol runs fine on the core i7 computers and successfully builds wonderful models from my SeMet SAD data. Thank you for Phenix and for any advice you can give. Best wishes, Nick Silvaggi P.S. I apologize if this message appears twice - I sent it once from an address that was not registered with the BB and decided to re-send from a registered address. Autobuild after MR: ----------------------------------------------------------------------------------------------------- Standard build in parallel starting with refine.pdb_1 This is the first try at building this model Try: 1 building 1 model Try: 2 building 1 model Try: 3 building 1 model Running up to 3 jobs in parallel... with total of 3 jobs Splitting work into 3 jobs and running with 3 processors using sh background=True in /home/nicks/projects/r51/phenix_mr/AutoBuild_run_2_/TEMP0 Starting job 1...Log will be: /home/nicks/projects/r51/phenix_mr/AutoBuild_run_2_/TEMP0/RUN_FILE_1.log Starting job 2...Log will be: /home/nicks/projects/r51/phenix_mr/AutoBuild_run_2_/TEMP0/RUN_FILE_2.log Starting job 3...Log will be: /home/nicks/projects/r51/phenix_mr/AutoBuild_run_2_/TEMP0/RUN_FILE_3.log Collecting all runs now *************ERROR ENDING ******************* ******************************************************************************** Failed to carry out AutoBuild_build_cycle: [Errno 5] Input/output error ******************************************************************************** *************ERROR ENDING ******************* Calculating SA composite omit map: ----------------------------------------------------------------------------------------------------- Starting job 34...Log will be: /home/nicks/projects/r51/phenix_mr/AutoBuild_run_3_/TEMP0/RUN_FILE_34.log Starting job 35...Log will be: /home/nicks/projects/r51/phenix_mr/AutoBuild_run_3_/TEMP0/RUN_FILE_35.log Starting job 36...Log will be: /home/nicks/projects/r51/phenix_mr/AutoBuild_run_3_/TEMP0/RUN_FILE_36.log Collecting all runs now *************ERROR ENDING ******************* ******************************************************************************** Failed to carry out AutoBuild_parallel_omit: [Errno 5] Input/output error ******************************************************************************** *************ERROR ENDING ******************* Nicholas R Silvaggi, Ph.D. University of Wisconsin, Milwaukee Department of Chemistry and Biochemistry 3210 North Cramer St Milwaukee, WI 53211 Office: Chemistry 272B Phone: 414-229-2647 silvaggi(a)uwm.edu

1 0

parallel structure comparison
by Kendall Nettles 14 Apr '11

14 Apr '11

Dear PHENIX developers, In this weeks build, when I launch coot from the structure comparison results, no models or maps are loaded. It works in the last major release. Also, in the major release version the ligands are stripped out. Is this a bug? If not, I would like to request that the ligands be left in the structures. Thanks so much Kendall

2 1

The individual site refinement for Low resolution data (~3.8A)
by Zhang yu 14 Apr '11

14 Apr '11

Dear phenixer, Phenix suggested "Individual site" refinement is only suitable for high or midi resolution data (<3.5A). What is the alternate strategy to refine the atom coordinates of low resolution data? Yu

5 17

occupancy refinement when alternate conformers clash
by Kendall Nettles 13 Apr '11

13 Apr '11

Hi, I am refining an RNA structure where on set of base pairs has alternate conformers, and the B alternates from the paired bases clash. That is, the structure should be viewed as A alternate in one base paired with B alternate in the other base. The refinement blew both B alts out of the density to prevent their clash. How can I fix this? Thanks for your help, Kendall Nettles

2 1

O-linked glycosylation refinement in Phenix
by Qiang Chen 12 Apr '11

12 Apr '11

Dear all, I'm refining a structure which has both N-linked and O-linked glycosylation. I use Phenix to do the refinement. It works well for the N-linked NAG. I defined the link as the following: apply_cif_link { data_link = "NAG-ASN" residue_selection_1 = "chain A and resname NAG and resid 701" residue_selection_2 = "chain A and resname ASN and resid 518" } But it doesn't work when I defined the fucose-threonine link as the similar way: apply_cif_link { data_link = "FUC-THR" residue_selection_1 = "chain A and resname FUC and resid 801" residue_selection_2 = "chain A and resname THR and resid 596" } The error message is "missing CIF link: data_link_FUC-THR". Does anyone have the cif link file of FUC-THR? Thanks a lot! Qiang The information in this e-mail is intended only for the person to whom it is addressed. If you believe this e-mail was sent to you in error and the e-mail contains patient information, please contact the Partners Compliance HelpLine at http://www.partners.org/complianceline . If the e-mail was sent to you in error but does not contain patient information, please contact the sender and properly dispose of the e-mail.

4 5

(Anisotropic ADP vs. isotropic ADP) coupling with occupancy refinement
by Young-Jin Cho 12 Apr '11

12 Apr '11

Dear phenix developer and users, I am just wondering while I was taking a small molecule X-ray course. I used to run refinement with anisotropic ADP when it is high resolution about better than 1.5 or 1.6 A. My question is when you do this with occupancy refinement, does phenix.refine take care the alternative conformations with aniostropic ADP (I guess so) even at low occupancy? Due to lack of electrons (scattering factor), I experienced that small molecule refinement applies Anisotropic ADP for the atoms (about more than 2 e- in occupancy) while applying isotropic ADP for the less occupied atoms (in an alternative conformation). Would it help if we apply this method to our macro molecular system when you have residues that have low occupied alternative residues? Is so, could we have something like script or else either to modify a pdb file or to make this in consideration during refinement? Any information or comments are welcome. Thanks in advance, Young-Jin

2 1

help me to install phenix on ubuntu10.4
by 王瑞 11 Apr '11

11 Apr '11

Excuse me,I'm a freshman ,could anyone tell me how to install phenix on ubuntu10.4? I have download phenix-installer-1.7-650- intel-linux-2.6-x86_64-fc12.tar.gz. when I install it in terminal,there was a error notice:"./install:no such file or directory"!Could anyone help me? Many thanks!

2 1

phenix don't start
by albert 10 Apr '11

10 Apr '11

Dear: I installed Phenix-1.7-650 in my SUSE Linux Enterprise 10 SP3 X64 system from source code. When I source /soft/phenix-1.7/phenix-1.7-650/phenix_env and try to run it by typing phenix in terminal, it is said: Traceback (most recent call last): File "/soft/phenix-1.7/phenix-1.7-650/phenix/wxGUI2/command_line/main.py", line 13, in <module> run(sys.argv[1:]) File "/soft/phenix-1.7/phenix-1.7-650/phenix/wxGUI2/command_line/main.py", line 7, in run from wxGUI2 import App File "/soft/phenix-1.7/phenix-1.7-650/phenix/wxGUI2/App.py", line 11, in <module> import wx ImportError: No module named wx The python in my OS is 2.4 with wxpython installed. I don't know why it cannot start. could you please give me some advices? THX

2 1

secondary structure restraints
by Panne Daniel 08 Apr '11

08 Apr '11

I am having problems with the secondary structure restraints in phenix.refine (1.6.4-486). I wrote a ss restraints file and use this as an input file. In the log file I can see that the file is read but the atom selections used by phenix.refine do not correspond to those set in the restraints file. I have turned off automatic ss detection. Any advice? Thanks, Daniel

3 3

difficulties with the refinement of a Cadmium containing pdb file
by Vellieux Frederic 08 Apr '11

08 Apr '11

Hi all, I have some difficulties refining (what happens is an error status, i.e. no refinement taking place) of a coordinate file that contains Cadmium ions. The person here who has already refined a structure containing Cd ions using Phenix is not here at the moment (otherwise I would have asked him). In case you wonder: the CD is shifted 1 position to the left as is required for metal ions in a PDB. Phenix choked. So I ran ReadySet on the input pdb file and gave the output file to Phenix for the next attempt. Still failed, and the error message I get is: Fatal problems interpreting PDB file: Number of atoms with unknown nonbonded energy type symbols: 2 Please edit the PDB file to resolve the problem and/or supply a Cif file with matching restraints definitions, along with apply_cif_modification and apply_cif_link parameter definitions if necessary. Also note that phenix_ready_set and phenix.elbow are available for creating restraints definitions (CIF files). Running phenix.elbow (phenix.elbow input_file.pdb --do-all) does not seem to do anything - no output file seems to be generated. Hence I am quite at loss here. Ready Set (from the gui) generated a file called name_of_the_coordinate_file.metal.edits , which is used by phenix for the next attempt but which does not seem to solve the problem. Thanks for any advice, Fred.

2 3

Re: [phenixbb] How to define link in phenix (rna3p) ?
by Ralf W. Grosse-Kunstleve 07 Apr '11

07 Apr '11

> Hi. I came across with the following post when I was dealing with a > similar issue (adding a modified nucleotide to RNA). I could create > .cif using elbow/reel as well, but I was wondering where to define > the link type between residues in phenix. You can use the apply_cif_link parameters. > For RNA, I typically see "rna3p', but with my newly created residue > linked to normal oligo, the link is formed, but recognized as Link ID > 'p'. Also, the newly created molecule is recognized as 'peptide', > although I defined _chem_comp.group in .cif file as RNA. The _chem_comp.group value is ignored by phenix.refine. The classifications are obtained via analysis of the atom names and the bonding in the monomer definitions. >------------------------------------------ >Chain: "Q1" > Number of atoms: 74 > Number of conformers: 1 > Conformer: "" > Number of residues, atoms: 3, 74 > Classifications: {'RNA': 2, 'peptide': 1} > Modifications used: {'5*END': 1, 'rna3p': 2} > Link IDs: {'p': 1, 'rna3p': 1} >------------------------------------------ It is strange that it thinks there is a peptide but still applies two rna/dna links. Could you send me (not the list) the part from the pdb file with the Q1 chain? Ralf

1 0

CUA Ligand issue
by Bradley Hintze 07 Apr '11

07 Apr '11

Hello, I am trying to refine PDB 1YEW using the default phenix settings. $phenix.refine 1yew.pdb 1yew.mtz yields: ''' Sorry: Fatal problems interpreting PDB file: Number of atoms with unknown nonbonded energy type symbols: 6 Please edit the PDB file to resolve the problems and/or supply a CIF file with matching restraint definitions, along with apply_cif_modification and apply_cif_link parameter definitions if necessary. ''' $ phenix.ready_set 1YEW_H.pdb yields: ''' Sorry: Atoms in the input are supposed to be bound CU1(CUA) - CU2(CUA) = 2.58 Fix the input file and retry ''' I assume this means that the Cu atoms in CUA are at an inappropriate distance. Is there an easy way to find out what the appropriate distance is? Can I fix this in coot? Thanks, Bradley

2 1

How to define link in phenix (rna3p) ?
by Maehigashi, Tatsuya 07 Apr '11

07 Apr '11

Hi. I came across with the following post when I was dealing with a similar issue (adding a modified nucleotide to RNA). I could create .cif using elbow/reel as well, but I was wondering where to define the link type between residues in phenix. For RNA, I typically see "rna3p', but with my newly created residue linked to normal oligo, the link is formed, but recognized as Link ID 'p'. Also, the newly created molecule is recognized as 'peptide', although I defined _chem_comp.group in .cif file as RNA. ------------------------------------------ Chain: "Q1" Number of atoms: 74 Number of conformers: 1 Conformer: "" Number of residues, atoms: 3, 74 Classifications: {'RNA': 2, 'peptide': 1} Modifications used: {'5*END': 1, 'rna3p': 2} Link IDs: {'p': 1, 'rna3p': 1} ------------------------------------------ Any input will be appreciated. Sincerely, "T" On Wed, Nov 17, 2010 at 7:43 AM, Gerwald Jogl <Gerwald_Jogl at brown.edu<http://phenix-online.org/mailman/listinfo/phenixbb>> wrote: > Hi All, > > I am refining several modified nucleotides in ribosomal RNA > (N2-methyl-guanosine and similar). While the restraints generation with > phenix.reel work quite well, I am unable to get the phosphate oxygen > restraints applied properly. Phenix.refine recognizes the links but for some > reason the phosphate groups come out quite distorted... > > For example, angles after refinement: > O3 (1) - P (2, modified) - OP2 (2) = 80 deg and > O3 (1) - P (2, modified) - OP1 (2) = 127 deg > > > Here is part of the phenix.refine log on reading the coordinates: > > Chain: "A" > Number of atoms: 32526 > Number of conformers: 1 > Conformer: "" > Number of residues, atoms: 1513, 32526 > Unusual residues: {'MA6%rna3p': 2, '5MC%rna3p': 3, '5MC%rna2p': 1, > '7M > G%rna3p': 1, '2MG%rna3p': 1, 'UR3%rna2p': 1, 'M2G%rna3p': 1, '4OC%rna3p': 1} > Classifications: {'undetermined': 11, 'RNA': 1502} > Modifications used: {'5*END': 1, 'rna3p': 1318, 'rna2p': 195} > Link IDs: {'rna3p': 1318, 'rna2p': 194} > Chain breaks: 1 > > I am at a loss about how to fix this as all links and restraints seem to be > in place. > > Many thanks for advice, > Gerwald ________________________________ This e-mail message (including any attachments) is for the sole use of the intended recipient(s) and may contain confidential and privileged information. If the reader of this message is not the intended recipient, you are hereby notified that any dissemination, distribution or copying of this message (including any attachments) is strictly prohibited. If you have received this message in error, please contact the sender by reply e-mail message and destroy all copies of the original message (including attachments).

1 0