- phenixbb - phenix-online.org

occupancy refinement
by Nicole Bischoff 30 Jan '09

30 Jan '09

Dear all, in my electron density I can see negative difference density (fo-fc) for the carboxyl groups of some aspartates or glutamates. Assuming that this is a result of partial decarboxylation due to radiation damage, I wanted to refine the occupancies of these carboxyl groups. To this end, I set the occupancy of these atoms to 0.99 in Coot so that they would be picked up by occupancy refinement in phenix.refine (using version 1.4-3). Phenix however just put the value back to one. It looks like the occupancy refinement in phenix is meant for alternative conformations where the sum of the occupancies of the different conformations is 1. How can I make phenix refine partially decarboxylated acidic side chains? It does not seem appropriate to me to make alternative conformations, one with and one without the carboxyl group. That would look funny at my resolution (2.7 Angstrom) and might not fix the problem. Any hints are appreciated. Thank you, Nicole ---- Nicole Bischoff, M. Sc. PhD student in the Gruetter lab Department of Biochemistry University of Zurich Winterthurerstrasse 190 8057 Zurich Switzerland +41 44 635 5556 nbischoff(a)bioc.unizh.ch

3 3

TLS question
by Jeffrey Chao 27 Jan '09

27 Jan '09

hello all- i apologize if this question has previously been posted, but i was not able to track down an answer by google. i am refining a 2.75A structure that has three copies of the protein in the asu with phenix 1.4-3 and using NCS with each chain as an NCS group. after molecular replacement, i set all of the b-factors to the wilson b-factor and then performed four different refinements: grouped b-factor (1 group per residue) : Rwork 22.6 Rfree 27.0 individual b-factor : Rwork 21.4 Rfree 26.6 grouped b + tls ( : Rwork 21.9 Rfree 25.2 individual + tls : Rwork 21.2 Rfree 24.8 i defined the tls groups with the following: adp { tls = chain A tls = chain B tls = chain C } after grouped b-factor the atoms of each residue all have the same b- factor as i expected. after grouped b-factor+tls, however, each atom of a residue has a unique b-factor, which i did not expect. i realize that this confusion is likely caused by me not fully understanding how tls refinement is implemented. if someone could explain this to me, i will be able to sleep better. thank you for your time, jeff Jeffrey Chao, PhD 1300 Morris Park Ave ASB-Golding 601 Bronx NY 10461 jchao(a)aecom.yu.edu 718-430-8597

2 1

rmsd: phenix vs refmac
by Frank von Delft 26 Jan '09

26 Jan '09

Hi I thought that a while ago, I read somewhere that rmsd(bonds) had not been calculated the same way in phenix as in refmac. I believe it was in connection with automated weights, and that rmsd(bonds) always much lower in phenix. I understood too that this had been adjusted since to match refmac. But we're still seeing some much lower rmsd(bond) from phenix: in particular, when we take a refmac model and plug it straight into phenix, the rmsd for cycle zero is already much lower than what it was at the end of refmac. Could someone enlighten...? (Yup, latest version of phenix.) Thanks! phx

5 5

bond and angle rmsds: phenix.refine, Refmac, hydrogens, ...
by Pavel Afonine 25 Jan '09

25 Jan '09

For all who asked about observed differences in bond rmsd reported by phenix.refine and Refmac, H atoms included in statistics calculation and more... 1) First, I wanted to check bond rmsd reported by phenix.refine. I made up a PDB file containing just one residue (SER, no H atoms) and I wrote a crystal-clear and simple Python script that computes bond rmsd for this residue. Here are the numbers: Bond length rmsd: Script: 0.148 phenix.pdbtools: 0.148 phenix.refine: 0.148 Refmac_5.5.0070: 0.152 Both programs are using the Monomer Library. Since both of these very different ways of computing bond rmsd: this simple script and phenix.refine, give the identical values, I tend to think that the PHENIX implementation is actually correct. I have no idea why Refmac gives a different value and I think it's now up to Refmac developers (or any other enthusiast) to find it out or show that I'm wrong in my calculations. It well could also be me who's incorrectly using Refmac - today was my second time I tried it (the first was 3 or 4 years ago) -:) As for the bond angles rmsd, I didn't bother to include it into the script, but for the rest of tools in question, the numbers are: phenix.pdbtools: 16.628 phenix.refine: 16.628 Refmac_5.5.0070: 16.628 so all results seem to be in-sync. 2) Some people reported that the geometry statistics for the exact same structure with and without hydrogens seems to be reported different in phenix.refine. To find out, I did the same test as above, where I used the same one residue model with all hydrogens added. Unfortunately, I failed to make Refmac accept a PDB file containing H atoms so I gave up after 20 minutes of struggle. Here are the numbers: Bond length rmsd: Script: 0.155 phenix.pdbtools: 0.155 phenix.refine: 0.155 So, it seems like the hydrogens are included... I have to fix it. Well at least those who are doing refinement using neutron data or both X-ray and neutron data, were not concerned, since in that case I presume the H and D atoms should be included. Anyway, I'm fixing it right now, I'm making a corresponding regression test (to make sure the fix stays there forever) and the fixed version will be available in the next release. By the way, the new tool, phenix.model_vs_data is and was smart about this: it automatically excludes H and D atoms if X-ray data is used. And this proves it: phenix.model_vs_data: 0.148 which is the exact same number as in "1)", where H were not present. Thanks to all who raised this issue (Frank von Delft, Jianghai Zhu, Karolina Michalska and others)! Pavel. PS> All files used and the script are here: http://cci.lbl.gov/~afonine/rmsd_issue/ Use "phenix.python run.py" to run the script.

1 0

IAS causes phenix.refine to die
by Anna Gardberg 25 Jan '09

25 Jan '09

Dear phenix list, My 0.75 A refinement has some problems in its higher resolution bins, so I thought I might try adding in IAS scatterers. However, the phenix.refine (version 2008_09_13_0905) crashed when I set "ias = True". I hope someone here can help me figure out the problem. I include snippets of output below, as the entire error message was very lengthy. Am I misusing the command, is there a problem with my system setup, or what? Regards, Anna Gardberg ----------Building IAS model---------- >>> Adding IAS.......... Total number of covalent bonds = 895 show_stack(1): /home/a1g/downloads/cci_apps_sources/mmtbx/ias/__init__.py(438) _ _init__ show_stack(2): /home/a1g/downloads/cci_apps_sources/mmtbx/ias/__init__.py(371) s et_peaks show_stack(3): /home/a1g/downloads/cci_apps_sources/mmtbx/ias/__init__.py(616) _ _init__ show_stack(4): /home/a1g/downloads/cci_apps_sources/mmtbx/model.py(494) add_ias show_stack(5): /home/a1g/downloads/cci_apps_sources/phenix/phenix/refinement/str ategies.py(409) refinement_machine show_stack(6): /home/a1g/downloads/cci_apps_sources/phenix/phenix/refinement/dri ver.py(1082) run <snip pages of output> /home/a1g/downloads/cci_apps_build/lib/cctbx_maptbx_ext.so(boost::python::obje cts::make_holder<4>::apply<boost::python::objects::value_holder<cctbx::maptbx::o ne_gaussian_peak_approximation>, boost::mpl::vector4<scitbx::af::const_ref<doubl e, scitbx::af::trivial_accessor> const&, scitbx::af::const_ref<double, scitbx::a f::trivial_accessor> const&, bool const&, bool const&> >::execute(_object*, scit bx::af::const_ref<double, scitbx::af::trivial_accessor> const&, scitbx::af::cons t_ref<double, scitbx::af::trivial_accessor> const&, bool const&, bool const&)+0x 75c) [0x2a9d79598c] /home/a1g/downloads/cci_apps_build/lib/cctbx_maptbx_ext.so(cctbx::maptbx::find _gaussian_parameters::find_gaussian_parameters(scitbx::af::const_ref<double, sci tbx::af::trivial_accessor> const&, scitbx::af::const_ref<double, scitbx::af::tri vial_accessor> const&, double const&, double const&, double)+0x3fe) [0x2a9d7827f e] /lib64/tls/libm.so.6 [0x32ee506263] /lib64/tls/libc.so.6 [0x32ee22e300] Floating-point error (Python and libc call stacks above) [1]+ Done gedit newscale-ias.eff

3 2

Unmasking cavities
by Frank von Delft 24 Jan '09

24 Jan '09

Hi After being through phenix.refine, I see in my hydrophobic core a big space (a few atoms wide) that is filled with strong negative difference density. I suspect the culprit is the bulk solvent mask, which is defined too tightly. The online manual mentions three parameters, but not what they do. solvent_radius, shrink_truncation_radius, grid_step_factor What *exactly* do they do? (I thought I'd elicit a contribution for the online docs this way :) Cheers phx

6 11

geometry of disulfide and NAG
by Li-Zhi Mi 23 Jan '09

23 Jan '09

Hi, I am using phenix to refine a low resolution structure. The problem I had is that the geometry of disulfide and NAG in the refined structure is not right. I did use a very tight geometry restrains. Both wxc_scale and wxu_scale are set to 0.05. In the refined structures, the sulfur of one cys is very close to the Cb of its disulfide partner. For NAG part, C2-N2-C7-O7-C8 should be in the same plane. But they are not. Is there a way to adjust these violations? Thanks in advance for the helps, Li-Zhi Mi,

7 9

HL A B C D
by Raja Dey 22 Jan '09

22 Jan '09

Hello, If I want to use F, SigF and HLA, B, C and D in the phenix refinement, should I use the same command like phenix.refine data.mtz model.pdb, where data.mtz contains FREER, F, SIGF, HLA, HLB, HLC and HLD? How do I check that phenix.refine used the HL coefficients? Looking forward to hearing from you... Thanks... Raja

2 1

Disappearing side chains after autobuilding
by Katya Heldwein 22 Jan '09

22 Jan '09

Hi, I ran phenix.autobuild with default setting (rebuild_in_place), and noticed that some side chains were gone in the resulting overall_best.pdb. That is, let's say, I started with Phe as residue 10. After rebuilding, the Phe10 is still called Phe in the output file, but the side chain atoms are gone, leaving residue 10 as Ala. This happened for several side chains, and I am at a loss for why this happened. Is this some kind of a bug? Katya -- Ekaterina Heldwein, Ph.D. Assistant Professor Department of Molecular Biology and Microbiology Tufts University School of Medicine 136 Harrison Ave Boston, MA 02111 phone: 617-636-0858 fax: 617-636-0337 e-mail: katya.heldwein(a)tufts.edu

2 1

missing map coefficients
by Anne Mølgaard 16 Jan '09

16 Jan '09

I am new to Phenix, so there might be an easy explanation for this problem, though I have not been able to find a solution in the documentation. I am using phenix version 1.3b-rc4, intel-linux-2.6. I have solved a structure using phenix.automr, followed by phenix.autobuild (run from command line using "phenix.autobuild data=fromXDS_unique.mtz model=AutoMR_run_2_/MR.1.pdb seq_file=sequence.fas"). The program runs succesfully, and I get the output files I should expect. However, the overall_best_refine_map_coeffs.mtz does not contain the map coefficients 2FOFCWT PH2FOFCWT for the 2mFo-DFc map and FOFC and PHFOFC for the Fo-Fc difference map. mtzdmp gives me the following columns: 1 ASC 0 28 0 100.00 10.6 10.6 28.28 1.85 H H 2 NONE 0 35 0 100.00 13.3 13.3 28.28 1.85 H K 3 NONE 0 71 0 100.00 26.6 26.6 28.28 1.85 H L 4 NONE 1.0 1888.8 0 100.00 148.71 148.71 28.28 1.85 F FP 5 BOTH 1.0 1.0 0 100.00 1.00 1.00 28.28 1.85 Q SIGFP 6 NONE -180.0 180.0 0 100.00 2.08 90.23 28.28 1.85 P PHIM 7 BOTH 1.000 1.000 0 100.00 1.000 1.000 28.28 1.85 W FOMM 8 NONE .0 19.0 0 100.00 9.50 9.50 28.28 1.85 I FreeR_flag 9 BOTH .0 .0 0 100.00 .00 .00 28.28 1.85 F FC What should I do to get these map coefficients? Thanks for your help! Best regards, Anne

2 2