- phenixbb - phenix-online.org

phenix server upgrade
by Ralf W. Grosse-Kunstleve 01 Dec '09

01 Dec '09

Over the weekend we've moved the phenix web server and mailing list to new hardware. If you notice any problems with the web content, the downloads, or the mailing list, please let us know. Ralf

1 0

Re: [phenixbb] sigma cutoff
by Kay Diederichs 30 Nov '09

30 Nov '09

> > > Hi, > > According to the Phenix manual, phenix.refine uses no sigma cutoff by > default. Since I haven't changed this manually, why do I see in my pdb > file MIN(FOBS/SIGMA_FOBS) 1.340? > > Thanks for help, > Karolina > > This one I can answer because I asked myself the same question recently. The answer is that Phenix does not impose any sigma cutoff, but rather your data are such that the minimum value of Fobs/Sigma_Fobs is 1.34 . You can verify this for yourself by going through your reflection file, line by line, and calculating Fobs/Sigma_Fobs for each reflection. You will find that all your reflections have a Fobs/Sigma_Fobs larger or equal to 1.34 . So it has nothing to with phenix; phenix just reports a property of your data. An easier way than doing it "manually" is using CCP4 sftools for calculating a new column "ratio", writing out the resulting MTZ file, and looking at this with mtzdmp. HTH, Kay -- Kay Diederichs http://strucbio.biologie.uni-konstanz.de email: Kay.Diederichs(a)uni-konstanz.de Tel +49 7531 88 4049 Fax 3183 Fachbereich Biologie, Universitaet Konstanz, Box M647, D-78457 Konstanz

2 1

Re: [phenixbb] Confused by the output from Phenix autobuild jobs
by Thomas C. Terwilliger 30 Nov '09

30 Nov '09

Hi Serah, You are right that the coefficients that you want for a difference map are not present in the overall_best_refine_map_coeffs.mtz file. The coefficients in that file (labelled FP PHIM) are the same as 2FOFCWT PH2FOFCWT from the refinement output, but there are no difference map coefficients. I will change this. In the meantime, I would do this after the end of autobuild: phenix.refine AutoBuild_run_1_/overall_best.pdb \ AutoBuild_run_1_/exptl_fobs_phases_freeR_flags.mtz and then you will have your map coefficients. >> I'm still wondering about the overall_best_*.* outputs - are these the >> really final output, or do we go for cycle_best_3.* series? The overall_best_xxx files are copied from the cycle_best_xxx files which are copied from the original files, all of which are in the TEMP0/ subdirectory. Thanks for pointing out that the difference coefficients are not present in the autobuild map file. I will fix the output to have the original refinement map file available. All the best, Tom T

1 1

alternate conformation
by subhash bihani 25 Nov '09

25 Nov '09

Hi, I am refining a structure with ligand in alternate conformations. I tried giving them A and B, it worked. but I feel that this doesnot make interactions between alternate conformations absolutely zero. Is there any command in phenix similar to Igroup statement in CNS which can be used for such cases. Regards subhash -- Subhash Bihani Scientific Officer Solid State Physics Division Bhabha Atomic Research Centre Trombay, Mumbai-85 Ph. 25594688 (o) alternate email: subhashbihani(a)yahoo.co.in

5 10

R-free rules
by Bryan W. Lepore 25 Nov '09

25 Nov '09

in the phenix.reflection_file_editor gui not sure i understand this: [x] Extend existing R-free array(s) ... [R-free flags generation] <--+ | | | | +-if this is checked, then do NOT check "Generate R-free flags" in here ... because "extend" is not the same as "generate", right? ... and thanks for the help before. -bryan

12 35

Age Concern edition of the Computing Commission and Teaching Commission newsletter
by Ralf W. Grosse-Kunstleve 24 Nov '09

24 Nov '09

Posted on behalf of Lachlan Cranswick: The Computing Commission and Teaching Commission has produced a joint November 2009 newsletter on the theme of "Age Concern", which deals with issue concerning "the slow march towards retirement of the major generation of crystallographic programmers", as further explained by David Watkin within the newsletter. http://www.iucr.org/resources/commissions/crystallographic-computing/newsle… or via http://www.iucr.org/resources/commissions/crystallographic-teaching/newslet… --------------------------- This issue's articles include : - the latest background and information from the UK based EPSRC funded Age Concern Project. - history and science within the Caltech based CRYRM single crystal suite; DIMS - Direct methods In Multidimensional Space software; Crunch 1.5 Direct Methods Software; the Glasgow based GX single crystal suite; the Italian based CAOS structure refinement software; Ton Spek's Platon; and a list of programs archived at Armel Le Bail's Crystallography Source Code Museum. - experience in conversion of large crystallographic Fortran-77 to C++; as well as the place of the Quick and Dirty Crystallographic Computer Program - Historical group photographs from some of the Crystallographic Computing and Teaching School of the 1970's. - Archived manuals and source code from the above mentioned software. - Teaching material of Isabella L. Karle, as she presented at the 1978 Erice School on Direct Methods for Solving Crystal Structures - A reprint of the The Analytical Theory of Point Systems (1923) by J. D. Bernal (1901-1971) Further Computing Commission newsletters are intended on theme of Age Concern : capturing the history and the crystallographic science behind past and present Small Molecule, Powder Diffraction, Protein and other crystallographic software. Submissions on this topic are welcome and encouraged. - Lachlan Cranswick ( lachlan.cranswick(a)nrc.gc.ca )

1 0

Confused by the output from Phenix autobuild jobs
by Serah Wangari 24 Nov '09

24 Nov '09

Dear Phenix users, I have run a number of Autobuild jobs using the Phenix GUI and in all cases, I'm not sure what output I should proceed with to the next step. - The first job was for rebuilding a loop section in the model (rebuild-in-place) - As far as I can remember, there was no mtz generated from this job, I had to generate map coefficients with the overall_best pdb and an old mtz using another refinement program. Did I miss something during the job input step, so that no mtz was generated, or how would I proceed in such a case? - The next job involved rebuilding the whole model (rebuild-in-place) - in this case, I took the overall_best pdb, but when I tried to open the overall_best_refine_map_coeffs.mtz in coot, I got this error that coot failed to find F/Phi columns in the mtz file. Then I checked in the TEMPO directory, and I managed to open a random mtz from one of the refinement cycles with coot, but I wasn't confident building into that, so again I had to generate map coefficients using overall_best pdb and an old mtz in another program. - The other job was rebuilding (rebuild-in-place) another model with composite maps generation - Here, in addition to the 12 mtz and pdb files from the 12 omit regions, I got resolve_composite_map.mtz I should assume as the only final output; and no final pdb. Again, this mtz file doesn't have the F/Phi columns as I got the same error as above in coot. I'm not sure how to proceed with this mtz...... I would really appreciate if someone could point me to the right direction or at least direct me to some literature where the output from Phenix autobuild jobs is explained - as I haven't come across such info yet. Thank you in advance, Serah --

3 2

Clean Jobs
by tea.pavkov＠uni-graz.at 24 Nov '09

24 Nov '09

Dear all, I would have a general question - how can I clean the Jobs (in the Job history list and all corresponding directories) in the new GUI? Thank you in advance! Best regards, Tea ----------------------- Dr. Tea Pavkov-Keller KF Uni Graz, ZMB Structural Biology Group Humboldtstr. 50/III 8010 Graz Austria

3 3

Problem installing on Leopard PowerPC G5
by Satpal Virdee 23 Nov '09

23 Nov '09

Hi, I am trying to install the latest 1.5-2 official release on an iMac PowerPC G5 running OSX 10.5.8. During installation I get the error: " Relocating Python framework. . . /usr/bin/install_name_tool: object: libwx_mac-2.8.0.5.0.dylib malformed object (unknown load command 20) Error: see /Users/Sat/Downloads/phenix-installer-1.5-2/build-binary/ mac-ppc-osx/cpemac406-18/log/binary.log for details " Inspection of the binary.log files reveals this: " libfreetype.6.3.10.dylib: /usr/local/phenix-1.5-2/build/mac-ppc-osx/ base/lib/libfreetype.6.3.10.dylib --> /usr/local/phenix-1.5-2/build/mac-ppc-osx/base/lib/libfreetype. 6.3.10.dylib libfreetype.6.dylib libfreetype.6.dylib: /usr/local/phenix-1.5-2/build/mac-ppc-osx/base/ lib/libfreetype.6.dylib --> /usr/local/phenix-1.5-2/build/mac-ppc-osx/base/lib/libfreetype. 6.dylib libfreetype.dylib libfreetype.dylib: /usr/local/phenix-1.5-2/build/mac-ppc-osx/base/lib/ libfreetype.dylib --> /usr/local/phenix-1.5-2/build/mac-ppc-osx/base/lib/libfreetype.dylib libwx_mac-2.8.0.5.0.dylib fix_mac_python: caught error - exiting. " Any help would be most appreciated. Regards, Satpal

2 1

baton build in coot
by Raja Dey 22 Nov '09

22 Nov '09

Dear Friends, How does CA ZONE > MAIN CHAIN under 'other modeling tools' in coot work? I was trying to rebuild main chain atoms for 3 residues(8,9,10) using CA ZONE > MAIN CHAIN. when I click on one atom in residue 8 and then one atom in residue 10, it failed with a message that says the zone is too short. I increased the zone to 5 residues, still same failed message. Also, is it possible to move an atom in coot? Thanking you in advance for your reply. Raja

2 1

phasing power value in Autosol/Autobuild
by vennila Natesan 19 Nov '09

19 Nov '09

Dear all... I need toknow the phasing power value for my solved structure. Where can i find it out in the output files of autosol and autobuild? I am using phenix1.14.4b and also i need to know the weighted mean phase errors vs resolution. Thanking you.

3 2

I vs F data input
by gerwald jogl 18 Nov '09

18 Nov '09

Hi All, A quick question about diffraction data input. I just noticed that phenix.refine uses Imean/SIGImean and converts these by default in my refinements. Is this the preferred way or would it be better to use F/SIGF (after CCP4)? Thanks, Gerwald

2 1

external ncs reference structures?
by Cathy Lawson 18 Nov '09

18 Nov '09

Dear Phenix developers, We are finishing up refinement of some multiprotein-DNA complex structures at 3-3.5 A resolution and I am wondering if it would be possible with phenix.refine to supply one or more external sets of coordinates from higher-resolution determinations to serve as the ncs references. If so what would be the best way to set this up? thanks, Cathy ********************************** Catherine Lawson, Ph.D. Rutgers University Chemistry & Chemical Biology 610 Taylor Rd Piscataway, NJ 08854 732 445 5494 cathy.lawson(a)rutgers.edu http://roma.rutgers.edu **********************************

2 1

Mac OSX Phenix GUI
by Miguel Ortiz Lombardia 18 Nov '09

18 Nov '09

Hi, This is a question/request to the phenix developers. Are there any plans to make the new phenix GUI to use an X11 python interface instead of the Aqua one? The reason I ask is that we use Mac OSX servers for calculations here and an X11 interface would be a lot more appropriate for this configuration. That way users could export the GUI to their local workstations. As it is implemented now this use of the GUI is impossible on Mac OSX. Before someone makes the suggestion: Remote Desktop simply allows to share the only Aqua desktop that the server can produce, that is, only one user can access the Aqua desktop, whether remotely or locally. We would be happy to test any X11 version you might be willing to develop. Thanks for your great software, which we can use anyway from the command line :-) Best regards, -- Miguel Architecture et Fonction des Macromolécules Biologiques (UMR6098) CNRS, Universités d'Aix-Marseille I & II Case 932, 163 Avenue de Luminy, 13288 Marseille cedex 9, France Tel: +33(0) 491 82 55 93 Fax: +33(0) 491 26 67 20 e-mail: miguel.ortiz-lombardia(a)afmb.univ-mrs.fr Web: http://www.pangea.org/mol/spip.php?rubrique2

2 1

Phaser output and error messages
by Ian Stokes-Rees 17 Nov '09

17 Nov '09

I'm having some discussion with a colleague about phaser output (we're using Phaser 2.1.4). We haven't been able to find any documentation which can clarify our situation, and I'm hoping someone on the list can help answer these questions. I should mention that I am relatively new to Phaser. 1. PHASER.sol files: Which "SOLU SET" does Phaser consider to be the best? The first or the last? Or the one with the highest LLG, wherever that may be? In our experience of running Phaser over several MTZ files with a range of models the best Phaser solution has always been the first, and this has had the highest LLG. Note: this is with "untuned" Phaser settings for identity, solvent fraction, or number of search models in ASU -- our goal is to do a first run with "generic" settings for these over a larger set of models, then (from TFZ and LLG scores) select a subset for which we will tune Phaser parameters and PDB search model variations. 2. If we are right that the first "SOLU SET" entry is indicative of the potential for the search model to form a good MR candidate, then is it the case that the first entry is the first Phaser solution that is computed? Or is the PHASER.sol file a sorted list output at the end of the run? From my reading of the documentation it is output in order of computation, and *for our purposes* (if my first statement in this question is correct) Phaser can stop after it outputs this first solution. Is there some way to tell Phaser to stop after the first solution is output? I realize that this doesn't sound like it makes sense (how could Phaser know to pick the best solution first, and even if it could, why would it ever continue past this point), however I ask because we have put a 10 minute timeout into our Phaser runs and we have many situations where we get a timeout but PHASER.sol has already been generated and the best LLG solutions are output first. It leaves me wondering why it didn't just stop on its own after outputting the first result instead of being aborted by our (external) timeout that terminates the process? 3. PHASER.sol files: For single domain search models, we usually get output of the form: SOLU SET RFZ=4.5 TFZ=5.2 PAK=0 LLG=14 LLG=14 SOLU 6DIM ENSE model1 EULER 242.049 45.040 326.088 FRAC -0.09425 0.50268 0.42575 however we see three variations: i) No LLG: SOLU SET RFZ=3.1 TFZ=5.0 PAK=0 SOLU 6DIM ENSE model2 EULER 59.983 69.335 319.701 FRAC -1.17131 -0.70030 0.23150 ii) One LLG: SOLU SET RFZ=3.7 TFZ=4.6 PAK=0 LLG=25 SOLU 6DIM ENSE model3 EULER 293.943 128.068 332.147 FRAC 0.06273 0.13175 0.25054 iii) Two LLG entries, but with different values: SOLU SET RFZ=3.8 TFZ=4.1 PAK=0 LLG=21 LLG=20 SOLU 6DIM ENSE model4 EULER 278.058 129.347 33.292 FRAC 0.28446 0.29011 -0.07986 4. Occasionally we get an error that we don't understand: FATAL RUNTIME ERROR: No scattering in pdbfile model1.pdb What does this mean? Is there a problem with the PDB file? We can't see anything obvious in the ones which produce this error. Thanks, Ian -- Ian Stokes-Rees, Research Associate SBGrid, Harvard Medical School http://sbgrid.org

2 2

Request: read definitions from stdin
by Joe Krahn 16 Nov '09

16 Nov '09

I want to be able to run PHENIX in a shell script, where stdin is accepted as a definition input file. It could be simply specified as the usual POSIX syntax where "-" defines stdin in place of an input file. If the front-end is based too heavily on file suffix mappings, you could add a flag "--stdin=SUFFIX" to define reading stdin as if it had the given SUFFIX. It would also be nice to avoid the duplication of output to both stdout and the .log file. Output can always be sent to /dev/null, but it would be nicer to have either a --quiet flag to suppress stdout, or a DEF option to disable the .log output, or both options to let the user decide. Most other computational programs can be run this way, so I assume that other people would find this useful, especially for running PHENIX from scripts. Hopefully, it is also a fairly easy feature to implement. Thanks, Joe Krahn

3 2

R-free rules
by Bryan W. Lepore 13 Nov '09

13 Nov '09

[ phenix.refine v1.5 rel. tag 2 intel linux 2.6 ] if i provide the test set .mtz (flag = "1") in refinement separately from the data .mtz, it looks like reflection_file_converter.py does not flag the work set as "0", if the work set lacks such a flag -- is that right? i.e. when providing the test set, it looks like all observed data must be in the test set, however the work set has flag "0". -bryan (some error output below) Sorry: R-free flags not compatible with F-obs array: missing flag for 22206 F-obs selected for refinement:

3 4

November 15, 2009 deadline- User proposal submission for Collaborative Crystallography at BCSB
by Banumathi Sankaran 12 Nov '09

12 Nov '09

Dear Users, The deadline for October/November/December 2009 Collaborative Crystallography proposals will be *November 15, 2009. * ------------------------------------------------------------------------ Through the Collaborative Crystallography Program (CC) at the Advanced Light Source (ALS), scientists can send protein crystals to Berkeley Center for Structural Biology (BCSB) staff researchers for data collection and analysis. The CC Program can provide a number of benefits to researchers: * Obtain high quality data and analysis through collaborating with expert beamline researchers; * Rapid turn around on projects; and * Reduced travel costs. To apply, please submit a proposal through the ALS General User proposal review process for beamtime allocation. Proposals are reviewed and ranked by the Proposal Study Panel, and beamtime is allocated accordingly. BCSB staff schedule the CC projects on Beamlines 5.0.1 and 5.0.2 to fit into the available resources. Only non-proprietary projects will be accepted. As a condition of participation, BCSB staff researchers who participate in data collection and/or analysis must be appropriately acknowledged - typically being included as authors on publications and in PDB depositions. Please consult the website for additional information at: http://bcsb.als.lbl.gov/wiki/index.php/Collaborative_Crystallography ----------------------------------------------------------------------------------------- How To Apply: To learn more, go to: http://www-als.lbl.gov/als/quickguide/independinvest.html To submit a proposal, go to: http://alsusweb.lbl.gov/. Scroll down to "*Structural Biology beamlines (includes protein SAXS)*" and click on "New Proposal." Enter your proposal information, with attention to the following details: * For question 3/First choice, select "5.0.1 (Monochromatic);" for question 3/Second choice, select "5.0.2 (MAD)." * Check box (yes) in response to question "Do you want collaborative crystallography (beamline 5.0.1 or 5.0.2 only) * In question 9, please describe a specific research project with a clear end point. In order to request CC time for October/November/December 2009 allocation period, proposals must be submitted by *November 15, 2009.* The deadline for CC proposals for the time period March/April 2010 will be January 15, 2010. Regards, Banumathi Sankaran

1 1

sulfate on a 2 fold axis
by David Cobessi 12 Nov '09

12 Nov '09

Hi, I am finishing to refine a structure using Phenix. A sulfate ion is located on a 2 fold axis. When I refine the structure with Refmac, the geometry of the ion is well conserved, while when I use Phenix the geometry is destroyed and the sulfate ion is like a cross!! Any idea to keep the good geometry of this ion in Phenix during the refinement? Thanks in advance, David -- David Cobessi Institut de Biologie Structurale 41, Rue Jules Horowitz 38027 Grenoble Cedex-1, France Tel:33(0)438789613 33(0)608164340 Fax:33(0)438785122

2 1

phenix.cif_as_mtz
by Bryan W. Lepore 12 Nov '09

12 Nov '09

[ v. 1.5 release tag 2 intex-linux-2.6 redhat-e4.8 ] `phenix.cif_as_mtz file.cif --symmetry=FILE.pdb --show-log --show-details-if-error` .. reports "No data extracted from input cif file.", nothing else. command works w/ other .cif/.pdb's, CRYST1 lines, etc. look comparable between them,... not sure where to look next. thanks for any suggestsions. -bryan

2 1

Nino Campobasso is out of the office.
by Nino.2.Campobasso＠gsk.com 11 Nov '09

11 Nov '09

I will be out of the office starting 11-Nov-2009 and will not return until 12-Nov-2009.

1 0

Failed to make map after phenix.refine
by Zhou, Tongqing (NIH/VRC) [E] 11 Nov '09

11 Nov '09

Hi, All I am having problem with the new 1.5.2. After few rounds of refinement, the program failed to write out maps. It ended with: "Sorry: Duplicate mtz_label_amplitudes:FOFCWT" Here is the portion in my param file which worked fine before: electron_density_maps { map_format = *xplor map_coefficients_format = mtz *phs suppress = None map { mtz_label_amplitudes = "2FOFCWT" mtz_label_phases = "PH2FOFCWT" likelihood_weighted = True obs_factor = 2 calc_factor = 1 } map { mtz_label_amplitudes = "FOFCWT" mtz_label_phases = "PHFOFCWT" likelihood_weighted = True obs_factor = 1 calc_factor = 1 } } Thanks in advance for your help. Tongqing

2 1

C-ter amidation of peptide
by Nate Nicely 10 Nov '09

10 Nov '09

Dear all, How best to handle a C-terminally amidated peptide in phenix? I defined the atom as 'N' and the resname as 'NH2' following a structure I found that also bears a C-terminally amidated peptide. However phenix.refine seems not to be making the linkage between the amide and the preceding residue. Putting the structure into phenix.elbow returns, "No unknown residues." Thanks, Nate N.

3 2

Glycan refinement
by Zhou, Tongqing (NIH/VRC) [E] 10 Nov '09

10 Nov '09

Dear All, I am refining a structure with a N-linked glycosylation site. There is density that goes beyond the first NAG, at least to MAN3. I know in Phenix there is data_link = NAG-ASN for the ASN-NAG connection, but I don't know how to do the "apply_cif_link' for the extensions such NAG-NAG, NAG-MAN, MAN-MAN...., any suggestions are appreciated. Thanks, Tongqing Tongqing Zhou, Ph.D. Staff Scientist Structural Biology Section Vaccine Research Center, NIAID/NIH Building 40, Room 4607B 40 Convent Drive, MSC3027 Bethesda, MD 20892 (301) 594-8710 (Tel) (301) 793-0794 (Cell) (301) 480-2658 (Fax) ****************************************************************** The information in this e-mail and any of its attachments is confidential and may contain sensitive information. It should not be used by anyone who is not the original intended recipient. If you have received this e-mail in error please inform the sender and delete it from your mailbox or any other storage devices. National Institute of Allergy and Infectious Diseases shall not accept liability for any statements made that are sender's own and not expressly made on behalf of the NIAID by one of its representatives. ******************************************************************

5 9

Problem reading ".edits" parameter file
by John Pascal 06 Nov '09

06 Nov '09

Hello, On a previous version of PHENIX, 1.4-153, I was able to read in a parameter file with geometry constraints. Upon loading the file, the program asked when to execute the commands, I said at run time, and everything worked fine. Using the same parameter file in PHENIX 1.5-2, I am no longer asked about when to execute the file, and when I proceed to run the refinement, the program quits stating an unexpected character. I can't find anything obviously wrong with the parameter file. Any suggestions? The previous version of PHENIX (the one that worked) is no longer located on my server, so I can't easily return to that version. Thanks. Best Regards, John -- John Pascal, PhD ph 215.503.4596 Assistant Professor fx 215.923.2117 Thomas Jefferson University Biochemistry & Molecular Biology 233 South 10th Street, BLSB 804 Philadelphia, Pennsylvania 19107 ---------------------------------------------------------------- This message was sent using IMP, the Internet Messaging Program.

2 1