- phenixbb - phenix-online.org

can phenix.refine output mtz with structure factors not just coeficients?
by amadrona＠uci.edu 02 Sep '10

02 Sep '10

Hi, I am going to deposit a structurein the pdb but noticed that all of the mtz files from phenix.refine only contain coeficients and not reflection info. Is it as simple as running another refinement with write_maps = true instead of false? Will this right out an mtz with structure factors instead of just coeficients. I am a newbie at crystalagraphy so I appolgize if this obvious. Ideally, I would like to create an mtz with structure factors as well as a map for CNS. Is there a way to do this without runing another refinment maybe using "maps.params". Thank you very much for the help! -Yarrow output { prefix = "CIN4_refine" serial = 12 serial_format = "%03d" title = None write_eff_file = True write_geo_file = True write_final_geo_file = False write_def_file = True export_final_f_model = mtz cns write_maps = false write_map_coefficients = True pickle_fmodel = False

3 3

atom selection in phenix.pdbtools
by Michael Hothorn 01 Sep '10

01 Sep '10

Hi, I want to calculate the mean b-value for all atoms of a specific ligand in my structure. I tried: phenix.pdbtools --show-adp-statistic model.pdb keep="chain A and residue 400" , but it still prints out all b-values including solvent etc. any ideas? thanks Michael

3 3

Re: [phenixbb] question on bad geometry
by Thomas C. Terwilliger 01 Sep '10

01 Sep '10

Hi Fengyun, Just to emphasize what Ed points out here, an essential step in completing a structure is to compare it, one residue, solvent molecule, or ligand at a time, with your electron density map, fixing any inconsistencies between the model and the map. There are many wonderful tools in all the crystallographic graphics programs to help you do this. All the best, Tom T >> On Tue, 2010-08-31 at 20:37 -0500, fn1(a)rice.edu wrote: >>> Does >>> that mean I need to manually adjust the model? >> >> Autobuilder is a fantastic tool which greatly simplifies the model >> building process, but you always must validate your structure by >> inspecting the electron density. I fully understand that (hopefully) it >> is what you plan to do and this is more of a reaction to the alarming >> number of obvious and easily correctable systematic errors one finds in >> structures deposited in the PDB these days. >> >> Thanks to Pavel for providing an excellent example to illustrate the >> point. One only needs to take a look at the model in context of >> electron density maps to confirm that autobuilder will fix some errors >> but also sometimes introduce some that require manual correction. >> >> Cheers, >> >> Ed. >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://phenix-online.org/mailman/listinfo/phenixbb >>

4 3

question on bad geometry
by fn1＠rice.edu 01 Sep '10

01 Sep '10

Hi everyone, I have a structure at 2.3 A refined in phenix to R/Rfree=0.23/0.28. However, the geometry is not that good, the followings are the output from phenix.model_vs_data, Stereochemistry statistics (mean, max, count): bonds : 0.0080 0.0547 3244 angles : 1.1669 12.1085 4374 dihedrals : 15.7241 85.6280 1202 chirality : 0.0706 0.2831 536 planarity : 0.0036 0.0246 573 non-bonded (min) : 2.1217 Ramachandran plot, number of: outliers : 17 (4.02 %) allowed : 37 (8.75 %) favored : 369 (87.23 %) Rotamer outliers : 30 (8.62 %) goal: < 1% Cbeta deviations >0.25A : 0 All-atom clashscore : 40.51 (steric overlaps >0.4A per 1000 atoms) I tried to use fix_rotamer=True, but it doesn't help at all. As to the rotamer outliers, the value is far higher than the ideal one. Does that mean I need to manually adjust the model? Or will autobuild help improve the geometry? Thanks! Fengyun

6 5

How to do a phenix.fmodel calculation incorporating the exponential scaling model only?
by zhangh1＠umbc.edu 31 Aug '10

31 Aug '10

Hi, I want to use phenix.fmodel to do a calculation incorporating the exponential scaling model only (Not the flat model). I am not sure whether I just need simply specify: mask { use_asu_masks = False } since the output maps are just identical no matter use_asu_masks = True or False. Should I do something else to enable exponential scaling model only? Thanks! Best Regards, Hailiang

2 2

0.00 Refined Isotropic ADPs
by Joseph Noel 31 Aug '10

31 Aug '10

Hi All. Me again! LOL I am sure this has come up and related to how the "B-factor" in the refined PDB's ADP column is calculated including the overall anisotropic B-factor, etc but is it unusual to have isotropic ADPs refine to 0.00? I have several nice structures with 0.00 isotropic ADPs. I am not refining TLS or anisotropic ADPs at all. Just the normal course of refinement. These are regions of a high degree of order and the maps are beautiful but just wandering what the 0.00 ultimately means particularly when I eventually submit to the Protein Data Bank. Joe P.S. Nat's tips on COOT and Phenix are working beautifully now! Thanks Nat! ___________________________________________________________ Joseph P. Noel, Ph.D. Investigator, Howard Hughes Medical Institute Professor, The Jack H. Skirball Center for Chemical Biology and Proteomics The Salk Institute for Biological Studies 10010 North Torrey Pines Road La Jolla, CA 92037 USA Phone: (858) 453-4100 extension 1442 Cell: (858) 349-4700 Fax: (858) 597-0855 E-mail: noel(a)salk.edu Web Site (Salk): http://www.salk.edu/faculty/faculty_details.php?id=37 Web Site (HHMI): http://hhmi.org/research/investigators/noel.html ___________________________________________________________

2 1

Re: [phenixbb] Phenix and COOT Together on a Mac (and Probe too?)
by Joseph Noel 31 Aug '10

31 Aug '10

Nat, One thing that seemed to solve all the problems at once was to completely remove any reference to the COOT path in the Phenix preferences! Is that weird? When I removed it completely only leaving the path to PYMOL then COOT is working just fine in conjunction with Phenix and even the Probe issues have disappeared. Prior to this, I was tried two different paths to COOT. 1. /sw64/bin/coot 2. /Applications/coot.app/Contents/MacOS/coot Neither behaved but if I don't specify anything it works. Joe ___________________________________________________________ Joseph P. Noel, Ph.D. Investigator, Howard Hughes Medical Institute Professor, The Jack H. Skirball Center for Chemical Biology and Proteomics The Salk Institute for Biological Studies 10010 North Torrey Pines Road La Jolla, CA 92037 USA Phone: (858) 453-4100 extension 1442 Cell: (858) 349-4700 Fax: (858) 597-0855 E-mail: noel(a)salk.edu Web Site (Salk): http://www.salk.edu/faculty/faculty_details.php?id=37 Web Site (HHMI): http://hhmi.org/research/investigators/noel.html ___________________________________________________________

2 1

Problem with low figure of merit
by sankaranarayanan srinivasan 31 Aug '10

31 Aug '10

Hi all, I am new to structure solving. I have a 2.9A dataset for a Seleno crystal of which i am trying to use Autosol to build an initial map. I used the settings for a thorough solution, instead of quick. All the other parameters were set to default settings. The solution that i got from autosol, gave me the following result Solution # 1 BAYES-CC: 23.7 +/- 30.9 Dataset #1 FOM: 0.27 Residues built: 238 Side-chains built: 139 Chains: 15 Overall model-map correlation: 0.65 R/R-free: 0.38/0.42. Score type: SKEW CORR_RMS Raw scores: 0.02 0.77 100x EST OF CC: 16.86 32.11 The density modification map matches, quite well with many of my side chains fitted. My question here is , IS THERE A WAY I CAN IMPROVE THE Figure Of Merit , BY CHANGING ANY PARAMETERS IN PHENIX-AUTOSOL. I HAVE A HIGH SOLVENT CONTENT IN MY CRYSTAL AND I DID NOT INPUT THAT IN MY INITIAL PARAMETERS. DOES PHENIX AUTOSOL CONSIDER SOLVENT FLATTENING FOR SUCH CRYSTALS? I would appreciate someone's help in this regard Thanks Shankar

2 1

Phenix and COOT Together on a Mac (and Probe too?)
by Joseph Noel 30 Aug '10

30 Aug '10

Hi All, I am having trouble getting COOT to link up with Phenix. I know that Bill Scott's fink version of COOT breaks the link but I have also tried the stand alone latest version for Mac OSX without success. I suspect there may be something additional I need to set up to get the two programs to talk to one and other. I can open COOT with the GUI link in Phenix but that's it. The two do not communicate and the "Phenix button" does not appear. I am running everything on a Mac laptop with an Intel chip and the latest snow leopard operating system. I can also run probe independently in COOT using the script suggestion for the .coot file from the MOLPROBITY web site. However, after several checks of bad contacts the system freezes at the level of XQuartz and I have to kill X-Windows to start over. I don't know if this is in someway linked to the problem getting COOT and Phenix to talk to one and other but I thought I would throw it out there. Finally, if I do run COOT by starting it from the Phenix GUI it also freezes after several rounds of building, etc just like the probes to check for clashes. If anyone has observed this behavior and has some suggestions I would much appreciate it! Its not a big deal as I can run things separately and work from a print out of Molprobity outliers but it would be sure nice to just click and correct. Sorry to keep going but I should say I have Bill Scott's fink'ed version of coot installed in as directed in /sw64 and the stand alone version installed along with CCP4 and Phenix in my Mac's /Applications directory where all my standard Mac apps appear. Thanks! Joe ___________________________________________________________ Joseph P. Noel, Ph.D. Investigator, Howard Hughes Medical Institute Professor, The Jack H. Skirball Center for Chemical Biology and Proteomics The Salk Institute for Biological Studies 10010 North Torrey Pines Road La Jolla, CA 92037 USA Phone: (858) 453-4100 extension 1442 Cell: (858) 349-4700 Fax: (858) 597-0855 E-mail: noel(a)salk.edu Web Site (Salk): http://www.salk.edu/faculty/faculty_details.php?id=37 Web Site (HHMI): http://hhmi.org/research/investigators/noel.html ___________________________________________________________

2 1

Can we untoggle molecular weight mean deviation check on ensembles in autoMR?
by James R. Thompson 30 Aug '10

30 Aug '10

Can folks untoggle the molecular weight mean deviation check in phenix.autoMR on the ensembles input? I want to use an ensemble of coordinates generated by the Mustangs server from a SSM alignment. Therefore, by design the 10 PDBs input to phenix.autoMR vary greatly in molecular weight from their mean value, but the whole is likely an excellent probe for molecular replacement. However, I am stymied by phenix.autoMR: ------------------------------------------------------------------------------------- FATAL RUNTIME ERROR: Molecular weight of /home/jrt04/Desktop/GB1STAS/MUSTANGServer_10_1.6_0_sculpt-coot-0.pdb (6610.56) deviates more than 20% from the mean (3977.08). This pdb file may contain domains or water not present in the other files. ------------------------------------------------------------------------------------- I understand why this check might be of utility to the unobservant. How can I disable this test? Many thanks, Jim -- James R. Thompson <thompson.james(a)mayo.edu>

1 0

stuck with phenix refine
by intekhab alam 29 Aug '10

29 Aug '10

Hi I am working with a twinned data set having R32 space group. when i tried to refine with phenix refine using the twin option it got stuck after some time. There is no error message showing up and the refinemnet donot progresses further. I have used the twin law suggested by Xtriage program in Phenix. Kindly help me out. Regards -- INTEKHAB ALAM LABORATORY OF STRUCTURAL BIOINFORMATICS KOREA UNIVERSITY, SEOUL

4 5

Re: [phenixbb] phenix.automr questions
by Kay Diederichs 29 Aug '10

29 Aug '10

Dear Randy, Nat, Tom, thank you for pointing out the existence of phenix.phaser to me, and sketching its usage. I think I'll make my way with that information; nevertheless it would be great if it were officially documented (Google returned nothing to me about phenix.phaser). BTW, I'm trying to put together some orthogonal Phenix documentation at http://strucbio.biologie.uni-konstanz.de/ccp4wiki/index.php/Phenix . I take the opportunity to state that although the URL refers to the CCP4 users' wiki, and obviously Phenix is not part of CCP4, the wiki covers a broad range of crystallographic information. Thus, I do think it is appropriate to collect Phenix-related information (from a user perspective) there, and would like to invite everybody to contribute to the Phenix entry (and others!). thanks, Kay

2 1

Re: [phenixbb] phenix.automr questions
by Kay Diederichs 29 Aug '10

29 Aug '10

Hi Tom, thanks for answering my questions! Yes, it would be good to have a flag to turn off the inverse hand because a) having the possibility to look at each spacegroup offers a very effective way of parallelization b) if I "grep" the spacegroup name in the logfile I find: dikay@turn10:-2010/phenix.automr% grep "P 4" mf92.log Space group: P 41 21 2 (No. 92) LIST OF ALL PLAUSIBLE SG: ['P 4 2 2', 'P 4 21 2', 'P 41 2 2', 'P 41 21 2', 'P 42 2 2', 'P 42 21 2', 'P 43 2 2', 'P 43 21 2'] Space group: P 41 21 2 (No. 92) NEW dictDATA {'raw_data': [' P 4abw 2nw', (101.71700286865234, 101.71700286865234, 344.74798583984375, 90.0, 90.0, 90.0), <scitbx_array_family_flex_ext.double object at 0xe554a10>, <scitbx_array_family_flex_ext.double object at 0xe554bd0>, <cctbx_array_family_flex_ext.miller_index object at 0xe4716c0>]} space_group : P 41 21 2 all_plausible_sg_list : ['P 4 2 2', 'P 4 21 2', 'P 41 2 2', 'P 41 21 2', 'P 42 2 2', 'P 42 21 2', 'P 43 2 2', 'P 43 21 2'] dictDATA: {'raw_data': [' P 4abw 2nw', (101.71700286865234, 101.71700286865234, 344.74798583984375, 90.0, 90.0, 90.0), <scitbx_array_family_flex_ext.double object at 0xe554a10>, <scitbx_array_family_flex_ext.double object at 0xe554bd0>, <cctbx_array_family_flex_ext.miller_index object at 0xe4716c0>]} HALL: P 4abw 2nw P 43 21 2 P 41 21 2 Space-Group Name (Number): P 41 21 2 (92) Space-Group Name: P 41 21 2 Space-Group Name (Number): P 41 21 2 (92) P 43 21 2 * P 41 21 2 Space-Group Name (Number): P 43 21 2 (96) Fast Translation Function Table: Space Group P 43 21 2 Space Group Name (Number): P 43 21 2 (96) Packing Table: Space Group P 43 21 2 Space-Group Name (Number): P 43 21 2 (96) Refinement Table (Unsorted): Space Group P 43 21 2 Refinement Table (Sorted): Space Group P 43 21 2 Best Space Group: P 43 21 2 P 43 21 2 P 41 21 2 * Space-Group Name (Number): P 41 21 2 (92) Fast Translation Function Table: Space Group P 41 21 2 Best Space Group so far: P 43 21 2 Space Group Name (Number): P 41 21 2 (92) Packing Table: Space Group P 41 21 2 Space-Group Name (Number): P 41 21 2 (92) Refinement Table (Unsorted): Space Group P 41 21 2 Refinement Table (Sorted): Space Group P 41 21 2 Space-Group Name (Number): P 43 21 2 (96) Best spacegroup from search with first component in permutation: P 43 21 2 Space-Group Name (Number): P 43 21 2 (96) Fast Translation Function Table: Space Group P 43 21 2 Space Group Name (Number): P 43 21 2 (96) Packing Table: Space Group P 43 21 2 Space-Group Name (Number): P 43 21 2 (96) Refinement Table (Unsorted): Space Group P 43 21 2 Refinement Table (Sorted): Space Group P 43 21 2 Best Space Group: P 43 21 2 Space-Group Name (Number): P 43 21 2 (96) Best spacegroup from search with first component in permutation: P 43 21 2 Space-Group Name (Number): P 43 21 2 (96) Fast Translation Function Table: Space Group P 43 21 2 Space Group Name (Number): P 43 21 2 (96) Packing Table: Space Group P 43 21 2 Space-Group Name (Number): P 43 21 2 (96) Refinement Table (Unsorted): Space Group P 43 21 2 Refinement Table (Sorted): Space Group P 43 21 2 Best Space Group: P 43 21 2 Space-Group Name (Number): P 43 21 2 (96) which shows that at a certain point the MTZ file's spacegroup is abandoned and the calculation finishes in the enantiomorph spacegroup, which has the drawback that I cannot really compare the best (complete) results of both enantiomorphs. And with one good and three bad search models I don't know if abandoning the given spacegroup early is definitively the right decision (but it's hard for me to tell exactly what criterion the decision was based upon - it may have been after an exhaustive search in #92 in which case my worry would not be justified). thanks, Kay > Message: 2 > Date: Sat, 28 Aug 2010 09:07:08 -0600 (MDT) > From: "Thomas C. Terwilliger"<terwilliger(a)lanl.gov> > To: "PHENIX user mailing list"<phenixbb(a)phenix-online.org> > Subject: Re: [phenixbb] phenix.automr questions > Message-ID:<14563.128.165.0.81.1283008028.squirrel(a)webmail.lanl.gov> > Content-Type: text/plain;charset=iso-8859-1 > > Hi Kay, > > I'll answer what I can; perhaps Randy the others may have more comments. > >>> Dear developers, >>> >>> I have four questions concerning phenix.automr (I'm using version 1.6.4): >>> >>> a) is there a way to make the program try all permutations for the >>> search order of my four search molecules, or do I have to write a script ? > > It would be a good idea for me to put that in. For now, yes, you have to > make a script to do that. > >>> >>> b) how can I completely switch off the search in the alternative >>> spacegroups - I am under the impression that although the "effective >>> parameters" are >>> all_plausible_sg_list = None >>> use_all_plausible_sg = False >>> the search is not only done in the spacegroup I specified (#92=P41212) >>> but also in #96=P43212 . I believe so because the output says: >>> ... >>> 2 alternative spacegroups >>> SpaceGroup will be determined by highest LLG in search for FIRST >>> ensemble >>> Search order: >>> P 43 21 2 >>> P 41 21 2 >>> ... >>> (I would rather run a second search "manually" in the other spacegroup!) > > Yes, it runs the inverse hand even if you say > "use_all_plausible_sg=False". However it should not run all the other > space groups (all_plausible_sg_list : ['P 4 2 2', 'P 4 21 2', 'P 41 2 2', > 'P 41 21 2', 'P 42 2 2', 'P 42 21 2', 'P 43 2 2', 'P 43 21 2']). I'll see > if we can put in a flag so it doesn't run the inverse hand if you don't > want it to. > > >>> >>> c) I have two occurrences of seq_file ("all.seq" and "all.dat"; the >>> latter is symlinked to the former). Am I right that I could just leave >>> out one of them? I did this after looking at >>> http://www.phenix-online.org/documentation/automr.htm#anch70 but I >>> wonder which of the specifications of sequences and masses overrides the >>> other one? or is the first specification used for editing the search >>> models, and the second for getting the likelihood terms right? > > Looking at your parameters file, I see what you mean. If you supply two > different files, then yes, the one you specified at the top of your file > with "seq_file=all.dat" is going to be used for editing the file for > autobuild, and the one you specified in the scope "component" is going to > be used for getting the likelihood terms right. > > Normally (always?) these should be the same. So normally, just supply the > first one "input_seq_file=seq.dat") and it will be used for both purposes. > > Or if you have several components, only specify them in the "component" > scopes, and all of them will be put together into a single file and used > as the input to autobuild. > > >>> >>> d) what kind of editing is done? Is phenix.sculptor run internally? I >>> find {a,b,c,d}_edited.pdb but are they actually used? > > Only basic editing to make the model compatible with autobuild is done. No > sculpting is done inside automr for now. > > All the best, > Tom T

1 0

phenix.automr questions
by Kay Diederichs 29 Aug '10

29 Aug '10

Dear developers, I have four questions concerning phenix.automr (I'm using version 1.6.4): a) is there a way to make the program try all permutations for the search order of my four search molecules, or do I have to write a script ? b) how can I completely switch off the search in the alternative spacegroups - I am under the impression that although the "effective parameters" are all_plausible_sg_list = None use_all_plausible_sg = False the search is not only done in the spacegroup I specified (#92=P41212) but also in #96=P43212 . I believe so because the output says: ... 2 alternative spacegroups SpaceGroup will be determined by highest LLG in search for FIRST ensemble Search order: P 43 21 2 P 41 21 2 ... (I would rather run a second search "manually" in the other spacegroup!) c) I have two occurrences of seq_file ("all.seq" and "all.dat"; the latter is symlinked to the former). Am I right that I could just leave out one of them? I did this after looking at http://www.phenix-online.org/documentation/automr.htm#anch70 but I wonder which of the specifications of sequences and masses overrides the other one? or is the first specification used for editing the search models, and the second for getting the likelihood terms right? d) what kind of editing is done? Is phenix.sculptor run internally? I find {a,b,c,d}_edited.pdb but are they actually used? thank you, Kay P.S. I run with "phenix.automr mr.eff" where mr.eff is automr { data = "m879_2_rg92.mtz" seq_file = all.dat ensemble { ensembleID = "mol1" copies_to_find = 1 coords = a.pdb identity = 80 } ensemble { ensembleID = "mol2" copies_to_find = 1 coords = b.pdb identity = 45 } ensemble { ensembleID = "mol3" copies_to_find = 1 coords = c.pdb identity = 45 } ensemble { ensembleID = "mol4" copies_to_find = 1 coords = d.pdb identity = 45 } component { component_type = *protein nucleic_acid seq_file = all.seq component_copies = 1 } autobuild_variables{ n_cycle_rebuild_max = 1 } }

4 3

Phenix error running MIR-example
by Shankar Prasad Kanaujia 29 Aug '10

29 Aug '10

Dear Developers/Users, I am getting an error while running "rh-dehalogenase-mir" example provided in the phenix example. I am using phenix 1.6.4 version. It runs properly using the command "phenix.run_example rh-dehalogenase-mir", whereas, it gives error in GUI mode. I am attaching the picture of my input and error (in the right-top corner). Please let me know the reason of this error. -- Yours Sincerely, Shankar Prasad Kanaujia Research Student C/O - Dr. K. Sekar Bioinformatics Centre, SERC Indian Institute of Science, Bangalore-12 Phone - 9480258032 Office - 080-2293-3059

2 1

Postdoctoral position in Protein X-ray Crystallography
by cdas＠purdue.edu 29 Aug '10

29 Aug '10

Applications are invited for a postdoctoral research position in the Das Lab, Department of Chemistry, Purdue University, West Lafayette, IN, USA. The position requires the applicant to work on structural and biochemical characterization of ubiquitin hydrolases (deubiquitinases) and their complexes with substrates, cofactors and inhibitors (see Boudreaux et al, PNAS, 107, 9117 (2010)). The applicant is also expected to work on an ongoing collaborative project on method development related to imaging of protein crystals (of ubiquitin hydrolases and others) using second order non-linear optical methods. We are looking for highly motivated candidates with extensive experience in biochemical studies, crystallization and structure determination of proteins by x-ray crystallography. Specifically, the candidate is required to have hands-on experience with protein purification, crystallization, X-ray data collection, structure determination and structural analysis. Applicants should send their CV and the names and contact information of at least three references via email to: Chitta Das Assistant Professor Department of Chemistry Purdue University 560 Oval Drive West Lafayette IN, 47907 Email address: cdas(a)purdue.edu

1 0

Can Phenix refine B factors for several residues only?
by zhangh1＠umbc.edu 29 Aug '10

29 Aug '10

Hi, I want to refine B factors for several residues only (all the other B factors and all coordinates fixed, I know it sounds weird but there is a reason to try that). Is there anyway Phenix can do this? Thanks for any suggestions! Best Regards, Hailiang

4 11

Re: [phenixbb] advice for refinement with anomalous information
by John Rose 27 Aug '10

27 Aug '10

Hi, I am refining my first structure (Cd-SAD) with the Phenix (GUI) and would like to refine anomalous f' and f". When i check the box, a menu appears and I have no idea of how to define the Cd atoms. I looked at the help and FAQ and found this: refinement.refine.anomalous_scatterers { group { selection = name BR f_prime = 0 f_double_prime = 0 refine = *f_prime *f_double_prime So I thought the name should be input as CD in the gui but when I submitted the refinement I got an error. Any help here would be appreciated? Thanks, John Rose

4 3

file type not supported???
by Joseph Noel 24 Aug '10

24 Aug '10

Hi All, I am having a problem that should be totally simple but I am stuck trying to find out what is happening. I use Text Wrangler on my Mac as a text editor of files. I have used it many many times to edit PDB files. It seems now that some PDB files that I have opened and removed a few lines now show up as an unrecognized file type when trying to read them into Phenix or into COOT. However, when I upload them on the Molprobity site they are read fine and processed. However, even the downloaded file from Molprobity gives the same error. I know it must be something simple but stumped and hoping someone else has experienced the same problem and has a quick fix. Thanks! Joe ___________________________________________________________ Joseph P. Noel, Ph.D. Investigator, Howard Hughes Medical Institute Professor, The Jack H. Skirball Center for Chemical Biology and Proteomics The Salk Institute for Biological Studies 10010 North Torrey Pines Road La Jolla, CA 92037 USA Phone: (858) 453-4100 extension 1442 Cell: (858) 349-4700 Fax: (858) 597-0855 E-mail: noel(a)salk.edu Web Site (Salk): http://www.salk.edu/faculty/faculty_details.php?id=37 Web Site (HHMI): http://hhmi.org/research/investigators/noel.html ___________________________________________________________

2 1

Re: [phenixbb] phenix.ready_set hydrogens not compatable with coot library
by Joseph Noel 24 Aug '10

24 Aug '10

Hi All, Indeed even the latest COOT still relies on the old PDB format and hydrogens fly. Its probably not prudent to use hydrogens in refinement then remove them and add back each time. Kind of defeats the purpose of using them in the first place. Until COOT is updated I found a pretty easy way around things as long as you have access to the internet. Just use the MOLPROBITY site (http://molprobity.biochem.duke.edu/). In addition to it being a fantastic tool for really cleaning up structures, you can quickly upload your refined file from Phenix with hydrogens, then hit the button to downgrade to PDB v23. It creates the appropriate file with renamed hydrogens that work great in COOT. You can then re-upload the COOT rebuilt model and it will "upgrade" naming to the latest PDB version. Download that file and there you go for Phenix refinement. Sounds like a bit but very fast and easy to implement until COOT is upgraded. Joe ___________________________________________________________ Joseph P. Noel, Ph.D. Investigator, Howard Hughes Medical Institute Professor, The Jack H. Skirball Center for Chemical Biology and Proteomics The Salk Institute for Biological Studies 10010 North Torrey Pines Road La Jolla, CA 92037 USA Phone: (858) 453-4100 extension 1442 Cell: (858) 349-4700 Fax: (858) 597-0855 E-mail: noel(a)salk.edu Web Site (Salk): http://www.salk.edu/faculty/faculty_details.php?id=37 Web Site (HHMI): http://hhmi.org/research/investigators/noel.html ___________________________________________________________

5 4

Phenix Structure Soluiton & Refinement webinar
by Angela Criswell 23 Aug '10

23 Aug '10

Dear colleagues, I would like to draw your attention to an upcoming educational webinar to be presented by Paul Adams titled "Structure Solution and Refinement with Phenix" which is scheduled to occur Thursday, August 26th at 10:00 PDT (13:00 PM EST /17:00 UTC/GMT). You can find more information including a registration link at the following site: http://www.rigaku.com/protein/webinars.html Best regards, Angela

2 1

Palestrant, Deborah is out of the office.
by deborah.palestrant＠novartis.com 23 Aug '10

23 Aug '10

I will be out of the office starting 08/23/2010 and will not return until 08/27/2010. I will respond to your message when I return.

1 0

labelit.version bug?
by Ben Eisenbraun 23 Aug '10

23 Aug '10

labelit seems to be working fine, but one of my users pointed this out to me: $ which labelit.version /programs/i386-linux/phenix/1.6.4-486/phenix-1.6.4-486/build/intel-linux-2.6-x86_64/bin/labelit.version $ labelit.version Traceback (most recent call last): File "/programs/i386-linux/phenix/1.6.4-486/phenix-1.6.4-486/labelit/command_line/version.py", line 18, in <module> run() File "/programs/i386-linux/phenix/1.6.4-486/phenix-1.6.4-486/labelit/command_line/version.py", line 7, in run assert (os.path.dirname(build) == os.path.dirname(repos)) AssertionError -ben -- | Ben Eisenbraun | Software Sysadmin | | Structural Biology Grid | http://sbgrid.org | | Harvard Medical School | http://hms.harvard.edu |

2 1

phenix.ready_set hydrogens not compatable with coot library
by amadrona＠uci.edu 23 Aug '10

23 Aug '10

hello, After reading many posts on the phenix and coot boards I realize that the hydrogens produced by phenix.ready_set are not compatable with coot. Specifically, coot does not recognize HB3 and HD3 hydrogens. When I do realspace refinement they fly off into space. I have edited HB3 to read HB2 and HB2 to read HB1 and everything works well. I was wondering if there is a fix for this besides editing the pdb each time you run phenix.ready_set. Is there a way to add HB3 hydrogens to the phenix library (CCP4 library)? Thank you -Yarrow

6 10

MD with phenix.refine?
by Ed Pozharski 20 Aug '10

20 Aug '10

Silly question: is there some way to run what would essentially be an MD simulation with phenix.refine? I mean, turn off X-ray-term by setting tardy.target_weights.fix_wxc=0, and then run cartesian dynamics? Ed. -- "I'd jump in myself, if I weren't so good at whistling." Julian, King of Lemurs

2 1

Ringer: A program to detect molecular motions by systematic X-ray electron-density sampling
by Terry Lang 20 Aug '10

20 Aug '10

The Alber lab at UC Berkeley is pleased to release of the code for Ringer version 1.0 (http://ucxray.berkeley.edu/ringer.htm). Ringer is a program to detect molecular motions by systematic X-ray electron-density sampling. The aim of Ringer is to go beyond static structural snapshots of proteins by uncovering structural ensembles in X-ray electron density. This information can reveal not only which parts of proteins are flexible and which parts are rigid, but it also can define alternate conformations that may be important for function. Alternate conformations of binding sites also may afford additional targets for drug design. The Ringer method is described in Lang et al. /Protein Sci/. 2010 Jul; 19(7):1420-31 <http://www.ncbi.nlm.nih.gov/pubmed/20499387>. An application of Ringer, determining the structural underpinnings of the side chain dynamics critical for the function of the enzyme proline isomerase, was published in Fraser JS et al. /Nature/. 2009 Dec 3;462(7273):669-73 <http://www.ncbi.nlm.nih.gov/pubmed/19956261>. Ringer is freely available to academics. -- PaulaTherese Lang Postdoctoral Scholar Alber Lab, UC Berkeley http://ucxray.berkeley.edu/ringer.htm

1 0

hydrogens in riding model create steric clashes after refinement
by amadrona＠uci.edu 20 Aug '10

20 Aug '10

Hi, I am having a problem with hydrogens causing steric clashes after refinement. I added hydrogens to my pdb using: "phenix.ready_set model.pdb". I was given two sets of coordinates (my model and a pdb for Iron atoms) I included this new model in the next round of refinement but I did not include the pdb for the Iron atoms. I ran refinement with the folowing parmeter file: refinement.main.ordered_solvent = true #refinement.main.ordered_solvent.mode=filter_only refinement.refine.strategy = individual_sites+individual_adp #refinement.refine.adp.individual.isotropic = "chain W" refinement.output.prefix = CIN4_refine refinement.output.serial = 18 refinement.output.write_def_file = True resolution is 1.37A with close to 400,000 unique observations (obs/param is still about 4) Stats: Rwork Rfree bonds angles start: .1444 .1723 .016 1.617 end: .1371 .1657 .010 1.301 I also ran wxu and wxc optimization from the command line. I know that the defualt is "riding-model" so that the hydrogen atoms are not refined. However, I thought that the program would be able to avoid steric clashes. Is there an aditional command I need to turn on to account for this? I am currently fixing every hydrogen that clashes on an individual basis in coot. There are also a handfull of hydrogens that exist out in space that I decided to delete. Thank you very much for your help. -Yarrow

4 11

altloc non-bonded interactions
by Paul Smith 19 Aug '10

19 Aug '10

Is it possible to tell phenix.refine to ignore certain non-bonded interactions? I have a structure at 1.1A with some interesting sulfate molecules that refuse to refine well. I suspect two positions are alternately occupied by sulfate in some cases and two water molecules in others. I would like to refine with both (a sulfate and two waters) with each component at 1/2 occupancy. With alternate conformations of a given sidechain, the two alternate locations automatically ignore one another w.r.t. non-bonded interactions. However, different residues, it would seem, do not. Has anyone encountered a similar situation or do the phenix gurus have any input? As an alternative, can I define an alternate sulfate residue with several missing atoms? Thanks in advance, --Paul

3 2

Re: [phenixbb] stuck with phenix refine
by Ian Munro Slaymaker 19 Aug '10

19 Aug '10

I also have been having a problem with phenix.refine freezing with no error message from the GUI. It usually freezes at the first macro cycle. I am using version 1.6.4-486. -Ian

3 4

Constraining carbohydrate conformation
by Kenneth Verstraete 19 Aug '10

19 Aug '10

Dear phenixbb users, I wonder how I can constrain the conformation of a NAG-residue during refinement. I have a dataset to 4.2 Angstrom. After refinement some sugars end up in the boat conformation. I would like to run the refinement so that the Asn side chain and the Asn-NAG linkage is refined (using xyz or torsion angle dynamics) but that the NAG-residue itself is considered as a rigid group. Maybe one way of doing this is constraining the bonds and angles of the NAG residue by the defining them as "custom bonds and angles" with sigma=0? Is there a possibility in phenix to generate these constraints from an ideal NAG-residue automatically? Thank you very much for you advice, Kenneth Verstraete

3 2

clean up temporary files
by Jianghai Zhu 18 Aug '10

18 Aug '10

Hi, I realize that after a phenix run, especially autobuild and autosol, there are a lot of temporary files. Is there a way to clean up the disk? Even the .geo file can be removed after some time, because it is huge and can be easily regenerated. -- Jianghai

4 3

ubuntu 10.04 phenix installation
by vincent rao 18 Aug '10

18 Aug '10

Dear Phenix bb user, i'm trying to install the latest version of phenix (1.6.4-486) on my pc with Ubuntu 10.04 installed. I'm using the*phenix-installer-1.6.4-486-intel-linux-2.6-x86_64-fc13.tar * installer. The installer complains that the * bundle-1.6.4-486-intel-linux-2.6-x86_64.tar.gz* is missing in the phenix-installer-1.6.4-486/bundles folder, and requires *bundle-1.6.4-486-intel-linux-2.6.tar.gz *to be present in this folder. So i simply renamed *bundle-1.6.4-486-intel-linux-2.6-x86_64.tar.gz *to *bundle-1.6.4-486-intel-linux-2.6.tar.gz *and this allowed the installation to complete without any further errors. However, when trying to start Phenix, i get the following error: laptop:~/Documents/installers/phenix-installer-1.6.4-486> phenix Traceback (most recent call last): File "/usr/local/phenix-1.6.4-486/phenix/wxGUI2/command_line/main.py", line 13, in <module> run(sys.argv[1:]) File "/usr/local/phenix-1.6.4-486/phenix/wxGUI2/command_line/main.py", line 7, in run from wxGUI2 import App File "/usr/local/phenix-1.6.4-486/phenix/wxGUI2/App.py", line 6, in <module> import wx.lib.dialogs File "/usr/local/phenix-1.6.4-486/build/intel-linux-2.6/base/lib/python2.6/site-packages/wx-2.8-gtk2-ansi/wx/__init__.py", line 45, in <module> from wx._core import * File "/usr/local/phenix-1.6.4-486/build/intel-linux-2.6/base/lib/python2.6/site-packages/wx-2.8-gtk2-ansi/wx/_core.py", line 4, in <module> import _core_ ImportError: /usr/local/phenix-1.6.4-486/build/intel-linux-2.6/base/lib/python2.6/site-packages/wx-2.8-gtk2-ansi/wx/_core_.so: wrong ELF class: ELFCLASS64 Not sure if this is of any use in solving the problem, but i typed the following command phenix.wxpython -c "import sys; print sys.path" and got the following output; ['', '/usr/local/phenix-1.6.4-486', '/usr/local/phenix-1.6.4-486/cctbx_project', '/usr/local/phenix-1.6.4-486/reel', '/usr/local/phenix-1.6.4-486/phaser', '/usr/local/phenix-1.6.4-486/tntbx', '/usr/local/phenix-1.6.4-486/phenix', '/usr/local/phenix-1.6.4-486/elbow', '/usr/local/phenix-1.6.4-486/PyQuante', '/usr/local/phenix-1.6.4-486/docutils', '/usr/local/phenix-1.6.4-486/PyCifRW', '/usr/local/phenix-1.6.4-486/cctbx_project/boost_adaptbx', '/usr/local/phenix-1.6.4-486/cctbx_project/libtbx/pythonpath', '/usr/local/phenix-1.6.4-486/build/intel-linux-2.6/lib', '/usr/local/phenix-1.6.4-486/build/intel-linux-2.6/base/lib/python26.zip', '/usr/local/phenix-1.6.4-486/build/intel-linux-2.6/base/lib/python2.6', '/usr/local/phenix-1.6.4-486/build/intel-linux-2.6/base/lib/python2.6/plat-linux2', '/usr/local/phenix-1.6.4-486/build/intel-linux-2.6/base/lib/python2.6/lib-tk', '/usr/local/phenix-1.6.4-486/build/intel-linux-2.6/base/lib/python2.6/lib-old', '/usr/local/phenix-1.6.4-486/build/intel-linux-2.6/base/lib/python2.6/lib-dynload', '/usr/local/phenix-1.6.4-486/build/intel-linux-2.6/base/lib/python2.6/site-packages', '/usr/local/phenix-1.6.4-486/build/intel-linux-2.6/base/lib/python2.6/site-packages/Numeric', '/usr/local/phenix-1.6.4-486/build/intel-linux-2.6/base/lib/python2.6/site-packages/PIL', '/usr/local/phenix-1.6.4-486/build/intel-linux-2.6/base/lib/python2.6/site-packages/wx-2.8-gtk2-ansi'] thanks

4 3

geometry weight and rmsd bonds - angles
by Lionel Costenaro 17 Aug '10

17 Aug '10

Hello, I'm using phenix.refine for the first time and have some problem to get it reach reasonable rmsd bonds (0.012-0.015) and angles (1.5-1.8). I set the bonds rmsd_max and angles rmsd_max to these values and try different values for: wxc_scale (from default 0.5 downto 0.015), wxu_scale (from default 1 down to 0.03), wc (from default 1 up to 8), with and without optimising wxc wxu. Whatever the parameters, phenix never reached final rmsd better than 0.022 and 2.2. I refine the structure with refmac so far, but need to refine the occupancies (so phenix). Am I doing something wrong? Should I just use phenix to refine only the occupancy (I am actually doing it)? Any help or comment are welcome. Some information: space group: P21, 2 molecules (180 aa) / au, 510 waters, about 20% alternate conformations, ligand covalently link to the protein (correct dictionary and link parameters for refmac or phenix -add link to cif), resolution: 1.32, refine xyz for all, Baniso for all (refmac), ncs not used, R 13.6 Rfree 17.7 fom 0.913 rmsd 0.012 1.56 (last refmac) phenix version 1.6 289 (gui) Best regards, Lionel

8 20

HKL_base
by Ed Pozharski 16 Aug '10

16 Aug '10

MTZ-files generated by phenix (at least some of them) are missing the base dataset (HKL_base) which is a requirement since CCP4 5.0. Is this deliberate? Ed. -- "I'd jump in myself, if I weren't so good at whistling." Julian, King of Lemurs

3 6

AutoSol scaling error
by Yogesh Gupta 16 Aug '10

16 Aug '10

I got following error while running MIR in phenix: "Failed to carry out AutoSol_scale_and_analyze_MIR: list index out of range". Could some one tell me how to fix this? Thanks, Yogesh

2 1

Re: [phenixbb] Restraints file for a new amino acid residue
by Joseph Noel 16 Aug '10

16 Aug '10

Pavel, I sent you the associated pdb file as an attachment. The site in question is 169 in both monomer A and monomer B. I think one way around the problem for me is to create alternative conformations for the residues on either end of the modified cysteine. In that way I can refine each as a section of alternative conformations and then I think things should link up both in Phenix and COOT. COOT is particularly problematic. It only recognizes conformer A and not any additional conformers for a particular position if the name of the amino acid changes. I guess my case is not that common - ie a mixture of oxidation states of a regular amino acid. Joe ___________________________________________________________ Joseph P. Noel, Ph.D. Investigator, Howard Hughes Medical Institute Professor, The Jack H. Skirball Center for Chemical Biology and Proteomics The Salk Institute for Biological Studies 10010 North Torrey Pines Road La Jolla, CA 92037 USA Phone: (858) 453-4100 extension 1442 Cell: (858) 349-4700 Fax: (858) 597-0855 E-mail: noel(a)salk.edu Web Site (Salk): http://www.salk.edu/faculty/faculty_details.php?id=37 Web Site (HHMI): http://hhmi.org/research/investigators/noel.html ___________________________________________________________

1 0

Re: [phenixbb] Restraints file for a new amino acid residue
by Joseph Noel 13 Aug '10

13 Aug '10

Thanks all for the advice. Basically that is what I have been doing. It is actually an oxidized cysteine to sulfenic acid. I am refining the site as alternative conformers with the normal Cys, the single oxidation state to sulfenic acid (CSO) and the double to sulfinic acid (CSD). I know that CSO is in the phenix library but it adds hydrogens to the amino group every time I run elbow or ready set. I then I have to edit the PDB file to delete the hydrogen or Phenix will not refine. I also have problems in COOT. It will not let me move the residues torsional angles. I have checked other structures in the PDB and the density fits perfectly with the structures of many such oxidized systems but for some reason Phenix is giving me some problems. I was going to create a "new" residue and explicitly define the oxygen as charged rather then protonated since the pKa of sulfenic acid is quite low. Maybe the problems all stem from refining this particular position with three alternative "conformers", each conformer of which is chemically distinct. For now, I have CSO labeled as convention BCSO for conformer B. Cys is ACYS. Joe ___________________________________________________________ Joseph P. Noel, Ph.D. Investigator, Howard Hughes Medical Institute Professor, The Jack H. Skirball Center for Chemical Biology and Proteomics The Salk Institute for Biological Studies 10010 North Torrey Pines Road La Jolla, CA 92037 USA Phone: (858) 453-4100 extension 1442 Cell: (858) 349-4700 Fax: (858) 597-0855 E-mail: noel(a)salk.edu Web Site (Salk): http://www.salk.edu/faculty/faculty_details.php?id=37 Web Site (HHMI): http://hhmi.org/research/investigators/noel.html ___________________________________________________________

2 1

Autobuild error message
by bschmidt＠berkeley.edu 13 Aug '10

13 Aug '10

Hello all, I am trying to use phenix.autobuild from the command line to generate either composite omit maps or prime and switch maps. In either case, when I try, the program starts to run for a while before failing and giving me the following error message: "Failed to carry out AutoBuild_build_cycle: Sorry...the file AutoBuild_run_2_/TEMP0/refine_pdb_in.mtz is missing..." I'm not sure why that would be the case, or where that mtz is supposed to come from. Can anyone explain what might be the problem? Thanks, Bryan

2 2

Phaser SAD module - resolution question
by Peter Grey 13 Aug '10

13 Aug '10

Dear Phenix community, I am using Phaser SAD module with a poor MR model to locate the initial heavy-atom sub-structure. I have a 3A dataset but the anomalous signal is weak at the higher resolution shells. I would like to ask for your advice regarding the resolution I should use in Phaser. I thought that for the sake of just locating the sites I should use relatively low resolution (let's say 3.5 or even 4A) and then run Phaser a second time. In this second run I would use the highest-resolution and let Phaser only refine the sites (located in the first run) during one cycle and calculate the combined phases. Thank you in advance for your suggestions and advice, Peter

2 1

Where to find optimized wxc_scale/wxu_scale weighting factors after optimization
by amadrona＠uci.edu 13 Aug '10

13 Aug '10

Hello, I ran a refinement optimizing geometry and B factor weighting factors using optimize_wxc=true optimize_wxu=true. I was given these statistics in the pdb header: Rwork Rfree Bonds Angles start .1588 .1838 .007 1.235 end .1454 .1780 .016 1.601 (1.37 A resolution) 1. I would like to find the optimized weighting factors in the log file but cannot. Instead I find the following (pasted below). However, the scale here doesn't make sense to me because many our greater than 0.5. How do I find the optimized values so that I can use them in further rounds of refinement? I'm sorry if this has been covered in previous posts but I could not find an answer. 2. Can anyone explain the concept of optimizing the B-factor? I don't understand how changing it brings you closer or further from an ideal model as oppose to xray. (which is how I make sense of refining the geometry) Thank you -Yarow ================================ xyz refinement =============================== Start r_work= 0.1462 r_free= 0.1784 bonds= 0.015 angles= 1.53 r_work= 0.1492 r_free= 0.1792 bonds= 0.005 angles= 1.08 scale= 1.00 r_work= 0.1484 r_free= 0.1788 bonds= 0.006 angles= 1.14 scale= 1.40 r_work= 0.1480 r_free= 0.1787 bonds= 0.007 angles= 1.19 scale= 1.80 r_work= 0.1476 r_free= 0.1785 bonds= 0.008 angles= 1.23 scale= 2.20 r_work= 0.1473 r_free= 0.1785 bonds= 0.008 angles= 1.26 scale= 2.60 r_work= 0.1471 r_free= 0.1784 bonds= 0.009 angles= 1.29 scale= 3.00 r_work= 0.1469 r_free= 0.1784 bonds= 0.010 angles= 1.32 scale= 3.40 r_work= 0.1467 r_free= 0.1783 bonds= 0.010 angles= 1.34 scale= 3.80 r_work= 0.1465 r_free= 0.1782 bonds= 0.011 angles= 1.36 scale= 4.20 r_work= 0.1464 r_free= 0.1782 bonds= 0.011 angles= 1.39 scale= 4.60 r_work= 0.1463 r_free= 0.1782 bonds= 0.011 angles= 1.40 scale= 5.00 r_work= 0.1461 r_free= 0.1781 bonds= 0.012 angles= 1.45 scale= 6.00 r_work= 0.1458 r_free= 0.1781 bonds= 0.013 angles= 1.49 scale= 7.00 r_work= 0.1456 r_free= 0.1781 bonds= 0.014 angles= 1.53 scale= 8.00 r_work= 0.1454 r_free= 0.1781 bonds= 0.015 angles= 1.57 scale= 9.00 r_work= 0.1453 r_free= 0.1780 bonds= 0.016 angles= 1.60 scale= 10.00 r_work= 0.1495 r_free= 0.1794 bonds= 0.005 angles= 1.05 scale= 0.85 r_work= 0.1501 r_free= 0.1797 bonds= 0.005 angles= 1.02 scale= 0.70 r_work= 0.1507 r_free= 0.1800 bonds= 0.004 angles= 0.98 scale= 0.55 Final r_work= 0.1453 r_free= 0.1780 bonds= 0.016 angles= 1.60

4 6

Restraints file for a new amino acid residue
by Joseph Noel 13 Aug '10

13 Aug '10

Hi All, How does one create a new CIF file using Elbow for a non-natural amino acid not found yet in the PDB? I know how to carry it out for the isolated residue readily (but of course thinks it is isolated and doesn't include the necessary peptide bonding) but how does one ensure the correct backbone restraints for it are included since it is not isolated but part of the chain? Joe ___________________________________________________________ Joseph P. Noel, Ph.D. Investigator, Howard Hughes Medical Institute Professor, The Jack H. Skirball Center for Chemical Biology and Proteomics The Salk Institute for Biological Studies 10010 North Torrey Pines Road La Jolla, CA 92037 USA Phone: (858) 453-4100 extension 1442 Cell: (858) 349-4700 Fax: (858) 597-0855 E-mail: noel(a)salk.edu Web Site (Salk): http://www.salk.edu/faculty/faculty_details.php?id=37 Web Site (HHMI): http://hhmi.org/research/investigators/noel.html ___________________________________________________________

3 2

too many outliers
by intekhab alam 13 Aug '10

13 Aug '10

Hi there I have recently solved a structure at 2.71A resolution using CCP4 (Refmac5) with R and R free of 22 and 29 respectively. The model is nicely fitted in the map and everything looks except in some loop region where the map is not fitted well. I just finished modelling and checked ramachandran plot which is showing unreasonably high number of outliers (12%) with only 65% in the favoured region. I am thinking to use phenix_refine with simulated annaeling and using the option fix bad chain rotamers. What else i can do aprt from this. i will be highly thankful for the suggestions. Apart from this can anyone let me know how -- INTEKHAB ALAM LABORATORY OF STRUCTURAL BIOINFORMATICS KOREA UNIVERSITY, SEOUL

2 1

wxu for anisotropic ADP refinement
by Joseph Noel 12 Aug '10

12 Aug '10

Dear Phenix User Community, If one would like to try different weights for anisotropic ADP refinement (ie fixing wxu) should one be focusing on numbers like - wxu = 199.878301 - instead of the usual low numbers such as - wxu = 0.25 - found during isotropic ADP refinement? While optimize X-ray/ADP weight is useful in phenix, I am still getting anistropic ADP values that are a bit outside what one might expect at 1.5 A compared to the PDB. For the latter, I check things using Parvati (http://skuld.bmsc.washington.edu/parvati/). Thanks! Joe ___________________________________________________________ Joseph P. Noel, Ph.D. Investigator, Howard Hughes Medical Institute Professor, The Jack H. Skirball Center for Chemical Biology and Proteomics The Salk Institute for Biological Studies 10010 North Torrey Pines Road La Jolla, CA 92037 USA Phone: (858) 453-4100 extension 1442 Cell: (858) 349-4700 Fax: (858) 597-0855 E-mail: noel(a)salk.edu Web Site (Salk): http://www.salk.edu/faculty/faculty_details.php?id=37 Web Site (HHMI): http://hhmi.org/research/investigators/noel.html ___________________________________________________________

2 1

difference in estimated error between phenix and refmac
by Ling Qin 12 Aug '10

12 Aug '10

Dear Colleagues, I am refining a 2.15 A structure that was previously refined with refmac by someone else. Without changing anything in the refmac-refined model and using all the default phenix.refine settings, I am getting pretty much the same R and Rfree (from 19.3%/21.7% to 18.9%/21.8%), with an ML estimated error of 0.22 A. In the refmac-refined pdb file, the estimated error is much lower: 0.16 (based on R), 0.144 (based on Rfree), and 0.096 (based on ML). I am using phenix.refine version 1.6_289, but also tried 1.5_2 with the same result. Is there an underlying reason for this difference? Thank you very much for your insightful opinions in advance. Cheers, Ling

3 2

autosol
by Abd Ghani Abd Aziz 12 Aug '10

12 Aug '10

Dear All, Just been trying to run autosol for SAD data with the new phenix (1.6.4-486) on ubuntu, but I get an error message: *********************************************************** Failed to carry out AutoSol_scale_and_analyze_mad: None of the solve versions worked ********************************************************************** xtriage appears to run OK, and I get the same error with unmerged or merged scalepack data. The same files in 1.6-289 (admittedly on OSX) work fine. I wonder if there is a known problem with this version of phenix? Thanks for any suggestions, Ghani

3 4

Disappearing density?
by Kyle Dolan 11 Aug '10

11 Aug '10

Hi all, I attempted to improve the density around a helix and short linker region using an average kick map but it didn't make things look much better. I then tried creating a composite iterative-build omit map from phenix.autobuild, but in this case the electron density around the helix/linker in question largely disappeared! Is there another approach I can use that might improve my maps, or should I start questioning my data? The dataset has >95% completeness to 2.7 A so I am very curious as to why I am missing the backbone now, where before I was only missing the side chain densities. One thing I am wondering is if the problem could have something to do with symmetry; the part of the map I am looking at may be close to a two-fold crystallographic symmetry axis and I have been told that these regions tend to be very noisy and hard to model--comments? Kyle Kyle T. Dolan Department of Biochemistry and Molecular Biology The University of Chicago ktd(a)uchicago.edu

4 6

Statistics for pdb deposition
by Elwood A. Mullins 11 Aug '10

11 Aug '10

I apologize if this has already been addressed. I'm trying to deposit a model to the pdb database, but there are a few required values in the Refinement Statistics section that I'm not sure where to find. 1,2) Cut-off sigma (F), Number unique reflections (observed) Looking at the log file, it appears that the cut-off is 0, which results in 330 reflections being rejected. The pdb header, however, lists the number of reflections at a value corresponding to 352 reflections being rejected. Why are these numbers different and which should I report? 3) R-factor (all) R-factor (all) includes both the working and test sets in addition to reflections rejected by the sigma cut-off. Does PHENIX calculate this value? I know R-factor (observed), which includes the working and test sets but not rejected reflections, is calculated. Thanks, Elwood

2 1

elbow: more detailed website : 404
by Lepore, Bryan 11 Aug '10

11 Aug '10

this link : http://www.phenix-online.org/elbow/docs/index_image.html is 404 the link is found on this page : http://www.phenix-online.org/documentation/elbow.htm ... IS there a "more detailed website"? i saw the reel site. -bryan

2 1

Random seeds?
by Joseph Noel 10 Aug '10

10 Aug '10

Hi All, Okay, a fairly mundane question concerning random seeds in Phenix. Obviously they are random but any restrictions on how many digits phenix will accept? Any other restrictions? Are they (it) only used for coordinate refinement or do the random seeds get used for other things as well? I know about hydrogen building but wandering about other refinement protocols and the use of "new" random seeds. I realize some use different random seeds to create slightly different models after refinement. Is it useful to change the number from time to time when going through iterative rounds of re-building and coordinate refinement? What about ADP refinement? Sorry for the mundane questions! Joe ___________________________________________________________ Joseph P. Noel, Ph.D. Investigator, Howard Hughes Medical Institute Professor, The Jack H. Skirball Center for Chemical Biology and Proteomics The Salk Institute for Biological Studies 10010 North Torrey Pines Road La Jolla, CA 92037 USA Phone: (858) 453-4100 extension 1442 Cell: (858) 349-4700 Fax: (858) 597-0855 E-mail: noel(a)salk.edu Web Site (Salk): http://www.salk.edu/faculty/faculty_details.php?id=37 Web Site (HHMI): http://hhmi.org/research/investigators/noel.html ___________________________________________________________

3 2

Restraint creation for ligand
by Hermella Woldemdihin 10 Aug '10

10 Aug '10

Hi! I was trying to fit my ligands into my model and density. And I have successfully fitted 2 ligands (same molecule) but my 3rd ligand is different molecule. After fitting my 3rd ligand (with COOT) and merged it with my protein, II created a restraint file with: phenix.elbow my_model_coot-0..pdb --do-all --output=all_ligands It created all_ligand.cif, which I used for in my next refinement with Phenix (GUI). I supplied my density, model and the restraint file and run phenix.refine. Unfortunatate I'm getting this error message when it starts to run phenix.refine and opening COOT window: Sorry: Fatal problems interpreting PDB file: Number of atoms with unknown nonbonded energy type symbols: 1 Please edit the PDB file to resolve the problems and/or supply a CIF file with matching restraint definitions, along with apply_cif_modification and apply_cif_link parameter definitions if necessary (see phenix.refine documentation). Also note that phenix.elbow is available to create restraint definitions for unknown ligands. I have use the same procedure for fitting my first 2 ligands. and it worked. But this time it failed! What do you think I should do? Thanks Hermi

3 2