- phenixbb - phenix-online.org

completeness problem in phenix.refine
by Yarrow Madrona 19 Nov '13

19 Nov '13

Hi, I recently merged two data sets, one having a phase and the other Amplitudes. For some reason my completeness is always 50%. I had anomalous data but the completeness is around 100% after processing and when I do an mtz dump. I'm confused why I get a completeness of 50%. Thank you for any help you can provide. -Yarrow -- Yarrow Madrona Graduate Student Molecular Biology and Biochemistry Dept. University of California, Irvine Natural Sciences I, Rm 2403 Irvine, CA 92697

2 2

Error: "Error running main function (INPUT): LLG completion peak selection method not recognized"
by Egor Tchesnokov 18 Nov '13

18 Nov '13

Hello, I am trying to run Phaser/SAD or Phaser/MR_SAD (GUI using defaults), but I get the following error: "Error running main function (INPUT): LLG completion peak selection method not recognized". Thank you for your help Egor

2 1

MSE out density after phenix.refine
by Pavel Afonine 16 Nov '13

16 Nov '13

Hello, anyone experiencing problems related to MSE residues being pushed out density after phenix.refine run should update to next Phenix nightly build, which should fix this problem. Please let me know otherwise off-list. Pavel

1 0

MR_Rosetta prerefine question
by Heather L Condurso 12 Nov '13

12 Nov '13

Dear all, I am trying to use the rosetta prerefinement fuction of MR_rosetta to complete a homology model with missing stretches of amino acids. Everything runs fine and residues are being added to fill the gaps using the fragment files, but rosetta is also altering parts of the model which should already be well placed. Is there a way to tell Rosetta only to prerefine certain residues? I've been following the command line suggestion: phenix.mr_rosetta \ seq_file=seq.dat \ data=coords1.mtz \ search_models=coords1.pdb \ run_prerefine=True \ number_of_prerefine_models=1000 \ fragment_files = test3.gz \ fragment_files = test9.gz \ rescore_mr.relax=False \ rosetta_models=100 \ ncs_copies=2 \ space_group=p212121 \ use_all_plausible_sg=False \ nproc=200 \ group_run_command=qsub Thanks, Heather Condurso condurso(a)bc.edu

4 5

Re: [phenixbb] Phenix.Rosetta_refine on Ubuntu 12.04 LTS
by qwertgfdsa78 11 Nov '13

11 Nov '13

Hi All, I tried Jan's tutorial and finally I was able to install and ran phenix.rosetta with test set 1yjp successfully. However, when I tried my data, I always got *"terminate called after throwing an instance of 'std::bad_alloc'"* Google suggests maybe memory is not enough. I used VirutalBox and assign 8gb memory to it. My protein has ~3000 residues. I guess this will use more memory but how much RAM I need? Or the error is due to some other problem? Thank for any inputs. Lee > ------------------------------ > > Message: 3 > Date: Tue, 05 Nov 2013 14:26:00 +0100 > From: Jan Gebauer <jan.gebauer(a)uni-koeln.de> > To: PHENIX user mailing list <phenixbb(a)phenix-online.org> > Subject: Re: [phenixbb] Phenix.Rosetta_refine on Ubuntu 12.04 LTS > Message-ID: <5278F1E8.9020602(a)uni-koeln.de> > Content-Type: text/plain; charset="utf-8" > > Hi David, hi all, > > I now managed to at least get a test refinement working. What I mostly > changed was the Ubuntu version. > I have no idea why, but I managed to get it to work on Ubuntu 10.04. > > In principle I installed 10.04 freshly in a virtualbox and from their it > worked rather straight-forward. > > Here is how I proceeded in detail for those of you who are interested: > > 1) Setting up a Virtualbox 4.3.2 r90405 on Ubuntu 12.04.3. > 2) Installing the 64bit version of Ubuntu 10.04 in a virtualbox (4 > processors / 90% speed cap / 4GB Ram /VT-x enabled) > 3) Install VBoxAdditions > 4) Update ubuntu (sudo apt-get update) > 5) Installing g++ and zlib library (sudo apt-get install g++ zlib1g-dev) > 6) Installing binary version of phenix-dev-1525 to > /usr/local/science/phenix-dev-1525 (sudo ./install > --prefix=/usr/local/science) > 7) Unpack rosetta_2013_wk42_bundle to > /usr/local/science/rosetta_2013_wk42_bundle > 8) set rossetta variables by inserting them at the start of > /usr/local/science/pheninx-dev-1525/phenix_env.sh (just in line 3ff): > export PHENIX_ROSETTA_PATH=/usr/local/science/rosetta_2013wk42_bundle > export ROSETTA3_DB=/usr/local/science/rosetta_2013wk42_bundle/main/database > export > PATH=$PATH:/usr/local/science/rosetta_2013wk42_bundle/main/source/bin > 9) assign ownership to active user (sudo chown gebauer:gebauer > /usr/local/science/ -R) > 10) Compile it (phenix.build_phenix_interface nproc=4) > It should run without problems > > Testing: > 1) Create a new directory (mkdir test && cd test) > 2) Download testset (suggested by Nat: phenix.fetch_pdb --mtz 1yjp) > 3) Run rosetta_refine (phenix.rosetta_refine 1yjp.pdb 1yjp.mtz nproc=4) > In my case rFree dropped from 18.21 to 14.90. It's not a low resolution > dataset, but it show that it general works and lowRes would potentially > take to long for a mere testing > > My real data are still being processed so I have no idea how well it > will works... only thing I found out already: > 4.5 GB RAM for 4 processes calculating ~500Aa via rosetta seems not to > be enough. I thereby reduced my "real run" to only one processor... > > Regards, > Jan > > >

3 3

structure factor labels in phenix.refine
by Tanner, John J. 08 Nov '13

08 Nov '13

I have two questions about which structure factor amplitudes are refined against in phenix.refine. 1. Which structure factor amplitude is used in refinement when the labels in the .def file are the following? labels = "F-obs(+),SIGF-obs(+),F-obs(-),SIGF-obs(-)" 2. Is anomalous data used in the refinement when the labels are as follows? labels = "F,SIGF,DANO,SIGDANO,ISYM" I ask these questions because I am wondering which structure factor amplitudes to deposit, and I am trying to rationalize differences between R-factors from refinement and Generate Table 1. Thanks. Jack Tanner John J. Tanner Professor of Biochemistry and Chemistry University of Missouri-Columbia 125 Chemistry Building Columbia, MO 65211 Phone: 573-884-1280 Fax: 573-882-2754 Email: tannerjj(a)missouri.edu<mailto:[email protected]> http://faculty.missouri.edu/~tannerjj/tannergroup/tanner.html

4 4

few questions about resolving new structure through MR
by Zhihong Yu 07 Nov '13

07 Nov '13

Hi, all I'm a rookie in resolving a brand new structure. I have some questions for my current case and look forward to some suggestions. Now I’m working on a protein like this: N-ter(55aa)—domainA(110aa)—linker(30aa)—domainB(150aa)—C-ter(20aa), I got a diffraction data just to 3.5Å, and there is no complete homology structure in pdb bank, but only a homology structure (named as structureX later) for domainB with ~30% sequence identity, so I have some questions as following: 1. Is it possible to find a resolution through MR approach using structureX as a search model? Especially considering that the resolution is only 3.5Å. Currently I just tried once using phaser and refine the structure for 3 cycyles, I can get a R/Rfree of 0.45/0.55, and it looks like most of backbone in the structureX, especially those within helix or sheet, can be well described by 2Fo-Fc density. Is this primary result promising or not? 2. If it’s possible, what’s the general optimal procedure I should follow? Really thanks for any advice and suggestions! Zhihong

2 3

Issue with SeMet residues
by Mark Wilson 07 Nov '13

07 Nov '13

Dear PHENIX users and developers, I am refining a structure containing SeMet residues at atomic (1.05 Å) resolution. The model was initially refined in Refmac5 and behaved as expected (i.e. reasonable R factors and agreement with maps). I've used this final Refmac5 model as a starting point for refinement in PHENIX 1.8.4 with standard settings, as below. phenix.refine data.mtz \ refinement.input.xray_data.labels="F, SIGF" model.pdb \ main.use_experimental_phases=false \ refinement.main.number_of_macro_cycles=5 output.prefix="model_refine"> phenix_refine_1.log Strangely, the conformations of the MSE residues are changed quite significantly after refinement in PHENIX, departing both from electron density and from the prior model. Any help would be appreciated. Best regards, Mark Mark A. Wilson Associate Professor Department of Biochemistry/Redox Biology Center University of Nebraska N118 Beadle Center 1901 Vine Street Lincoln, NE 68588 (402) 472-3626 mwilson13(a)unl.edu

3 3

Issues using reference model with Autobuild
by Dan McNamara 07 Nov '13

07 Nov '13

Hello all, I'm hoping to troubleshoot a problem with autobuild and the use of a reference model. Autobuild has been always failing once it gets to the initial refinement step and attempts to add reference model restraints. I've included my autobuild parameterization. Refinement is where it dies and the contents of LAST.LOG and my error log file are below too. I'm not sure how to get the reference model file location properly interpreted. I even tried manually moving the reference PDB into the AutoBuild_run_X_/ and TEMP0/ folders as they're generated to no avail. Best, Dan phenix.autobuild \ data=/home/.../xyz.mtz \ refinement_file=/home/.../xyz.mtz \ input_refinement_labels="FP SIGFP FreeRflag" \ model=/home/.../xyz.pdb \ seq_file=/home/.../xyz.dat \ s_annealing=False \ rebuild_in_place=True \ n_cycle_rebuild_max=3 \ refine_eff_file=/home/.../xyz.eff \ refine_with_ncs=True \ model_building.reference_model=/home/.../def.pdb \ #I have also tried locations like ../def.pdb and def.pdb that give Errno 21 or Errno 2 below nproc=200 \ general.run_command="/home/.../xyz.sh" AutoBuild_run_1_/LAST.LOG: ================== Extract refinement strategy and selections ================= Refinement flags and selection counts: individual_sites = True (30132 atoms) torsion_angles = False (0 atoms) rigid_body = True (30120 atoms in 24 groups) individual_adp = False (iso = 0 aniso = 0) group_adp = True (30132 atoms in 4092 groups) tls = True (30120 atoms in 147 groups) occupancies = True (12 atoms) group_anomalous = False n_use = 30132 n_use_u_both = 30120 n_use_u_iso = 12 n_use_u_aniso = 0 n_grad_site = 0 n_grad_u_iso = 0 n_grad_u_aniso = 0 n_grad_occupancy = 0 n_grad_fp = 0 n_grad_fdp = 0 total number of scatterers = 30132 *** Adding Reference Model Restraints *** errors.txt: Sorry: ******************************************************************************** Failed to carry out AutoBuild_build_cycle: Initial phenix refinement of AutoBuild_run_1_/TEMP0/NO_LIG_NO_WAT_xyz_ed.pdb with aniso_refinement_data_PHX.mtz failed...quitting... /home/.../AutoBuild_run_1_/LAST.LOG may have additional information... Cannot open file for reading: "/" [Errno 21] Is a directory: '/' ******************************************************************************** Sorry: ******************************************************************************** Failed to carry out AutoBuild_build_cycle: Initial phenix refinement of AutoBuild_run_2_/TEMP0/NO_LIG_NO_WAT_xyz_ed.pdb with aniso_refinement_data_PHX.mtz failed...quitting... /home/.../AutoBuild_run_2_/LAST.LOG may have additional information... Cannot open file for reading: "." [Errno 21] Is a directory: '.' ******************************************************************************** Sorry: ******************************************************************************** Failed to carry out AutoBuild_build_cycle: Initial phenix refinement of AutoBuild_run_3_/TEMP0/NO_LIG_NO_WAT_xyz_ed.pdb with aniso_refinement_data_PHX.mtz failed...quitting... /home/.../AutoBuild_run_3_/LAST.LOG may have additional information... Cannot open file for reading: "d" [Errno 2] No such file or directory: 'd' ********************************************************************************

2 2

Sigmas of F+ and F- equal in _data file
by Tanner, John J. 07 Nov '13

07 Nov '13

I have a question about the _data.mtz file that phenix.refine outputs (version 1.8-1069). It appears that sigmaF+ equals sigmaF- for each reflection in the _data.mtz file, although they were unequal in the input .mtz file. You can compare the input mtz and the output _data.mtz file by viewing these links: http://faculty.missouri.edu/~tannerjj/input.jpeg http://faculty.missouri.edu/~tannerjj/samesigma.jpeg Thanks. Jack Tanner John J. Tanner Professor of Biochemistry and Chemistry University of Missouri-Columbia 125 Chemistry Building Columbia, MO 65211 Phone: 573-884-1280 Fax: 573-882-2754 Email: tannerjj(a)missouri.edu<mailto:[email protected]> http://faculty.missouri.edu/~tannerjj/tannergroup/tanner.html

2 1