- phenixbb - phenix-online.org

phenix.map_to_object for HA sites
by Michael Hothorn 16 Dec '10

16 Dec '10

Hi, I have a .pdb with heavy atoms that are dispersed all over the unit cell. I would like to move all of them (about 40) as close as possible to my refined structure. Could phenix.map_to_object do that for all heavy atoms at once, or would I need to do this in an atom by atom manner? Thanks Michael

2 1

Map filenames are not honored [dev-610]
by Schubert, Carsten [PRDUS] 15 Dec '10

15 Dec '10

Hi, I just noticed that v610 does not honor map file names defined in the .def file. Instead it constructs its own file name from the prefix, serial and map_type settings. That works fine for run-of-the-mill fofc and 2fofc maps, but causes a clash in map filenames when one also computes kicked maps or other maps with non-default settings beside the default fofc and 2fofc maps. Cheers, Carsten

2 1

phenix.find_tls_groups: a new tool for automated partitioning a model into TLS groups
by Joseph Noel 15 Dec '10

15 Dec '10

Pavel or Nat, Can you point me to the button in the GUI to use to run this? I am pretty much a GUI guy although not adverse to running things through the command like like the old days (yes, I was a command like geek when I started at Yale in 1990 and until I stopped (temporarily) doing anything worthwhile around 2001. Now that I am back in the lab doing work (but probably still not worthwhile) I have become addicted to GUIs and buttons. Blame my daughters! LOL Joe ___________________________________________________________ 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 ___________________________________________________________

2 2

phenix.find_tls_groups: a new tool for automated partitioning a model into TLS groups
by Pavel Afonine 14 Dec '10

14 Dec '10

PHENIX users: starting dev-610 (development version of PHENIX) there is a new tool available for completely automated partitioning a model into TLS groups: http://www.phenix-online.org/download/nightly_builds.cgi To run: phenix.find_tls_groups model.pdb or if you have a multiple CPU machine: phenix.find_tls_groups model.pdb nproc=N where N is the number of CPUs available (thanks Nat for parallelization!). There is no parameters that a user is supposed to tweak (except defining the number of CPUs, if desired). The result of running the above command are atom selections that define TLS groups. These atom selections are ready to use in phenix.refine. This is available from PHENIX GUI too, where automatically defined TLS groups can be readily visualized and checked on the graphics (thanks Nat!). The algorithm is fast. For example, for a GroEL structure (3668 residues, 26957 atoms, 7 chains) it takes only 135 seconds using 1 CPU, and 44 seconds using 10 CPUs. Analogous job takes 3630 seconds using TLSMD server. For a lysozime structure it takes 9.5 seconds with one CPU, and 2.5 seconds using 10 CPUs. The timing results may vary depending on the performance of your computer. There is ongoing work that will slightly improve phenix.find_tls_groups within the next few weeks / a month; however the current version is functional and can be tried now. An example of such improvements are analyzing (scoring) user-defined TLS groups (for example, TLS groups from PDB file header), automated combining cross-chain TLS groups (non-contiguous segments) that will be obtained through connectivity analysis, better handling non-protein chains, and more. Integration with phenix.refine is also planned. Any feedback is very much appreciated! Thanks, Pavel.

2 2

request for openmp compile comments
by Bryan Lepore 14 Dec '10

14 Dec '10

if anyone has specifically gotten automr to run on an openmp-compiled phenix installation with something like the specs i give below, i would appreciate any details you could offer. (in my hands, automr always crashes, refine works.) -Bryan any recent nightly release with --openmp gcc/g++ 4.4.5 (or thereabouts) 2.6.35-23-generic SMP i686 (or thereabouts) amd, intel, (or thereabouts) Ubuntu (recent - but not Fedora)

1 0

phenix.refine use_experimental_phases vs target=mlhl
by Francis E Reyes 13 Dec '10

13 Dec '10

Are these settings redundant? i.e. what happens when you have target=mlhl, use_experimental_phases=False? Thanks! 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

6 5

AutoSol/PhaserSAD with more than one anomalously-scattering atom produces error...
by Francis E Reyes 13 Dec '10

13 Dec '10

Hi all I'm getting an error... "'InputEP_SAD' object has no attribute 'addSCAT_TYPE'" from autosol. Facts: - phenix-1.6-289 - SAD experiment with two types of anomalous scatterers in the xtal, data is the peak wavelength of one of the anomalous scatterers. - sites are input from a previous phaser run (ha_x.pdb_formatted.pdb) - the atom positions are likely correct, but the atom types are all the same atom (I'd like phenix/phaser to go through the list and reannotate/refine based on the proper heavy atom identities) - additional atom types/f'/f'' defined in autosol.eff, Using the documentation example "SAD with more than one anomalously-scattering atom" (http://www.phenix-online.org/documentation/autosol.htm). Thanks! 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

changing test set flag
by Michael Hothorn 11 Dec '10

11 Dec '10

Dear all, I used phenix.refine to generate the test set flag for my dataset which it did with a value of 1. I now use a program that expects the test set flag in the mtz file to be 0. Is there a simple way to change the flag from 1 to 0 keeping the flagged reflections identical? Thanks! Michael

3 3

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