- phenixbb - phenix-online.org

test set count question
by Filip YABUKARSKI 08 Aug '12

08 Aug '12

Dear all, I am refining a structure with phenix.refine using a .mtz which contains F/SigF and DF/SigDF (CCP format from xdsconv) and i have something like 10000 unique reflections. The test set was created with the uniqueify script. If I understand correctly in the manual it is written that phenix.refine will separate the <F>/<SigF> into "F(+) SIGF(+) F(-) SIGF(-)), instead of keeping the original non-anomalous amplitudes". My question is about the "FREE R VALUE TEST SET COUNT" in the .pdb after phenix.refine. If I understand correctly the number for the test set will also be "doubled" because the <F> will be separated into F+ F-? Logically both F+/F- after separation are kept in the test set is this right? And an impair number for the test set in this case means that only the F+ OR F- was observed? Sorry if I ask a trivial question but I just want to be sure I am using the program in the right way Thanks in advance Filip

3 4

Re: [phenixbb] test set count question
by Kay Diederichs 08 Aug '12

08 Aug '12

Dear Filip, concerning XDSCONV: if FRIEDEL'S_LAW=FALSE, using filetype of CCP4, you get F,SigF,DF,SigDF,isym using filetype CCP4_F, you get F,SigF,F(+),SigF(+),F(-),SigF(-) using filetype CCP4_I, you get IMEAN,SIGIMEAN,I(+),SIGI(+),I(-),SIGI(-) so you can (and should) directly output from xdsconv what you need for a downstream program. A short explanation is at http://strucbio.biologie.uni-konstanz.de/xdswiki/index.php/Xdsconv Detailed explanation of each format is at http://homes.mpimf-heidelberg.mpg.de/~kabsch/xds/html_doc/xdsconv_formats.h… HTH, Kay

1 0

changing a project directory
by Kendall Nettles 08 Aug '12

08 Aug '12

Dear PHENIX developers, I want to migrate my PHENIX project directories. I moved one project directory as a test, and the GUI doesn't give an opportunity to relocate, which would be a nice feature. Is there a way to edit their location for the PHENIX GUI? Kendall

2 1

two problems about phenix.mr_rosetta
by chen c 08 Aug '12

08 Aug '12

Dear all: I've collected a data set for my crystal, which is a two-protein complex. However, normal molecular replacement method gives no solution and thus didn't work. Right now I'm struglling with the phenix.mr_rosetta method. However, I've got several problems when reviewing it, and thus looking forward to get some advice here. 1: There are such words *"NOTE 2: If your structure contains more than one chain or requires more than one homology model to represent the structure, then you need to use mr_model_preparation and phenix.automr to place your model. Then you can start phenix.mr_rosetta with already_placed=True**"* in phenix documentation page. My understanding for this sentences is: the solution derived from *automr* is treated as an original searching model for phenix.mr_rosetta, in which case, we should add the command line *"run_prerefine=True \ number_of_prerefine_models=1000"*. However, what is *mr_model_preparation * for? And if no solution is given by *automr*, what should I do? What's more important, why not endow the phenix.mr_rosetta program to be capable of dealing with data set from multi-chain protein complexes? Since I thought the idea of preliminary manipulation by the program *automr* would decrease the power of *phenix.mr_rosetta*. 2: Two .hhr files are produced. But as displayed in the phenix.mr_rosetta documentation page, only one .hhr file like *"hhr_files=bfr258e.hhr"* are listed in the command line. Why? Could somebody give me any advice on this, please? Thank you and best regards chen -- Cheng Chen, Ph.D. Candidate Laboratory of Structural Biology Life Science Building,Tsinghua University Beijing 100084 China -

1 0

refining multiple conformations
by Kendall Nettles 08 Aug '12

08 Aug '12

Hi, We have a structure with the ligand showing two overlapping conformers. When we refine it with both conformers separately, it is pretty clear that there are substantial differences in the protein as a result, for about a third of the protein chain. My question is, would it be better to try to define altloc for those specific regions, or would it be OK to refine with two entire alternate protein chains? There is also a second protein chain that shows only a single binding mode for the ligand. It's a 2.0 angstrom structure with 4000 atoms and 27,000 reflections. The yellow map goes with the green model in the attached pic. Also, do we want to let each amino acid have its own occupancy? or should one ligand copy and one chain all have the same occupancy? I'm leaning towards the latter since the differences should be directly tied to the ligand binding mode. Kendall Nettles [cid:[email protected]]

2 1

mr_rosetta problems
by Ian Slaymaker 07 Aug '12

07 Aug '12

I am trying to get mr_rosetta working on our cluster and am running into some problems I am not sure how to solve. The phenix_regression.wizards.test_command_line_rosetta_quick_tests tests went smoothly were completed. However, my first test run with all the defaults (not submitted to multiple processors) using my reflection file, HHR analysis file, sequence file, and the 2 fragment files failed with the below errors. I checked all the files, and it seems this is a problem with the fragment files but I am not sure exactly what they are supposed to look like. I thought it could be because our cluster doesnt have access to any outside networks...is this a problem with rosetta not being able to download pdb files? if so, is there a way to supply them internally? The aligned model looks somewhat what I expected it to (same beta sheet arrangement) but it seems to be incomplete. Any help would be very much appreciated! Splitting work into 100 jobs and running with 1 processors using sh background=None in /MR_ROSETTA_1/GROUP_OF_ROSETTA_REBUILD_1/RUN_1/REBUILD_IN_SETS_1 Starting job 1...Log will be: /MR_ROSETTA_1/GROUP_OF_ROSETTA_REBUILD_1/RUN_1/REBUILD_IN_SETS_1/RUN_FILE_1.log ******************* ERROR ENDING *************** mr_rosetta failed... ******************************************************************************** Sorry, a subprocess has failed... END OF LOG FILE /MR_ROSETTA_1/GROUP_OF_ROSETTA_REBUILD_1/RUN_1/REBUILD_IN_SETS_1/RUN_FILE_1.log : to match the actual PDB sequence: (protein sequence) New alignment line 4 for PDB: (protein sequence) Start position of target sequence in sequence file: 0 Aligned model for comparison: /MR_ROSETTA_1/GROUP_OF_ROSETTA_REBUILD_1/RUN_1/REBUILD_IN_SETS_1/RUN_1/WORK_1/aligned_model.pdb Log for rebuilding: /MR_ROSETTA_1/GROUP_OF_ROSETTA_REBUILD_1/RUN_1/REBUILD_IN_SETS_1/RUN_1/WORK_1/rebuild.log **************ERROR ENDING********************** Failed to run REBUILD...see log file in /MR_ROSETTA_1/GROUP_OF_ROSETTA_REBUILD_1/RUN_1/REBUILD_IN_SETS_1/RUN_1/WORK_1/rebuild.log The rebuild log had repeats of several errors that I am not sure how to address: core.conformation.Conformation: [ WARNING ] missing heavyatom: CD2 on residue PHE 3 (repeated several times) core.fragments.ConstantLengthFragSet: line too short: Skipping line ' (repeated several times) core.fragments.ConstantLengthFragSet: finished reading top 0 0mer fragments from file/dataFiles/aat000_03_05.200_v1_3 protocols.loops.loops_main: Frag libraries debug none 1

4 5

Ordered solvent distances
by Abhinav Kumar 03 Aug '12

03 Aug '12

Hi, The default minimum distance that a water gets placed near protein/water atoms using 'ordered_solvent=true' option in phenix.refine is 1.8. Isn't this too close? Also, if two waters can't get placed in the density due to distance criterion, can phenix place one water in dual conformations in that density? -- Thanks, Abhinav JCSG@SSRL, SLAC (650) 926-2992

3 4

Lost my job history, can it be recovered? (GUI)
by Andrew Waight 02 Aug '12

02 Aug '12

Hello All, It just so happened that my lab's server ran out of space (someone else uploading datasets) while my phenix rosetta MR job was running. Now the GUI no longer shows any job history in the job history tab (although on the title page the number of jobs is correctly displayed. I cannot find any information on how or if its possible to recover my job history in the GUI. Any tips? Thanks so much. Drew

2 1

July 2012 Computational Crystallography Newsletter
by Nigel Moriarty 31 Jul '12

31 Jul '12

Folks The latest issue of the Computational Crystallography Newsletter is available on-line at http://www.phenix-online.org/newsletter for download. The articles are of general interest to protein crystallographers and are listed below. Enjoy Nigel Moriarty Articles ----------- cctbx tools for transparent job execution on clusters On the analysis of residual density distributions on an absolute scale Short communications -------------------------------- DETAC: tools to detect alternative conformations by unrestrained refinement ERRASER, a powerful new system for correcting RNA models Fitting tips - Rotamer correction with backrub -- Nigel W. Moriarty, PhD Building 64R0246B, Physical Biosciences Division Lawrence Berkeley National Laboratory Berkeley, CA 94720-8235 Phone : 510-486-5709 Email : NWMoriarty(a)LBL.gov Fax : 510-486-5909 Web : CCI.LBL.gov

1 0

what installer would work with a cluster with Darwin 9.8.0 if any?
by Tatyana Sysoeva 27 Jul '12

27 Jul '12

Hi all, I am starting to use a new computer cluster and just found out by uname that it runs Darwin OS? It is totally new to me and I am wondering if Phenix would run on it? should I use just some Mac installer? I would really appreciate any advices! Thank you! Sincerely, Tanya

3 3