- 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 30 Oct '12

30 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 30 Oct '12

30 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 30 Oct '12

30 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 28 Oct '12

28 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 27 Oct '12

27 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 26 Oct '12

26 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 26 Oct '12

26 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

1.8.1 1168
by Colin Levy 25 Oct '12

25 Oct '12

Hi all, I have just updated to the latest version of Phenix 1.8.1 1168 but I now appear to have lost the ability to drag and drop files between different windows. Running on a mac 10.7.5 OSX. Thanks Colin Manchester Protein Structure Facility MIB 1.025 Tel 0161 306 5185 c.levy(a)manchester.ac.uk<mailto:[email protected]>

2 1

question du jour: monomer libs where?
by Patrick Loll 25 Oct '12

25 Oct '12

Where does PHENIX find the monomer libraries? Here's why I ask: I submit a refinement job for a peptide containing a mix of D- and L-amino acids. The program seems to have no problem with D-versions of Glu, Asn, or Phe (DGL, DAS, or DPN), but it flags the atoms in a D-ornithine (ORD) residue as having unknown nonbonded energy type symbols. Curiously, when I look in the relevant directories in ...Contents/phenix-1.8.1-1168/chem_data/chemical_components I can find cif files corresponding to ALL of these D-amino acids, including ORD (...chemical_components/o/data_ORD.cif). Files are also present for DGL, DAS, and DPN. However, these files don't contain all of the information I'm accustomed to seeing in monomer libraries (e.g., explicit bond lengths & angles), so I'm not sure if the same information is somehow encoded differently (via the idealized coordinates, perhaps?), or if I'm barking up the wrong tree. So: Are these the actual monomer libraries that phenix.refine consults? And in this case why doesn't phenix.refine recognize ORD? Or are the true monomer libraries hiding elsewhere? As always, thanks, Pat ps the directory structure I'm describing is under OS X 10.7.5 --------------------------------------------------------------------------------------- Patrick J. Loll, Ph. D. Professor of Biochemistry & Molecular Biology Director, Biochemistry Graduate Program Drexel University College of Medicine Room 10-102 New College Building 245 N. 15th St., Mailstop 497 Philadelphia, PA 19102-1192 USA (215) 762-7706 pjloll(a)gmail.com pat.loll(a)drexelmed.edu pat.loll(a)drexel.edu

4 7

some basic questions about phenix autobuild and phenix.refine
by Niks 24 Oct '12

24 Oct '12

Dear All, I have recently started using PHENIX so I am a beginer in PHENIX field. I am getting better statistics (Rfactor/Rfree) as compared to REFMAC so having strong inclination to continue with PHENIX but there are some discripancies I found when comparing with CCP4 programs.Therefore i I have some basic questions which I request you to throw some light and wisdom. 1) I have solved my structure using AutoMR and got two files MR.1.mtz and MR.1.pdb as output with high LLG of 8236.889. When going for AutoBuild afterwards it uses output pdb (MR.1.pdb) which is fine but uses MR.1.mtz as map coefficient. Shouldnot it uses original mtz from scala? 2) Same in Phenix.refine when going after Autobuild finishes and it asks for refining the model it uses same AutoMR output mtz (MR.1.mtz). As I was suggested by many that we should use original mtz coming from scala for Refinement to reduce model bias. shouldnot phenix.refine use Original mtz file here too? 3) After going for AutoBuild after AutoMR, even when I unselect "Place Waters" options in parameters for Refinement, I see output pdb file with waters added to it. I do not want waters to be added till I refine my model to the best , but Autobuild option comes out with waters added. Is there anyway I can have autobuilt pdb without waters added? 4) Another question that troubling me most that in Phenix.refine even if I use original mtz which clearly specifies resoultion range from 127.488-2.5A my refined pdb files comes out with Resolution range used for refinement as 37.8-2.5A. I have tried to input low resolution range as 127.488 but deafult is always 37.8 selected. My data may be very poor in lower resolutions than 37.8A but my completeness is also going down (65.5% as opposed to 79.8% in REFMAC). Can I select my refinement using my specified values of Resolution from mtz file coming out from SCALA? Sorry for the lengthy mail and many questions. Thanks and regards Nishant PhD Student, C/O Dr. C. G. Suresh, National Chemical Laboratory, Pune, India -- "The most difficult phase of life is not when No one understands you;It is when you don't understand yourself"

4 3

phenix GUI crash
by Sébastien Moniot 24 Oct '12

24 Oct '12

Hi, Little problem with the phenix GUI crashing ... we are running phenix via ssh login on a workstation and each time the REEL window or the graphical interface to make atom selection are called the whole GUI crashes. In the first case, only the following message is reported in the shell. --- segmentation fault ----- In the second case, --- (python2.7:29303) : Gtk-CRITICAL **: gtk_window_resize: assertion `width > 0' failed segmentation fault ---- Any idea how to fix this ? Best, Dr. Sebastien Moniot Department of Biochemistry University of Bayreuth Bayreuth, Germany

3 3

hydrogen nomenclature
by Patrick Loll 24 Oct '12

24 Oct '12

Hi, I'm trying to sort out some monomer libraries for a high-resolution refinement involving non-canonical amino acids, but now I've become side-tracked with a hydrogen nomenclature issue that seems to affect even canonical amino acids. I will use the example of L-leucine to illustrate: In the CCP4 monomer library LEU.cif, the two hydrogens on CB are named HB1 and HB2. However, if I look at the LEU entry on the RCSB PDB site, these two hydrogens are named HB2 and HB3. Looking at pdbe.org is not too helpful, since it lists those two hydrogens as HB2 and HB3, but includes alternate "PDB names" of HB1 and HB2. Similar situations pertain for other atoms and other residues. Presumably I need to stick with the nomenclature found in the monomer library in order for refinement to work properly; but will this come back to bite me later when I submit coordinates to the PDB? Any insights welcome. Thanks, Pat

4 3

Durandal clustering applied to homology models
by Quentin Delettre 23 Oct '12

23 Oct '12

Dear Phenix users, I am new to phenix and the Durandal tool, i would like to know if Durandal is an appropriate way to cluster many homology models of a protein. And if so, what is the syntax of the output ? The man page does not provide any information beside inputs commands. Thanks.

2 1

problem installing Phenix 1.8.1 on Fedora16
by Marianna Biadene 22 Oct '12

22 Oct '12

Dear all, I am trying to install Phenix 1.8.1-1168 on my Linux 64bit, Fedora 16. If I try to install it in the default directory (/usr/local) I have no problem, but when I try to install it in /usr/local/bin using the command written in the documentation (http://www.phenix-online.org/documentation/install.htm) ./install --prefix=</usr/local/bin> I got the following error message: -bash: syntax error near unexpected token `newline' I am not a super expert in Linux and I have tried to search in Internet for a solution but I am lost. Any idea of how could I solve this problem? Thanks a lot, Marianna -- Dr. Marianna Biadene Division of Biological Chemistry - Group Prof. Scheffzek Biocenter Medical University Innsbruck Innrain 80 / 4 6020 Innsbruck Austria Phone: +43 512 9003 70336

3 2

xyz and adp weight optimisation
by Alex Theodossis 19 Oct '12

19 Oct '12

Dear bb, I have made use of the target weight optimisation feature with great results. The question is, having optimised these weighting terms, how can I use the optimised values in subsequent runs without having to keep the optimisation option switched on? In the last round of refinement I get: Best XYZ weight: 9.000 Best ADP weight: 0.695 But the output .eff file does not seem to have saved these values for use in subsequent refinements: target_weights { optimize_xyz_weight = True optimize_adp_weight = True wxc_scale = 0.5 wxu_scale = 1 wc = 1 wu = 1 fix_wxc = None fix_wxu = None Is it appropriate to enter the optimised values in the fix_wxc and fix_wxu fields? Thanks in advance

2 1

Re: [phenixbb] Metal Coordination and symmetry operation
by Pavel Afonine 19 Oct '12

19 Oct '12

Hi Alex, thanks for the file - I will have a closer look. Did you mean the restraint did not work or it worked but you wanted something else? (I'm just trying to understand if this is a Phenix problem that we need to fix). Pavel On 10/18/12 10:30 PM, Alex Theodossis wrote: > Hi Pavel, > > thanks for the info. With it I was able to correctly define the > symmetry-related links as per this example: > > refinement.geometry_restraints.edits { > bond { > action = *add > atom_selection_1 = name ZN and chain C and resname ZN and resseq > 1 > atom_selection_2 = name NE2 and chain A and resname HIS and altid > A and resseq 39 > symmetry_operation = -x+1,y+1/2,-z+1/2 > distance_ideal = 2.030000 > sigma = 0.050 > } > > Unfortunately, this didn't have the desired result. Instead, the > addition of riding hydrogens fixed the problem. > I've also attached the relevant fragment of the refined model as > requested. > > Cheers, > > Alex > > > > On 17 October 2012 13:04, Pavel Afonine <pafonine(a)lbl.gov > <mailto:[email protected]>> wrote: > > Hi Alex, > > Nigel is the best person to reply, so he may provide more > comments. My understanding is that phenix.metal_coordination > command (or whatever the equivalent click in the GUI is) should > account for symmetry related atoms. If not you can fix it > manually. First create edits file using phenix.metal_coordination > and then add symmetry related bonds as described here: > > http://www.phenix-online.org/documentation/refinement.htm#anch353 > > Also, could you please send me the PDB file or its portion > containing relevant to your question atoms and I will have a > closer look? > > Pavel > > > On 10/16/12 6:33 PM, Alex Theodossis wrote: >> >> Dear bb, >> >> I am trying to refine a Zn ion coordinated by symmetry >> related residues. To avoid groups being pushed away. I have >> use phenix.metal_coordination to >> generate the necessary link edits. However only non-symmetry >> related interactions are accounted for. Is there a relevant >> option I can switch on? >> Alternatively, does anyone have an example of the syntax >> required to define the necessary symmetry operations in the >> elbow.edits file? >> >> Thanks in advance. >> >> >> >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org <mailto:[email protected]> >> http://phenix-online.org/mailman/listinfo/phenixbb > >

1 0

Re: [phenixbb] Metal Coordination and symmetry operation
by Alex Theodossis 19 Oct '12

19 Oct '12

> > Dear bb, > > I am trying to refine a Zn ion coordinated by symmetry related residues. > To avoid groups being pushed away. I have use phenix.metal_coordination to > generate the necessary link edits. However only non-symmetry related > interactions are accounted for. Is there a relevant option I can switch on? > Alternatively, does anyone have an example of the syntax required to > define the necessary symmetry operations in the elbow.edits file? > > Thanks in advance. >

2 2

Re: [phenixbb] Table 1 stastics
by Nathaniel Echols 19 Oct '12

19 Oct '12

On Thu, Oct 18, 2012 at 8:45 AM, Katherine Sippel <katherine.sippel(a)gmail.com> wrote: > Part of the problem could lie in whether he input the original mtz as > intensities or not. In older versions of phenix some intensities were > dropped if they were below the cutoff and particularly for high res data the > completeness would go down significantly. I know you have implemented the > French and Wilson conversion in the newer versions but if he comes from a > lab like my old one and is refining off of intensities that could be the > culprit. This is a good point, but it depends on the version of Phenix, the type of input data, and how the data file was processed. I checked my code for Table 1 and if the input data are intensities, the negative observations should *not* be discarded. The same goes for phenix.refine, among other programs. However, if amplitudes were used as input for Table 1, and had previously been filtered through a program that does not perform the French & Wilson treatment, they may be missing some data (most of which will probably be at high resolution). One caveat: the reflection file editor in Phenix has an option to convert intensities to amplitudes, but this will not run F&W, which is perhaps a bug... -Nat

1 0

Table 1 stastics
by Andrew Waight 18 Oct '12

18 Oct '12

Hi All, 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. From Table 1 Completeness (%) 98.62 (90.57) Mean I/sigma(I) 7.39 (1.59) Wilson B-factor 59.78 R-sym R-factor 0.2006 (0.3275) R-free 0.2241 (0.3739) From XSCALE SUBSET OF INTENSITY DATA WITH SIGNAL/NOISE >= -3.0 AS FUNCTION OF RESOLUTION RESOLUTION NUMBER OF REFLECTIONS COMPLETENESS R-FACTOR R-FACTOR COMPARED I/SIGMA R-meas CC(1/2) Anomal SigAno Nano LIMIT OBSERVED UNIQUE POSSIBLE OF DATA observed expected Corr 20.00 313 69 72 95.8% 3.6% 4.6% 313 30.88 4.3% 99.9* 8 0.725 54 10.00 2174 471 471 100.0% 3.9% 4.8% 2174 28.71 4.4% 99.8* 34* 0.972 382 6.00 8195 1983 1992 99.5% 4.7% 5.3% 8160 20.70 5.4% 99.8* 12* 0.913 1299 5.00 8055 1818 1820 99.9% 7.0% 6.9% 8051 17.73 7.9% 99.5* 7 0.892 1337 4.50 7207 1634 1636 99.9% 6.0% 6.2% 7207 20.06 6.8% 99.6* -1 0.824 1163 4.00 9603 2479 2521 98.3% 6.8% 6.7% 9536 17.05 7.9% 99.4* -2 0.819 1369 3.60 13641 3141 3156 99.5% 9.6% 9.5% 13621 13.71 11.0% 99.1* -2 0.799 2073 3.20 22112 4948 4951 99.9% 16.5% 16.6% 22083 8.82 18.8% 98.1* 4 0.810 3383 3.00 14155 3550 3579 99.2% 27.8% 27.7% 14058 5.45 32.1% 94.6* -1 0.742 1875 2.80 19676 4584 4611 99.4% 44.9% 46.2% 19604 3.71 51.2% 88.5* 3 0.740 2750 2.60 26749 6192 6202 99.8% 76.3% 78.4% 26686 2.20 87.2% 71.8* 1 0.717 3598 2.40 33999 8339 8388 99.4% 130.9% 132.7% 33810 1.26 150.6% 48.5* 1 0.671 4170 2.30 22347 5332 5348 99.7% 198.4% 206.4% 22276 0.85 227.7% 31.5* 0 0.623 2685 total 188226 44540 44747 99.5% 9.0% 9.4% 187579 7.09 10.3% 99.8* 2 0.755 26138 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. Thanks, Drew

2 2

close contact and torsion angle
by hchxia＠ujs.edu.cn 13 Oct '12

13 Oct '12

Hi, There are some close contacts (hydrogen bonds) in my refined model. Can I just ignore it, delete it or what? I added hydrogen for refinement since the resolution is 1.47 A. Should I remove all these added hydrogens and refine again? Second question is there are several outilers of torsion angles (6 from Ramachandran plot in coot ). But 16 are outside expected Ramachandran regions from the validation report from PDB bank(ADIT) . What is the cause? And should I re-refine to correct them? Or just go ahead to deposit it? Thanks! chuan

4 4

how to continue an aborted phenix.refine
by hchxia＠ujs.edu.cn 11 Oct '12

11 Oct '12

Hi! I aborted a run of phenix.refine due to linux problem. Is there a way to continue from the aborted point to finish the run? or I have to start from 0 to rerun it ? Thanks a lot! chuan

2 1

OS X 10.8 Mountain Lion
by Alexander Batyuk 10 Oct '12

10 Oct '12

Dear all, Could anybody confirm that Phenix does or does not run on OS X 10.8 Mountain Lion, please? Best wishes, Alex -- Alex Batyuk The Plueckthun Lab www.bioc.uzh.ch/plueckthun

3 5

labelit and latest mosflm
by "Weiergräber, Oliver H." 05 Oct '12

05 Oct '12

Hello, the labelit.index executable provided with phenix 1.8.1 has a problem with the latest mosflm (7.0.9, as included in ccp4 6.3.0): "MosflmVersionError: mosflm versions after 7.0.8 not supported" Is there any reason to believe that mosflm 7.0.9 has really become incompatible with labelit? If not, can the version check be overridden as a temporary fix? Best, Oliver ================================================ PD Dr. Oliver H. Weiergräber Institute of Complex Systems ICS-6: Structural Biochemistry Tel.: +49 2461 61-2028 Fax: +49 2461 61-9540 ================================================ ------------------------------------------------------------------------------------------------ ------------------------------------------------------------------------------------------------ Forschungszentrum Juelich GmbH 52425 Juelich Sitz der Gesellschaft: Juelich Eingetragen im Handelsregister des Amtsgerichts Dueren Nr. HR B 3498 Vorsitzender des Aufsichtsrats: MinDir Dr. Karl Eugen Huthmacher Geschaeftsfuehrung: Prof. Dr. Achim Bachem (Vorsitzender), Karsten Beneke (stellv. Vorsitzender), Prof. Dr.-Ing. Harald Bolt, Prof. Dr. Sebastian M. Schmidt ------------------------------------------------------------------------------------------------ ------------------------------------------------------------------------------------------------ Kennen Sie schon unsere app? http://www.fz-juelich.de/app

2 1

automr
by Shya Biswas 03 Oct '12

03 Oct '12

Hi all, I was wondering if anyone knows how to change the packing clash distance in Automr gui. Thanks, Shya

2 1

changes from 1.8.1-1161 to official 1.8.1?
by Eric Williams 02 Oct '12

02 Oct '12

Where any major bugs fixed between the 1161 nightly build and the official 1.8.1 release? I have a large batch job running that uses 1161. If the squashed bugs are big enough, it might be worth stopping the job to install the new release. Otherwise, I'll just let it run to completion. Thanks. :) Eric

2 1

Phenix version 1.8.1 released
by Paul Adams 02 Oct '12

02 Oct '12

The Phenix developers are pleased to announce that version 1.8.1 of Phenix is now available. Binary installers for Linux, Mac OSX, and Windows platforms are available at the download site: http://phenix-online.org/download/ Some of the new features in this version are: Graphical interface: ==================== - New GUIs: phenix.find_alt_orig_sym_mate General: ======== - Multiple bug fixes - New commands: phenix.erraser, phenix.merging_statistics - Added geometry restraints for over 400 non-standard amino acid entries Refinement: =========== - Overall B-cart now applied to individual atomic B-factors (restores pre-1.8 behavior) - Improvements to torsion NCS rotamer correction - Real-space refinement (including rotamer fitting) for twinned data New command: ============ - phenix.erraser: - New method from Fang-Chieh Chou and Rhiju Das (Stanford) combining Rosetta nucleic acid rebuilding with phenix.refine (requires separate installation of Rosetta 3.4 or more recent) Molecular Replacement: ====================== - Phaser updated to version 2.5.2: - Scattering from HETATM records that represent modified amino acids (e.g. MSE/SeMet) included so that maps no longer show holes for these HETATM records in the model; occupancies will not be reset. - Packing test was incorrectly applied to reoriented molecules where there was PTNCS Experimental Phasing and Model Building: ======================================== - AutoSol: - Scale of derivative in SIR (for RIP phasing) can be set with keyword "derscale" - AutoBuild: - Unassigned chains are now given the segid of "UNK" and are normally given chain ID="U" Ligands: ======== - eLBOW: - New option, --multiple-planes, to generate planes around double bonds - Saturated rings now have dihedral restraints to maintain chair conformation - REEL: - View relationships between ligands using PubChem fingerprint and tanimoto score. Miscellaneous: ============== - phenix.real_space_correlation: - Simplified command-line behavior - phenix.merging_statistics (new command): - Calculate statistics for unmerged data in a variety of formats, including <I/sigmaI>, R-merge, R-meas, R-pim, CC1/2 - References: Diederichs K & Karplus PA (1997) Nature Structural Biology 4:269-275 Weiss MS (2001) J Appl Cryst 34:130-135. - phenix.kinemage: - Improvements to handling of alternate conformations. - phenix.combine_models: - Combine two models, taking the best parts of each to create a new model For a full list of changes see: http://www.phenix-online.org/documentation/CHANGES Please note that this publication should be used to cite use of Phenix: PHENIX: a comprehensive Python-based system for macromolecular structure solution. P. D. Adams, P. V. Afonine, G. Bunkóczi, V. B. Chen, I. W. Davis, N. Echols, J. J. Headd, L.-W. Hung, G. J. Kapral, R. W. Grosse-Kunstleve, A. J. McCoy, N. W. Moriarty, R. Oeffner, R. J. Read, D. C. Richardson, J. S. Richardson, T. C. Terwilliger and P. H. Zwart. Acta Cryst. D66, 213-221 (2010). Full documentation is available here: http://www.phenix-online.org/documentation/ There is a Phenix bulletin board: http://www.phenix-online.org/mailman/listinfo/phenixbb/ Please consult the installer README file or online documentation for installation instructions. Direct questions and problem reports to the bulletin board or: help(a)phenix-online.org and bugs(a)phenix-online.org Commercial users interested in obtaining access to Phenix should visit the Phenix website for information about the Phenix Industrial Consortium. The development of Phenix is principally funded by the National Institute of General Medical Sciences (NIH) under grant P01-GM063210. We also acknowledge the generous support of the members of the Phenix Industrial Consortium. -- Paul Adams Deputy Division Director, Physical Biosciences Division, Lawrence Berkeley Lab Division Deputy for Biosciences, Advanced Light Source, Lawrence Berkeley Lab Adjunct Professor, Department of Bioengineering, U.C. Berkeley Vice President for Technology, the Joint BioEnergy Institute Laboratory Research Manager, ENIGMA Science Focus Area Building 64, Room 248 Tel: 1-510-486-4225, Fax: 1-510-486-5909 http://cci.lbl.gov/paul Lawrence Berkeley Laboratory 1 Cyclotron Road BLDG 64R0121 Berkeley, CA 94720, USA. Executive Assistant: Louise Benvenue [ LBenvenue(a)lbl.gov ][ 1-510-495-2506 ] --

1 0

point-n-click installer is bjorking OS X Dock
by Scott Classen 27 Sep '12

27 Sep '12

Hello, After running the OS X point-n-click installer for the 1.8.1-1160 nightly build, the dock becomes unresponsive and I can no longer command+tab to switch between applications. I've reproduced this on 2 Macs. One running 10.8.2 and the other running 10.7.5 Relaunching Finder does not fix the proton dock/app switching Logging out and then logging back in fixes the problem. Scott

2 5

den refinement with multiple reference files
by Andreas Förster 27 Sep '12

27 Sep '12

Dear all, I've started DEN refinement on a structure obtained by MR from data optimistically going to 3.8A. There are several chains of two kinds in the structure. High-resolution structures exist for both. How do I use two reference files for DEN refinement, and how do I assign which chain corresponds to which reference? (Currently, I'm using reference_file = None.) Thanks. Andreas -- Andreas Förster, Research Associate Paul Freemont & Xiaodong Zhang Labs Department of Biochemistry, Imperial College London http://www.msf.bio.ic.ac.uk

2 2

Phenix cut out density labels choice
by Mikalai Lapkouski 21 Sep '12

21 Sep '12

Hello, I have a following issue with running Phenix.cutout density either on command line or in GUI. I have multiple columns in the mtz.file and I want to choose map coefficients FOFCWT_d1-8,PHFOFCWT_d1-8 for the difference map to be cut around pdb selection. I run a command: phenix.cut_out_density model.pdb my.mtz cutout_fix_position=True cutout_model_radius=10 cutout_type=model \ map_coeff_labels=FOFCWT_d1-8,PHFOFCWT_d1-8 Problem: The output cutout.mtz has FWT and PHWT (which I did not have at all in the input) and not my difference map coefficients, and obviously the map is not my difference map. In GUI is the same problem, I choose columns I need but it outputs FWT and PHWT again. Please, let me know what I am doing wrong. Thank you very much! Mikalai

3 2

Combine map labels
by Mikalai Lapkouski 21 Sep '12

21 Sep '12

Hello, I have a following question. I have generated several Fo-Fc sa omit maps (FOFCWT, PHFOFCWT) for different portions of a DNA molecule. Now that I have each omit map as a separate file to display the whole DNA omit map I have to load them in coot one by one. Is there a way of how to combine all this maps into one file and have a single name for label columns so that the whole map will open at once? I managed to put them all in one file but column labels are separate: FOFCWT_1, PHFOFCWT_1 FOFCWT_2, PHFOFCWT_2...... Can they all go under FOFCWT, PHFOFCWT? Thank you very much for any suggestion Mikalai

5 5

Specifying union of two chains in pdb_atom_selection
by Francis E Reyes 21 Sep '12

21 Sep '12

Hi all, I issue the following command to pdb_atom_selection but I get 0 atoms selected at the end? $ phenix.pdb_atom_selection refine_to_submit.pdb "chain X and chain B" refine_to_submit.pdb Monomer Library directory: "/usr/local/phenix-1.8-1069/chem_data/mon_lib" Total number of atoms: 3817 Number of models: 1 Model: "" Number of chains: 6 Chain: "B" Number of atoms: 1800 Number of conformers: 1 Conformer: "" Number of residues, atoms: 84, 1800 Classifications: {'RNA': 84} Modifications used: {'rna3p_pyr': 26, 'rna2p_pur': 5, 'rna3p_pur': 44, 'rna2p_pyr': 9} Link IDs: {'rna3p': 70, 'rna2p': 13} Chain: "X" Number of atoms: 1800 Number of conformers: 1 Conformer: "" Number of residues, atoms: 84, 1800 Classifications: {'RNA': 84} Modifications used: {'rna3p_pyr': 26, 'rna2p_pur': 5, 'rna3p_pur': 44, 'rna2p_pyr': 9} Link IDs: {'rna3p': 70, 'rna2p': 13} chain X and chain B 0 atoms selected Thanks! F --------------------------------------------- Francis E. Reyes PhD 215 UCB University of Colorado at Boulder

2 1

Scaling for Isomorphous Difference Maps
by Michael Thompson 20 Sep '12

20 Sep '12

Hello All, I'm currently using the GUI to calculate some Fo-Fo isomorphous difference maps. From looking at the Phenix documentation, it hasn't become apparent to me if the tool that performs this calculation automatically scales the two sets of amplitudes to one another, or if the input mtz files need to be scaled together elsewhere before running the map utility. Any info would be appreciated! Thanks, Mike -- Michael C. Thompson Graduate Student Biochemistry & Molecular Biology Division Department of Chemistry & Biochemistry University of California, Los Angeles miket(a)chem.ucla.edu

3 2

Two versions of TFZ, molecular replacement
by 赵岩 20 Sep '12

20 Sep '12

Hi everyone, I am trying to solve a structure though molecular replacement and I have obtained a suspected solution. The header of the solution as follows, REMARK TITLE [no title set] REMARK Log-Likelihood Gain: 238.884 REMARK RFZ=10.9 TFZ=13.9 PAK=2 LLG=235 TFZ==5.0 LLG=239 TFZ==5.0 REMARK ENSEMBLE 1 EULER 15.53 0.10 344.54 FRAC -0.499 -1.000 -0.499 At first sight, the solution has no problem(TFZ>8), but the second and third TFZ are very low. For a correct solution, the later TFZs always higher than the first TFZ. So I don't konw how to understand the later TFZs and what extent I can refer to them. Thanks a lot! --

5 4

advice for speeding up cif_as_mtz
by Eric Williams 20 Sep '12

20 Sep '12

I'm currently running cif_as_mtz and model_vs_data on the whole PDB as part of my dissertation. Time has become a critical concern, and cif_as_mtz (dev-1141) is running very, very slowly (1-4 minutes per entry). Here's how I'm calling it: phenix.cif_as_mtz <sf_cif_fn> --ignore-bad-sigmas --merge --remove-systematic-absences --map-to-asu --symmetry=<pdb_model_fn> --output-file-name=<sf_mtz_fn> I'm hoping that perhaps someone could help me determine which command line argument(s) could be removed to provide maximum speed-up. Thanks. :) Eric

4 8

B-iso vs. B-aniso
by Yuri 19 Sep '12

19 Sep '12

Dear community, The protein model I am refining has 400 amino acids (3320 atoms). Some real quick calculations tell me that to properly refine it anisotropically, I would need 119,520 observations. Given my unit-cell dimension and space-group it is equivalent to about a 1.24 A complete data set. However, I have had a couple of cases where anisotropic B-factor refinement significantly improved R-work and R-free, while maintaining a reasonable R-gap, for lower resolution models (1.4-1.5 A, around 70,000 reflections). What is the proper way of modelling the B-factors? Any thoughts/opinions/call-outs from the community are welcome. Cheers, -- Yuri Pompeu

4 3

Monitoring MR_ROSETTA progress on a cluster
by Dmitry Rodionov 18 Sep '12

18 Sep '12

Hi Tom, How are you doing? Sorry to bother you again. I got my hands on a login to local cluster (MOAB) and started MR_rosetta there as well. Problem is, I have hard time figuring out what it is doing at this moment and if everything is under control. All I see in the queue are the original submission and a process it spawned less then a minute later. While a day ago there were 30 processes (that's how much I can have). All mr_rosetta.log files end with a list of methods to be run. MR_ROSETTA_1/GROUP_OF_PLACE_MODEL_1/LIST_OF_JOBS_RUNNING says: /sb/home/drodionov/MR_ROSETTA_1/GROUP_OF_PLACE_MODEL_1/RUN_1 /sb/home/drodionov/MR_ROSETTA_1/GROUP_OF_PLACE_MODEL_1/RUN_FILE_1.log MISSING FINAL COUNT: 0 But the file happily exists : GROUP_OF_PLACE_MODEL_1/RUN_FILE_1.log .... Starting job 1000...Log will be: /sb/home/drodionov/MR_ROSETTA_1/GROUP_OF_PLACE_MODEL_1/RUN_1/REBUILD_IN_SETS_1/RUN_FILE_1000.log And that job successfully finished a while ago. Is everything ok with the run? Am I just being impatient? Thanks! Best regards, Dmitry

1 0

Custom restraints editor
by Kristof Van Hecke 17 Sep '12

17 Sep '12

Dear, When I try to define additional 'planarity' restraints via the 'Custom restraints editor', I get the following error (see below): (I'm using the Phenix nightly build dev-1148 on a Mac Book Pro, OS X 10.5.8, but had the same problem with the official release version 1.8) I also tried to add these additional restraints to my restraints .cif file, but then rmsd values increase to unreasonably high values. (I just want to keep 2 chlorine atoms on a phenyl ring better in the same plane) Any ideas please..? Thank you very much Kristof Van Hecke --------------------- untimeError : Syntax error: missing closing quote (python_object, line 1) Traceback (most recent call last): File "/Applications/phenix-dev-1148/phenix/wxGUI2/Programs/Refine/ Dialogs.py", line 296, in OnUpdate new_phil = self.GetPhilObject() File "/Applications/phenix-dev-1148/cctbx_project/wxtbx/ custom_restraints.py", line 646, in GetPhilObject new_phil = self._index.master_phil.format(python_object=self._params) File "/Applications/phenix-dev-1148/cctbx_project/libtbx/phil/ __init__.py", line 1714, in format result.append(object.format(sub_python_object)) File "/Applications/phenix-dev-1148/cctbx_project/libtbx/phil/ __init__.py", line 1714, in format result.append(object.format(sub_python_object)) File "/Applications/phenix-dev-1148/cctbx_project/libtbx/phil/ __init__.py", line 1714, in format result.append(object.format(sub_python_object)) File "/Applications/phenix-dev-1148/cctbx_project/libtbx/phil/ __init__.py", line 1727, in format result.append(object.format(sub_python_object_i)) File "/Applications/phenix-dev-1148/cctbx_project/libtbx/phil/ __init__.py", line 1727, in format result.append(object.format(sub_python_object_i)) File "/Applications/phenix-dev-1148/cctbx_project/libtbx/phil/ __init__.py", line 1094, in format words = type_as_words(python_object=python_object, master=self) File "/Applications/phenix-dev-1148/cctbx_project/libtbx/phil/ __init__.py", line 151, in as_words source_info="python_object") File "/Applications/phenix-dev-1148/cctbx_project/libtbx/phil/ __init__.py", line 60, in tokenize_value_literal tokenizer.settings(contiguous_word_characters="")])) File "/Applications/phenix-dev-1148/cctbx_project/libtbx/phil/ tokenizer.py", line 249, in next where_str(O.source_info, char_iter.line_number))) RuntimeError: Syntax error: missing closing quote (python_object, line 1) Platform info: __FILE__ = /net/patchnose/scratch1/phenix/phenix-dev-1148/ cctbx_project/boost_adaptbx/meta_ext.cpp __DATE__ = Sep 3 2012 __TIME__ = 06:05:56 __i386__ __APPLE_CC__ = 5367 __GNUC__ = 4 __GNUC_MINOR__ = 0 __GNUC_PATCHLEVEL__ = 1 boost::python::cxxabi_cxa_demangle_is_broken(): false __GXX_WEAK__ = 1 __VERSION__ = 4.0.1 (Apple Computer, Inc. build 5367) FE_INEXACT = 32 FE_DIVBYZERO = 4 FE_UNDERFLOW = 16 FE_OVERFLOW = 8 FE_INVALID = 1 FE_ALL_EXCEPT = 61 __SSE2__ = 1 BOOST_LIB_VERSION = 1_51 PY_VERSION = 2.7.3 PYTHON_API_VERSION = 1013 sizeof(bool) = 1 sizeof(short) = 2 sizeof(int) = 4 sizeof(long) = 4 sizeof(std::size_t) = 4 sizeof(void*) = 4 sizeof(long long) = 8 sizeof(float) = 4 sizeof(double) = 8 sizeof(long double) = 16 sizeof(boost::int32_t) = 4 sizeof(boost::uint32_t) = 4 sizeof(wchar_t) = 4 __PTRDIFF_TYPE__ Py_USING_UNICODE sizeof(PY_UNICODE_TYPE) = 2 BOOST_ADAPTBX_TYPE_ID_SIZE_T_EQ_UNSIGNED_LONG boost/TAG: 79991 os.name: posix sys.platform: darwin sys.byteorder: little platform.platform(): Darwin-9.8.0-i386-32bit platform.architecture(): ('32bit', '') floating_point_exceptions.division_by_zero_trapped: False floating_point_exceptions.invalid_trapped: False floating_point_exceptions.overflow_trapped: False number of processors: 2 Memory total: 4,294,967,296 Memory free: unknown import thread: OK "hello" = ['h', 'e', 'l', 'l', 'o'] u"hello" = ['h', '\x00', 'e', '\x00', 'l', '\x00', 'l', '\x00', 'o', '\x00'] u"\u00C5" = ['\xc5', '\x00'] as utf-8 = ['\xc3', '\x85'] LATIN CAPITAL LETTER A WITH RING ABOVE = √Ö

2 2

Sa_Omit_map_DNA_low_resolution
by Mikalai Lapkouski 16 Sep '12

16 Sep '12

Hello, I want to generate a simulated annealing omit map around DNA molecule in my protein-DNA complex for publication. My resolution is very low - 5A. 1. I deleted the whole DNA (947 atoms of DNA against 7543 for protein part) from pdb and ran phenix.refine with simulated_annealing=True. But resulted Fo-Fc in the omitted region is very weak. Is it because of the low resolution and because big part of the model was omitted? 2. Next, I wanted to generate a composit_sa_omit map as in this case DNA would be deleted in small portions and the resulted map would give me a map I can use for publication: phenix.autobuild data=inp.mtz model=whole.pdb (contains DNA) maps_only=True (to skip rebuilding) composite_omit_type=sa_omit omit_selection=DNA.pdb (contains only DNA coordinates). But the job runs too long and still it seems like it is doing rebuilding. Is this a correct command and correct action to create a sa_omit map for the DNA in my complex? Thank you a lot for the help, Nick

3 3

Re: [phenixbb] Two versions of TFZ, molecular replacement
by yzhao 14 Sep '12

14 Sep '12

Hi, Some information may be needed. The model used to do molecular replacement is not an original homolog PDB file. I feed the homolog structure and sequence alignment to MR_rosetta, and I obtained an overall_best.pdb from AutoBuild(diffraction resolution is about 3. 5A). So I use this overall_best.pdb to run Phaser, and achieved a solution as I just mentioned. the overall_best.pdb from AutoBuild has no obvious regular secondary structure relative to the homolog structure, but their overall appearances are very similar. Before I run Phaser, I also checked native Patterson, there is no large off- origin peak. Additionally, the self rotation function and matthews_coef indicate there is a monomer in ASU. Thanks a lot! -- Yan Zhao Randy Read 写:Hi, There's an explanation of the TFZ== entry and how to interpret the job history on our web page:http://www.phaser.cimr.cam.ac.uk/index.php/ Molecular_Replacement#Has_Phaser_Solved_It.3F Basically, the TFZ=13.9 means that, in the translation search using the orientation that came from the rotation function, the Z-score of this peak was 13.9. However, we found that the TFZ was too sensitive to random errors in the orientation, so Phaser now computes what the TFZ score would have been if the orientation in the translation search had been the optimised one from the rigid- body refinement. (We call this the TFZ-equivalent, denoted as TFZ==.) So after the refinement (giving LLG=235), a random sampling of translations shows that the TFZ score would have been 5.0 for the refined orientation. The next two entries are the LLG after the final refinement at full resolution (239) and the TFZ score for that refined solution (again 5.0). What you see is very unusual behaviour. The TFZ -equivalent is almost always higher than the original TFZ, and I can't recall a case where it was significantly lower. One possible explanation occurs to me. If your crystal has translational NCS ( indicated by a large off-origin peak in the native Patterson) and you're searching for a dimer, a translation search with the dimer axis parallel to a crystallographic 2-fold would give a very sudden increase in the LLG when the dimer is in the right place relative to the crystallographic 2-fold to generate a model that obeys the translational NCS. Then, if the refinement changed the orientation of the dimer so that its axis was less parallel to the crystallographic 2-fold, that would reduce the modulation of the calculated data and reduce the TFZ. If these are the circumstances of your job, then you should be able to fix it by making sure that you have a recent version of Phenix, then running a job looking for 2 copies of the monomer. In that case, Phaser will take proper account of the translational NCS and there won't be any strange artefacts. If those aren't the circumstances, then I'm very puzzled and would appreciate receiving the data, model and other information needed to reproduce your job ( off-line) to check out what is going on! Best wishes, Randy Read On 14 Sep 2012, at 20:13, 赵岩 wrote:Hi everyone, I am trying to solve a structure though molecular replacement and I have obtained a suspected solution. The header of the solution as follows, REMARK TITLE [no title set] REMARK Log-Likelihood Gain: 238.884 REMARK RFZ=10.9 TFZ=13.9 PAK=2 LLG=235 TFZ==5.0 LLG=239 TFZ==5.0 REMARK ENSEMBLE 1 EULER 15.53 0.10 344.54 FRAC -0.499 -1.000 -0.499 At first sight, the solution has no problem(TFZ>8), but the second and third TFZ are very low. For a correct solution, the later TFZs always higher than the first TFZ. So I don't konw how to understand the later TFZs and what extent I can refer to them. Thanks a lot! -- _______________________________________________ phenixbb mailing list phenixbb(a)phenix-online.org http://phenix-online.org/mailman/listinfo/phenixbb ------Randy J. ReadDepartment of Haematology, University of CambridgeCambridge Institute for Medical Research Tel: + 44 1223 336500Wellcome Trust/MRC Building Fax: + 44 1223 336827Hills Road E-mail: rjr27(a)cam.ac.ukCambridge CB2 0XY, U.K. www-structmed.cimr.cam.ac.uk

2 1

Anisotropic scale matrix in phenix.refine
by Alex Taylor 13 Sep '12

13 Sep '12

2 1

Ligand geometry, Ideal vs. Observed
by Yuri 13 Sep '12

13 Sep '12

Hi everyone, I am trying to show that a ligand underwent catalysis during a soaking experiment. One of the things I would like to show is the geometry of the ligand, bond angles/lengths, dihedrals, etc... One of my models has a hi-res of 1.18A and the ligand density is really clear and complete. What is the best way to refine the ligand unrestrained and then generate measurements? Also, the idea is to finally compare to ideal geometry. How should I generate these values (any softwares in mind)? ANy idea is welcome. Thanks a lot -- Yuri Pompeu

2 1

Bond/Angle restraints for multiple copies of ligands/cofactors
by Eckhard Hofmann 12 Sep '12

12 Sep '12

Hi all, if the xray data reach real atomic resolution (say <1A), one can nicely refine cofactors w/o restraints. In intermediate cases, the existence of multiple copies allows to use the "SAME" option in SHELX refinement to restrain identical molecules to each other, somewhat like NCS, but on bonds and angles. Is there a possibility to do something similar with PHENIX? I have structures with up to 16 or more copies in the ASU (and several copies per chain), and 1.1A resolution. This would give me a very stable way to access unbiased information about stereochemistry. (Using NCS is of course not an powerful idea, as the torsion angles might vary between the different copies). Thanks for feedback, Eckhard -- Eckhard Hofmann <eckhard.hofmann(a)bph.ruhr-uni-bochum.de> Ruhr-Uni Bochum AG Proteinkristallographie, LS Biophysik, ND04/316 44780 Bochum Tel: +49-(0)234/32-24463, Sekr. -24461, FAX: -14762

2 1

Advice on final refinement of twinned R3 structure: DETWIN and ML-target vs twin_lsq_sq target
by Bjørn Panyella Pedersen 11 Sep '12

11 Sep '12

Dear list I would like some advice on how to do the final refinement of a twinned structure: I have a 3.0A structure (sg R3) which is twinned. I have collected ~50 datasets and the estimated twin-fraction varies from 0.2 to 0.45 based on phenix.xtriage analysis (britton plot etc). The structure could only be solved by experimental phasing using _detwinned_ datasets (detwinned by the ccp4 program DETWIN). Refinement of the model was initially done using my highest resolution _detwinned_ dataset (3.0A, britton_alpha:0.21) and a ML target in phenix.refine. Lately I have switched to th twin_lsq_f target and refined using the 'unmodified'/twinned dataset. 1. I have tried systematically to detwin the dataset with different twin-fractions and then refine using the ML-target and my best refinements comes using an alpha of 0.2 (see attached figure). This fits the Britton plot etc. 2. xtriage estimated the twin-fraction to by ~0.21, but using twin_lsq_f, phenix.refine now gets alpha to be 0.31. This seems like a large increase to me? 3. Using the manually detwinned dataset and a ML-target gives better R-factors than using the twin_lsq_f target (see attached figure). Bottom line: How should I report/deposit this structure: A. Should I use the manually detwinned data and an ML-target to get a model that explains the data best (lowest R-factors). And deposit my structure with modified (i.e DETWINNED) structure factors. I fear that people would react strongly to the deposition of modified F_obs in the pdb. I could deposit both twinned and detwinned structure-factores but this might also confuse? OR B. Should I use the twin_lsq_f target and accept that my deposited model is not as good as I know it could be. And deposit my structure with the real (i.e TWINNED) structure factors. OR C. Something I'm not aware of means that refinement using the ML-target is not allowed on detwinned datasets, causing me to have artificially low R-factors. I hope I'm not missing something basic here.:) Thank for any advice, pointers and ideas! -Bjørn -- Bjørn Panyella Pedersen, PhD Macromolecular Structure Group Dept. of Biochemistry and Biophysics University of California, San Francisco MC2240, 600 - 16th Street S414 San Francisco, CA 94158-2517 Phone: +1 415-476-3937 E-mail: bjopp(a)msg.ucsf.edu http://www.msg.ucsf.edu

2 1

error message on dev-1145 during den session
by fulvio saccoccia 08 Sep '12

08 Sep '12

Dear Phenix users, here you have an error occurred during a DEN refinement cycle on a twinned structure. Please, help me in order to overcome the problem. This is the command line I used to run phenix: phenix.den_refine data.mtz model.pdb twin_law=h,-k,-l output.prefix=DEN10 --overwrite main.ncs=true set_b_iso=50 convert_to_isotropic=true xray_data.high_resolution=3.6 nproc=4 secondary_structure_restraints=true optimize_xyz_weight=true optimize_adp_weight=true reference_model.file=3ZVJ-coot1.pdb use_geometry_restraints=true den.params fix_rotamers=true den.params is attached I was running phenix dev-1145 on a 64bit Debian machine. This is the error message: Running DEN torsion optimization on 4 processors... Running torsion simulated annealing ...trying gamma 0.0, weight 100.0 ...trying gamma 0.2, weight 100.0 ...trying gamma 0.0, weight 300.0 ...trying gamma 0.2, weight 300.0 Traceback (most recent call last): File "/usr/local/phenix-dev-1145/build/intel-linux-2.6-x86_64/../../phenix/phenix/command_line/den_refine.py", line 21, in <module> den_refine.run(command_name="phenix.den_refine", args=sys.argv[1:]) File "/usr/local/phenix-dev-1145/phenix/phenix/refinement/den_refine.py", line 113, in run customized_master_params=master_phil) File "/usr/local/phenix-dev-1145/phenix/phenix/refinement/command_line.py", line 146, in run call_back_handler=call_back_handler) File "/usr/local/phenix-dev-1145/phenix/phenix/refinement/driver.py", line 1532, in run ncs_manager = self.ncs_manager) File "/usr/local/phenix-dev-1145/phenix/phenix/refinement/strategies.py", line 830, in refinement_machine log=log) File "/usr/local/phenix-dev-1145/cctbx_project/mmtbx/den/refinement.py", line 88, in __init__ "serious error; please contact bugs(a)phenix-online.org.") % so) RuntimeError: DEN weight optimization failed: DEN cycle 1 rfree at start of SA cycle: 0.4642 k_rep = 1.00 tardy_tree summary: number of vertices: 16570 number of edges: 16855 find cluster loops: 0 repeats number of fixed vertex lists: 0 number of fixed vertices: 0 number of clusters: 7857 merge clusters with multiple connections: 1 pass number of hinge edges: 7857 number of loop edges: 0 number of loop edge bendings: 0 number of fixed hinges: 0 CAUGHT EXCEPTION: (argument #0) ARGUMENT LEADING TO EXCEPTION: (0.0, 100.0) Traceback (most recent call last): File "/usr/local/phenix-dev-1145/cctbx_project/libtbx/easy_mp.py", line 155, in __call__ result = O.func(arg) File "/usr/local/phenix-dev-1145/cctbx_project/mmtbx/den/refinement.py", line 227, in try_den_weight_torsion call_back_after_step=False) File "/usr/local/phenix-dev-1145/cctbx_project/mmtbx/refinement/tardy.py", line 228, in run reduced_geo_manager=reduced_geo_manager) File "/usr/local/phenix-dev-1145/cctbx_project/mmtbx/refinement/tardy.py", line 86, in __init__ O.weights = target_weights.xyz_weights_result AttributeError: 'NoneType' object has no attribute 'xyz_weights_result' This is a serious error; please contact bugs(a)phenix-online.org. Thanks in advance Fulvio Saccoccia PhD student Dept. of Biochemical Sciences Sapienza University of Rome +390649910556

2 6

Fwd: mr_rosetta and Condor
by Dmitry Rodionov 08 Sep '12

08 Sep '12

re: > Submitting job(s) > ERROR: No such directory: /Volumes/condor/Sht1/MR_ROSETTA_5/GROUP_OF_PLACE_MODEL_1/"/Volumes/condor/Sht1/MR_ROSETTA_5/GROUP_OF_PLACE_MODEL_1" 1)It appears that the submission problem comes from the Condor job description file: Universe = vanilla Error = err.log Output = out.log Requirements = Memory > 256 copy_to_spool = False Initialdir = "/Volumes/condor/Och1/MR_ROSETTA_32/GROUP_OF_PLACE_MODEL_1" Executable = "/Volumes/condor/Och1/MR_ROSETTA_32/GROUP_OF_PLACE_MODEL_1/RUN_FILE_10.sh" Queue I am able to submit the job manually if i remove the double quotes from initialdir and executable lines. Sounds like a very simple problem. Does anyone know how can this be fixed? 2) shortly after submission, the modified (no quotes) job stops with the following: RUN_FILE_1.log: Rebuilding the model /Volumes/condor/Och1/MR_ROSETTA_32/WORK_1/global_1/2vk9_1_chain.pdb with rosetta ================================================================================ Splitting generation of rebuilt models into 1000 jobs with 1 structures generated per job ================================================================================ ******************* ERROR ENDING *************** 'NoneType' object has no attribute 'endswith' ******************* ERROR ENDING *************** Oddly enough, if I run the script locally it runs without this problem. I understand that the error is Python. Does anybody know where it might be coming from? 2)Another thing is that Condor's err.log and out.log are created in the same location and under the same name for all jobs which "may lead to bizarre behaviour" according to the Condor's manual. Might not be a bad idea to name the files something like run_xx_err.log or stick them into the RUN_xx folder for the sake of safety. Thanks! Dmitry

1 0

restraining covalent links between symmetry related molecules
by Dave Lodowski 07 Sep '12

07 Sep '12

Dear all, I am attempting to add a bond restraint that connects a serine in one molecule to a carbon located in its symmetry mate in the adjacent unit cell. Atom selection 1: chain 'A' and resid ' 230 ' and name ' C ' and altloc 'A' Atom selection 2: chain 'A' and resid ' 179 ' and name ' OG ' and altloc 'A' and (*insert symmetry keyword here*) 'x-1,y,z' (altloc B in both cases is not covalently attached) I cannot seem to find a complete description of the atom selection syntax in the Phenix documentation to specify which symmetry mate. and the gui/ phenix graphics does not have any way to generate symmetry mates. What is the keyword that specifies a symmetry mate? Thanks, Dave Lodowski

2 1

Autosol: MIRAS
by ash.k＠aol.com 07 Sep '12

07 Sep '12

Hi all, I am tryingto solve a structure through experimental phasing using AUTOSOL. I have acouple of heavy atom derivative datasets (Hg, La, Eu, Cd) and also a SeMetdata. Unfortunately all the datasets are of low resolution (3.7-4.2A) and thereare possibly 4-8 molecules in the asu. MIR, SAD and MAD alone did not give any convincing solution. However,MIRAS, with a combination of few heavy atom datasets and the anomalous datafrom SeMet crystals, gave a somewhat satisfactory map. But the heavy atom sitepicked by AUTOSOL list only one of the heavy atoms i.e. Lanthanum. In another set of run, the solution of which was not convincing, the heavy atom substructure had only Hg. There are 6Met out of 200 residues in one molecule and mass spec results show that Seincorporation is 100%. Now, mydoubt is that why does the heavy atom substructure contain only La and how canI get the substructure involving Se from this solution (or the datasets used)? Se location is going to help me a lot forfinding a starting point to assign side chains. Anysuggestion would be greatly appreciated. Thanks AK

4 8