- phenixbb - phenix-online.org

Help needs: SAD-phasing for Cobalt ion
by YoungJin 10 Dec '10

10 Dec '10

Dear helpers in phenix, I have anomalous data (Co; wavelength: 1.609A) that I want to use its Phase information from SAD-phasing method (autosol). As I have a native data (from nomal wavelength ~1A) with higher resolution (~2A) than SAD data (~2.5A), I wanted to apply phase information to this native data for further refinement: My goals were 1. I want to confirm where is my heavy metal is; 2. I want to refine further using native data as has higher resolution. At first, when I did autosol, it didn't give me perfect result. So far what I tired is as following: 'SAD peak' from imosflm generated mtz file (scaled), 'Phases' from first refmac generated mtz file from the native data, and 'MIR native' from the native data mtz file (imosflm and scaled). FOM: 0.186; overall score was 10.75 +/- 12.33. FYI. I input 2 sites and result suggested 3 sites. Each molecule has one binding site of metal ion. I was using old Mac version:1.6.4-486 Anyway, I simply passed the result over Autobuild. It ran and gave me better statistics as following: R/Rfree=0.2734/0.3312; CC 0.76. However (I have two molecules/ASU) one of models was not great. But this time, anomalous density map well phased(?!) to 2Fo-Fc map and I could see - it was before only one heavy metal I clearly could see- 2 heavy metal signal from 2 molecules and those positions are well matched with native density maps. So I continued phenix.refine using Autobuild made 'expt_fobs_phase_freeR_flags_1.mtz'. The question is the statistics from anomalous data itself (it has slight dispcrepancy between anomalous map and 2Fo-Fc/Fo-Fc map due to phase problem) or even those from native data is better than my current model which came from native data set (SAD-phased as described above) with higher resolution. I think one of main reasons is one of molecule has a lid that does not perfectly fit into current density map. However, in terms of heavy metal position and the active site is better from current map (Autobuild generated model). FYI: current model has R/Rfree=0.21/0.27(best was 0.21/0.26); 0.19/0.24 (for native data only used model) .... Am I too concerning for statistics or is there something I am missing? Now my questions are as following plus any suggestions are welcome: 1. For SAD-phasing, could you tell me which is best way to do so with my case (one data from the edge of heavy metal (lower resolution) and native data (higher resolution) or where is the point I can convince. Also I tried to use DANO/SIGDANO column (from hkl2000 and ccp4) but seems like it (Autosol) doesn't take them into account. Could you tell me how can I improve phase by using anomalous data? 2. In case of Co, does anyone have experience about coordination numbers? I thought it should have 6 (or 4) coordination numbers. The question what I have now is about the numbers and their lengths ... I though it should be within 2A, but not true for my current models (And somehow distorted octahedron or else with 2.07-3.28A). I also searched PDB files, and found that many deviations from 1.98 to 2.82A (depending on corresponding atoms and resolutions: 3KOI, 3KON, 2DD5, and 2XJM). I would appreciate if you can give me any tips regarding this issue. Many thanks in advance, Young-Jin

2 1

phenix.phaser ep_sad mode
by Bryan Lepore 09 Dec '10

09 Dec '10

dev-601 i found an ep_sad script on the phaser wiki^[1] that works ok^[2] for old versions. for the dev-601 phaser, i think it comes down to me not understanding this line : LLGComplete CRYStal insulin COMPLETE ON SCATtering ELEMent S i tried jumbling the commands/arguments around^[3], leading to syntax errors for each attempt (output excerpt below). if this "should work", i'll gladly return to the manual - or is there something unusual between versions, am i missing something, etc.... -Bryan [1]http://www.phaser.cimr.cam.ac.uk/index.php/Experimental_Phasing#Combined_MR_and_SAD_Phasing [2] seem not to write ATOM records in the .pdb when sites are in fact found. [3]http://www.phaser.cimr.cam.ac.uk/index.php/Keywords#LLGCOMPLETE output excerpt : MODE EP_SAD TITLe sad phasing insulin with best AutoMR solution HKLIn truncate.mtz CRYStal insulin DATAset sad LABIn F+=F(+) SIG+=SIGF(+) F-=F(-) SIG-=SIGF(-) ! LLGComplete CRYStal insulin COMPLETE ON SCATtering ELEMent S LLGComplete CRYStal  <B><FONT COLOR="#FF8800"> $TEXT:Warning: $$ Baubles Markup $$ ------------------------------------------------------------------------------------- SYNTAX ERROR: Use (optionally) COMPLETE TYPE SCATTERER ELEMENT CLUSTER SIGMA NCYC CLASH ANOM REAL MAPS or METHOD -------------------------------------------------------------------------------------

3 5

phaser ep_sad zscore
by Bryan Lepore 09 Dec '10

09 Dec '10

dev-601 in phaser ep_sad mode, i tried ZSCORE SOLVED 0 to get more sites output, and i can see a peak of about Z=4 in the stdout, but no peaks are written to the .pdb am i applying ZSCORE wrong here? -Bryan

3 3

reflection outliers detection in phenix.refine
by guanrj＠foxmail.com 08 Dec '10

08 Dec '10

Hi All I have a question regarding reflection outliers detection in phenix.refine: in my input reflection file, there are 24432 reflections, but only 22769 of them were used in refinement. I checked the log file and saw 1663 outliers. (below if from the log) basic_wilson_outliers = 14 extreme_wilson_outliers = 18 beamstop_shadow_outliers = 0 model_based_outliers = 1657 total = 1663 What is the "model_based_outliers"? What causes so many reflections rejected and not used in the refinement? Thanks Rongjin

3 3

Multistart simulated annealing refinement
by Joseph Noel 08 Dec '10

08 Dec '10

Hi All, I'm sure this has already been addressed but is there a way through the PHENIX GUI to carry out multistart simulated annealing refinement either with different starting seeds (velocities) or temperatures automatically or does one have to start a fresh run each time? Thanks! Joe P.S. Similar procedure as described in the PNAS paper: "Multistart simulated annealing refinement of the crystal structure of the 70S ribosome" ___________________________________________________________ Joseph P. Noel, Ph.D. Investigator, Howard Hughes Medical Institute Professor, The Jack H. Skirball Center for Chemical Biology and Proteomics The Salk Institute for Biological Studies 10010 North Torrey Pines Road La Jolla, CA 92037 USA Phone: (858) 453-4100 extension 1442 Cell: (858) 349-4700 Fax: (858) 597-0855 E-mail: noel(a)salk.edu Web Site (Salk): http://www.salk.edu/faculty/faculty_details.php?id=37 Web Site (HHMI): http://hhmi.org/research/investigators/noel.html ___________________________________________________________

3 2

phenix.fmodel or other utility: calculate HL's as well?
by Francis E Reyes 08 Dec '10

08 Dec '10

I love the phenix cmd line utilities, however is there one that'll calculate HLA/HLB/HLC/HLD given FModel/Phimodel? My expectation is that phenix.fmodel would do this with an option toggle or something. F --------------------------------------------- Francis E. Reyes M.Sc. 215 UCB University of Colorado at Boulder gpg --keyserver pgp.mit.edu --recv-keys 67BA8D5D 8AE2 F2F4 90F7 9640 28BC 686F 78FD 6669 67BA 8D5D

2 1

No clash evaluation between symmetry mates in phenix.refine
by Keitaro Yamashita 08 Dec '10

08 Dec '10

Dear Phenix developers, I have a question about phenix.refine. I'm refining complex structure against X-ray diffraction data with packing disorder. (Some domains overlap with their symmetry mates (4-fold), so their occupancies are set to 0.25) So I want phenix.refine to ignore clashes between symmetry mates. Can phenix.refine exclude any interaction between symmetry mates from geometric term in target function: E = E(x-ray) + E(geometry)? Is there any option? Thank you very much in advance, K. Yamashita

3 4

switch off anti-bump restrain
by Agnieszka 08 Dec '10

08 Dec '10

Hi, I have double conformation of my residue. There is a water molecule that interacts with conformation A. In the same time the water is very close to conformation B. During the refinement the water is pushed out. I try to switch off the site refinement of the water and some atoms of conformation b but it didn't work out for me. Do you know how can I switch off the anti-bump restrain of this two guys? Thank you Agnieszka Kiliszek

2 1

ncs restraints
by Timothy Springer 08 Dec '10

08 Dec '10

Hi, I am using NCS restraints for 4 copies of a molecule at 3.05 Å. We are using the standard weights: coordinate_sigma = 0.05 b_factor_weight = 10 and the 4 molecules look identical. What are suggested values for looser restraints? One of my colleagues said that varying them did not make a difference- ncs was still very strict. Is there a way to make them less strict? Thanks, Timothy A. Springer, Ph.D.

4 4

phaser version
by Bryan Lepore 08 Dec '10

08 Dec '10

using dev-596 from source, i see `phenix.phaser --version` is Phaser-2.2.4 (svn ) that's true for dev-601 (Intel binary) also. the phaser wiki says 2.3 is the latest, and to compile it from the phenix distribution. so... wondering how to get phaser 2.3 (before asking questions about ep_sad...) -Bryan

4 3

anyone with a scalepack file with R3 or R32 using rhombohedral axes?
by Ralf W. Grosse-Kunstleve 08 Dec '10

08 Dec '10

I'm looking for a little help from the phenix user community. Does anyone have a scalepack output file with data in space group R3 or R32 using rhombohedral axes, i.e. unit cell parameters a=b=c, alpha=beta=gamma? What is the space group symbol in the file? I'm only interested to know this for files that are the direct output of scalepack. I'd like to verify that this part of the scalepack documentation is actually true: http://www.hkl-xray.com/hkl_web1/hkl/manual_online.pdf Page 139: 146 R3 R3 in hexagonal setting 446 H3 R3 in primitive setting 155 R32 R32 in hexagonal setting 455 H32 R32 in primitive setting If this is actually true, it is the exact opposite of what's used by the PDB. Thanks in advance for feedback. Ralf

2 1

modelling of a novel protein
by intekhab alam 07 Dec '10

07 Dec '10

Hi I am trying to solve the structure of RNA polymerase in complex with its protein primer. The size of protein primer is 70 amino acids which i complexed with the polymerase and collected a data of the complex crysta at 2.3A. The space group is C2 and after Molecular replacemnet with full length polymerase in auto MR i got 3 subunits in ASU. The parameters of MR solution RFZ=71.5 TFZ=41.2 PAK=0 LLG=8812 suggests that the solution is correct. I used autobuild for model building which gave me the model of the polymerase as well as some extra density that i am considering as the other protein. After several round of refinemnt i have R and R free of 27and 32 respectively. My queries: 1. I want to know ,how can i improve my R as well as R free. 2. How can i model other protein as the density is weak and not complete, how can i improve the map of the other protein so that i can b ebale to trace the c-alpha backbone, i was able to model only 10 amino acids so far. Please guide. Regards -- INTEKHAB ALAM LABORATORY OF STRUCTURAL BIOINFORMATICS KOREA UNIVERSITY, SEOUL

2 1

distorted (collapsed) zinc-finger motifs
by Hasan Demirci 07 Dec '10

07 Dec '10

Hi, I am in the middle of refining a 3.5 A (low resolution) ribonucleo-protein complex structure. Two of the proteins in this complex contain "zinc-finger" motif. These zinc-finger motifs instead of forming a cysteine-zinc interactions, cysteines are forming di or in these case tetra-sulfide bridges. Initially, zinc is centered but after refinement 4 cysteines come out collapsed and forming covalent bonds to each other and also distorted and zinc is not at the center anymore. Below is the phenix.refine zn_link.param file I am using. Many thanks for advice, Hasan Demirci ---------------------------------------------------------------------------- refinement.pdb_interpretation { disulfide_distance_cutoff = 0.1 } refinement.pdb_interpretation.apply_cif_link { data_link = ZN-CYS residue_selection_1 = chain D and resname ZN and resid 306 residue_selection_2 = chain D and resname CYS and resid 9 } refinement.pdb_interpretation.apply_cif_link { data_link = ZN-CYS residue_selection_1 = chain D and resname ZN and resid 306 residue_selection_2 = chain D and resname CYS and resid 12 } refinement.pdb_interpretation.apply_cif_link { data_link = ZN-CYS residue_selection_1 = chain D and resname ZN and resid 306 residue_selection_2 = chain D and resname CYS and resid 26 } refinement.pdb_interpretation.apply_cif_link { data_link = ZN-CYS residue_selection_1 = chain D and resname ZN and resid 306 residue_selection_2 = chain D and resname CYS and resid 31 } refinement.pdb_interpretation.apply_cif_link { data_link = ZN-CYS residue_selection_1 = chain N and resname ZN and resid 307 residue_selection_2 = chain N and resname CYS and resid 24 } refinement.pdb_interpretation.apply_cif_link { data_link = ZN-CYS residue_selection_1 = chain N and resname ZN and resid 307 residue_selection_2 = chain N and resname CYS and resid 27 } refinement.pdb_interpretation.apply_cif_link { data_link = ZN-CYS residue_selection_1 = chain N and resname ZN and resid 307 residue_selection_2 = chain N and resname CYS and resid 40 } refinement.pdb_interpretation.apply_cif_link { data_link = ZN-CYS residue_selection_1 = chain N and resname ZN and resid 307 residue_selection_2 = chain N and resname CYS and resid 43 } # data_link_ZN-CYS # loop_ _chem_link_bond.link_id _chem_link_bond.atom_1_comp_id _chem_link_bond.atom_id_1 _chem_link_bond.atom_2_comp_id _chem_link_bond.atom_id_2 _chem_link_bond.type _chem_link_bond.value_dist _chem_link_bond.value_dist_esd ZN-CYS 1 ZN 2 SG single 2.340 0.020 loop_ _chem_link_angle.link_id _chem_link_angle.atom_1_comp_id _chem_link_angle.atom_id_1 _chem_link_angle.atom_2_comp_id _chem_link_angle.atom_id_2 _chem_link_angle.atom_3_comp_id _chem_link_angle.atom_id_3 _chem_link_angle.value_angle _chem_link_angle.value_angle_esd ZN-CYS 1 ZN 2 SG 2 CB 109.000 3.000 ---------------------------------------------------------------------------------------------------- -- Hasan DeMirci, Ph.D. Postdoctoral Fellow Department of Molecular Biology, Cellular Biology & Biochemistry Brown University 185 Meeting Street Providence, RI 02912 (401) 863-3652 lab (SFH) (401) 863-6124 lab (ship st) (401) 226-7852 cell Hasan_DeMirci(a)Brown.edu demircha99(a)gmail.com

2 2

Re: [phenixbb] Using Reference Structure in Phenix Refine
by John Pascal 05 Dec '10

05 Dec '10

Got it! Thank you, Pavel and Jeff. Best, 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.

1 0

Using Reference Structure in Phenix Refine
by John Pascal 04 Dec '10

04 Dec '10

Hello Everyone, I am refining against low resolution data, and I would like to use a higher resolution model that is available as a reference structure. There are multiple copies of the molecule in the asymmetric unit, each with a different chain ID. Does the reference PDB file need to contain multiple copies, such that there is a match for each chain (e.g. RESIDs and SEGIDs)? Or does phenix figure out that the reference structure refers to several copies of the molecule, and you only need a single copy in the reference molecule PDB file? Or will NCS restraints take care of this? Thanks. -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.

3 2

Can PHENIX write the peak list to a PDB file?
by Joe Krahn 03 Dec '10

03 Dec '10

PHENIX has an option to find positive and negative peaks, but only lists them in the log file along with the residue they are closest to. Is it possible to write out the actual peak coordinates, such as to a PDB file? Thanks, Joe Krahn

5 7

heavy atom search within phaser molecular replacement
by Bryan Lepore 03 Dec '10

03 Dec '10

i know that phaser will find heavy atoms in autosol, and that "MR-SAD" can be run as two separate steps, such that with some scripting, one could in principle cycle through a bunch of heavy atom searches using MR phases. however, it is not clear to me whether phaser can be run such that e.g. the refined MR solutions can have their anomalous Fourier's searched for peaks in a single MR job. seems possible in principle - before i go trying this, any comments would be appreciated.... i.e. run heavy atom searches as part of molecular replacement, for convenience. hope that makes sense. -Bryan

1 0

Issue with phenix.superpose_maps
by Stephen Graham 03 Dec '10

03 Dec '10

Hi all, I'm trying to use phenix.superpose_maps through the GUI to get superposed maps covering the same region of NCS-related chains A and B for a figure. I have used PDB Tools to split the model into two files with only chain A or chain B. However, when I run superpose_maps specifying the two individual chains for the input PDB files but the same mtz file for the map coefficients I only get one output map (I think it overwrites the output file...). Selected input from the generated .eff file is as follows: input_files { map_1 = "/home/stepheng/blah/Native_refine_6_map_coeffs.mtz" labels_1 = 2FOFCWT,PH2FOFCWT map_2 = "/home/stepheng/blah/Native_refine_6_map_coeffs.mtz" labels_2 = 2FOFCWT,PH2FOFCWT pdb_1 = "/home/stepheng/blah/Native_refine_6_chainA.pdb" d_min_1 = None pdb_2 = "/home/stepheng/blah/Native_refine_6_chainB.pdb" d_min_2 = None output_dir = "/home/stepheng/blah" } And the end of the output is as follows: saving CCP4 map to Native_refine_6_map_coeffs_superposed.ccp4 Wrote Native_refine_6_chainA_superposed.pdb Wrote Native_refine_6_map_coeffs_superposed.ccp4 saving CCP4 map to Native_refine_6_map_coeffs_superposed.ccp4 Wrote Native_refine_6_chainB_superposed.pdb Wrote Native_refine_6_map_coeffs_superposed.ccp4 A bug? Stephen P.S. Incidentally, is there any reason you write both maps in a P1 box rather than putting map 2 on top of map 1 in the position of the first (reference) model? The latter would be more convenient as you can then re-use your carefully chosen pymol viewing matrices for all figures... -- Dr Stephen Graham 1851 Research Fellow Cambridge Institute for Medical Research Wellcome Trust/MRC Building Addenbrooke's Hospital, Hills Road Cambridge, CB2 0XY, UK Phone: +44 1223 762 638

2 1

question about xtriage
by crystallogrphy 02 Dec '10

02 Dec '10

Hi, I am checking the space group of my data set using phenix.xtriage. However, the program terminated in the following scripts. Does anyone knows why and what should I do? Thanks! phenix.xtriage B13_native.sca Traceback (most recent call last): File "/usr/local/phenix-1.6.4-486/cctbx_project/mmtbx/command_line/xtriage.py", line 20, in <module> xtriage.run(args=sys.argv[1:]) File "/usr/local/phenix-1.6.4-486/cctbx_project/mmtbx/scaling/xtriage.py", line 866, in run output_file = open( params.scaling.input.parameters.reporting.log ,'w') IOError: [Errno 13] Permission denied: 'logfile.log'

2 1

force field of PHENIX.refine ?
by chene＠umbi.umd.edu 02 Dec '10

02 Dec '10

Dear authors of PHENIX, I just want to know what force field is used by phenix.refine, especially for simulated annealing, CHARMM, AMBER, or what else? Thanks. Chen Chen Chen, PhD, Research Associate Center for Advanced Research in Biotechnology (CARB), University of Maryland Biotechnology Institute (UMBI), 9600 Gudelsky Drive, Rockville, MD 20850 240-314-6122

3 2

distance constrain in phenix
by r n 02 Dec '10

02 Dec '10

Hi I would like to do distance restrain refinement for metal coordination. How can I do metal - protein strict distance restrain refinement in phenix? Thanks ram

3 4

Hydrogen occupancies
by Florian Sauer 01 Dec '10

01 Dec '10

Dear all, I'm currently refining a small protein against 0.96A data using phenix.refine (1.6-486) from the GUI. Refinement without hydrogens yielded Rwork/free of 12.4/14.2. I then added H-atoms with phenix.reduce and switched to individual refinement of hydrogens which dropped the R-factors to 11.3/13.4. After this refinement the occupancies of some hydrogens dropped to 0 which makes me a little concerned. See e.g.: ATOM 752 H LYS A 45 -7.443 2.613 9.537 1.00 4.73 H ATOM 753 HA LYS A 45 -9.707 1.128 9.401 0.55 2.42 H ATOM 754 HB2 LYS A 45 -8.476 2.797 7.435 1.00 10.58 H ATOM 755 HB3 LYS A 45 -9.795 1.948 7.166 0.95 8.70 H ATOM 756 HG2 LYS A 45 -10.840 3.210 8.889 0.00 33.83 H ATOM 757 HG3 LYS A 45 -9.604 4.206 8.775 0.00 33.83 H ATOM 758 HD2 LYS A 45 -11.227 3.615 6.563 0.00 62.10 H ATOM 759 HD3 LYS A 45 -11.508 4.867 7.507 0.00 62.10 H ATOM 760 HE2 LYS A 45 -9.229 5.527 6.985 1.00 83.62 H ATOM 761 HE3 LYS A 45 -9.348 4.488 5.786 0.00 83.86 H ATOM 762 HZ1 LYS A 45 -10.909 5.758 4.853 1.00 99.63 H ATOM 763 HZ2 LYS A 45 -11.131 6.534 6.042 0.35100.44 H ATOM 764 HZ3 LYS A 45 -9.917 6.715 5.257 0.39100.08 H I would like to ask whether it is possible in phenix to group the occupancies to keep them at the same value as the rest of the residue but still refine xyz and isotropic adp of the individual hydrogens? Thank you in advance for your answers, Florian

2 1

wander wall interaction
by subhash bihani 01 Dec '10

01 Dec '10

hi all I am working on a enzyme-ligand complex structure. For some analysis I need to remove wander wall interactions between protein and ligand. How can I do that in Phenix? Something similar to igroup statement in CNS? Thanks in adwance -- 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

4 3

phaser and threads
by Bryan Lepore 01 Dec '10

01 Dec '10

[ dev-595 ] [ GNU/Linux 2.6.31-20-generic #58-Ubuntu SMP four-processors ] [ compiled from source with --no-gui --no-apps --openmp ] i noticed phaser splits the rotation onto 2 threads - just wondering why that might be, since i used nproc=4 on the command line.... oh, hmmm... i seem to have just got this at the packing function : |python: malloc.c:4591: _int_malloc: Assertion `(unsigned long)(size) >= (unsigned long)(nb)' failed. -Bryan

2 3

building from source : svn_checkout
by Bryan Lepore 01 Dec '10

01 Dec '10

[ dev-596 from source using --no-gui --no-apps --openmp ] README-DEV mentions ./bin/svn_checkout ... i cannot find this (but i do find bundle_update and recompile). any pointers appreciated. -Bryan

2 1

Refinement of a complex with a symmetric inhibitor lying on the axis of symmetry
by Ina Lindemann 01 Dec '10

01 Dec '10

Hi, I try to refine a crystal structure of a protein - inhibitor complex. The protein is a homodimer and the two molecules have a crystallographic symmetry. The space group is P6122. The symmetric inhibitor binds symmetrically to the protein so that I refine the half of it. The core of the inhibitor is a five-membered ring with a nitrogen lying directly on the axis of symmetry. I achieved to refine this by setting the occupancy to 1.0. And now, here is my problem: Two C atoms of this ring system are bound across the symmetry axis - how could I include this in the refinement? Is there any possibility to generate restraints to the other part of the inhibitor (of the symmetry equivalent)?? Thanks for your help in advance, Ina Ina Lindemann Philipps-Universität Marburg Pharmazeutische Chemie AG Klebe Marbacher Weg 6 35032 Marburg Tel.: 06421/2825908

2 1

autosol : unmerged data with input_phase_file
by Bryan Lepore 30 Nov '10

30 Nov '10

dev-595 i gather that autosol will not take unmerged data when the input_phase_file originates from merged data. e.g. in "method #2 for MRSAD". -Bryan

2 2

dev-590 compile from source : -lGLU, -lGL
by Bryan Lepore 30 Nov '10

30 Nov '10

[ dev-590 ] [ gcc 4.4.4 ] [ GNU/Linux 32-bit ] [ Ubuntu 10.10-ish ] * OpenGL/GLUT is of unstable nature* got an error compiling from source (see after sig). i tried using --openmp first, then --openmp --nogui in the same build directory and the error occurs both times. just trying this out, not critical - was interested to know if there's an easy workaround anyone might think of (besides touching GL/GLUT - shudder ). -Bryan /usr/bin/ld: cannot find -lGLU /usr/bin/ld: cannot find -lGL collect2: ld returned 1 exit status scons: *** [lib/gltbx_fonts_ext.so] Error 1

2 5

sculptor : "Wrong alignment format:"
by Bryan Lepore 29 Nov '10

29 Nov '10

(starting a new thread?) trouble with alignment format. tried this : fubar.fas : Sorry: Unknown file(s): fubar.fas fubar.fasta : Sorry: Wrong alignment format: fubar.fasta i compared to one that worked yesterday and its not obvious : first sequence is the target, the pdb sequence is in there (somewhere), i changed a ">sp" to ">gi".... if anyone knows what i'm missing here,... -Bryan

3 13

International Conference on Structural Genomics May 10-14, 2011, Toronto -- registration now open
by Tom Terwilliger 29 Nov '10

29 Nov '10

Dear Colleagues, We are pleased to announce that registration is now open for the International Conference on Structural Genomics 2011, which will be held in Toronto, Canada on May 10-14, 2011. This meeting is the 6th in this series of biennial meetings of the International Structural Genomics Organization. The meeting is designed to serve as a forum to discuss the most recent developments in structural genomics, structural/chemical biology, and their impact on research in biology, medicine and disease. A substantial number of short talks will be selected from submitted poster abstracts. Also, a limited number of student travel fellowships will be available for travel to this meeting. Please see http://www.icsg2011.org for further details. We hope to see you there! Best regards, Cheryl Arrowsmith ICSG 2011 Organizer Ted Baker Stephen Burley Dino Moras Joel Sussman Shigeyuki Yokoyama Tom Terwilliger ISGO Executive Committee

1 0

occupancy refinement
by Dalibor Milic 25 Nov '10

25 Nov '10

Dear all! How to refine occupancies of mutually dependent structural fragments found in alternate conformations, but not consecutive in a sequence? E.g. I want to refine occupancies of two alternate conformations of residues 30 and 451 (see below the corresponding parts of a PDB file). Because of the nature of the structural disorder, an occupancy for altLoc A should be strictly the same in both fragments - the same is valid for altLoc B. Of course, the sum of occupancies for altLoc A and B should add up to 1. Thanks for your suggestions, Dalibor Milic ... ATOM 376 N AGLU A 30 23.241 86.975 31.961 0.52 17.13 N ATOM 377 N BGLU A 30 18.971 88.677 32.817 0.48 9.05 N ATOM 378 CA AGLU A 30 23.240 85.507 31.907 0.52 18.89 C ATOM 379 CA BGLU A 30 19.108 87.234 32.991 0.48 10.00 C ATOM 380 CB AGLU A 30 21.848 84.946 32.256 0.52 18.78 C ATOM 381 CB BGLU A 30 17.727 86.562 33.046 0.48 11.01 C ATOM 382 CG AGLU A 30 21.247 85.562 33.524 0.52 26.38 C ATOM 383 CG BGLU A 30 16.935 87.037 34.275 0.48 19.26 C ATOM 384 CD AGLU A 30 20.264 84.659 34.277 0.52 33.36 C ATOM 385 CD BGLU A 30 15.694 86.223 34.542 0.48 26.60 C ATOM 386 OE1AGLU A 30 19.437 83.968 33.620 0.52 31.85 O ATOM 387 OE1BGLU A 30 15.401 85.288 33.762 0.48 33.05 O ATOM 388 OE2AGLU A 30 20.310 84.678 35.535 0.52 29.12 O ATOM 389 OE2BGLU A 30 15.007 86.527 35.538 0.48 30.95 O ATOM 390 C AGLU A 30 23.750 84.946 30.575 0.52 18.22 C ATOM 391 C BGLU A 30 19.976 86.652 31.877 0.48 10.17 C ATOM 392 O AGLU A 30 24.319 83.853 30.534 0.52 19.40 O ATOM 393 O BGLU A 30 20.621 85.606 32.059 0.48 10.67 O ... ATOM 4627 N AALA A 451 31.047 82.837 22.335 0.52 18.30 N ATOM 4628 N BALA A 451 29.107 84.432 26.354 0.48 9.81 N ATOM 4629 CA AALA A 451 31.958 83.502 23.284 0.52 18.73 C ATOM 4630 CA BALA A 451 29.655 85.424 27.274 0.48 9.73 C ATOM 4631 CB AALA A 451 32.880 84.486 22.541 0.52 20.38 C ATOM 4632 CB BALA A 451 30.469 86.474 26.478 0.48 11.60 C ATOM 4633 C AALA A 451 31.176 84.206 24.406 0.52 18.79 C ATOM 4634 C BALA A 451 28.556 86.128 28.089 0.48 6.84 C ATOM 4635 O AALA A 451 30.019 84.593 24.197 0.52 18.25 O ATOM 4636 O BALA A 451 27.441 86.357 27.604 0.48 6.81 O ... -- Dalibor Milic, Dr. Sc. Laboratory of General and Inorganic Chemistry Department of Chemistry Faculty of Science University of Zagreb Horvatovac 102a HR-10000 Zagreb Croatia phone: +385 1 460 6377 fax: +385 1 460 6341 e-mail: dmilic(a)chem.pmf.hr

3 8

Re: [phenixbb] occupancy refinement
by Dalibor Milic 25 Nov '10

25 Nov '10

Just an additional comment to my previous e-mail: the refined final value of occupancies in the expression occupancy(AGLU A 30) = occupancy(BGLU A 30) = occupancy(AALA A 451) = occupancy(BALA A 451) were not equal to 0.50 and I got occupancy(AGLU A 30) + occupancy(BGLU A 30) = occupancy(AALA A 451) + occupancy(BALA A 451) < 1.00. Regards, Dalibor -- Dalibor Milic Laboratory of General and Inorganic Chemistry Department of Chemistry Faculty of Science University of Zagreb Horvatovac 102a HR-10000 Zagreb Croatia phone: +385 1 460 6377 fax: +385 1 460 6341 e-mail: dmilic(a)chem.pmf.hr

1 0

same origin
by r n 24 Nov '10

24 Nov '10

Hi all I wanted to compare the two models from two data set. I want the MR replacement in same origin to compare their them with respect to the map and model. I tried reforigin in ccp4 , it gave me some error , that I am not aware or not clear. I did tried overlay the two model and rfree is going up. Most likely different hands, so Is there any way to compare the two model and data and modify them in Phenix either reindex or shift the model origin to same in all cases. Thought phenix.xmanip could help, not seen any option to do it command line or in the GUI. I did tried def file, but phenix.xmanip params.def , it complains me that unknown file or keyword? any help will be appreciated. Thanks ram

2 2

sculptor
by Bryan Lepore 22 Nov '10

22 Nov '10

dev-572 I am checking out sculptor, looks pretty cool, but I provided what i thought was a clustal format alignment to sculptor using --verbose, and it reports "Sorry: Wrong alignment format:" that's it. sounds like user error, but i'm not sure where to proceed next. i tried mafft and clustalw2 output for this. a sample of a bit of the top of the .aln is below. -Bryan p.s: shouldn't readseq convert it to .pir? that'd be code 14 PIR/CODATA? CLUSTAL format alignment by MAFFT L-INS-i (v6.833b) gi|16132042|ref --ANVRLQVEGLSGQLEKNVRAQLSTIESDEVTPDRRFRARVDDAIREGLKALGYYQPTI gi|157163695|re --ANVRLQVEGLSGQLEKNVRAQLSTIESDEVTPDRRFRARVDDAIREGLKALGYYQPTI

4 11

sculptor : alternate conformations
by Bryan Lepore 22 Nov '10

22 Nov '10

dev-586 i'll probably try this anyways, but what is the advice on alternate conformations with sculptor? -Bryan

2 1

Re: [phenixbb] pid's from output
by Thomas C. Terwilliger 20 Nov '10

20 Nov '10

Hi Bryan, I added output such as the following to all autosol/autobuild/ligandfit/automr log files: HOST = sigma.lbl.gov HOSTTYPE = x86_64-linux USER = terwill PID = 27371 This will be available starting with tomorrow's PHENIX build. Thanks for suggesting this! -Tom T any way to get pid's of a given phenix process from their output (or any other tricks)? e.g. tried lsof, etc., dug through some title.dat, correlating date.. this would be helpful e.g. when a bunch of autobuilds start after an automr with an expired pid. -Bryan _______________________________________________ phenixbb mailing list phenixbb(a)phenix-online.org http://phenix-online.org/mailman/listinfo/phenixbb

1 0

phenix refinement question
by Katarina Moravcevic 20 Nov '10

20 Nov '10

Hi all, I have only started using PHENIX and have made an interesting observation. If I run my model through PHENIX refine (individual sites, individual ADPs, occupancies, optimized X-ray/stereochemistry weight, optimized X-ray/ADP weight options selected) I obtain R/Rfree 0.15/0.18 and very nice geometry values (data is to 1.7A). However, if I now take this refined model and run it through Refmac (without any refinement) I obtain R/Rfree of 0.19/0.21. Can anyone explain why is there such a big difference? Will this be an issue during PDB submission? Any and all comments are appreciated. Thank you all in advance Katarina M.

6 10

phenix.elbow with gamess
by William Scott 19 Nov '10

19 Nov '10

The command phenix.elbow --games --smiles="CCO" now fails (1.61 and most recent nightly), but works in an older version of phenix (1.6). Do I need to change input syntax? (I've tried both with older and more recent 32 and 64-bit gamess). Thanks. Bill William G. Scott Contact info: http://chemistry.ucsc.edu/~wgscott/ PS: Here is the error output 0:07 WriteGAMESSInput WriteGAMESSInput WriteGAMESSInput WriteGAMESSInput Write 0 -3.5160526786 -17.580263393 1 4.03396746237 20.1698373118 2 -2.978694357 -14.893471785 3 2.92412711501 14.620635575 4 -3.25561740601 -16.2780870301 5 3.2398773487 16.1993867435 [-3.5160526785963961, 4.0339674623667952, -2.9786943570030093, 2.9241271150050951, -3.2556174060111247, 3.2398773486970533] Traceback (most recent call last): File "/usr/local/xtal/phenix-1.6.1-357/elbow/elbow/command_line/builder.py", line 1691, in <module> run() File "/usr/local/xtal/phenix-1.6.1-357/elbow/elbow/command_line/builder.py", line 1583, in run call_back_function=call_back_function, File "/usr/local/xtal/phenix-1.6.1-357/elbow/elbow/command_line/builder.py", line 874, in post_process call_back_function=call_back_function, File "/usr/local/xtal/phenix-1.6.1-357/elbow/elbow/utilities/optimisation_manager.py", line 514, in optimisation_manager function(*args, **kwds) File "/usr/local/xtal/phenix-1.6.1-357/elbow/elbow/quantum/GAMESSMixins.py", line 403, in WriteGAMESSInput points, outl = get_density_grid(self) File "/usr/local/xtal/phenix-1.6.1-357/elbow/elbow/quantum/GAMESSMixins.py", line 190, in get_density_grid assert 0 AssertionError

2 1

pid's from output
by Bryan Lepore 19 Nov '10

19 Nov '10

any way to get pid's of a given phenix process from their output (or any other tricks)? e.g. tried lsof, etc., dug through some title.dat, correlating date.. this would be helpful e.g. when a bunch of autobuilds start after an automr with an expired pid. -Bryan

1 0

relationship between B factors and Koff
by Sebastiano Pasqualato 19 Nov '10

19 Nov '10

Hi all, I have a crystallographical/biochemical problem, and maybe some of you guys can help me out. We have recently crystallized a protein:protein complex, whose Kd has been measured being ca. 10 uM (both by fluorescence polarization and surface plasmon resonance). Despite the 'decent' affinity, we couldn't purify an homogeneous complex in size exclusion chromatography, even mixing the protein at concentrations up to 80-100 uM each. We explained this behavior by assuming that extremely high Kon/Koff values combine to give this 10 uM affinity, and the high Koff value would account for the dissociation going on during size exclusion chromatography. We have partial evidence for this from the SPR curves, although we haven't actually measured the Kon/Koff values. We eventually managed to solve the crystal structure of the complex by mixing the two proteins (we had to add an excess of one of them to get good diffraction data). Once solved the structure (which makes perfect biological sense and has been validated), we get mean B factors for one of the component (the larger) much lower than those of the other component (the smaller one, which we had in excess). We're talking about 48 Å^2 vs. 75 Å^2. I was wondering if anybody has had some similar cases, or has any hint on the possible relationship it might (or might not) exist between high a Koff value and high B factors (a relationship we are tempted to draw). Thanks in advance, best regards, ciao s -- Sebastiano Pasqualato, PhD IFOM-IEO Campus Dipartimento di Oncologia Sperimentale Istituto Europeo di Oncologia via Adamello, 16 20139 - Milano Italy tel +39 02 9437 5094 fax +39 02 9437 5990

1 0

automr/phaser : b-factors
by Bryan Lepore 18 Nov '10

18 Nov '10

does (or can) phaser use the search probe B-factors in any way? -Bryan

2 3

CRYST card being changed
by Nino Campobasso 18 Nov '10

18 Nov '10

Phenix developers, When I run a refinement, it appears that the CRYST card is not being written as expected by programs like pymol and coot. I'll start with a file that has C 1 2 1 and phenix.refine will write out C2 . If I attempt to view symmetry related molecule within coot or pymol, these programs complain that there is no symmetry op. If "fix" the spacegroup in the CRYST line, pymol and coot are happy and generate the related molecules. Best, Nino

2 1

autobuild : rigid?
by Bryan Lepore 18 Nov '10

18 Nov '10

dev-572 my understanding of the autobuild<->refine keyword table on the website is that there's no rigid body available through autobuild - is that right? -Bryan

6 6

reflection file converter
by r n 17 Nov '10

17 Nov '10

Hi I did tried reflection file converter in phenix. I do have the anomalous data, like to keep all anomalous difference and also merged F and FreeR for refinement. Phenix GUI give me an warning message that "File are not suitable for amplitudes change them as F or leave them as Intensities" Also have hard time to convert them in command line, Is there any auto script that can work similar to scaelpack2mtz , which keeps all I+, I- and F , SigF, Rfree for refinement. Converted using Phenix.refine , that keeps separate I+, I- and Rfree+, Rfree-. Thanks ram

2 1

link restraints for modified nucleotides in RNA
by Gerwald Jogl 17 Nov '10

17 Nov '10

Hi All, I am refining several modified nucleotides in ribosomal RNA (N2-methyl-guanosine and similar). While the restraints generation with phenix.reel work quite well, I am unable to get the phosphate oxygen restraints applied properly. Phenix.refine recognizes the links but for some reason the phosphate groups come out quite distorted... For example, angles after refinement: O3 (1) - P (2, modified) - OP2 (2) = 80 deg and O3 (1) - P (2, modified) - OP1 (2) = 127 deg Here is part of the phenix.refine log on reading the coordinates: Chain: "A" Number of atoms: 32526 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1513, 32526 Unusual residues: {'MA6%rna3p': 2, '5MC%rna3p': 3, '5MC%rna2p': 1, '7M G%rna3p': 1, '2MG%rna3p': 1, 'UR3%rna2p': 1, 'M2G%rna3p': 1, '4OC%rna3p': 1} Classifications: {'undetermined': 11, 'RNA': 1502} Modifications used: {'5*END': 1, 'rna3p': 1318, 'rna2p': 195} Link IDs: {'rna3p': 1318, 'rna2p': 194} Chain breaks: 1 I am at a loss about how to fix this as all links and restraints seem to be in place. Many thanks for advice, Gerwald

2 1

phaser and automr
by Bryan Lepore 16 Nov '10

16 Nov '10

IIUC, phaser keywords cannot be provided to phenix.automr - is that right? i.e. if that's the case, then i ought to run phenix.phaser. -bryan

3 2

automr : [1] MAP_COEFFS [2] output pdb's of ensemble
by Bryan Lepore 16 Nov '10

16 Nov '10

dev-572 [1] in automr, i do not understand why only the last MR.MAP_COEFFS.n.mtz is written- e.g. i don't get 1-9 if 10 models are requested as output. (MR.[1-n].mtz are written ok). [2] when using an ensemble, i gather that the output is *hardwired* for the best rmsd model - i.e. there are no keywords to write out all parts of the ensemble. (i know i can run pdbset for them, but anyways...) -Bryan

2 2

automr : *refinement* resolution limit keyword?
by Bryan Lepore 16 Nov '10

16 Nov '10

[ dev-572 ] i'm wondering if users can set the *refinement* resolution limits in automr? the webpage suggests this, but it wasn't clear to me. i can set the search resolution fine, however by default as indicated, all data is used - and the only other keyword i see that would at all be relevant is the building resolution. -Bryan

3 3

refine : annealing temperature
by Bryan Lepore 13 Nov '10

13 Nov '10

just wondering if the annealing temperature in phenix is absolutely the same as in cns. -Bryan

2 1

Question about restraints
by Schubert, Carsten [PRDUS] 12 Nov '10

12 Nov '10

Hi, I was wondering what the default values for restraints in reel/elbow represent in terms of 'tightness' for lack of a better word. I am trying to tighten up the geometry of a badly behaving ligand and want to avoid over tightening. What is the consensus out there what kind of ESDs represent a 'loose' 'medium' and 'tight' restraint for bond, angle, dihedral, planar and chiral restraints. Thanks Carsten

3 6