- phenixbb - phenix-online.org

Reproducing a run in phenix.hyss
by Jeffrey, Philip D. 19 Aug '12

19 Aug '12

I've got an interesting potential S-SAD solution that crops up rarely in a script I was using to trawl parameter space in phenix.hyss. Since I was not optimistic of the outcome I did not save all output files. Hyss appears to use a random number generator and allows the seed to be specified on the command line, so is there any way to figure out what random number seed it used based on the information written to stdout while the program is running ? Merely running the command again does not trivially recreate the solution and egrep -i 'rand|seed' LOG only pulls out: Nigel Moriarty, Reetal Pai, Randy Read, Jane Richardson, David Richardson, time_log: p1_recycling random_omit: 990 1.49 0 0.00151 time_log: random_omit: 891 0.62 0 0.000696 Cheers Phil Jeffrey Princeton

2 1

phenix.refine crash
by Maxime Cuypers 18 Aug '12

18 Aug '12

> > Hi, > > i am refining a neutron structure and i keep having this error: > > > File > "/usr/local/phenix-dev-1129/build/intel-linux-2.6/../../phenix/phenix/command_line/refine.py", > line 11, in <module> > command_line.run(command_name="phenix.refine", args=sys.argv[1:]) > File > "/usr/local/phenix-dev-1129/phenix/phenix/refinement/command_line.py", line > 146, in run > call_back_handler=call_back_handler) > File "/usr/local/phenix-dev-1129/phenix/phenix/refinement/driver.py", > line 1533, in run > ncs_manager = self.ncs_manager) > File > "/usr/local/phenix-dev-1129/phenix/phenix/refinement/strategies.py", line > 1109, in refinement_machine > h_params = h_params) > File > "/usr/local/phenix-dev-1129/cctbx_project/mmtbx/refinement/adp_refinement.py", > line 212, in __init__ > nproc = nproc) > File > "/usr/local/phenix-dev-1129/cctbx_project/mmtbx/refinement/adp_refinement.py", > line 315, in __init__ > result = self.try_weight(weight, print_stats=True) > File > "/usr/local/phenix-dev-1129/cctbx_project/mmtbx/refinement/adp_refinement.py", > line 401, in try_weight > minimized = self.minimize() > File > "/usr/local/phenix-dev-1129/cctbx_project/mmtbx/refinement/adp_refinement.py", > line 470, in minimize > self.model.set_refine_individual_adp() > File "/usr/local/phenix-dev-1129/cctbx_project/mmtbx/model.py", line > 1231, in set_refine_individual_adp > iselection = selection_iso.iselection()) > RuntimeError: cctbx Internal Error: > /net/longnose/scratch2/phenix/phenix-dev-1129/cctbx_project/cctbx/xray/scatterer_flags.h(417): > CCTBX_ASSERT(f.use_u_iso()) failure. > > i see the last line is already a known issue as it was already posted in > phenixbb in 2010 but was perhaps solved offline. > > thanks in advance for any help (as i am really stuck) > > -- > Dr. Maxime Cuypers > Institut Laue Langevin (ILL), Deuteration lab > Rue Jules Horowitz, 6, BP 156 > 38042 Grenoble cedex 9 > France > > >

7 8

Real-space and twin refinement
by Haim Rozenberg 13 Aug '12

13 Aug '12

Hello, We've just downloaded and install the new dev-1129 release of phenix hoping to be able to run RS refinement on our twin structures. According to the change log it should be possible with this version. ------------------------------------------------------------------------ Version 1.8.1 and nightly builds (TBA) ====================================== - phenix.refine: - overall B-cart now applied to individual atomic B-factors (restores pre-1.8 behavior) - improvements to torsion NCS rotamer correction - real-space refinement for twinned structures ------------------------------------------------------------------------ Unfortunately, while trying to run - through the GUI - phenix dev-1129 with our twin setting and the "Real space" option, the usual message popped up: "PHENIX error: Real-space refinement is not available for twinned structures." FYI, it runs just fine without the "Real-space" option. Did we have false hopes? Did we misunderstood? Is that a bug? All the best, Haim ============================================ Haim Rozenberg, PhD. Department of Structural Biology Kimmelman Building Weizmann Institute of Science 76100 Rehovot, Israel Phone: +972-8-9344289 FAX:+972-8-9344136 E-Mail: haim.rozenberg(a)weizmann.ac.il<mailto:[email protected]> ============================================ "Quis custodiet ipsos custodes?" D.I. Iuvenalis

2 1

difficult phasing problem
by Ian Slaymaker 09 Aug '12

09 Aug '12

I am trying to solve the structure of a ~30kd monomeric protein for which we have 2.5A diffraction data, and a search model with ~65 % similarity (i have tried, homology models of this and others, as well as using sculpter to build several variations of core search models). I have tried all my possible space groups and a huge range of parameters, and mr_rosetta to rebuild some of my better (but still bad) MR solutions. My LLG is persistently negative (<-1000), which may be due to translational noncrystallographic symmetry which was identified by xtriage. We have unfortunatly been unable to phase this data with selmet SAD, so I am looking for some other MR things to try. Any suggestions for working with tNCS (if this is really the problem) would be much appreciated.

2 1

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