- phenixbb - phenix-online.org

Re: [phenixbb] special positions
by o.h.weiergraeber＠fz-juelich.de 02 Jul '08

02 Jul '08

Adding to my previous post, I have found another strange issue (possibly a bug?) concerning special positions: In my case, the "special" atom (a cation on a three-fold in sg I23), while being considered during refinement, seems to be completely omitted (or heavily downweighted) for map calculations! The 2fo-fc map shows a "hole" at this position, while there is huge difference density in the fo-fc map (largely irrespective of whether or not the occupancy is reduced in the input.) The same pdb file, used in cns just for map calculation, gives a beautiful density at the cation position, while the rest of the map is more or less identical. At the moment I cannot think of any error in the input since refinement of the atom does take place, with the proper restrictions due to the three-fold axis. But something is going seriously wrong. Any ideas? Oliver ------------------------------- Dr. Oliver H. Weiergraeber Institute of Neurosciences and Biophysics Molecular Biophysics Research Centre Juelich D-52425 Juelich Germany Phone: +49-2461-612028 Fax: +49-2461-612020 ------------------------------- ----- Original Message ----- From: o.h.weiergraeber(a)fz-juelich.de Date: Wednesday, July 2, 2008 2:55 pm Subject: [phenixbb] special positions > Hello, > > a very simple question (although I could not find *anything* about > this in the documentation!): > What does phenix.refine expect as occupancy for an atom on a > special position: reduced occupancy (like refmac) or full occupancy > (like cns_solve)? > > Oliver > > > > ------------------------------- > Dr. Oliver H. Weiergraeber > Institute of Neurosciences and Biophysics > Molecular Biophysics > Research Centre Juelich > D-52425 Juelich > Germany > Phone: +49-2461-612028 > Fax: +49-2461-612020 > ------------------------------- > > > > ------------------------------------------------------------------- > ------------------------------------------------------------------- > Forschungszentrum Juelich GmbH > 52425 Juelich > > Sitz der Gesellschaft: Juelich > Eingetragen im Handelsregister des Amtsgerichts Dueren Nr. HR B 3498 > Vorsitzende des Aufsichtsrats: MinDir'in Baerbel Brumme-Bothe > Geschaeftsfuehrung: Prof. Dr. Achim Bachem (Vorsitzender), > Dr. Ulrich Krafft (stellv. Vorsitzender), Prof. Dr. Harald Bolt, > Dr. Sebastian M. Schmidt > ------------------------------------------------------------------- > ------------------------------------------------------------------- > _______________________________________________ > phenixbb mailing list > phenixbb(a)phenix-online.org > http://www.phenix-online.org/mailman/listinfo/phenixbb > ------------------------------------------------------------------- ------------------------------------------------------------------- Forschungszentrum Juelich GmbH 52425 Juelich Sitz der Gesellschaft: Juelich Eingetragen im Handelsregister des Amtsgerichts Dueren Nr. HR B 3498 Vorsitzende des Aufsichtsrats: MinDir'in Baerbel Brumme-Bothe Geschaeftsfuehrung: Prof. Dr. Achim Bachem (Vorsitzender), Dr. Ulrich Krafft (stellv. Vorsitzender), Prof. Dr. Harald Bolt, Dr. Sebastian M. Schmidt ------------------------------------------------------------------- -------------------------------------------------------------------

2 1

Fwd: [COOT] Using Coot with PHENIX
by Stephen Graham 02 Jul '08

02 Jul '08

Below is an email from Kevin Cowtan (posted on the COOT bb) that may be of interest to people who use COOT for building. Stephen ---------- Forwarded message ---------- From: Kevin Cowtan <cowtan(a)ysbl.york.ac.uk> Date: Jul 2, 2008 9:31 AM Subject: [COOT] Using Coot with PHENIX To: COOT(a)jiscmail.ac.uk There are some issues which may be affecting some people using Coot with PHENIX. If you try and use the two together in the wrong way, then you may get unexpected results. (Problems hashed out with help from Dirk). phenix.refine produces two different sets of output files. I quote from the documentation: At the end of refinement the program generates: ... two maps: likelihood weighted mFo-DFc and 2mFo-DFc. These are in ASCII X-PLOR format. A reflection file with map coefficients is also generated for use in Coot or XtalView (e.g. lysozyme_refine_001_map_coeffs.mtz); There are some CNS maps, which may be used in conventional graphics programs, and there is an MTZ file, which may be used in Coot or XtalView. The difference is important. Coot and XtalView have a fundamentally different concept of what a map is to other graphics programs. As a result they require different information. Phenix.refine recognises this and provides appropriate information for each type of program. In the last year we incorporated the ability to read CNS files into coot. As a result, some users have started reading the CNS map files into coot instead of the MTZ file. Unfortunately, this does not always work the way you would expect, because the map files has not been generated with Coot in mind. This may lead to misleading sigma levels, or, if the map is modified or masked in any way before giving it to Coot, much more serious artifacts. The solution is simple: When using Coot, use the MTZ file which the phenix developers have provided for that purpose, *not* the CNS map files. In addition to giving the right answers, there are some other benefits: - The MTZ is much smaller, in case you want to archive it. - Using auto-open mtz, you get both maps at once, and the difference map is automatically displayed as such, saving mouse clicks. - The spacegroup and cell information are in the file, so it works even if you don't open the PDB file first. .. You can stop reading here, the rest is detail for those who are interested: Coot and XtalView use a different representation of electron density, a continuous crystal representation. Most other graphics programs use a bounded box representation. The continuous crystal representation means that the electron density is everywhere, it is omnipresent. "If I scroll up to the heavens, density is there; if I place my ligand in the depths, density is there. If I drag over to the far side of the solvent, even there the density will guide my fitting." (*) To generate a continuous crystal representation, you need a complete and consistent description of the electron density over the asymmetric unit (any asymmetric unit - it doesn't matter where in space). The bounded box representation need the electron density over the volume of the molecule. It is not always possible to get the asymmetric unit representation from the bounded box representation, because in creating the bounded box some information may have been lost. However, you can always get it from an MTZ file, or from a map file on the condition that the map file is appropriately bounded. Kevin (*) Psalm 139 v8-10 -- Dr Stephen Graham Nuffield Medical Fellow Division of Structural Biology Wellcome Trust Centre for Human Genetics Roosevelt Drive Oxford OX3 7BN United Kingdom Phone: +44 1865 287 549

1 0

distance restraints
by o.h.weiergraeber＠fz-juelich.de 02 Jul '08

02 Jul '08

Hello, I am about to define distance restraints for a metal coordination site using a refinement.geometry_restraints.edits entry. Unfortunately, two of the four ligand side chains come from a symmetry-related molecule! What is the syntax to properly refer to a symm in such a case? Oliver ------------------------------- Dr. Oliver H. Weiergraeber Institute of Neurosciences and Biophysics Molecular Biophysics Research Centre Juelich D-52425 Juelich Germany Phone: +49-2461-612028 Fax: +49-2461-612020 ------------------------------- ------------------------------------------------------------------- ------------------------------------------------------------------- Forschungszentrum Juelich GmbH 52425 Juelich Sitz der Gesellschaft: Juelich Eingetragen im Handelsregister des Amtsgerichts Dueren Nr. HR B 3498 Vorsitzende des Aufsichtsrats: MinDir'in Baerbel Brumme-Bothe Geschaeftsfuehrung: Prof. Dr. Achim Bachem (Vorsitzender), Dr. Ulrich Krafft (stellv. Vorsitzender), Prof. Dr. Harald Bolt, Dr. Sebastian M. Schmidt ------------------------------------------------------------------- -------------------------------------------------------------------

1 0

tlsanal
by john kryst 02 Jul '08

02 Jul '08

Hi all !! I have coupls of questions regarding TLS in phenix. 1. Has anyone tried using program tlsviewer with phenix refined pdb. It is not able to detect TLS groups unless you modify your residue range in the pdb to refmac convention. Is there any other program available to visualize the TLS movements. 2. I would like to use tlsanal to see how much is the contibution from TLS to the ADP. Is there any utility availble in TLS to do it. If i need to use tlsanal.. what should i modify in the pdb. thanks regards John

2 1

Postdoctoral Position - Chemical Imaging of Plant Biomass with Micro- and Nano-Raman Spectroscopy
by Paul Adams 01 Jul '08

01 Jul '08

Postdoctoral Researcher Chemical Imaging of Plant Biomass with Micro- and Nano-Raman Spectroscopy Job ID: UCB-013 Institution: Energy Biosciences Institute, Berkeley California Date Opened: 5/10/2008 Position Description: Development of computational tools for the analysis of chemical imaging data. Computational construction of three- dimensional maps of plant cell wall material from Raman spectroscopy measurements. Position Responsibilities: Development of algorithms to indentify chemical species from Raman spectroscopy measurements. Processing of data to generate two-dimensional images of chemical components in plant cell wall material. Use of tomographic methods to construct three-dimensional images. Development of maximum likelihood methods to extract information from weak spectroscopic data. Position Qualifications: Strong computational skills with experience in their practical application to imaging or complex data analysis. Preference will be given to candidates with experience in the analysis of spectroscopic data or imaging of biological samples. Experience with programming languages: C/C++ and scripting languages: Python/ Perl. Demonstrated publication record in the computational analysis of data and algorithm development. Ability to work in a team environment and good communication skills. Ph.D. in physical chemistry, biology, physics, computer science or related field. See the program description for more details: http://energybiosciencesinstitute.org/index.php?option=com_content&task=vie… How to Apply: Submit a Letter of Interest, CV, and three (3) contacts for references via email to: Paul Adams Lawrence Berkeley Lab, BLDG 64R0121 1 Cyclotron Road Berkeley, CA 94720 510-486-4225 pdadams(a)berkeley.edu The University of California, Berkeley is an Equal Opportunity Employer committed to diversity in its staff, faculty, and student body, and invites all qualified people to apply, including minorities and women, veterans and individuals with disabilities. -- Paul Adams Adjunct Professor, Department of Bioengineering, U.C. Berkeley Tel: 510-486-4225, Fax: 510-486-5909 http://cci.lbl.gov/paul --

1 0

Fwd: Re: help with refinement in Phenix
by hosur＠barc.gov.in 30 Jun '08

30 Jun '08

Hi PHENIX developers, I am forwarding herewith my correspondence with Dr. Pavel Afonine regarding a reinement issue in PHENIX. I would appreciate if other developers can help me with the solution before Pavel comes back on June 25th. Thanks. The requirement is to be able to refine a model in which a piece of electron density is interpreted as partial occupancy protein residue and a set of water molecules. Hosur M.V. ----- Forwarded message from Pavel Afonine <pafonine(a)lbl.gov> ----- Date: Sat, 14 Jun 2008 21:57:36 -0700 From: Pavel Afonine <pafonine(a)lbl.gov> Reply-To: Pavel Afonine <pafonine(a)lbl.gov> Subject: Re: help with refinement in Phenix To: hosur(a)barc.gov.in Hi Hosur, I'm out of town and not regular at email, so the delay with the answer,, sorry for this. Actually, it is the best to address your questions to phenixbb (http://www.phenix-online.org/mailman/listinfo/phenixbb) so all the PHENIX developers get it and you don't hang up with the delay if a particular person is temporarily absent. I will be back in office on June 25. Regarding your question: > In my case, the set of 4 water molecules has to have the same occupancy, and > this occupancy with the occupancy of a side chain must make up unity. That is > the sum of two occupancies must be unity. > This particular feature is currently not available in phenix.refine. I will add this option once I'm back in office. > If I understand correctly, if i identify each of the four water molecules by a > different altloc label, then the sum of 5 occupancies will be adjusted to one, > with the possibility of each water occupancy converging to different value. > Yes, this is expected behavior. > Further, the final occupancy would be close to 0.25 of the required value. The > complication, in my opinion, is because the set of four water molecules is not > like set of atoms that make one residue which can be assigned a single altloc > label. I hope I have been able to convey what I mean. Let me have the suitable > PDB, which I can try out. > Yes, this is very clear for me and like I said above I will add this feature to PHENIX as soon as I get back to office. All the best, Pavel. ----- End forwarded message ----- -------------------------------------------------

2 1

refinement with several twins / modulated data / pseudo-translational symmetry
by iulek 29 Jun '08

29 Jun '08

Dear all, I have a strongly modulated dataset as can be observed in the truncate output below. Specially, the modulation is higher for L odd/even and at lower resolution: !!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Analysis of mean intensity by parity for reflection classes For each class, Mn(I/sig(I)) is given for even and odd parity with respect to the condition, eg group 1: h even & odd; group 7 h+k+l even & odd; group 8 h+k=2n & h+l=2n & k+l=2n or not Range Min_S Dmax Nref 1 2 3 4 5 6 7 8 h k l h+k h+l k+l h+k+l h+k,h+l,k+l 1 0.00018 20.85 519 13.2 8.9 13.9 8.2 22.3 2.4 13.3 8.9 12.5 9.6 12.9 9.0 12.6 9.5 16.5 9.2 2 0.00230 15.04 994 15.9 16.0 15.3 16.5 30.1 2.3 15.7 16.2 15.1 16.7 14.9 16.8 15.1 16.7 13.5 16.6 3 0.00442 12.36 1289 17.2 16.0 17.0 16.2 33.8 3.7 17.2 16.0 14.9 18.0 14.7 18.2 14.9 18.0 13.1 17.4 4 0.00655 10.74 1568 20.7 18.9 20.3 19.2 35.5 4.9 20.4 19.2 18.7 20.7 18.3 21.0 18.5 20.9 17.7 20.3 5 0.00867 9.63 1759 20.3 15.8 20.1 16.0 33.2 4.9 20.1 15.9 18.7 17.3 18.6 17.4 18.8 17.2 21.6 16.9 6 0.01079 8.80 2042 19.6 15.4 19.3 15.7 30.5 5.8 19.5 15.5 17.6 17.3 17.5 17.4 17.6 17.3 19.9 16.7 7 0.01291 8.15 2206 15.9 12.0 15.8 12.1 22.6 5.0 15.5 12.4 14.3 13.5 14.1 13.7 14.0 13.8 16.2 13.2 8 0.01504 7.63 2393 11.8 9.0 11.5 9.3 18.6 3.6 11.7 9.1 10.5 10.2 10.3 10.3 10.4 10.3 12.0 9.9 9 0.01716 7.20 2526 9.8 6.5 9.8 6.5 13.8 3.2 9.8 6.6 9.0 7.3 9.0 7.3 9.0 7.4 11.6 7.1 10 0.01928 6.83 2747 11.1 7.5 11.0 7.6 15.7 3.0 10.8 7.8 10.0 8.5 9.9 8.6 9.8 8.6 12.3 8.3 11 0.02141 6.52 2867 10.3 7.3 10.2 7.3 14.4 3.0 10.3 7.3 9.4 8.1 9.4 8.1 9.4 8.2 11.7 7.9 12 0.02353 6.24 3122 8.6 5.8 8.5 5.8 11.6 3.0 8.4 5.9 7.3 6.9 7.3 6.9 7.2 7.0 8.8 6.6 13 0.02565 6.00 3124 7.9 5.7 7.7 5.9 11.1 2.9 7.7 5.9 6.8 6.7 6.8 6.8 6.8 6.7 7.8 6.5 14 0.02777 5.78 3346 9.8 6.7 9.7 6.7 12.7 3.5 9.6 6.8 8.8 7.7 8.7 7.8 8.6 7.9 10.7 7.5 15 0.02990 5.59 3428 9.0 5.6 8.9 5.7 11.2 3.1 8.8 5.8 7.9 6.7 7.8 6.8 7.7 6.9 9.9 6.4 16 0.03202 5.41 3633 8.3 5.5 8.2 5.5 10.9 3.5 8.2 5.6 7.1 6.6 7.1 6.6 7.1 6.6 8.8 6.2 17 0.03414 5.25 3658 9.6 6.2 9.4 6.3 12.4 3.7 9.4 6.4 8.4 7.3 8.4 7.3 8.2 7.5 10.6 7.0 18 0.03627 5.10 3696 10.2 6.4 10.1 6.5 12.5 4.3 9.8 6.7 8.8 7.8 8.7 7.9 8.6 7.9 10.8 7.4 19 0.03839 4.97 3955 10.1 6.9 9.9 7.0 12.6 4.5 9.8 7.1 8.8 8.2 8.6 8.3 8.6 8.3 10.4 7.9 20 0.04051 4.84 3902 10.9 7.2 10.9 7.2 13.4 4.6 10.7 7.4 9.7 8.4 9.7 8.4 9.7 8.4 12.1 8.0 21 0.04263 4.73 4207 11.4 6.8 11.3 6.9 13.7 4.8 11.1 7.0 9.5 8.5 9.4 8.6 9.5 8.6 12.0 8.1 22 0.04476 4.62 4169 9.7 6.8 9.6 6.9 12.4 4.6 9.5 7.0 8.3 8.1 8.3 8.2 8.2 8.2 9.7 7.8 23 0.04688 4.52 4100 11.1 6.7 11.0 6.8 12.6 5.5 10.8 7.0 9.0 8.8 8.9 8.8 8.9 8.9 11.0 8.2 24 0.04900 4.42 4605 10.8 6.5 10.7 6.5 12.3 4.9 10.5 6.8 9.5 7.8 9.4 7.9 9.3 7.9 12.2 7.5 25 0.05113 4.33 4449 9.6 6.1 9.5 6.2 11.1 4.6 9.4 6.3 8.3 7.4 8.2 7.5 8.2 7.5 10.2 7.0 26 0.05325 4.25 4528 9.1 5.8 9.1 5.8 10.0 4.8 8.8 6.0 7.6 7.2 7.5 7.3 7.5 7.3 9.2 6.8 27 0.05537 4.17 4782 8.6 5.3 8.6 5.3 9.7 4.6 8.5 5.4 7.5 6.5 7.4 6.5 7.4 6.5 9.5 6.1 28 0.05749 4.10 4774 8.3 5.2 8.2 5.3 9.6 3.9 8.0 5.4 7.0 6.4 7.0 6.4 6.9 6.5 8.6 6.1 29 0.05962 4.02 4837 8.1 5.6 8.1 5.7 9.7 4.4 7.9 5.8 7.2 6.6 7.2 6.6 7.2 6.6 8.6 6.4 30 0.06174 3.96 5066 8.2 4.9 8.2 4.9 8.9 4.2 8.0 5.1 6.9 6.2 6.8 6.2 6.7 6.3 8.6 5.8 31 0.06386 3.89 4954 7.4 4.7 7.3 4.8 7.9 4.1 7.2 4.9 6.2 5.9 6.1 6.0 6.2 5.9 7.5 5.6 32 0.06599 3.83 5052 7.6 4.6 7.6 4.6 8.2 4.0 7.4 4.8 6.4 5.8 6.4 5.8 6.3 5.9 8.0 5.5 33 0.06811 3.77 5159 6.9 4.4 6.8 4.5 7.2 4.0 6.6 4.6 5.8 5.4 5.7 5.5 5.6 5.6 6.9 5.2 34 0.07023 3.72 5226 6.7 4.1 6.7 4.2 7.2 3.8 6.6 4.2 5.6 5.2 5.6 5.2 5.6 5.2 7.0 4.9 35 0.07236 3.66 5359 6.6 4.1 6.6 4.2 6.7 4.1 6.4 4.4 5.4 5.3 5.4 5.3 5.3 5.4 6.6 5.0 36 0.07448 3.61 5293 6.4 3.9 6.3 3.9 6.9 3.5 6.2 4.1 5.4 4.8 5.4 4.9 5.4 4.8 6.7 4.6 37 0.07660 3.56 5633 6.7 4.4 6.7 4.4 7.4 3.7 6.6 4.5 5.8 5.3 5.7 5.4 5.7 5.4 7.0 5.1 38 0.07872 3.52 5479 5.9 3.8 5.9 3.8 6.2 3.6 5.8 3.9 5.0 4.7 5.0 4.7 5.0 4.7 6.1 4.5 39 0.08085 3.47 5578 5.2 3.3 5.2 3.3 5.4 3.1 5.0 3.5 4.4 4.0 4.4 4.0 4.3 4.1 5.4 3.8 40 0.08297 3.43 5580 5.3 3.3 5.3 3.3 5.3 3.3 5.2 3.4 4.6 4.0 4.6 4.1 4.6 4.1 5.8 3.8 41 0.08509 3.39 5936 5.3 3.3 5.3 3.2 5.4 3.2 5.1 3.4 4.5 4.0 4.5 4.0 4.4 4.1 5.5 3.8 42 0.08722 3.35 5783 4.7 3.0 4.7 3.0 4.5 3.3 4.6 3.1 3.9 3.8 4.0 3.8 3.9 3.8 4.8 3.6 43 0.08934 3.31 5889 4.4 3.0 4.4 2.9 4.5 2.9 4.3 3.1 3.8 3.6 3.8 3.5 3.7 3.6 4.5 3.4 44 0.09146 3.27 5966 4.4 2.8 4.5 2.8 4.5 2.8 4.3 2.9 3.8 3.5 3.8 3.5 3.8 3.5 4.7 3.3 45 0.09358 3.23 6145 4.0 2.6 4.0 2.6 4.1 2.5 3.9 2.8 3.5 3.2 3.5 3.2 3.5 3.2 4.2 3.0 46 0.09571 3.20 6202 3.6 2.3 3.7 2.3 3.6 2.3 3.6 2.4 3.2 2.8 3.2 2.7 3.2 2.7 4.0 2.6 47 0.09783 3.16 6042 3.7 2.3 3.8 2.2 3.5 2.5 3.7 2.3 3.1 2.9 3.1 2.9 3.1 2.9 4.0 2.7 48 0.09995 3.13 6198 3.5 2.2 3.5 2.1 3.3 2.4 3.4 2.2 2.9 2.8 2.9 2.8 2.9 2.8 3.6 2.6 49 0.10208 3.10 6421 3.3 2.1 3.3 2.1 3.1 2.3 3.2 2.1 2.7 2.6 2.7 2.6 2.7 2.6 3.4 2.4 50 0.10420 3.07 6403 3.2 2.2 3.3 2.1 3.1 2.3 3.2 2.2 2.8 2.6 2.8 2.6 2.8 2.6 3.3 2.5 51 0.10632 3.04 6357 2.9 2.0 3.0 1.9 2.8 2.1 2.9 2.0 2.6 2.4 2.6 2.3 2.6 2.4 3.1 2.2 52 0.10844 3.01 6389 2.6 1.9 2.6 1.9 2.6 2.0 2.6 1.9 2.4 2.2 2.4 2.2 2.3 2.2 2.8 2.1 53 0.11057 2.98 6580 2.8 1.9 2.8 1.9 2.7 2.0 2.8 2.0 2.5 2.3 2.4 2.3 2.4 2.3 2.9 2.2 54 0.11269 2.95 6663 2.8 1.9 2.8 1.8 2.7 1.9 2.7 1.9 2.5 2.1 2.5 2.1 2.5 2.1 3.1 2.0 55 0.11481 2.92 6582 2.6 1.9 2.7 1.8 2.6 1.9 2.6 1.9 2.4 2.1 2.4 2.1 2.4 2.1 2.9 2.0 56 0.11694 2.90 6572 2.6 1.8 2.6 1.8 2.5 1.9 2.6 1.8 2.3 2.1 2.3 2.1 2.3 2.1 2.7 2.0 57 0.11906 2.87 6820 2.4 1.7 2.4 1.7 2.3 1.8 2.4 1.8 2.1 2.0 2.1 2.0 2.1 2.0 2.5 1.9 58 0.12118 2.85 6970 2.3 1.6 2.3 1.6 2.2 1.8 2.3 1.6 2.0 1.9 2.0 1.9 2.0 2.0 2.3 1.8 59 0.12330 2.82 6717 2.4 1.7 2.5 1.6 2.2 1.8 2.4 1.7 2.0 2.0 2.0 2.0 2.0 2.0 2.4 1.9 60 0.12543 2.80 5677 2.1 1.6 2.1 1.6 2.0 1.7 2.0 1.6 1.8 1.8 1.8 1.8 1.8 1.8 2.0 1.8 Totals: 277916 6.4 4.3 6.4 4.3 7.6 3.2 6.3 4.4 5.5 5.2 5.5 5.2 5.5 5.2 6.6 4.9 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! As you can imagine, I have pseudo-translation peaks , the strongest at 48 % height (of the origin peak) at 0.0 0.0 0.5, other two at ~ 30 % height at 0.0 0.5 0.0 and 0.5 0.0 0.0 (and so at 0.5 0.0 0.5, etc.). The cell is large, 229 229 72 90 90 120, data processed at several space group options with good statistics; the highest symmetry is P622. With much struggle, some solutions from molecular replacement came out. Basically, the monomers (~274 residues each) makes "strings" (pack of one monomer over the other in one dimension), then these strings seem to be able to pack themselves in parallel (within them there must be an up/down feature) and parallel to the shorter cell axis. Given the cell dimensions, it seems most probable that I have 48 mol/au, although I cannot discard 36 mol/au (or even something between, depending upon the symmetry). I attach a picture of the supposed pack in case of 48 mol/au, viewed down to the 120 angle. The packing has made me practically rule out space groups P622 and P312, yet the need for a 2-fold axis passing through both 0, 1, 0 and 1, 0 0. Whenever I process the data in a lower space group, xtriage suggests strongly twinning, and close to 50 % in almost any operation. Also, whatever space group I choose for refinement, the results are all some how alike, R-free stucks (with any twinning operation that I can use to refine) at about 40/41 %; non-twinned refinement stucks R-free around 48 %. The investigated space groups include P321, P6, P3, P2, C222, C2 and P1 (besides, of course, at the beginning, P622 and P312), with several combinations of molecules/au (I there might be screw axes, although no systematic absences). With P1, I have the lowest reason between R and R-free (because R is higher). I must add that I have already tested some refinement alternatives, such as bulk solvent models, fix these, etc., but of course not all alternatives. Much information has been produced; one other highlight I think is the fact that most of the time atomic B-factors are refined to the ground, I suppose because of the abnormal intensity distribution/data modulation. I have a couple of questions and other information which I prefer to number since from now. a) Is there a source of information for refinement of such modulated data? Or previous cases reported in the literature? b) I have not tried yet P1 with 48 molecules. I run out of chain ID's, what does phenix accept besides the alphabet and 0-9? c) For P3, with any of the 3 twin operators the result is practically the same. To my knowledge, only shelxl can deal, currently, with more that one twinned domain, but these can only be produced in a iterative fashion. So, what might be the other way round? Is this development expected for phenix? For P1, I would have 11 pseudo-merohedral twin operations... I think the kind of packing is quite prone to twinning. d) Concerning what is raised in c, how necessary is to combine the 3 possible twin operators with the 3 reindexing possibilities.... e) How can I convince phenix to use the suggested twin operation 1/2*h-1/2*k,3/2*h+1/2*k,l in C222 (and some others with the * symbol in C2)? Much more information is available; I think I could pick the main ones. Thanks for any help. Jorge

2 2

Error when specifying anomalous scatterers in autosol.
by Francis E Reyes 27 Jun '08

27 Jun '08

------------------------------------------------------------ Starting AutoSol with the command: phenix.autosol data=ScalAveragedSG_76_30.ref seq_file=seq_from_pdb.dat \ refine_eff_file_list=./refinement/non_default.eff sites=2 atom_type=I \ f_prime=-0.582314 f_double_prime=6.84299 res_phase=3.4 resolution_build=3.0 \ resolution=3.0 wavelength=1.5418 ha_sites_file=./ha_sites.pdb Sending output to AutoSol_run_7_/AutoSol_run_7_1.log ******************************************************************************** Failed to carry out AutoSol_build: Refinement failed...perhaps something wrong with input refinement file or data file or the column labels for them? Error message from phenix.refine: Empty atom selection: refinement.refine.anomalous_scatterers.group.selection="name I" ******************************************************************************** Cuts out in the first round of refinement (after initial Build_1.pdb) . I guess I should not have specified anomalous scatterers (at least to autosol). Does it refine anomalous automagically in the autobuild when doing a SAD via autosol? How do I recover from this error? (continue with the autobuild without anomalous scattering it now seems). Thanks FR --------------------------------------------- Francis 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 3

Twinnig
by Michael Gorman 27 Jun '08

27 Jun '08

Dear all, After over 15 years solving structures I have (knowingly) encountered my first twinned data set. Short history.. Processed synchrotron data with HKL2000 SG P 42 21 2 ( Cell 105.485, 105.485, 62.445, 90, 90, 90) Rmerge 6% to 2A Solved MR using Phaser ( 1 mol in ASU) and built using COOT Refined using Phenix to a 'final' Rwork 29.4% Rfree 31.7% The maps look good but not great. Spurious peaks in difference maps. Couldn't get R's down any further. Xtriage suggested 'Overmerging pseudo-symmetric or twinned data' Re-processed in P2 ( cell 105.299, 62.385, 105.548, 90, 90.021, 90) Two molecules in ASU Xtriage suggested the presence of 3 pseudo-merohedral twin operators (one 4-fold and two 2-folds) Refined with twin law h,-k,-l (one of the pseudo 2-folds) using Phenix Rwork 26% Rfree 28% (no further building carried out at the moment) Wondered if it was a higher space group so.... re-processed in P222 (again Rmerge 6% to 2A) and did MR with Phaser in all the alternative space groups. The highest LL-gain was for P 21 21 21 (cell 62.377 105.297 105.546 90 90 90) I re-processed the data in P 21 21 21, all the systematic absences looked 'absent'. Xtriage suggested the presence of only one 4-fold Pseudo-merohedral twin law h,l,-k, R obs 5.1% Simulated annealing and ADP refinement in Phenix with twin law h,l,-k gave Rwork 33% Rfree 40%. Here are my questions... 1) Should I go for P2 due to the low R-factors and build from my good model ( beta angle of 90.021 seems unlikely ) 2) Should I try other possible SG ( HKL2000 suggested, primitive tetragonal, primitive orthorhombic, C centred orthorhombic, and monoclinic). 3) I'm concerned that my free-R reflections will be screwed-up since I have pseudo-merohedral twinning and some Rwork reflections may be the same as Rfree ones. I have read that Phenix is rather clever and chooses the highest observed symmetry to ovoid this problem but I can't see it being reported in the log-file that, in the case of P2, pseudo P42 21 2 is used for Rfree calculations. All help is much appreciated Mike -- Michael Gorman PhD St. Vincent's Institute of Medical Research 9 Princes Street Fitzroy Victoria 3065 Australia Tel : (03) 9288 2519

2 1

Re: [phenixbb] Rfree increases during ordered solvent filtering
by Jennifer L Ekstrom 25 Jun '08

25 Jun '08

Hi Pavel, Thanks. I don't believe that is the cause of this refinement anomaly. If I run the same job without the filter_only option, I am still seeing a very large jump in R/Rfree at the ordered solvent location/analysis/update step (before 2_xyz). ordered_solvent=true 0 : 0.1564 0.1709 2.905699e+00 3.076326e+00 1_bss: 0.1145 0.1431 2.837441e+00 3.027552e+00 1_xyz: 0.1159 0.1433 2.846219e+00 3.029651e+00 1_adp: 0.1159 0.1433 2.846219e+00 3.029651e+00 1_occ: 0.1155 0.1429 2.844726e+00 3.028603e+00 2_bss: 0.1154 0.1427 2.844702e+00 3.028393e+00 2_xyz: 0.1313 0.1555 2.950732e+00 3.111378e+00 2_adp: 0.1307 0.1543 2.945031e+00 3.101981e+00 2_occ: 0.1306 0.1542 2.945015e+00 3.101954e+00 3_bss: 0.1301 0.1538 2.944930e+00 3.101580e+00 3_xyz: 0.1312 0.1565 2.951306e+00 3.117549e+00 3_adp: 0.1306 0.1546 2.940355e+00 3.100448e+00 3_occ: 0.1305 0.1546 2.940347e+00 3.100439e+00 3_bss: 0.1292 0.1534 2.940856e+00 3.100060e+00 However, if I run this job with no ordered solvent manipulation, the R/Rfree are stable at their previous lower value. ordered_solvent=false 0 : 0.1564 0.1709 2.905699e+00 3.076326e+00 1_bss: 0.1145 0.1431 2.837441e+00 3.027552e+00 1_xyz: 0.1159 0.1433 2.846219e+00 3.029651e+00 1_adp: 0.1159 0.1433 2.846219e+00 3.029651e+00 1_occ: 0.1155 0.1429 2.844726e+00 3.028603e+00 2_bss: 0.1154 0.1427 2.844702e+00 3.028393e+00 2_xyz: 0.1151 0.1427 2.842538e+00 3.027628e+00 2_adp: 0.1151 0.1427 2.842538e+00 3.027628e+00 2_occ: 0.1151 0.1427 2.842433e+00 3.027654e+00 3_bss: 0.1151 0.1426 2.842430e+00 3.027578e+00 3_xyz: 0.1156 0.1428 2.846107e+00 3.029162e+00 3_adp: 0.1156 0.1428 2.846106e+00 3.029160e+00 3_occ: 0.1156 0.1429 2.846122e+00 3.029225e+00 3_bss: 0.1156 0.1429 2.846122e+00 3.029225e+00 I wouldn't expect to see a 1% jump in Rfree after minor tweaking of a few dozen waters. This is not a small structure. There must be something strange in the refinement protocol during the ordered water treatment. I had thought it could be a weighting issue, but subsequent refinement with optimization of wxc and wxu does not bring the R/Rfree down to the original range. Thanks, Jennifer (Running version 2008_05_03_2330) ................................................... Pavel Afonine writes: > Hi Jennifer, > > you are using "ordered_solvent.mode=filter_only". This means that the new > waters are not added but the existing ones are filtered (removed) at final > step based on selection criteria (B-factors, distances, H-bonding, > etc...). Do you really want this? If not, then do not use > "ordered_solvent.mode=filter_only" (the default is "auto" which is good in > most of cases). > > Please let me know if you have any questions! > > Pavel. > > > On 6/24/2008 2:42 PM, Jennifer L Ekstrom wrote: >> I am running phenix.refine (version 2008_05_03_2330) with the >> command-line options ordered_solvent=true >> ordered_solvent.mode=filter_only >> >> During the final step (3_bss), my R and Rfree increase significantly. >> Any thoughts on why this jump would occcur? Per the notations in the >> logfile, only 30 waters were removed in the final model. >> >> Thanks, Jennifer >> >> >> R-factors, x-ray target values and norm of gradient of x-ray target >> stage r-work r-free xray_target_w xray_target_t >> 0 : 0.1564 0.1709 2.905699e+00 3.076326e+00 >> 1_bss: 0.1145 0.1431 2.837441e+00 3.027552e+00 >> 1_xyz: 0.1159 0.1433 2.846219e+00 3.029651e+00 >> 1_adp: 0.1159 0.1433 2.846219e+00 3.029651e+00 >> 1_occ: 0.1155 0.1429 2.844726e+00 3.028603e+00 >> 2_bss: 0.1154 0.1427 2.844702e+00 3.028393e+00 >> 2_xyz: 0.1151 0.1427 2.842538e+00 3.027628e+00 >> 2_adp: 0.1151 0.1427 2.842538e+00 3.027628e+00 >> 2_occ: 0.1151 0.1427 2.842433e+00 3.027654e+00 >> 3_bss: 0.1151 0.1426 2.842430e+00 3.027578e+00 >> 3_xyz: 0.1156 0.1428 2.846107e+00 3.029162e+00 >> 3_adp: 0.1156 0.1428 2.846106e+00 3.029160e+00 >> 3_occ: 0.1156 0.1429 2.846122e+00 3.029225e+00 >> 3_bss: 0.1306 0.1512 2.943092e+00 3.079441e+00 >> >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://www.phenix-online.org/mailman/listinfo/phenixbb >> >

2 1