- phenixbb - phenix-online.org

Rwork Rfree Bonds Angles
by Hexi Zhang 01 Nov '12

01 Nov '12

Dear All, During Phenix refinement, it request the Rwork and Rfree graduately decrease (at least no increase). But there is no request that the Bonds (should be RMS?) and the Angles (should also be RMS?) should graduately decrease. Am I right? I am looking forward to getting your reply. Hexi

2 1

"H" atom refining occupancy with previous altconf residue in sequence
by Alexander Scouras 01 Nov '12

01 Nov '12

After a series of residues (1+) with alternate conformations, the H atom of the next residue seems to lock its occupancy to that of the final conformation of the previous residue. Hence I'm getting several H's with 0.1-0.5 occupancy on single conformation residues, where the entire rest of the residue has occupancy 1. This happens after every altconf series, and showed up in two separate structures I've worked on this week. Within a series of 2+ residues with altconfs, the H's refine correctly with their altconf groups. So it isn't just that all subsequent residue's H's are getting overridden. This behavior seems to have started with 1.8.1-1168, or at least previous refinements of the same structure didn't do this. 64 bit build on Retina MBP running OS X 10.8.2. I couldn't google any reference to this behavior, so assuming it's a new bug I can send along source files. Below I've pasted a couple residues from one of my PDBs. I added a ***** to mark the offending H with the wrong occupancy. -Alex Scouras ATOM 668 N CMET A 29 -17.970 6.464 59.720 0.21 10.72 N ANISOU 668 N CMET A 29 1641 1300 1133 -328 -379 -50 N ATOM 669 CA CMET A 29 -17.578 5.189 60.316 0.21 9.05 C ANISOU 669 CA CMET A 29 1545 882 1012 -248 -350 -140 C ATOM 670 C CMET A 29 -17.248 5.317 61.797 0.21 7.83 C ANISOU 670 C CMET A 29 1363 666 947 173 -134 166 C ATOM 671 O CMET A 29 -17.496 4.395 62.578 0.21 7.92 O ANISOU 671 O CMET A 29 1187 751 1071 371 -145 338 O ATOM 672 CB CMET A 29 -16.376 4.609 59.573 0.21 9.09 C ANISOU 672 CB CMET A 29 1561 834 1059 -243 -283 -203 C ATOM 673 CG CMET A 29 -16.681 4.210 58.145 0.21 9.05 C ANISOU 673 CG CMET A 29 1529 756 1154 -28 -331 -235 C ATOM 674 SD CMET A 29 -18.049 3.042 58.009 0.21 13.05 S ANISOU 674 SD CMET A 29 2356 1325 1279 -218 -284 -50 S ATOM 675 CE CMET A 29 -17.505 1.757 59.124 0.21 12.06 C ANISOU 675 CE CMET A 29 1866 1485 1232 -630 -667 420 C ATOM 676 H CMET A 29 -17.333 6.878 59.318 0.21 12.87 H ATOM 677 HA CMET A 29 -18.320 4.557 60.216 0.21 10.86 H ATOM 678 HB2CMET A 29 -15.671 5.275 59.551 0.21 10.91 H ATOM 679 HB3CMET A 29 -16.069 3.819 60.044 0.21 10.91 H ATOM 680 HG2CMET A 29 -16.916 5.004 57.640 0.21 10.86 H ATOM 681 HG3CMET A 29 -15.894 3.794 57.761 0.21 10.86 H ATOM 682 HE1CMET A 29 -18.164 1.059 59.140 0.21 14.48 H ATOM 683 HE2CMET A 29 -16.666 1.409 58.813 0.21 14.48 H ATOM 684 HE3CMET A 29 -17.401 2.130 60.002 0.21 14.48 H ATOM 685 N MET A 30 -16.667 6.445 62.187 1.00 6.77 N ANISOU 685 N MET A 30 1361 564 647 -30 -56 18 N ATOM 686 CA MET A 30 -16.416 6.699 63.606 1.00 7.66 C ANISOU 686 CA MET A 30 1362 946 601 -42 -173 2 C ATOM 687 C MET A 30 -17.759 6.724 64.367 1.00 7.62 C ANISOU 687 C MET A 30 1310 944 641 -9 -24 -22 C ATOM 688 O MET A 30 -17.852 6.212 65.486 1.00 9.05 O ANISOU 688 O MET A 30 1571 1151 718 84 98 231 O ATOM 689 CB MET A 30 -15.620 7.987 63.790 1.00 7.65 C ANISOU 689 CB MET A 30 1299 894 714 26 -120 -166 C ATOM 690 CG MET A 30 -14.229 7.965 63.231 1.00 8.62 C ANISOU 690 CG MET A 30 1464 1072 741 31 -163 -10 C ATOM 691 SD MET A 30 -13.130 6.745 64.006 1.00 9.04 S ANISOU 691 SD MET A 30 1555 961 919 105 -180 -94 S ATOM 692 CE MET A 30 -12.812 7.527 65.563 1.00 9.77 C ANISOU 692 CE MET A 30 1540 1284 889 150 -190 -62 C ATOM 693 H MET A 30 -16.395 7.062 61.653 0.21 8.12 H ***** ATOM 694 HA MET A 30 -15.872 5.965 63.962 1.00 9.19 H ATOM 695 HB2 MET A 30 -16.099 8.709 63.352 1.00 9.18 H ATOM 696 HB3 MET A 30 -15.550 8.173 64.739 1.00 9.18 H ATOM 697 HG2 MET A 30 -14.278 7.758 62.284 1.00 10.35 H ATOM 698 HG3 MET A 30 -13.831 8.840 63.354 1.00 10.35 H ATOM 699 HE1 MET A 30 -12.228 6.968 66.079 1.00 11.73 H ATOM 700 HE2 MET A 30 -12.397 8.378 65.406 1.00 11.73 H ATOM 701 HE3 MET A 30 -13.644 7.648 66.026 1.00 11.73 H

4 10

Re: [phenixbb] Troublesome molecular replacement
by Yuri 29 Oct '12

29 Oct '12

Nat, Just a few simple thoughts: --1) Are you certain about the space group, and did you try letting --Phaser pick one for you? (It sounds like you did, but I can't tell --from your description.) I searched for a solution in all possible space groups in P 4. I thought I was certain of the space group, but maybe not... --2) How much refinement did you run to get an R-free of 42$, and did --you already try running AutoBuild with rebuild-in-place turned off? --Refinement alone may not have the radius of convergence you need. True, I only did a couple of rounds of rebuild/refine --3) What is the sequence identity? My model is a point-mutant of the crystal. --4) Are you using the latest version of Phenix? Phaser has had TNCS --support for almost a year now but I keep hearing from people who are --using ancient versions... I am using the latest official PHENIX built. --- -- Yuri Pompeu

2 1

Troublesome molecular replacement
by Yuri 29 Oct '12

29 Oct '12

Dear community, I recently collected on a crystal that diffracted pretty well to 2.0 A. It indexes and scales fairly well in P 4 21 2 with cell dimensions 100, 100, 200. This would predict 4 molecules/ASU. I found a MR solution in P 41 21 2 but it does not refine very well Rfree =42%. Xtriage detects strong peaks at 0,0,0.25 for Translational NCS. Any suggestions on how to go about this one? Thanks -- Yuri Pompeu

2 1

autobuilding homomers
by Ingemar Andre 29 Oct '12

29 Oct '12

I have a question about the proper procedure for building homomers. My test case is a trimeric coiled-coil with 3-fold symmetry (trimer in assymetric unit). My input model is a trimer with three chains (A,B,C) with identical sequence. If I run phenix autobuild (rebuild_in_place=False) with this model I get a structure with decent R/Rfree. However, there are some issues: 1) Helix 1 consists of chain A in the output structure. There is a jump in the sequence number for chain A at the end. However, disregarding this chain break the full sequence is correct for helix 1. 2) Helix 2 consists of chain B + C in the output structure. The chain B has the correct sequence for helix 2. Chain C has partly the right sequence for helix 2. 3) Helix 3 consists of chain C in the output structure. Chain C has partly the correct sequence for helix 3. Lots jumps in the sequence numbering. GIven that I know its has three subunits with identical sequences I gather it should be possible to have autobuild generate structures with three chains with identical sequences. I have tested a few things. With and without giving it the sequence file corresponding to the sequence of a monomer. With and without enforcing ncs based on the input model (find_ncs=False input_ncs_file=find_ncs.ncs_spec). I am never able to build sequence-symmetric structures with rebuild_in_place=False for any homomer I have tried to build. I have limited experience with phenix so it is likely I am missing something obvious. In any case, I would appreciate your input on this. Best regards, Ingemar

3 3

Parallella
by Morten Groftehauge 27 Oct '12

27 Oct '12

Hi guys, So the Parallella project seems like it is going forward (with the smaller processor first off) http://www.kickstarter.com/projects/adapteva/parallella-a-supercomputer-for… I figure that phenix would have problems with low amount of RAM but are there other reasons that Phenix and Parallella might not play well together? Supposedly ships in May, btw. Cheers, Morten -- Morten K Grøftehauge, PhD Pohl Group Durham University

2 1

HYSS and Xtriage general questions
by Heather Condurso 26 Oct '12

26 Oct '12

Dear all, I am currently solving a structure using anomalous data from a Se-Met derivative. I collected a SAD set recently and am trying to use HYSS to find the sites. I am still not 100% sure which space group it is, but I am fairly confident that it is P43212/P41212. I have scaled the data with HKL and gave the resulting sca files to Xtriage. The recommended resolution cut off from Xtriage is 4.0 with the optimistic cut off at 3.6A. The anomalous signal is 0.15. I am hoping this data is good enough, but am also trying to get better data. In the meantime I am chugging along with the data I have. I have tried the fast version of HYSS with all of the higher P4 space groups varying the resolution and number of sites. I often find that for a given sca file the resolution given is important with 4.0A giving 1 site and 3.6A giving 24 out of 24 but with a cc of ~0.2. However, when I look at the sites and maps it seems a bit odd that all of sites lie in nearly a plane. So here are! a few questions: 1. My cell is rather large with dimensions of 132, 132, 195. Does expected CC for a good solution differ based on the size of the cell or the total number of expected sites? I expect 4 chains each with 6 sites. 2. Is it likely that I will find a solution with the anomalous signal only going to 4.0A? Are there other parameters that may clue me in as to whether or not it may be possible with this data? My data is also anisotropic. Does that play a role in finding a solution? I have other data sets, but in those the mosaicity varies outside the ideal range. In this set the mosaicity varies between 0.2 and 0.5 when I integrate the images without fixing the mosaicity. 3. I've also tried to use the full search, but I'm not sure if it is working or not. I am using the GUI with version 1.7.3-928 on a Mac running OS 10.6.8. I just realized now that this is not the latest version and I can try it on another machine that I know is running 1.8.1-1168. None the less, the run status window only shows that it is running, but the log display only shows the initial search parameters, not what it is currently doing. I also tried using different rescoring and the log display is blank. Does this mean it isn't doing anything? Or should I be more patient? Some of these have been running for days.... I appreciate any advice you may have to help me along! Heather Condurso UF-Department of Chemistry condurso(a)ufl.edu

2 1

adding some user-defined graphical objects in a PDB file (and displaying them)
by Francois Berenger 26 Oct '12

26 Oct '12

Hello, For some new project, I'd like to be able to generate things and store them in PDB format. For example, a triangle, a line segment, a square, a cube, a sphere, an arrow, etc. Being able to change the color and "line width" would be nice. Is there some official recommended way of doing this? Is there some software able to read and display such graphical annotations of PDB files? I'll also need the format description in that case. I want to be able to process a PDB file and store the result of my processing in the same PDB file as some kind of annotation. My current way of doing this is to discretize my objects as H atoms in some other output PDB file, but that's just a temporary workaround. My current search got me this: http://www.cgl.ucsf.edu/eccc1/ So, maybe there is some support for what I am looking for into Chimera. Thanks a lot for your suggestions, Francois.

1 0

Re: [phenixbb] Table 1 stastics
by Kay Diederichs 25 Oct '12

25 Oct '12

Dear all, I just want to point out that one possible source of discrepancy between I/sigma reported by a data processing program, and I/sigma reported by a program that uses the relevant cctbx routine is the following: that cctbx routine calculates the variance of merged data as max("internal variance", "external variance") which is different from what the data processing programs do (they calculate the "internal variance" - I hope I didn't get it the wrong way round). SHELXC and SHELXL also calculate the variance of the merged data like the data processing programs do. I dislike this "feature" of cctbx but I do not know if phenix.table_one actually uses this routine for this purpose. thanks, Kay P.S. you can find the discussion over at the cctbx mailing list; see e.g. http://phenix-online.org/pipermail/cctbxbb/2012-September/000530.html and other posts. On 10/18/2012 09:08 PM, phenixbb-request(a)phenix-online.org wrote: > Date: Thu, 18 Oct 2012 08:07:44 -0700 > From: Nathaniel Echols<nechols(a)lbl.gov> > To: PHENIX user mailing list<phenixbb(a)phenix-online.org> > Subject: Re: [phenixbb] Table 1 stastics > Message-ID: > <CALeAa1NDDy7=ZCUqThgun+k-Si7-Gf920ROvuYwpPgsShAhHVQ(a)mail.gmail.com> > Content-Type: text/plain; charset=ISO-8859-1 > > On Thu, Oct 18, 2012 at 7:53 AM, Andrew Waight<waight(a)msg.ucsf.edu> wrote: >> Thanks for Phenix! Anyway quick question, I tried out the "generate >> Table1" utility for fun and noticed that my I/SigI and completeness do not >> match what is found in my XSCALE file. >> ... >> How does Phenix compute these parameters? Obviously my I/sigI has gotten >> much better (1.59 versus 0.85) but the scaling statistics should be the >> correct ones. > > It computes these statistics directly from the reflection data you > provide, not the log file. (The reasoning being that since these are > the data that Phenix actually uses, the statistics calculated from > them are more accurate and relevant than whatever the scaling program > thinks.) My guess is that some of the internal processing accounts > for the difference - could you please send me the data file? > > FYI, the plan for the future is to recalculate all of these statistics > from unmerged data rather than rely on parsing data processing > logfiles, which obey no standard format and are subject to change > without notice. > > -Nat > -- Kay Diederichs http://strucbio.biologie.uni-konstanz.de email: Kay.Diederichs(a)uni-konstanz.de Tel +49 7531 88 4049 Fax 3183 Fachbereich Biologie, Universität Konstanz, Box M647, D-78457 Konstanz This e-mail is digitally signed. If your e-mail client does not have the necessary capabilities, just ignore the attached signature "smime.p7s".

4 5

Problem with elbow
by Subhani Bandara 25 Oct '12

25 Oct '12

Dear all, I am using the PHENIX version 1.8.1-1168 and need to get the pdb file and the cif file for the metal-chelate complex Gd-DTPA for refinement of a structure. I got the pdb file and the cif file for a dimer of the complex. Then I edited those to contain only one metal chelator complex and used as inputs in elbow. I have attached the pdb file and the cif file I used as inputs. I gave cif file as chemical file type with pdb file as initial geometry and ran elbow. Then I used pdb file only or cif file only as inputs separately. Every time it ended with the error “Failed to determine the bonding of a fragment of the molecule” Also I tried using smile string C(CN(CC(=O)[O-])CC(=O)[O-])N(CCN(CC(=O)[O-])CC(=O)[O-])CC(=O)[O-].[Gd+3]and pdb file as inputs. Then it runs without an error, but it outputs cif and pdb files containing only Gd. How can I get the cif file for the metal-chelate complex? I want to use the metal-chelate complex as one entity in the refinement, so that it will have the correct geometry of the chelator. Thanks Subhani

2 1