- phenixbb - phenix-online.org

phenix version 1.4-3: problem with bulk solvent scaling
by Oliv Eidam 31 Jan '09

31 Jan '09

Hi, I am currently refining a good 1.6 Angstroem data set using Phenix. Thereby I noticed a significant difference of R-free (>3%) depending on whether I used phenix version 1.3-rc6 (installed on my Mac laptop) or phenix version 1.4-3 (installed on my linux work station). Rfree differs by more than 3% (!) although I am using the same starting model, the same reflection file and the same refinement protocol (shown at the end of this email with detailed R and R-free factors for the 3 steps). In my opinion, my phenix 1.4-3 installation has a problem in the "bulk solvent modeling and scaling" steps. Indeed, I find in the log file that the "overall anisotropic scale matrix" shows all values set to zero. Analogous the "maximum likelihood estimate for coordinate error". This is true in all steps where an "overall anisotropic scale matrix" is calculated. Here is one example during 1.4-3 refinement: ====================== bulk solvent modeling and scaling ====================== |--(resolution: 1.61 - 29.59 A; n_refl. = 99535)------------------------------| | | | r_work= 0.3419 r_free= 0.3358 ksol= 0.39 Bsol= 19.12 scale= 1.247 | | | | overall anisotropic scale matrix (Cartesian basis; B11,B22,B33,B12,B13,B23):| | (0.00,0.00,0.00,0.00,0.00,0.00); trace/3= 0.00 | | | | maximum likelihood estimate for coordinate error: -0.00 A | | x-ray target function (ml) for work reflections: 5.981207 | |-----------------------------------------------------------------------------| The same example during 1.3-rc6 refinement: ===================== bulk solvent correction and scaling ===================== |-macro cycle = 1-(resolution: 1.61 - 29.59 A; n_refl. = 99537)---------------| | | | r_work= 0.3153 r_free= 0.3064 ksol= 0.40 Bsol= 43.79 scale= 1.087 | | | | overall anisotropic scale matrix (Cartesian basis; B11,B22,B33,B12,B13,B23):| | (-4.01,-5.36,-8.20,-0.00,2.74,0.00); trace/3= -5.86 | | | | maximum likelihood estimate for coordinate error: 0.33 A | | x-ray target function (ml) for work reflections: 5.946583 | |-----------------------------------------------------------------------------| I am asking now: Did something go wrong during my Phenix installation? I already re-installed Phenix (no errors): same result. Do I have to change something in my input files? I should probably mention that "bulk_solvent_and_scale = True" in my input .def files. Any help appreciated. Many thanks, Oliv -------------------------------------------------------------------------------- Refinement protocol and development of R and Rfree during each refinement round: -------------------------------------------------------------------------------- 1)Rigid body refinement 2) Simulated annealing (5000K) 3) Ordered solvent (plus annealing @ 2000K) phenix version 1.3-rc6: 1) R=31.5, Rf=30.6 2) R=23.9, Rf=25.8 3) R=19.8, Rf=22.0 phenix version 1.4-3: 1) R=31.6, Rf=30.6 2) R=27.8, Rf=29.8 3) R=23.7, Rf=25.2 -- Oliv Eidam, Ph.D. Postdoctoral fellow University of California, San Francisco Dept. of Pharmaceutical Chemistry 1700 4th Street, Byers Hall North, Room 501 San Francisco, CA 94158 - Box 2550 Phone: 415-514-4253 Fax : 415-514-4260 Email: eidamo(a)blur.compbio.ucsf.edu

3 6

joint xn refinement question
by Walter Novak 31 Jan '09

31 Jan '09

Hi, I am new to the list, but I didn't see this issue in the list archive. I am doing a joint refinement using x-ray and neutron data (each set is ~2.0 Å). I added H positions using reduce and have been switching to Ds as the density suggests for exchangeable positions. My problem is that I cannot seem to get the H (or D) positions to refine into the neutron density. Perhaps I need an additional command, so I am attaching my refinement script. I can provide a few pictures from coot to demonstrate if helpful. Thanks! Wally phenix.refine prot_joint4_r8_001.pdb \ neutron_data.file_name=prot_best_mcH_D2O_2_data.mtz \ xray_data.file_name=prot_best_mcH_D2O_2_data.mtz \ output.prefix=prot_joint4_r9 \ strategy=rigid_body+individual_sites+individual_adp \ simulated_annealing=False \ input.xray_data.labels=F-obs input.neutron_data.labels=F-obs-neutron \ input.xray_data.r_free_flags.label=R-free-flags \ input.xray_data.r_free_flags.file_name=prot_best_mcH_D2O_2_data.mtz \ input.xray_data.r_free_flags.test_flag_value=1 \ input.neutron_data.r_free_flags.label=R-free-flags-neutron \ input.neutron_data.r_free_flags.file_name=prot_best_mcH_D2O_2_data.mtz \ input.neutron_data.r_free_flags.test_flag_value=1 \ main.number_of_macro_cycles=5 \ restraints_edits.params \ refinement.hydrogens.refine=individual \ find_and_add_hydrogens=True \ switch_to_isotropic_high_res_limit=2.5 \ optimize_wxc=true \ optimize_wxu=true \ --overwrite Walter R.P. Novak, Ph.D. Postdoctoral Fellow Rosenstiel Basic Medical Research Center Brandeis University 415 South St. MS 029 Waltham, MA 02454-9110 Phone: (781) 736-4944 Fax: (781) 736-2405

2 1

Site-wide defaults for phenix.refine parameters
by Stephen Graham 30 Jan '09

30 Jan '09

Hi all, Apologies if I've missed this in the documentation, but is there any way to set 'site-wide' defaults for, say, phenix.refine? That is, I'd like to modify the default parameters for phenix.refine so that occupancy refinement is turned off by default, the output.serial_format is different, geo files aren't created, always use F/SIGF not I/SIGI, etc. Cheers, Stephen P.S. I know I can change this by hand in my .def file and then propagate the .def file through subsequent rounds of refinement for a given project. It would just be nice to have my default personal preference parameters set (~/.phenix?) so I don't need to change them each time I start a new project... -- 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

4 5

occupancy refinement
by Nicole Bischoff 30 Jan '09

30 Jan '09

Dear all, in my electron density I can see negative difference density (fo-fc) for the carboxyl groups of some aspartates or glutamates. Assuming that this is a result of partial decarboxylation due to radiation damage, I wanted to refine the occupancies of these carboxyl groups. To this end, I set the occupancy of these atoms to 0.99 in Coot so that they would be picked up by occupancy refinement in phenix.refine (using version 1.4-3). Phenix however just put the value back to one. It looks like the occupancy refinement in phenix is meant for alternative conformations where the sum of the occupancies of the different conformations is 1. How can I make phenix refine partially decarboxylated acidic side chains? It does not seem appropriate to me to make alternative conformations, one with and one without the carboxyl group. That would look funny at my resolution (2.7 Angstrom) and might not fix the problem. Any hints are appreciated. Thank you, Nicole ---- Nicole Bischoff, M. Sc. PhD student in the Gruetter lab Department of Biochemistry University of Zurich Winterthurerstrasse 190 8057 Zurich Switzerland +41 44 635 5556 nbischoff(a)bioc.unizh.ch

3 3

TLS question
by Jeffrey Chao 27 Jan '09

27 Jan '09

hello all- i apologize if this question has previously been posted, but i was not able to track down an answer by google. i am refining a 2.75A structure that has three copies of the protein in the asu with phenix 1.4-3 and using NCS with each chain as an NCS group. after molecular replacement, i set all of the b-factors to the wilson b-factor and then performed four different refinements: grouped b-factor (1 group per residue) : Rwork 22.6 Rfree 27.0 individual b-factor : Rwork 21.4 Rfree 26.6 grouped b + tls ( : Rwork 21.9 Rfree 25.2 individual + tls : Rwork 21.2 Rfree 24.8 i defined the tls groups with the following: adp { tls = chain A tls = chain B tls = chain C } after grouped b-factor the atoms of each residue all have the same b- factor as i expected. after grouped b-factor+tls, however, each atom of a residue has a unique b-factor, which i did not expect. i realize that this confusion is likely caused by me not fully understanding how tls refinement is implemented. if someone could explain this to me, i will be able to sleep better. thank you for your time, jeff Jeffrey Chao, PhD 1300 Morris Park Ave ASB-Golding 601 Bronx NY 10461 jchao(a)aecom.yu.edu 718-430-8597

2 1

rmsd: phenix vs refmac
by Frank von Delft 26 Jan '09

26 Jan '09

Hi I thought that a while ago, I read somewhere that rmsd(bonds) had not been calculated the same way in phenix as in refmac. I believe it was in connection with automated weights, and that rmsd(bonds) always much lower in phenix. I understood too that this had been adjusted since to match refmac. But we're still seeing some much lower rmsd(bond) from phenix: in particular, when we take a refmac model and plug it straight into phenix, the rmsd for cycle zero is already much lower than what it was at the end of refmac. Could someone enlighten...? (Yup, latest version of phenix.) Thanks! phx

5 5

bond and angle rmsds: phenix.refine, Refmac, hydrogens, ...
by Pavel Afonine 25 Jan '09

25 Jan '09

For all who asked about observed differences in bond rmsd reported by phenix.refine and Refmac, H atoms included in statistics calculation and more... 1) First, I wanted to check bond rmsd reported by phenix.refine. I made up a PDB file containing just one residue (SER, no H atoms) and I wrote a crystal-clear and simple Python script that computes bond rmsd for this residue. Here are the numbers: Bond length rmsd: Script: 0.148 phenix.pdbtools: 0.148 phenix.refine: 0.148 Refmac_5.5.0070: 0.152 Both programs are using the Monomer Library. Since both of these very different ways of computing bond rmsd: this simple script and phenix.refine, give the identical values, I tend to think that the PHENIX implementation is actually correct. I have no idea why Refmac gives a different value and I think it's now up to Refmac developers (or any other enthusiast) to find it out or show that I'm wrong in my calculations. It well could also be me who's incorrectly using Refmac - today was my second time I tried it (the first was 3 or 4 years ago) -:) As for the bond angles rmsd, I didn't bother to include it into the script, but for the rest of tools in question, the numbers are: phenix.pdbtools: 16.628 phenix.refine: 16.628 Refmac_5.5.0070: 16.628 so all results seem to be in-sync. 2) Some people reported that the geometry statistics for the exact same structure with and without hydrogens seems to be reported different in phenix.refine. To find out, I did the same test as above, where I used the same one residue model with all hydrogens added. Unfortunately, I failed to make Refmac accept a PDB file containing H atoms so I gave up after 20 minutes of struggle. Here are the numbers: Bond length rmsd: Script: 0.155 phenix.pdbtools: 0.155 phenix.refine: 0.155 So, it seems like the hydrogens are included... I have to fix it. Well at least those who are doing refinement using neutron data or both X-ray and neutron data, were not concerned, since in that case I presume the H and D atoms should be included. Anyway, I'm fixing it right now, I'm making a corresponding regression test (to make sure the fix stays there forever) and the fixed version will be available in the next release. By the way, the new tool, phenix.model_vs_data is and was smart about this: it automatically excludes H and D atoms if X-ray data is used. And this proves it: phenix.model_vs_data: 0.148 which is the exact same number as in "1)", where H were not present. Thanks to all who raised this issue (Frank von Delft, Jianghai Zhu, Karolina Michalska and others)! Pavel. PS> All files used and the script are here: http://cci.lbl.gov/~afonine/rmsd_issue/ Use "phenix.python run.py" to run the script.

1 0

IAS causes phenix.refine to die
by Anna Gardberg 25 Jan '09

25 Jan '09

Dear phenix list, My 0.75 A refinement has some problems in its higher resolution bins, so I thought I might try adding in IAS scatterers. However, the phenix.refine (version 2008_09_13_0905) crashed when I set "ias = True". I hope someone here can help me figure out the problem. I include snippets of output below, as the entire error message was very lengthy. Am I misusing the command, is there a problem with my system setup, or what? Regards, Anna Gardberg ----------Building IAS model---------- >>> Adding IAS.......... Total number of covalent bonds = 895 show_stack(1): /home/a1g/downloads/cci_apps_sources/mmtbx/ias/__init__.py(438) _ _init__ show_stack(2): /home/a1g/downloads/cci_apps_sources/mmtbx/ias/__init__.py(371) s et_peaks show_stack(3): /home/a1g/downloads/cci_apps_sources/mmtbx/ias/__init__.py(616) _ _init__ show_stack(4): /home/a1g/downloads/cci_apps_sources/mmtbx/model.py(494) add_ias show_stack(5): /home/a1g/downloads/cci_apps_sources/phenix/phenix/refinement/str ategies.py(409) refinement_machine show_stack(6): /home/a1g/downloads/cci_apps_sources/phenix/phenix/refinement/dri ver.py(1082) run <snip pages of output> /home/a1g/downloads/cci_apps_build/lib/cctbx_maptbx_ext.so(boost::python::obje cts::make_holder<4>::apply<boost::python::objects::value_holder<cctbx::maptbx::o ne_gaussian_peak_approximation>, boost::mpl::vector4<scitbx::af::const_ref<doubl e, scitbx::af::trivial_accessor> const&, scitbx::af::const_ref<double, scitbx::a f::trivial_accessor> const&, bool const&, bool const&> >::execute(_object*, scit bx::af::const_ref<double, scitbx::af::trivial_accessor> const&, scitbx::af::cons t_ref<double, scitbx::af::trivial_accessor> const&, bool const&, bool const&)+0x 75c) [0x2a9d79598c] /home/a1g/downloads/cci_apps_build/lib/cctbx_maptbx_ext.so(cctbx::maptbx::find _gaussian_parameters::find_gaussian_parameters(scitbx::af::const_ref<double, sci tbx::af::trivial_accessor> const&, scitbx::af::const_ref<double, scitbx::af::tri vial_accessor> const&, double const&, double const&, double)+0x3fe) [0x2a9d7827f e] /lib64/tls/libm.so.6 [0x32ee506263] /lib64/tls/libc.so.6 [0x32ee22e300] Floating-point error (Python and libc call stacks above) [1]+ Done gedit newscale-ias.eff

3 2

Unmasking cavities
by Frank von Delft 24 Jan '09

24 Jan '09

Hi After being through phenix.refine, I see in my hydrophobic core a big space (a few atoms wide) that is filled with strong negative difference density. I suspect the culprit is the bulk solvent mask, which is defined too tightly. The online manual mentions three parameters, but not what they do. solvent_radius, shrink_truncation_radius, grid_step_factor What *exactly* do they do? (I thought I'd elicit a contribution for the online docs this way :) Cheers phx

6 11

geometry of disulfide and NAG
by Li-Zhi Mi 23 Jan '09

23 Jan '09

Hi, I am using phenix to refine a low resolution structure. The problem I had is that the geometry of disulfide and NAG in the refined structure is not right. I did use a very tight geometry restrains. Both wxc_scale and wxu_scale are set to 0.05. In the refined structures, the sulfur of one cys is very close to the Cb of its disulfide partner. For NAG part, C2-N2-C7-O7-C8 should be in the same plane. But they are not. Is there a way to adjust these violations? Thanks in advance for the helps, Li-Zhi Mi,

7 9

HL A B C D
by Raja Dey 22 Jan '09

22 Jan '09

Hello, If I want to use F, SigF and HLA, B, C and D in the phenix refinement, should I use the same command like phenix.refine data.mtz model.pdb, where data.mtz contains FREER, F, SIGF, HLA, HLB, HLC and HLD? How do I check that phenix.refine used the HL coefficients? Looking forward to hearing from you... Thanks... Raja

2 1

Disappearing side chains after autobuilding
by Katya Heldwein 22 Jan '09

22 Jan '09

Hi, I ran phenix.autobuild with default setting (rebuild_in_place), and noticed that some side chains were gone in the resulting overall_best.pdb. That is, let's say, I started with Phe as residue 10. After rebuilding, the Phe10 is still called Phe in the output file, but the side chain atoms are gone, leaving residue 10 as Ala. This happened for several side chains, and I am at a loss for why this happened. Is this some kind of a bug? Katya -- Ekaterina Heldwein, Ph.D. Assistant Professor Department of Molecular Biology and Microbiology Tufts University School of Medicine 136 Harrison Ave Boston, MA 02111 phone: 617-636-0858 fax: 617-636-0337 e-mail: katya.heldwein(a)tufts.edu

2 1

missing map coefficients
by Anne Mølgaard 16 Jan '09

16 Jan '09

I am new to Phenix, so there might be an easy explanation for this problem, though I have not been able to find a solution in the documentation. I am using phenix version 1.3b-rc4, intel-linux-2.6. I have solved a structure using phenix.automr, followed by phenix.autobuild (run from command line using "phenix.autobuild data=fromXDS_unique.mtz model=AutoMR_run_2_/MR.1.pdb seq_file=sequence.fas"). The program runs succesfully, and I get the output files I should expect. However, the overall_best_refine_map_coeffs.mtz does not contain the map coefficients 2FOFCWT PH2FOFCWT for the 2mFo-DFc map and FOFC and PHFOFC for the Fo-Fc difference map. mtzdmp gives me the following columns: 1 ASC 0 28 0 100.00 10.6 10.6 28.28 1.85 H H 2 NONE 0 35 0 100.00 13.3 13.3 28.28 1.85 H K 3 NONE 0 71 0 100.00 26.6 26.6 28.28 1.85 H L 4 NONE 1.0 1888.8 0 100.00 148.71 148.71 28.28 1.85 F FP 5 BOTH 1.0 1.0 0 100.00 1.00 1.00 28.28 1.85 Q SIGFP 6 NONE -180.0 180.0 0 100.00 2.08 90.23 28.28 1.85 P PHIM 7 BOTH 1.000 1.000 0 100.00 1.000 1.000 28.28 1.85 W FOMM 8 NONE .0 19.0 0 100.00 9.50 9.50 28.28 1.85 I FreeR_flag 9 BOTH .0 .0 0 100.00 .00 .00 28.28 1.85 F FC What should I do to get these map coefficients? Thanks for your help! Best regards, Anne

2 2

Will phenix.refine read CNS *.cv files?
by Scott Classen 15 Jan '09

15 Jan '09

Hi all, I'm embarking on a collaboration with someone who is providing me with CNS cross validated data files. They look something like this. NREFlection= 92086 ANOMalous=FALSe { equiv. to HERMitian=TRUE} DECLare NAME=FOBS DOMAin=RECIprocal TYPE=REAL END DECLare NAME=SIGMA DOMAin=RECIprocal TYPE=REAL END DECLare NAME=IOBS DOMAin=RECIprocal TYPE=REAL END DECLare NAME=SIGI DOMAin=RECIprocal TYPE=REAL END DECLare NAME=TEST DOMAin=RECIprocal TYPE=INTE END INDE 3 0 0 FOBS= 100.142 SIGMA= 1.478 IOBS= 10028.400 SIGI= 293.800 TEST= 0 INDE 4 0 0 FOBS= 158.996 SIGMA= 2.329 IOBS= 25279.600 SIGI= 735.200 TEST= 0 INDE 5 0 0 FOBS= 492.449 SIGMA= 17.677 IOBS= 0.242506E+06 SIGI= 17098.000 TEST= 0 INDE 7 0 0 FOBS= 246.492 SIGMA= 3.127 IOBS= 60758.400 SIGI= 1531.800 TEST= 0 INDE 8 0 0 FOBS= 32.790 SIGMA= 0.549 IOBS= 1075.200 SIGI= 35.700 TEST= 1 INDE 9 0 0 FOBS= 426.866 SIGMA= 8.886 IOBS= 0.182215E+06 SIGI= 7507.000 TEST= 0 INDE 10 0 0 FOBS= 411.347 SIGMA= 8.597 IOBS= 0.169206E+06 SIGI= 6999.000 TEST= 0 INDE 11 0 0 FOBS= 87.275 SIGMA= 1.346 IOBS= 7617.000 SIGI= 233.200 TEST= 0 When I try to run phenix.refine --dry_run my_model.pdb my_CNS_cv_file.cv I get the following error: Processing inputs. This may take a minute or two. Sorry: Unknown file format: my_CNS_cv_file.cv Any ideas? Thanks, Scott ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Scott Classen, Ph.D. SIBYLS Beamline 12.3.1 http://bl1231.als.lbl.gov Advanced Light Source Lawrence Berkeley National Laboratory 1 Cyclotron Rd MS6R2100 Berkeley, CA 94720 C) 510.206.4418 O) 510.495.2697 Beamline) 510.495.2134 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

2 3

Re: [phenixbb] Will phenix.refine read CNS *.cv files?
by Scott Classen 15 Jan '09

15 Jan '09

Hi Nick, I think I tried that at some point, but the problem lies in the IOBS. Some are regular float fortran format: INDE 3 0 0 FOBS= 100.142 SIGMA= 1.478 IOBS= 10028.400 SIGI= 293.800 TEST= 0 while other IOBS are large and for some reason CNS decided to put them in scientific notation such as INDE 5 0 0 FOBS= 492.449 SIGMA= 17.677 IOBS= 0.242506E+06 SIGI= 17098.000 TEST= 0 As far as I can tell this requires two sets of fortran format statements for f2mtz to work properly. format '(6X,3F5.0,6X,F10.3,7X,F10.3,6X,E15.3,/,24X,F10.3,6X,F10.0)' or format '(6X,3F5.0,6X,F10.3,7X,F10.3,6X,F10.3,/,24X,F10.3,6X,F10.0)' I'm playing around with with grep to separate the original cv file into reflections with IOBS in F10.3 format and those with IOBS in E15.3 format, convert with f2mtz, then CAD them together, but this is bloody cumbersome. I was hoping phenix.refine would just automagically figure it all out. Scott On Jan 14, 2009, at 5:35 PM, Noinaj, Nicholas (NIH/NIDDK) [F] wrote: > Scott, > > <not sure if anyone has remedied your problem yet> > > I don't know if you can read cns file directly into phenix.refine. > but > here is another option....convert the cns file to *.mtz using CCP4. > If > you use the CCP4i GUI --> Reflection Data utilities --> Convert > to/modify/extend MTZ. You can also retain the original flagged > reflections here also. You can then use the *.mtz file in > phenix.refine. > > If you need further help, just drop me a quick email. I can also > try it > here if you supply me with the data file, but this may be too 'top > secret' (lol). Best of Luck! > > > > > > Cheers, > Nick > > > > -----Original Message----- > From: Scott Classen [mailto:[email protected]] > Sent: Wednesday, January 14, 2009 8:06 PM > To: PHENIX user mailing list > Subject: [phenixbb] Will phenix.refine read CNS *.cv files? > > Hi all, > I'm embarking on a collaboration with someone who is providing me with > CNS cross validated data files. They look something like this. > > > NREFlection= 92086 > ANOMalous=FALSe { equiv. to HERMitian=TRUE} > DECLare NAME=FOBS DOMAin=RECIprocal TYPE=REAL END > DECLare NAME=SIGMA DOMAin=RECIprocal TYPE=REAL END > DECLare NAME=IOBS DOMAin=RECIprocal TYPE=REAL END > DECLare NAME=SIGI DOMAin=RECIprocal TYPE=REAL END > DECLare NAME=TEST DOMAin=RECIprocal TYPE=INTE END > INDE 3 0 0 FOBS= 100.142 SIGMA= 1.478 IOBS= 10028.400 > SIGI= 293.800 TEST= 0 > INDE 4 0 0 FOBS= 158.996 SIGMA= 2.329 IOBS= 25279.600 > SIGI= 735.200 TEST= 0 > INDE 5 0 0 FOBS= 492.449 SIGMA= 17.677 IOBS= > 0.242506E+06 > SIGI= 17098.000 TEST= 0 > INDE 7 0 0 FOBS= 246.492 SIGMA= 3.127 IOBS= 60758.400 > SIGI= 1531.800 TEST= 0 > INDE 8 0 0 FOBS= 32.790 SIGMA= 0.549 IOBS= 1075.200 > SIGI= 35.700 TEST= 1 > INDE 9 0 0 FOBS= 426.866 SIGMA= 8.886 IOBS= > 0.182215E+06 > SIGI= 7507.000 TEST= 0 > INDE 10 0 0 FOBS= 411.347 SIGMA= 8.597 IOBS= > 0.169206E+06 > SIGI= 6999.000 TEST= 0 > INDE 11 0 0 FOBS= 87.275 SIGMA= 1.346 IOBS= 7617.000 > SIGI= 233.200 TEST= 0 > > When I try to run phenix.refine --dry_run my_model.pdb > my_CNS_cv_file.cv > > I get the following error: > > Processing inputs. This may take a minute or two. > Sorry: Unknown file format: my_CNS_cv_file.cv > > > Any ideas? > > Thanks, > Scott > > > ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ > Scott Classen, Ph.D. > SIBYLS Beamline 12.3.1 > http://bl1231.als.lbl.gov > Advanced Light Source > Lawrence Berkeley National Laboratory > 1 Cyclotron Rd > MS6R2100 > Berkeley, CA 94720 > C) 510.206.4418 > O) 510.495.2697 > Beamline) 510.495.2134 > ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ > > > > > _______________________________________________ > phenixbb mailing list > phenixbb(a)phenix-online.org > http://www.phenix-online.org/mailman/listinfo/phenixbb ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Scott Classen, Ph.D. SIBYLS Beamline 12.3.1 http://bl1231.als.lbl.gov Advanced Light Source Lawrence Berkeley National Laboratory 1 Cyclotron Rd MS6R2100 Berkeley, CA 94720 C) 510.206.4418 O) 510.495.2697 Beamline) 510.495.2134 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

1 0

f_prime and f_double_prime in MIRAS/SIRAS
by Leo Sazanov 15 Jan '09

15 Jan '09

Hi, how can one enter f_prime and f_double_prime values for SIRAS/MIRAS experiments in phenix.autosol? Could not find any examples of that in manual, only for MAD. Just putting them in as f_prime=... and f_double_prime=... does not work, it then assumes MAD experiment. Thanks! -- Dr. Leonid A. Sazanov Research group leader Medical Research Council Dunn Human Nutrition Unit Wellcome Trust / MRC Building Hills Road Cambridge CB2 0XY WEB: www.mrc-dunn.cam.ac.uk Tel: +44-1223-252910 Fax: +44-1223-252915

3 3

Jan 15th deadline - User proposal submission for Collaborative Crystallography at BCSB
by Banumathi Sankaran 14 Jan '09

14 Jan '09

Dear Users, /Deadline for March/April 2009 proposals will be Jan 15, 2009. / Through the Collaborative Crystallography Pilot Program (CC) at the Advanced Light Source, scientists can send protein crystals to Berkeley Center for Structural biology (BCSB) staff researchers for data collection and analysis. The CC Pilot Program can provide an number of benefits to researchers: * Obtain high quality data and analysis through collaborating with expert beamline researchers; * Rapid turn around on projects; * Reduced travel costs. Thumbnail Description: ====================== The Collaborative Crystallography program will be piloted on beamlines 5.0.1 and 5.0.2 for one year; if successful it will be implemented on a permanent basis. CC proposals will go through the regular ALS General User proposal review process for beamtime allocation. Proposals will be reviewed and ranked by the Proposal Study Panel, and beamtime will be allocated accordingly. BCSB staff will schedule the projects to fit into the available resources. Only non-proprietary projects will be accepted. As a condition of participation, BCSB staff researchers who participate in data collection and/or analysis must be appropriately acknowledged - typically being included as authors on publications and in PDB depositions. Please consult the website for additional information at: http://bcsb.als.lbl.gov/wiki/index.php/Collaborative_Crystallography ============= How To Apply: To submit a proposal, go to: http://www-als.lbl.gov/als/quickguide/independinvest.html Click on "protein crystallography" to begin proposal. For question 3 select "Collaborative PX 501/502." For question 9, please describe a specific research project with a clear end point. In order to request CC time for March/April 2009 allocation period, proposals must be submitted by January 15, 2009. The deadline for CC proposals for the time period May/June 2009 will be March 15, 2009. Regards, Banumathi Sankaran

1 0

Fwd: Refinement using new Phenix causes electron density/model mismatch
by Miles Pufall 13 Jan '09

13 Jan '09

The bb bounced my email back because the enclosures were too big - let me know if you'd like them. Thanks - Miles Begin forwarded message: > From: Miles Pufall <mpufall(a)cmp.ucsf.edu> > Date: January 12, 2009 12:16:06 PM PST > To: phenixbb(a)phenix-online.org > Subject: Refinement using new Phenix causes electron density/model > mismatch > > Hi Guys - > > I am having a very puzzling problem using the new Phenix (1.4-3). I > have about 20 structures that I am preparing to deposit with an > accepted manuscript, and decided to do a final refine using the > phenix.refine_gui not only to see if I got better numbers, but to > use the integrated validation tools, which I love. Unfortunately, > often, but not always, when I refine the pdb against the > *_refine_data.mtz the statistics are a bit worse, and when I open > the resulting *.map_coeffs.mtz in Coot it appears to be shifted with > tons of negative electron density. Needless to say, this is making > my nervous, so I've tried a few things: > > 1) Tried refining the same .pdb and .mtz files in Refmac - I seem > to get the same shifting. > 2) Try using the old .def file left over from the final run (just > to make sure I didn't screw up some setting) - and I get the same > shifting. > 3) Try refining against the original .sca file - and I get the same > shifting. > 4) Try molecular replacement of the refined .pdb against the > original .sca file - and sometimes, but not always, get a shift. > > I am nervous and concerned - can you help me figure out what > happened? I am enclosing an example - > > refine_001.pdb = pdb I refined in new phenix > refine_001map_coeffs.mtz = mtz output from new phenix refinement > (lots of red) > oste44a_output.sca = original scale file from HKL2000 > MR.1_refine_data.mtz = data from which I am refining > sgk_ost_temp_H_refine_004.pdb > > Thanks - > > Miles > > Miles Pufall > Postdoctoral Scholar > Yamamoto Lab > UC San Francisco > Cellular and Molecular Pharmacology > Mail Stop 2280 > 600 16th Street, Genentech Hall S-574 > San Francisco, California 94158-2517 > (415)476-4480 > > Miles Pufall Postdoctoral Scholar Yamamoto Lab UC San Francisco Cellular and Molecular Pharmacology Mail Stop 2280 600 16th Street, Genentech Hall S-574 San Francisco, California 94158-2517 (415)476-4480

2 1

Deadline for ALS PX beamtime this Thursday, Jan 15.
by Peter Zwart 13 Jan '09

13 Jan '09

Dear All, The deadline for general user proposals for PX beamtime at the ALS is coming up soon (Jan 15)! At the BCSB, we have 5 beamlines for which you can apply for beamtime. At the Berkeley Center for Structural Biology (BCSB) we have recently completed a number of upgrades: Sector 8: Recent improvements include: * A high-accuracy microdiffractometer on BL 8.2.1, allowing for centering and viewing of very small crystals * Robotic automounters which can increase the screening throughput at the beamlines, and allows for the implementation of remote data collection. * A new large format (315mm active surface area) CCD detector on BL 8.2.1, facilitating data collection on large unit cell crystals Sector 5: Beamlines (5.0.1 and 5.0.3) now generate 1.5x10^11 photons per second at 0.97 Angstrom. This wavelength is slightly above the Se K edge making them ideal for performing Se SAD experiments. Beamline 5.0.2 (the tunable beamline) generates a peak flux of 8x10^11 photons per second, and the energy range of 5.0.2 has been extended to 17 keV, enabling the routine use of shorter wavelengths and anomalously scattering elements such as bromine. Remote data collection for general users on sector 5 is currently being tested for a limited number of users. Please contact beamline staff for details. Please visit http://bcsb.lbl.gov/ for more details about the Center and its beamlines. If you'd like to apply for next round of beamtime (March-April) at the Advanced Light Source, please submit a General User proposal by Januari 15, 2009. To find out more, click on: http://www-als.lbl.gov/als/quickguide/independinvest.html We invite you to submit a proposal at: http://alsusweb.lbl.gov/4DCGI/WEB_GetForm/PXProposalEntry.shtml/Initialize If you have any questions or would like to request open beamtime, please e-mail BCSBBeamtime(a)lbl.gov. (Please note that executed user agreements must be received by LBNL prior to beamtime. Proprietary fees, if applicable, must be received by LBNL at least five working days prior to scheduled beamtime.) ----------------------------------------------------------------- P.H. Zwart Beamline Scientist Berkeley Center for Structural Biology Lawrence Berkeley National Laboratories 1 Cyclotron Road, Berkeley, CA-94703, USA Cell: 510 289 9246 BCSB: http://bcsb.als.lbl.gov PHENIX: http://www.phenix-online.org CCTBX: http://cctbx.sf.net -----------------------------------------------------------------

1 0

Re: [phenixbb] TLS in phenix.refine_gui
by Pavel Afonine 11 Jan '09

11 Jan '09

Hi Miles, > Thanks for the info. I see where I can load a parameter file in the > gui (File => Open parameter file) but I don't think I have it with the > right extension (.eff or .param). I'll give it a whirl tomorrow. Oh good, I'm glad to hear that this exists! I'm not sure if it relies on the file extension (it wouldn't be right if it does). I think we should have this option (to load a parameters file) right next to the prompt for PDB, data and CIF files. Best, Pavel.

1 0

Scratch last help request
by Miles Pufall 11 Jan '09

11 Jan '09

As it turns out there were some problems with my .pdb file - missing chain label somehow. Once I fixed it everything works - sorry for the trouble. Miles Miles Pufall Postdoctoral Scholar Yamamoto Lab UC San Francisco Cellular and Molecular Pharmacology Mail Stop 2280 600 16th Street, Genentech Hall S-574 San Francisco, California 94158-2517 (415)476-4480

1 0

Error running refine_gui
by Miles Pufall 10 Jan '09

10 Jan '09

Hi Guys - I was running the phenix_gui ad had just started a refinement run with tls, adp, occupancy, individual sites, simulated annealing, update waters and add hydrogen when the run stop and this error popped up: RuntimeError : mean_weighted() argument is an empty array Traceback: Traceback (most recent call last): File "/usr/local/phenix-1.4-3/phenix/wxGUI2/Refine.py", line 259, in OnExcept raise event.data RuntimeError: mean_weighted() argument is an empty array CHeers - M Miles Pufall Postdoctoral Scholar Yamamoto Lab UC San Francisco Cellular and Molecular Pharmacology Mail Stop 2280 600 16th Street, Genentech Hall S-574 San Francisco, California 94158-2517 (415)476-4480

1 0

Postdoc Crystallographer Position - Berkeley
by Paul Adams 10 Jan '09

10 Jan '09

Lawrence Berkeley Laboratory Postdoctoral Researcher Job ID: 22539 Division: Physical Biosciences Date Opened: 1/5/2009 Summary: A postdoctoral position is immediately available to study structure and function of macromolecular complexes and enzymes using X-ray crystallographic methods. This position is joint between the Center for Protein Folding Machinery (www.proteinfoldingcenter.org) and the Joint BioEnergy Institute (www.jbei.org). Robotic hardware for performing crystallization trials and imaging results are available. Crystals will be characterized and data collected using the beamline resources of the Berkeley Center for Structural Biology at the Advanced Light Source. The biological systems to be studied are: - Type-II chaperonins responsible for the refolding of unfolded proteins in archeal and mammalian cells. - Glycosyl transferases responsible for the synthesis of hemicellulose in plants. - Novel glycosyl hydrolases that are able to breakdown cellulosic material for biofuels production. Duties: Using robotic hardware to perform screens for crystallization conditions. Optimization of crystallization conditions. Biophysical characterization of protein samples using standard techniques, such as dynamic light scattering. Characterization of crystals using the Berkeley Center for Structural Biology beamlines at the Advanced Light Source. Collection and analysis of diffraction data, model building and structure refinement. Extensive email and verbal interaction with other researchers. Qualifications: Ph.D. or equivalent in a scientific discipline, preferably structural biology, biology, or chemistry. Demonstrated experience of scientific research using crystallographic methods. Experience with crystallization methods. Solid interpersonal skills and the ability to work in a team environment are critical. Ability to communicate with a broad range of researchers. Note that this is a one-year appointment with renewal based on performance and continued funding. Lawrence Berkeley National Laboratory is a world leader in science and engineering research, with 11 Nobel Prize recipients over the past 75 years, and 59 present members of the National Academy of Sciences. LBNL conducts unclassified research across a wide range of scientific disciplines and hosts four national user facilities. AA/EEO employer committed to the development of a safe and diverse workforce. Learn more at http://www.lbl.gov. To apply visit http://cjo.lbl.gov/ and search for the job number 22539 or visit: http://jobs.lbl.gov/LBNLCareers/details.asp?jid=22539&p=1 -- Paul Adams Deputy Division Director, Physical Biosciences Division, Lawrence Berkeley Lab Adjunct Professor, Department of Bioengineering, U.C. Berkeley Vice President for Technology, the Joint BioEnergy Institute Head, Berkeley Center for Structural Biology Building 64, Room 248 Tel: 1-510-486-4225, Fax: 1-510-486-5909 http://cci.lbl.gov/paul Lawrence Berkeley Laboratory 1 Cyclotron Road BLDG 64R0121 Berkeley, CA 94720, USA. Administrator: Kristine Sawyer [ KKSawyer(a)lbl.gov ][ 1-510-486-4151 ] --

1 0

phenix v1.4-3 and NX
by Jesper Lykkegaard Karlsen 10 Jan '09

10 Jan '09

Hi, After upgrading our workstations to phenix v1.4-3, the user started having problems connecting to the machines through nxclient. The login ended up with a black screen on the client and absolutely no error messages on neither client nor server. The problem seemed to come from setting up the global environment for phenix (phenix_env.sh). I narrowed the problematic statement down to this: /if [ "$PHENIX_MTYPE" = "unknown-unknown" ]; then exit 1 fi/ Although "$PHENIX_MTYPE" = "intel-linux-2.6-x86_64" when logged in, somehow the "nx user" gets "$PHENIX_MTYPE" = "unknown-unknown" when it logs in, through it's authorization key, to establish the users session. To make a already to long story short. The problem can be fixed by commenting out the problematic statement above. And one question to the developers: Could this be avoided in the future? Perhaps give a warning instead off exiting completly! Sincerely Jesper -- Jesper L. Karlsen Scientific Computing Centre for Structural Biology Department of Molecular Biology University of Aarhus Gustav Wieds Vej 10C 8000 Aarhus C Phone: +45 8942 5268 E-mail: jelka(a)bioxray.au.dk

2 1

Re: [phenixbb] phenix v1.4-3 and NX
by Jesper Lykkegaard Karlsen 10 Jan '09

10 Jan '09

Hi Paul, > what happens if you change the statement to say: > > exit 0 > > instead of exit 1? No did not help, still get black screen on the nxclient. Only out-commenting/deleting helps. Cheers, Jesper > > On Jan 9, 2009, at 3:52 AM, Jesper Lykkegaard Karlsen wrote: > >> Hi, >> >> After upgrading our workstations to phenix v1.4-3, the user started >> having problems connecting to the machines through nxclient. The >> login ended up with a black screen on the client and absolutely no >> error messages on neither client nor server. >> >> The problem seemed to come from setting up the global environment for >> phenix (phenix_env.sh). I narrowed the problematic statement down to >> this: >> >> if [ "$PHENIX_MTYPE" = "unknown-unknown" ]; then >> exit 1 >> fi >> >> Although "$PHENIX_MTYPE" = "intel-linux-2.6-x86_64" when logged in, >> somehow the "nx user" gets "$PHENIX_MTYPE" = "unknown-unknown" when >> it logs in, through it's authorization key, to establish the users >> session. >> >> To make a already to long story short. The problem can be fixed by >> commenting out the problematic statement above. >> >> And one question to the developers: >> >> Could this be avoided in the future? >> Perhaps give a warning instead off exiting completly! >> >> Sincerely >> Jesper >> >> -- Jesper L. Karlsen Scientific Computing Centre for Structural >> Biology Department of Molecular Biology University of Aarhus Gustav >> Wieds Vej 10C 8000 Aarhus C Phone: +45 8942 5268 E-mail: >> jelka(a)bioxray.au.dk_______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://www.phenix-online.org/mailman/listinfo/phenixbb >

1 0

fixing a partial solution of mol. rep.
by chern 07 Jan '09

07 Jan '09

I am doing a molecular replacement. There are two types of protein molecules in the crystal. One type is forming a hexameric core; there are 6 molecules of the second type that I need to locate. I found a solution for the hexameric core, but the other 6 molecules are difficult to find. 1. Can I fix the solution for the core, so that I don't have to find it every time. The other inconvenience with that is that the origine is different for every new run. If I fix the hexamer, then I fix the origine. It's important to me to compare different runs of MR. 2. Some of the six other molecules that MR finds are not the closest to the core. How can I request to output the closest molecules? Maia

3 2

TRIED resolve ...NEED A BIGGER VERSION... message
by Ben Eisenbraun 05 Jan '09

05 Jan '09

Howdy Phenixians, One of my users is running phenix.autobuild, and while it seems to run okay, the output contains this warning: TRIED resolve ...NEED A BIGGER VERSION... Looking through phenix/phenix/autosol/run_resolve.py, it seems like this is an informational message, but I'm not sure, so I figured I'd ask. Is this actually an error? This is phenix-1.3-final on OS X 10.5. Thanks. -ben -- Ben Eisenbraun Structural Biology Grid Harvard Medical School http://sbgrid.org http://hms.harvard.edu

3 3

question about xn joint refinement
by crystallogrphy 03 Jan '09

03 Jan '09

Deall all: I have a question about xray and neutron joint refinement. When I start phenix joint refinement like this, phenix gave the wrong message: phenix.refine dhfr-mtx-rt-xray.hkl neutron_data.file_name=newdhfr_un_change.hkl xn_easy_080326a.pdb input.xray_data.labels=FOBSx input.neutron_data.labels=FOBSn xray_data.r_free_flags.generate=True neutron_data.r_free_flags.generate=True refine_default.params ===================== joint neutron/xray refinement start ===================== Traceback (most recent call last): call_back_after_monitor_collect=call_back_after_monitor_collect) File "/usr/local/phenix-1.3-rc4/phenix/phenix/refinement/driver.py", line 1031, in run call_back_after_monitor_collect = call_back_after_monitor_collect) File "/usr/local/phenix-1.3-rc4/phenix/phenix/refinement/strategies.py", line 184, in refinement_machine * assert not neutron_refinement AssertionError* I guess in my default parameters, something is wrong, particularly the bold word I highlighted. But I do not know how to fix. Does any one have a clue?\ Thanks and Happy New Year! main { bulk_solvent_and_scale = True simulated_annealing = False ordered_solvent = False ncs = False ias = False number_of_macro_cycles = 3 max_number_of_iterations = 25 use_form_factor_weights = False tan_u_iso = False use_convergence_test = False target = *ml mlhl ml_sad ls min_number_of_test_set_reflections_for_max_likelihood_target = 50 max_number_of_resolution_bins = 30 reference_xray_structure = None use_experimental_phases = None compute_optimal_errors = False random_seed = 2679941 *scattering_table = wk1995 it1992 *n_gaussian neutron* use_normalized_geometry_target = True target_weights_only = False use_f_model_scaled = False max_d_min = 0.25 fake_f_obs = False -- Qun Wan Department of Pharmacology, Case Western Reserve University 10900 Euclid Avenue, Cleveland, OH 44106 email: qunwan1(a)gmail.com lab phone: 216-368-3337

3 2

Error
by Miles Pufall 28 Dec '08

28 Dec '08

Hi - It appears there was a problem with the NCS flag being on, and when I took it off I received this error. Let me know if you'd like data to recreate this. Cheers - Miles RuntimeError : OpenGL: invalid operation Traceback: Traceback (most recent call last): File "/usr/local/phenix-1.4-3/phenix/wxGUI2/Refine.py", line 188, in run self.check_params_and_run() File "/usr/local/phenix-1.4-3/phenix/wxGUI2/Refine.py", line 229, in check_params_and_run self.model.add_hydrogens() File "/usr/local/phenix-1.4-3/phenix/wxGUI2/Model.py", line 369, in add_hydrogens self.update_from_interface() File "/usr/local/phenix-1.4-3/phenix/wxGUI2/Model.py", line 181, in update_from_interface self.update_viewer() File "/usr/local/phenix-1.4-3/phenix/wxGUI2/Model.py", line 317, in update_viewer self .viewer_frame.glview.update_structure_from_files(self.get_all_files()) File "/usr/local/phenix-1.4-3/phenix/wxGUI2/GLViewer.py", line 166, in update_structure_from_files self.update_structure(pdb_hierarchy, atomic_bonds, redraw=True) File "/usr/local/phenix-1.4-3/phenix/wxGUI2/GLViewer.py", line 180, in update_structure self.OnUpdate(recenter=True) File "/usr/local/phenix-1.4-3/phenix/wxGUI2/GLViewer.py", line 208, in OnUpdate self.update_view(True, True) File "/usr/local/phenix-1.4-3/phenix/wxGUI2/GLViewer.py", line 566, in update_view self.OnRedraw() File "/usr/local/phenix-1.4-3/gltbx/wx_viewer.py", line 473, in OnRedraw self.OnRedrawGL(event) File "/usr/local/phenix-1.4-3/phenix/wxGUI2/GLViewer.py", line 118, in OnRedrawGL wx_viewer.wxGLWindow.OnRedrawGL(self, event) File "/usr/local/phenix-1.4-3/gltbx/wx_viewer.py", line 476, in OnRedrawGL gltbx.util.handle_error() RuntimeError: OpenGL: invalid operation Miles Pufall Postdoctoral Scholar Yamamoto Lab UC San Francisco Cellular and Molecular Pharmacology Mail Stop 2280 600 16th Street, Genentech Hall S-574 San Francisco, California 94158-2517 (415)476-4480

2 1

PHENIX 1.4-3 Available
by Paul Adams 24 Dec '08

24 Dec '08

= = = = = ======================================================================== PHENIX 1.4-3 Development Release = = = = = ======================================================================== The PHENIX developers are pleased to announce that the 1.4-3 development version of PHENIX is now available. Binary installers for Linux, and Mac OSX platforms are available at the download site: http://phenix-online.org/download/ This version of PHENIX includes the following new features: - A new GUI for phenix.xtriage: - To start the GUI version type phenix.xtriage_gui - A new GUI for phenix.refine: - To start the GUI version type phenix.refine_gui - Incorporates a subset of Molprobity validation tools - Documentation for the new GUIs: - See the help menu or button after starting the GUI - Structure refinement (phenix.refine) - Default is to fill missing reflections in map calculations with DFc - Default is to write phase combined maps if experimental phases are used in refinement - Improved capabilities for selection of atoms for occupancy refinement - Xplor maps no longer written by default - Numerous minor bug fixes - Multi-crystal averaging (phenix.multi_crystal_average) - Fast secondary structure building (phenix.find_helices_strands) - Find all ligands in a map (phenix.find_all_ligands) - Minor bug fixes in phenix.elbow and phenix.reel Full documentation is available here: http://www.phenix-online.org/documentation/ There is a PHENIX bulletin board: http://www.phenix-online.org/mailman/listinfo/phenixbb/ Please consult the installer README file or online documentation for installation instructions. Direct questions and problem reports to the bulletin board or: help(a)phenix-online.org and bugs(a)phenix-online.org Commercial users interested in obtaining access to PHENIX should visit the PHENIX website for information about the PHENIX Industrial Consortium. The development of PHENIX is principally funded by the National Institute of General Medical Sciences (NIH) under grant P01-GM063210. We also acknowledge the generous support of the members of the PHENIX Industrial Consortium. = = = = = ======================================================================== PHENIX 1.4-3 Development Release = = = = = ======================================================================== -- Paul Adams Deputy Division Director, Physical Biosciences Division, Lawrence Berkeley Lab Adjunct Professor, Department of Bioengineering, U.C. Berkeley Vice President for Technology, the Joint BioEnergy Institute Head, Berkeley Center for Structural Biology Building 64, Room 248 Tel: 510-486-4225, Fax: 510-486-5909 http://cci.lbl.gov/paul Lawrence Berkeley Laboratory 1 Cyclotron Road BLDG 64R0121 Berkeley, CA 94720, USA. --

1 0

Selection of structure factors for DM
by Graeme Winter 18 Dec '08

18 Dec '08

Hi Folks, If I am using phenix.autosol for my phasing and dm, is there a way I can specify the "native" FP to use for phase extension? What I would like to be able to do is use a named wavelength for this rather than a synthetic data set composed from all input data... I used to do this with resolve by cadding in the chosen columns (F, SIGF) after solve and then telling resolve what columns to use. Would be cute to be able to not do this! Am running from the command line interface. Thanks, Graeme

2 2

Re: [phenixbb] phenixbb Digest, Vol 37, Issue 3
by CNB 18 Dec '08

18 Dec '08

Thanks for the reply. I think I didn't explain myself properly, of course I want to do bss, but we have noticed that previous versions of phenix (1.24 I think) did bss once per macro cycle and now with the newer versions, we get the same plus a final bbs cycle at the end of the whole process -shown in red- Is that final cycle introduced by default in the newer versions? Thanks again. Cesar. REMARK 3_bss: 1.222 0.008 0.075 19.562 0.005 4.107 1.2876e-01 REMARK 3_xyz: 1.021 0.005 0.063 19.049 0.004 4.111 1.0306e-01 REMARK 3_adp: 1.021 0.005 0.063 19.049 0.004 4.111 1.0306e-01 REMARK 3_bss: 1.021 0.005 0.063 19.049 0.004 4.111 1.0306e-01 On Dec 3, 2008, at 9:00 PM, phenixbb-request(a)phenix-online.org wrote: > Send phenixbb mailing list submissions to > phenixbb(a)phenix-online.org > > To subscribe or unsubscribe via the World Wide Web, visit > http://www.phenix-online.org/mailman/listinfo/phenixbb > or, via email, send a message with subject or body 'help' to > phenixbb-request(a)phenix-online.org > > You can reach the person managing the list at > phenixbb-owner(a)phenix-online.org > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of phenixbb digest..." > > > Today's Topics: > > 1. Bulk solvent (csantiag(a)cnb.csic.es) > 2. Re: Bulk solvent (junfeng liu) > 3. Re: Bulk solvent (Pavel Afonine) > 4. Re: Bulk solvent (Pavel Afonine) > 5. Re: Question regarding phenix.elbow (Nigel W Moriarty) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 3 Dec 2008 10:15:30 +0100 (CET) > From: csantiag(a)cnb.csic.es > Subject: [phenixbb] Bulk solvent > To: phenixbb(a)phenix-online.org > Message-ID: > <49732.150.244.83.180.1228295730.squirrel(a)webmail.cnb.csic.es> > Content-Type: text/plain;charset=iso-8859-1 > > Hi all, > I'm refining a protein and after running phenix.refine I always get > a bulk > solvent refinement cycle at the end of the process. Is there any > way to > get rid of the last step? > > ---------------------------------------------------------------------- > -- > REMARK stage angl bond chir dihe plan repu > geom_target > REMARK 0 : 1.007 0.027 0.237 18.975 0.004 4.107 > 2.7953e-01 > REMARK 1_bss: 1.007 0.027 0.237 18.975 0.004 4.107 > 2.7953e-01 > REMARK 1_xyz: 1.844 0.015 0.087 19.502 0.005 4.106 > 1.9824e-01 > REMARK 1_adp: 1.844 0.015 0.087 19.502 0.005 4.106 > 1.9824e-01 > REMARK 2_bss: 1.844 0.015 0.087 19.502 0.005 4.106 > 1.9824e-01 > REMARK 2_xyz: 1.222 0.008 0.075 19.562 0.005 4.107 > 1.2876e-01 > REMARK 2_adp: 1.222 0.008 0.075 19.562 0.005 4.107 > 1.2876e-01 > REMARK 3_bss: 1.222 0.008 0.075 19.562 0.005 4.107 > 1.2876e-01 > REMARK 3_xyz: 1.021 0.005 0.063 19.049 0.004 4.111 > 1.0306e-01 > REMARK 3_adp: 1.021 0.005 0.063 19.049 0.004 4.111 > 1.0306e-01 > REMARK 3_bss: 1.021 0.005 0.063 19.049 0.004 4.111 > 1.0306e-01 > REMARK > ---------------------------------------------------------------------- > -- > > > C?sar Santiago > Dept of Macromolecules Structure > Centro Nacional de Biotecnolog?a, CSIC > Campus Cantoblanco > 28049 Madrid > Spain > > csantiag(a)cnb.csic.es > > > > > ------------------------------ > > Message: 2 > Date: Wed, 03 Dec 2008 09:36:30 +0000 > From: junfeng liu <jliu(a)nimr.mrc.ac.uk> > Subject: Re: [phenixbb] Bulk solvent > To: PHENIX user mailing list <phenixbb(a)phenix-online.org> > Message-ID: <4936531E.5030505(a)nimr.mrc.ac.uk> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > Add > > bulk_solvent_and_scale=False > when you run it again. > > > > csantiag(a)cnb.csic.es wrote: >> Hi all, >> I'm refining a protein and after running phenix.refine I always >> get a bulk >> solvent refinement cycle at the end of the process. Is there any >> way to >> get rid of the last step? >> >> --------------------------------------------------------------------- >> --- >> REMARK stage angl bond chir dihe plan repu >> geom_target >> REMARK 0 : 1.007 0.027 0.237 18.975 0.004 4.107 >> 2.7953e-01 >> REMARK 1_bss: 1.007 0.027 0.237 18.975 0.004 4.107 >> 2.7953e-01 >> REMARK 1_xyz: 1.844 0.015 0.087 19.502 0.005 4.106 >> 1.9824e-01 >> REMARK 1_adp: 1.844 0.015 0.087 19.502 0.005 4.106 >> 1.9824e-01 >> REMARK 2_bss: 1.844 0.015 0.087 19.502 0.005 4.106 >> 1.9824e-01 >> REMARK 2_xyz: 1.222 0.008 0.075 19.562 0.005 4.107 >> 1.2876e-01 >> REMARK 2_adp: 1.222 0.008 0.075 19.562 0.005 4.107 >> 1.2876e-01 >> REMARK 3_bss: 1.222 0.008 0.075 19.562 0.005 4.107 >> 1.2876e-01 >> REMARK 3_xyz: 1.021 0.005 0.063 19.049 0.004 4.111 >> 1.0306e-01 >> REMARK 3_adp: 1.021 0.005 0.063 19.049 0.004 4.111 >> 1.0306e-01 >> REMARK 3_bss: 1.021 0.005 0.063 19.049 0.004 4.111 >> 1.0306e-01 >> REMARK >> --------------------------------------------------------------------- >> --- >> >> >> C?sar Santiago >> Dept of Macromolecules Structure >> Centro Nacional de Biotecnolog?a, CSIC >> Campus Cantoblanco >> 28049 Madrid >> Spain >> >> csantiag(a)cnb.csic.es >> >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://www.phenix-online.org/mailman/listinfo/phenixbb >> >> > > > > ------------------------------ > > Message: 3 > Date: Wed, 03 Dec 2008 08:09:57 -0800 > From: Pavel Afonine <pafonine(a)lbl.gov> > Subject: Re: [phenixbb] Bulk solvent > To: PHENIX user mailing list <phenixbb(a)phenix-online.org> > Cc: csantiag(a)cnb.csic.es > Message-ID: <4936AF55.8010508(a)lbl.gov> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > Hi C?sar, > > this step is important. Bulk solvent modeling is based on mask > calculation and since during refinement of coordinates the atoms shift > the mask gets invalidated. This is why we do a final bulk solvent > correction step at the very last step of refinement just before the > final statistics is reported. So I'm puzzled about why you want to > turn > it off. If it causes any problem, it's much better to find out the > primary reason for it rather than turn bulk solvent correction. Could > you please give more details? > > If you do "main.bulk_solvent_and_scale=false" phenix.refine will > not do > bulk solvent correction and anisotropic scaling at all which is bad. > > Pavel. > > On 12/3/2008 1:15 AM, csantiag(a)cnb.csic.es wrote: >> Hi all, >> I'm refining a protein and after running phenix.refine I always >> get a bulk >> solvent refinement cycle at the end of the process. Is there any >> way to >> get rid of the last step? >> >> --------------------------------------------------------------------- >> --- >> REMARK stage angl bond chir dihe plan repu >> geom_target >> REMARK 0 : 1.007 0.027 0.237 18.975 0.004 4.107 >> 2.7953e-01 >> REMARK 1_bss: 1.007 0.027 0.237 18.975 0.004 4.107 >> 2.7953e-01 >> REMARK 1_xyz: 1.844 0.015 0.087 19.502 0.005 4.106 >> 1.9824e-01 >> REMARK 1_adp: 1.844 0.015 0.087 19.502 0.005 4.106 >> 1.9824e-01 >> REMARK 2_bss: 1.844 0.015 0.087 19.502 0.005 4.106 >> 1.9824e-01 >> REMARK 2_xyz: 1.222 0.008 0.075 19.562 0.005 4.107 >> 1.2876e-01 >> REMARK 2_adp: 1.222 0.008 0.075 19.562 0.005 4.107 >> 1.2876e-01 >> REMARK 3_bss: 1.222 0.008 0.075 19.562 0.005 4.107 >> 1.2876e-01 >> REMARK 3_xyz: 1.021 0.005 0.063 19.049 0.004 4.111 >> 1.0306e-01 >> REMARK 3_adp: 1.021 0.005 0.063 19.049 0.004 4.111 >> 1.0306e-01 >> REMARK 3_bss: 1.021 0.005 0.063 19.049 0.004 4.111 >> 1.0306e-01 >> REMARK >> --------------------------------------------------------------------- >> --- >> >> >> C?sar Santiago >> Dept of Macromolecules Structure >> Centro Nacional de Biotecnolog?a, CSIC >> Campus Cantoblanco >> 28049 Madrid >> Spain >> >> csantiag(a)cnb.csic.es >> >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://www.phenix-online.org/mailman/listinfo/phenixbb >> > > > ------------------------------ > > Message: 4 > Date: Wed, 03 Dec 2008 08:13:21 -0800 > From: Pavel Afonine <pafonine(a)lbl.gov> > Subject: Re: [phenixbb] Bulk solvent > To: PHENIX user mailing list <phenixbb(a)phenix-online.org> > Message-ID: <4936B021.7010401(a)lbl.gov> > Content-Type: text/plain; charset="iso-8859-1" > > Yes, this will turn bulk solvent correction and anisotropic scaling > off > completely. In some cases it may result in R-factors increase by > ~5% or > similar. So it may only make sense in exotic cases or if you do > numerical experiments (simulated data etc..). > > Pavel. > > > On 12/3/2008 1:36 AM, junfeng liu wrote: >> Add >> >> bulk_solvent_and_scale=False >> when you run it again. >> >> >> >> csantiag(a)cnb.csic.es wrote: >> >>> Hi all, >>> I'm refining a protein and after running phenix.refine I always >>> get a bulk >>> solvent refinement cycle at the end of the process. Is there any >>> way to >>> get rid of the last step? >>> >>> -------------------------------------------------------------------- >>> ---- >>> REMARK stage angl bond chir dihe plan repu >>> geom_target >>> REMARK 0 : 1.007 0.027 0.237 18.975 0.004 4.107 >>> 2.7953e-01 >>> REMARK 1_bss: 1.007 0.027 0.237 18.975 0.004 4.107 >>> 2.7953e-01 >>> REMARK 1_xyz: 1.844 0.015 0.087 19.502 0.005 4.106 >>> 1.9824e-01 >>> REMARK 1_adp: 1.844 0.015 0.087 19.502 0.005 4.106 >>> 1.9824e-01 >>> REMARK 2_bss: 1.844 0.015 0.087 19.502 0.005 4.106 >>> 1.9824e-01 >>> REMARK 2_xyz: 1.222 0.008 0.075 19.562 0.005 4.107 >>> 1.2876e-01 >>> REMARK 2_adp: 1.222 0.008 0.075 19.562 0.005 4.107 >>> 1.2876e-01 >>> REMARK 3_bss: 1.222 0.008 0.075 19.562 0.005 4.107 >>> 1.2876e-01 >>> REMARK 3_xyz: 1.021 0.005 0.063 19.049 0.004 4.111 >>> 1.0306e-01 >>> REMARK 3_adp: 1.021 0.005 0.063 19.049 0.004 4.111 >>> 1.0306e-01 >>> REMARK 3_bss: 1.021 0.005 0.063 19.049 0.004 4.111 >>> 1.0306e-01 >>> REMARK >>> -------------------------------------------------------------------- >>> ---- >>> >>> >>> C?sar Santiago >>> Dept of Macromolecules Structure >>> Centro Nacional de Biotecnolog?a, CSIC >>> Campus Cantoblanco >>> 28049 Madrid >>> Spain >>> >>> csantiag(a)cnb.csic.es >>> >>> >>> _______________________________________________ >>> phenixbb mailing list >>> phenixbb(a)phenix-online.org >>> http://www.phenix-online.org/mailman/listinfo/phenixbb >>> >>> >>> >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://www.phenix-online.org/mailman/listinfo/phenixbb >> >

2 1

B-group and TLS
by Engin Ozkan 17 Dec '08

17 Dec '08

Dear phenixbb people, I am performing group ADP refinement and TLS at the same time. At the end, the total B factors in the pdb are not identical within residue groups. ATOM 7 N ILE A 11 23.863 6.894 77.827 1.00151.81 A N ANISOU 7 N ILE A 11 20044 20245 17393 -796 -2448 -568 A N ATOM 8 CA ILE A 11 23.647 6.312 76.502 1.00149.93 A C ANISOU 8 CA ILE A 11 19407 19827 17732 -699 -2577 -541 A C ATOM 9 CB ILE A 11 24.451 7.060 75.423 1.00142.17 A C ANISOU 9 CB ILE A 11 18258 18757 17003 -502 -2834 -817 A C ATOM 10 CG2 ILE A 11 24.819 6.135 74.260 1.00140.21 A C ANISOU 10 CG2 ILE A 11 17760 18301 17213 -445 -3136 -730 A C ATOM 11 CG1 ILE A 11 23.648 8.266 74.929 1.00141.20 A C ANISOU 11 CG1 ILE A 11 17949 18698 17002 -369 -2502 -1103 A C ATOM 12 CD1 ILE A 11 22.184 7.945 74.630 1.00147.23 A C ANISOU 12 CD1 ILE A 11 18421 19495 18025 -392 -2156 -1014 A C ATOM 13 C ILE A 11 24.000 4.831 76.450 1.00152.71 A C ANISOU 13 C ILE A 11 19772 20013 18238 -798 -2860 -201 A C ATOM 14 O ILE A 11 24.895 4.377 77.158 1.00153.29 A O ANISOU 14 O ILE A 11 20114 20074 18055 -857 -3135 -9 A O I would like to know if the refinement was truly grouped B-factor refinement - I am suspecting the residual B-factors might be grouped, but I haven't checked for that. Could anyone, maybe Pavel, confirm that is the case? Is there anyplace in the log files I could see the number of B groups? Thanks, Engin

2 1

Re: [phenixbb] Completeness issue by phenix
by Young-Jin Cho 17 Dec '08

17 Dec '08

Thanks Ralf, You were right. I found that my .log file has Anomalous flag:True. As it was pointed out, my mtz file was came from ccp4 where it uses anomalous data by default. Another small wondering is about differently automatic(?) resolution pick up: my refmac assigned resolution range between 1.48 and 54.23 whereas my phenix.refine set it between 1.363 and 31.041 by default. Also while the refmac provides both entire and high resolution bin completeness, the phenix only gives a whole completeness. Is there any way I can investigate high resolution bin completeness? and why I have different resolution ranges by default without specifically assigning them? Thanks again. YoungJin ----- Original Message ----- From: "Ralf W. Grosse-Kunstleve" <rwgk(a)cci.lbl.gov> To: phenixbb(a)phenix-online.org Sent: Monday, December 15, 2008 12:59:12 PM GMT -05:00 US/Canada Eastern Subject: Re: [phenixbb] Completeness issue by phenix Hi YoungJin, I suspect phenix.refine decided your mtz file contains anomalous data when you actually have non-anomalous data. Could you look in the phenix.refine output for something like this: ================================== X-ray data ================================= F-obs: 1yjp.mtz:FOBS_X,SIGFOBS_X R-free flags: 1yjp.mtz:R-free-flags Miller array info: 1yjp.mtz:FOBS_X,SIGFOBS_X Observation type: xray.amplitude Type of data: double, size=495 Type of sigmas: double, size=495 Number of Miller indices: 495 Anomalous flag: False ^^^^^ in particular here If you see "Anomalous flag: True" try adding this to the phenix.refine command line: xray_data.force_anomalous_flag_to_be_equal_to=False Let me know if this doesn't help. If it does help, I'd be interested to know how the .mtz file was generated. Ralf _______________________________________________ phenixbb mailing list phenixbb(a)phenix-online.org http://www.phenix-online.org/mailman/listinfo/phenixbb

2 1

Re: [phenixbb] Completeness issue by phenix
by Ralf W. Grosse-Kunstleve 16 Dec '08

16 Dec '08

Hi Rupert, Thank you very much for the files! I found a simple, universal trick to "do the right thing", not just in phenix.refine. This will be in future releases. Ralf

1 0

Completeness issue by phenix
by Young-Jin Cho 16 Dec '08

16 Dec '08

Hi, I was happy for using phenix but am struggling with a strange result on completeness calculation. One typical data set I used kept showing low completeness, so I used the same mtz, and pdb file on Refmac: it turned out has more than 98% whereas phenix.refine gave me about 46% completeness. I have slightly lower value of R but higher Rfree from phenix calculation without adding H atoms. I don't know if this is a bug or a problem on my data set (hopefully not). thanks in advance, YoungJin

4 4

phenix GUI takes longer time to comeup
by shivesh kumar 15 Dec '08

15 Dec '08

Dear all The GUI of phenix installed in Fedora 8 takes longer time, sometimes 30mins. Can anyone help in this regard, what is the problem and how it can be sorted out? Thanx in advance Shivesh

1 0

Re: [phenixbb] individual anisotropic adps
by Pavel Afonine 14 Dec '08

14 Dec '08

Currently phenix.refine does not support refinement of individual anisotropic ADP and TLS for the same set of atoms (see my previous email). However, in phenix.refine you can refine TLS and individual anisotropic ADP for non-overlapping set of atoms. Pavel. On 12/11/2008 8:34 AM, gerwald jogl wrote: > Thanks, that's roughly what I thought. How about tls and individual > anisotropic together? same strategy, just observe Rfree?

3 2

individual anisotropic adps
by gerwald jogl 12 Dec '08

12 Dec '08

Hi All, I was wondering if there is information/recommendations out there for the resolution required to refine individual anisotropic adps. I think the recommendation for shelx is at least 1.5 A resolution. In addition, is phenix set up and would it make sense to combine tls with individiual anisotropic adps? All comments welcome, Gerwald

4 3

Keeping the same origin
by Mo Wong 12 Dec '08

12 Dec '08

Hi all, I would like to compare maps from different phasing runs, so it would make my life much easier if I could keep the same origin. I recall this could be done in solve/resolve, but can't seem to find it mentioned in the phenix documentation. Any suggestions? Thanks!

2 1

resolution bins
by Daniele de Sanctis 11 Dec '08

11 Dec '08

Hi, I was wondering how phenix.refine decide the number of resolution bins in which it shows the refinement statistics (R-free, phase error...). In the case I'm working on it prints out only three resolution bins. Thanks d -- Daniele de Sanctis, PhD Macromolecular Crystallography Group ESRF, Grenoble, France Tel 33 (0)4 76 88 2869

2 1

Question regarding phenix.elbow
by Young-Jin Cho 04 Dec '08

04 Dec '08

Hi I would like to know how I can add H atoms to my ligand, ADP in the protein complex. I was trying to add H atoms to one of my ligand, ADP. When I types "phenix.elbow --final-geometry=test10bcoot_h_refine_001.pdb --residue=ADP --output=model_h", it went on but failed with the following error message: Traceback (most recent call last): File "/home2/yjcho/phenix-1.3-final/elbow/elbow/command_line/builder.py", line 1201, in <module> run() File "/home2/yjcho/phenix-1.3-final/elbow/elbow/command_line/builder.py", line 1108, in run molecule = post_process(molecule, options) File "/home2/yjcho/phenix-1.3-final/elbow/elbow/command_line/builder.py", line 577, in post_process options, File "/home2/yjcho/phenix-1.3-final/elbow/elbow/utilities/optimisation_manager.py", line 331, in optimisation_manager pre_opt)) File "/home2/yjcho/phenix-1.3-final/elbow/elbow/chemistry/IOMixins.py", line 857, in WritePDB original_order=original_order, File "/home2/yjcho/phenix-1.3-final/elbow/elbow/chemistry/IOMixins.py", line 965, in WritePDB2String self.AdjustHydrogenNames() File "/home2/yjcho/phenix-1.3-final/elbow/elbow/chemistry/IOMixins.py", line 940, in AdjustHydrogenNames if test_name[3]==" ": test_name[3]="0" TypeError: 'str' object does not support item assignment FYI. ADP is names as ADP as a residue name in the pdb file but has different chain names. By the way, I should mention that I used phenix.reduce first for adding H atoms but it was not able to add sugar H's in ADP molecule albeit other residues were fine. Thanks in advance for your help. Best, YoungJin

3 3

Bulk solvent
by csantiag＠cnb.csic.es 04 Dec '08

04 Dec '08

Hi all, I'm refining a protein and after running phenix.refine I always get a bulk solvent refinement cycle at the end of the process. Is there any way to get rid of the last step? ------------------------------------------------------------------------ REMARK stage angl bond chir dihe plan repu geom_target REMARK 0 : 1.007 0.027 0.237 18.975 0.004 4.107 2.7953e-01 REMARK 1_bss: 1.007 0.027 0.237 18.975 0.004 4.107 2.7953e-01 REMARK 1_xyz: 1.844 0.015 0.087 19.502 0.005 4.106 1.9824e-01 REMARK 1_adp: 1.844 0.015 0.087 19.502 0.005 4.106 1.9824e-01 REMARK 2_bss: 1.844 0.015 0.087 19.502 0.005 4.106 1.9824e-01 REMARK 2_xyz: 1.222 0.008 0.075 19.562 0.005 4.107 1.2876e-01 REMARK 2_adp: 1.222 0.008 0.075 19.562 0.005 4.107 1.2876e-01 REMARK 3_bss: 1.222 0.008 0.075 19.562 0.005 4.107 1.2876e-01 REMARK 3_xyz: 1.021 0.005 0.063 19.049 0.004 4.111 1.0306e-01 REMARK 3_adp: 1.021 0.005 0.063 19.049 0.004 4.111 1.0306e-01 REMARK 3_bss: 1.021 0.005 0.063 19.049 0.004 4.111 1.0306e-01 REMARK ------------------------------------------------------------------------ César Santiago Dept of Macromolecules Structure Centro Nacional de Biotecnología, CSIC Campus Cantoblanco 28049 Madrid Spain csantiag(a)cnb.csic.es

3 3

anomalous scatters refinement
by Daniele de Sanctis 03 Dec '08

03 Dec '08

Hi, I'm making some refinement tests with anomalous data and I wanted to try to refine the anomalous contribution. When would it be the best stage when to introduce it? Thanks D -- Daniele de Sanctis, PhD Macromolecular Crystallography Group ESRF, Grenoble, France Tel 33 (0)4 76 88 2869

2 1

phenix.autobuild- multiple models
by meindert lamers 30 Nov '08

30 Nov '08

Hi all, I am a bit puzzled by an error message I get when trying to run phenix.autobuild. (I am using the latest/final version of phenix: PHENIX VERSION: 1.3 of 11-08-2008) phenix.autobuild data=01.mtz model=3mei_coot.pdb multiple_models=True include_input_model=True \ multiple_models_number=1 n_cycle_rebuild_max=5 This gives me: **************************************** AutoBuild Input failed Sorry, multiple_models is only available with rebuild_in_place **************************************** However, when I run the same commands using an older version of phenix (PHENIX VERSION: 1.3 of 23-03-2008), all runs fine. Any ideas ? Meindert

2 1

R/Rfree from phenix.autobuild
by Graeme Winter 29 Nov '08

29 Nov '08

Hi Folks, Quick question about R/Rfree from phenix - I ran the autobuild and it kindly built pretty much all of the visible residues, and it gave me the following output: SOLUTION CYCLE R RFREE BUILT PLACED 1 1 0.22 0.23 169 169 2 2 0.21 0.23 170 170 3 4 0.22 0.23 170 170 Now, I'm not one to complain about good results, but I did feel that the Rfree looked a little close to the R factor - usually I would expect a difference of maybe 4-5%. Is there a reason for this that I am missing? The model seemed to fit very nicely, closely corresponded to the PDB entry and the waters were pretty good. There is no NCS (one copy in the ASU) and the phases are from SeMet mad data from Autosol. Resolution is ~1.6A, if that is any help. Just noticed this and thought it was interesting. I am used to Arp/wArp for automated tracing, which typically gives a bigger difference... Many thanks, Graeme

3 2

0 for wxc_scale and wxu_scale?
by Michaela Kramer, Biochemisches Inst. 27 Nov '08

27 Nov '08

I am refining a structure with a resolution of 2.1 Angstroem using Phenix. The spacegroup is P 31 2 1 (No. 152) with a twin fraction of 0.37 and twin law -h,-k,l. In the refinement procedure the NCS (4 chains) constraints are on, using individual sites, ADP and occupancy refinement. In the last refinement step, I use the option to optimize the values for wxu_scale and wxc_scale. These values are optimizes to 0 for wxc_scale and close to 0 for wxu_scale (0.2) When I tried this at earlier refinement stages wxu_scale was also here optimized in some refinement macrocycles to 0. Interestingly I obtain better R values (0.1948/0.2184), than without the optimization by using the default values (0.1391/0.2438 Rwork/Rfree). In addition to that the statistics are better: clashscore 14.41 (compared to 25.73 without wxc_scale & wxu_scale optimization), rotamer outliers 2.33% (compared o 4.66%) and Ramachandran favoured (Molprobity). How is it possible that wxu_scale and wxc_scale can be optimized to 0 by having better statistics? The total refinement target is defined as Etotal= wxc_scale x wxc x Exray + wc x Egeom As I understand it, the crystallographic refinement target (Exray) would now not be considered in the refinement procedure. Would this refinement procedure be valid? And how is wxc (output of target weights: x-ray data) defined? Thank you very much for your help Michaela Here is the output belonging to the optimization: Refinement Macrocycle 1 of 3 ========================== Target weights: x-ray data ========================= |-----------------------------------------------------------------------------| | XYZ refinement: T = Eexperimental * wxc * wxc_scale + Echem * wc | | wxc = 1.649809 wxc_scale = 0.500 wc = 1.000 | | angle between x-ray and geometry gradient vectors: 93.128 (deg) | | | | ADP refinement: T = Eexperimental * wxu * wxu_scale + Eadp * wu | | wxc = 0.888942 wxc_scale = 1.000 wc = 1.000 | | angle between Xray and ADP gradient vectors: 85.407 (deg) | |-----------------------------------------------------------------------------| ================================ xyz refinement =============================== Start r_free = 0.2118 scale= 0.0000 total_weight= 0.0000 r_free= 0.2181 0 scale= 1.0000 total_weight= 1.6498 r_free= 0.2336 0 scale= 1.5000 total_weight= 2.4747 r_free= 0.2379 1 scale= 2.0000 total_weight= 3.2996 r_free= 0.2412 2 scale= 0.6667 total_weight= 1.0999 r_free= 0.2297 0 scale= 0.5000 total_weight= 0.8249 r_free= 0.2245 1 scale= 0.4000 total_weight= 0.6599 r_free= 0.2226 2 Best r_free = 0.2181 ================================ ADP refinement =============================== ----------Individual ADP refinement---------- Start r_free = 0.2181 scale= 0.0000 total_weight= 0.0000 r_free= 0.2247 0 scale= 1.0000 total_weight= 0.8889 r_free= 0.2262 0 scale= 1.5000 total_weight= 1.3334 r_free= 0.2293 1 scale= 2.0000 total_weight= 1.7779 r_free= 0.2310 2 scale= 0.6667 total_weight= 0.5926 r_free= 0.2236 0 scale= 0.5000 total_weight= 0.4445 r_free= 0.2221 0 scale= 0.4000 total_weight= 0.3556 r_free= 0.2208 0 scale= 0.3333 total_weight= 0.2963 r_free= 0.2203 0 scale= 0.2857 total_weight= 0.2540 r_free= 0.2195 0 scale= 0.2500 total_weight= 0.2222 r_free= 0.2191 0 scale= 0.2222 total_weight= 0.1975 r_free= 0.2186 0 scale= 0.2000 total_weight= 0.1778 r_free= 0.2183 0 Best r_free = 0.2183 Refinement Macrocycle 2 of 3 ========================== Target weights: x-ray data ========================= |-----------------------------------------------------------------------------| | XYZ refinement: T = Eexperimental * wxc * wxc_scale + Echem * wc | | wxc = 1.737535 wxc_scale = 0.500 wc = 1.000 | | angle between x-ray and geometry gradient vectors: 89.612 (deg) | | | | ADP refinement: T = Eexperimental * wxu * wxu_scale + Eadp * wu | | wxc = 0.228155 wxc_scale = 1.000 wc = 1.000 | | angle between Xray and ADP gradient vectors: 155.392 (deg) | |-----------------------------------------------------------------------------| ================================ xyz refinement =============================== Start r_free = 0.2179 scale= 0.0000 total_weight= 0.0000 r_free= 0.2223 0 scale= 1.0000 total_weight= 1.7375 r_free= 0.2272 0 scale= 1.5000 total_weight= 2.6063 r_free= 0.2310 1 scale= 2.0000 total_weight= 3.4751 r_free= 0.2332 2 scale= 0.6667 total_weight= 1.1584 r_free= 0.2267 0 scale= 0.5000 total_weight= 0.8688 r_free= 0.2258 1 scale= 0.4000 total_weight= 0.6950 r_free= 0.2269 2 Best r_free = 0.2223 ================================ ADP refinement =============================== ----------Individual ADP refinement---------- Start r_free = 0.2223 scale= 0.0000 total_weight= 0.0000 r_free= 0.2228 0 scale= 1.0000 total_weight= 0.2282 r_free= 0.2235 0 scale= 1.5000 total_weight= 0.3422 r_free= 0.2256 1 scale= 2.0000 total_weight= 0.4563 r_free= 0.2275 2 scale= 0.6667 total_weight= 0.1521 r_free= 0.2219 0 scale= 0.5000 total_weight= 0.1141 r_free= 0.2213 0 scale= 0.4000 total_weight= 0.0913 r_free= 0.2207 0 scale= 0.3333 total_weight= 0.0761 r_free= 0.2204 0 scale= 0.2857 total_weight= 0.0652 r_free= 0.2201 0 scale= 0.2500 total_weight= 0.0570 r_free= 0.2197 0 scale= 0.2222 total_weight= 0.0507 r_free= 0.2195 0 scale= 0.2000 total_weight= 0.0456 r_free= 0.2194 0 Best r_free = 0.2194 Macrocycle 2 of 3 ========================== Target weights: x-ray data ========================= |-----------------------------------------------------------------------------| | XYZ refinement: T = Eexperimental * wxc * wxc_scale + Echem * wc | | wxc = 2.157676 wxc_scale = 0.500 wc = 1.000 | | angle between x-ray and geometry gradient vectors: 93.409 (deg) | | | | ADP refinement: T = Eexperimental * wxu * wxu_scale + Eadp * wu | | wxc = 0.228399 wxc_scale = 1.000 wc = 1.000 | | angle between Xray and ADP gradient vectors: 140.390 (deg) | |-----------------------------------------------------------------------------| ================================ xyz refinement =============================== Start r_free = 0.2163 scale= 0.0000 total_weight= 0.0000 r_free= 0.2189 0 scale= 1.0000 total_weight= 2.1577 r_free= 0.2299 0 scale= 1.5000 total_weight= 3.2365 r_free= 0.2340 1 scale= 2.0000 total_weight= 4.3154 r_free= 0.2341 2 scale= 0.6667 total_weight= 1.4385 r_free= 0.2282 0 scale= 0.5000 total_weight= 1.0788 r_free= 0.2255 1 scale= 0.4000 total_weight= 0.8631 r_free= 0.2249 2 Best r_free = 0.2189 ================================ ADP refinement =============================== ----------Individual ADP refinement---------- Start r_free = 0.2189 scale= 0.0000 total_weight= 0.0000 r_free= 0.2218 0 scale= 1.0000 total_weight= 0.2284 r_free= 0.2231 0 scale= 1.5000 total_weight= 0.3426 r_free= 0.2250 1 scale= 2.0000 total_weight= 0.4568 r_free= 0.2272 2 scale= 0.6667 total_weight= 0.1523 r_free= 0.2215 0 scale= 0.5000 total_weight= 0.1142 r_free= 0.2209 0 scale= 0.4000 total_weight= 0.0914 r_free= 0.2201 0 scale= 0.3333 total_weight= 0.0761 r_free= 0.2197 0 scale= 0.2857 total_weight= 0.0653 r_free= 0.2195 0 scale= 0.2500 total_weight= 0.0571 r_free= 0.2193 0 scale= 0.2222 total_weight= 0.0508 r_free= 0.2191 0 scale= 0.2000 total_weight= 0.0457 r_free= 0.2190 0 Best r_free = 0.2190 Michaela Kramer PhD student Institute of Biochemistry Winterthurerstrasse 190 8057 Zuerich Switzerland

7 6

trouble in phenix.refine
by Huiying Li 27 Nov '08

27 Nov '08

I was running phenix.refine in tls mode with the following command line: phenix.refine mydata.mtz mymodel.pdb strategy=tls tls_group_selection.param mylibrary.cif against a 2.1A data set, P21211 SG, a homodimer in AU, totally 10 tls groups. The following errors occured and program stopped. I do not know how to fix the problem. Traceback (most recent call last): File "/usr/local/phenix-1.3-final/phenix/phenix/command_line/refine.py", line 5, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/usr/local/phenix-1.3-final/phenix/phenix/refinement/command_line.py", line 68, in run log=log) .... ... .... File "/usr/local/phenix-1.3-final/mmtbx/mmtbx/monomer_library/pdb_interpretation.py", line 970, in get_lib_link atom_dicts[angle.atom_3_comp_id].get(angle.atom_id_3, None)] IndexError: list index out of range Your help is highly appreciated. Huiying -- Huiying Li, Ph. D Department of Molecular Biology and Biochemistry Natural Sciences I, Rm 2443 University of California at Irvine Irvine, CA 92697, USA Tel: 949-824-4322(or -1953); Fax: 949-824-3280 email: hli(a)uci.edu

2 2

twin axis for p21
by HUA YUAN 26 Nov '08

26 Nov '08

Dear Colleagues, I'd like to share with you with a twinning puzzle that's been troubling me for a while. I have a ~2.2 A data merged in p212121 space group with unit cell: 41.8300 117.4490 134.4840 90.0000 90.0000 90.0000 Translational NCS was detected but structure solution was able to be obtained by phased molecular replacement as we have a Au derivative data set. Obtaining experimental phase or MR alone didn't work. However, with this model and data set, the refinement didn't work (with R-free ~40% or higher if with more cycles). Different twinning tests have been carried out and no sign of twinning. All possible p222 space groups have been tried but no luck or even worse. Considering that translational symmetry and perfect twinning may make things deceiving, the data was merged into p21 and a twin-refinement was performed in phenix.refine. After 3 cycles of rigid-body and group ADP, a R-free of ~30% was obtained and the output map looks OK. It looks like we got it. But then things get interesting when we tested other two axis as the p21 unique axis. With all the three choices plus corresponding twin-refinement, we could get a similar R-free (30%). Besides, the three maps all look OK except for one missing some side chains. Then the structure was solve in p1 and Xtriage was used with the hope that one axis would stand out. However, the RvsR for all three axis are almost the same. Did I miss anything here? Of course, we could pick up the one axis we like best and get the structure solved ^---^. However, this puzzle would eat me alive for a very long time. Best, Hua

2 1

solvent flipping
by Mark Collins 25 Nov '08

25 Nov '08

Hi Is there a solvent flipping (like solomon, ccp4) rather than flattening option in phenix? Mark

4 3