- phenixbb - phenix-online.org

sculptor : alternate conformations
by Bryan Lepore 22 Nov '10

22 Nov '10

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

2 1

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

20 Nov '10

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

1 0

phenix refinement question
by Katarina Moravcevic 20 Nov '10

20 Nov '10

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

6 10

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

19 Nov '10

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

2 1

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

19 Nov '10

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

1 0

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

19 Nov '10

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

1 0

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

19 Nov '10

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

2 3

CRYST card being changed
by Nino Campobasso 18 Nov '10

18 Nov '10

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

2 1

autobuild : rigid?
by Bryan Lepore 18 Nov '10

18 Nov '10

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

6 6

reflection file converter
by r n 17 Nov '10

17 Nov '10

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

2 1

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

17 Nov '10

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

2 1

phaser and automr
by Bryan Lepore 16 Nov '10

16 Nov '10

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

3 2

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

16 Nov '10

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

2 2

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

16 Nov '10

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

3 3

refine : annealing temperature
by Bryan Lepore 13 Nov '10

13 Nov '10

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

2 1

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

13 Nov '10

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

3 6

automr : ensemble on command line
by Bryan Lepore 12 Nov '10

12 Nov '10

dev-572 i keep getting this error : The number of entries in identity_list ([80.0, 80.0]) needs to be the same as the number of entries in ensemble_1.coords (1) on the command line, i specify the ensemble as : coords='MR_probe_5_match_to_9.pdb MR_probe_9.pdb' identity='80 80' it gets picked up like this : ... ensembleID = "ensemble_1" copies_to_find = None coords = "MR_probe_5.pdb MR_probe_9.pdb" identity = 80 80 ... i tried using quotes or not, and rmsd. i'm using "\" for the shell - what am i missing here? any clues to look for? -Bryan

3 5

Two Post-Doc opportunities in cell biology and biochemistry in Milan, Italy
by Sebastiano Pasqualato 11 Nov '10

11 Nov '10

Hi all, I'm sending this on behalf of a colleague, so please reply to giorgio.scita(at)ifom-ieo-campus.it if interested. PS: to those concerned, apologies for the double post. -- Sebastiano Pasqualato, PhD IFOM-IEO Campus Dipartimento di Oncologia Sperimentale Istituto Europeo di Oncologia via Adamello, 16 20139 - Milano Italy tel +39 02 9437 5094 fax +39 02 9437 5990

1 0

phenix.autosol : scala unmerged ok?
by Bryan Lepore 10 Nov '10

10 Nov '10

[ phenix-dev-572 ] i gave autosol unmerged I, SIGI, M/ISYM from scala for a derivative and a native. from the docs, it appears this is automated. but it appears that solve thinks its scalepack_merge, and other things (LAST.LOG) what keywords should i consider here? should i just use scalepack format? -bryan

2 3

autosol/phaser : weak site search and sites
by Bryan Lepore 10 Nov '10

10 Nov '10

dev-572 is there a keyword to disable phaser's weak site search in an autosol run? it seems that using sites=1, max_single_sites=1, sites_file=one_site.pdb isn't doing it. -Bryan

2 1

autosol : space group
by Bryan Lepore 10 Nov '10

10 Nov '10

dev-572 i gave P 2 21 21 data to autosol - seems OK until the .pdb's get written out - they have P 21 21 2 instead. will this be of any consequence - or should i use a standard setting? -Bryan

4 5

autosol : crystal info
by Bryan Lepore 09 Nov '10

09 Nov '10

dev-572 i gather that a .params for isomorphous replacement gets one cell in the crystal_info scope? -bryan

2 1

C and C++
by Hena Dutta 09 Nov '10

09 Nov '10

Dear Members, Sorry for unrelated questions. I thought some of you might have information could be helpful for me. Is there any free website for learning C and C++? Thanks... Hena

2 2

geometry
by Jianghai Zhu 09 Nov '10

09 Nov '10

Hi all, I am refining a 2.75 A structure. The wxc_scale is 0.02, optimize_wxc is false. But the geometry I got doesn't make sense at all. Any explanations? ------------------------------------------------------------------------ stage angl bond chir dihe plan repu geom_target 0 : 0.617 0.004 0.039 14.275 0.003 4.123 5.8886e-02 1_bss: 0.617 0.004 0.039 14.275 0.003 4.123 5.8886e-02 1_xyz: 3.391 0.059 0.209 16.068 0.010 4.114 1.3055e+00 1_adp: 3.391 0.059 0.209 16.068 0.010 4.114 1.3055e+00 2_bss: 3.391 0.059 0.209 16.068 0.010 4.114 1.3055e+00 2_xyz: 2.435 0.032 0.180 15.961 0.007 4.114 5.4206e-01 2_adp: 2.435 0.032 0.180 15.961 0.007 4.114 5.4206e-01 3_bss: 2.435 0.032 0.180 15.961 0.007 4.114 5.4206e-01 3_xyz: 2.139 0.025 0.191 15.477 0.007 4.114 3.7907e-01 3_adp: 2.139 0.025 0.191 15.477 0.007 4.114 3.7907e-01 ------------------------------------------------------------------------ Thanks. -- Jianghai

2 2

Re: [phenixbb] question on the missing fobs filled maps
by Alexandre OURJOUMTSEV 09 Nov '10

09 Nov '10

Hello, Continuing and confirmating Pavel's last comments: yes, we saw a number of practical cases when non-uniformly distributed missed data, even a small percent, caused great map distortions (some example are listed in our manuscript recently submitted to J.Appl.Cryst). To track this situation, recently we have developed a small python-based program that starting from your MTZ file searches for connected regions of such unmeasured reflections, gives their characteristics and visualizes the regions on the user's request. Presence of such regions may make your maps ugly even when the phases are perfect. For an official release when ready the program will be available at the Web site of the institute; for a time being if you want to test the current version please send me a mail to sacha(a)igbmc.fr With best regards, Sacha Urzhumtsev -----Message d'origine----- De : phenixbb-bounces(a)phenix-online.org [mailto:[email protected]] De la part de Pavel Afonine Hi Fengyun, > I am interesting in that at what completeness of the dataset, will one > use the missing fobs filled map for model building confidently without > too much bias included? I'm not aware of any systematic study on this matter, although at some point I reviewed the available literature. There are numerous examples of how the data incompleteness distorts the map, and literature that discusses this. Interestingly, sometimes much smaller amount of systematically missing reflections, such as plane or cone in reciprocal space, may have much drastic effect than a larger amount of randomly missing data (if you are "lucky" enough it can mask entire structural domain).

2 1

moving user home directories
by Ben Eisenbraun 09 Nov '10

09 Nov '10

Hi PHENIXians, Our local sysadmins moved everyone's home directory from /nfs/home/foo to /nfs/userdocs/home/foo, and it appears to have broken most people's PHENIX projects. Is there a straightforward way to fix this? Poking around in the ~/.phenix/project_db.pkl files, I can see the old paths, but I'm not exactly sure how to rewrite. My brute force attempt with sed was unsuccssful. Thanks. -ben -- | Ben Eisenbraun | Software Sysadmin | | Structural Biology Grid | http://sbgrid.org | | Harvard Medical School | http://hms.harvard.edu |

10 15

Nov 15, 2010 deadline- User proposal submission for Collaborative Crystallography at BCSB
by Banumathi Sankaran 08 Nov '10

08 Nov '10

> Dear Users, > > The deadline for Jan/Feb 2010 Collaborative > Crystallography proposals will be *Nov 15, 2010. * > > ---------------------------------------------------------------------------------------- > Through the Collaborative Crystallography Program (CC) at the > Advanced Light Source (ALS), scientists can send protein crystals > to Berkeley Center for Structural Biology (BCSB) staff researchers > for data collection and analysis. The CC Program can provide a > number of benefits to researchers: > > * Obtain high quality data and analysis through collaborating with > expert beamline > researchers; > * Rapid turn around on projects; and > * Reduced travel costs. > > To apply, please submit a proposal through the ALS General User > proposal review process for beamtime allocation. Proposals are > reviewed and ranked by the Proposal Study Panel, and beamtime is > allocated accordingly. BCSB staff schedule the CC projects on > Beamlines 5.0.1 and 5.0.2 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 learn more, go to: > > http://www.als.lbl.gov/als/quickguide/becomealsuser.html > > To submit a proposal, go to: http://alsusweb.lbl.gov/. Scroll > down to "*Structural Biology beamlines (includes protein SAXS)*" > and click on "New Proposal." Enter your proposal information, with > attention to the following details: > > * For question 3/First choice, select "5.0.1 > (Monochromatic);" for question 3/Second choice, select "5.0.2 (MAD)." > * Check box (yes) in response to question (7) "Do you want collaborative > crystallography > (beamline 5.0.1 or 5.0.2 only) > * In question 4, please describe a specific research project with a clear > end point. > > In order to request CC time for Jan/Feb 2010 allocation period, > proposals must be submitted by *Nov 15, 2010.* > > The deadline for CC proposals for the time period Mar/April 2010 will be > Jan 15, 2010. > > Regards, > Banumathi Sankaran >

1 1

Fixing a ligand by a symmetry axis
by Simon Kolstoe 04 Nov '10

04 Nov '10

Dear Phenix bb, I have a palindromic ligand bound to my protein that crosses a symmetry axis. Although the head-group that is in close contact with the protein refines beautifully (in phenix..refine version 1.6.4 run from the command line on a mac), the atom closest to the symmetry axis get knocked away from the equivalent atom that it binds to on the other side of the symmetry axis. Is there any way of stopping phenix.refine from refining the atoms closest to the symmetry axis, or should I just manually position the wayward atom in the appropriate orientation prior to my structure deposition? Thanks, Simon

3 5

More trained monkeys please (request for improvements to auto-NCS)
by Kendall Nettles 03 Nov '10

03 Nov '10

Dear PHENIX developers, I am writing to ask about the automated NCS selection procedure. I have a structure with 2 homodimers in the ASU, and there are significant differences in the two monomers that make up teach dimer. When I run automr/autobuild I am very pleased with the result overall, but it does not choose the best NCS groups- it groups all chains to the A chain, instead of one dimer to the other. The result is a helix that is obviously out of place (see attached image). I wonder if this could be fixed by implementing a throw away round of refinement, in order to identify optimal NCS groups? Though it looks like there are only a few areas of significant difference, the structures end up better by grouping one dimer against the other, rather than just the regions that are obviously different. This brings me to my second question. Might it be worthwhile to implement an NCS optimization procedure, which would compare different NCS groups, similarly to the optimize solvents, weights, etc? Kendall Nettles[cid:[email protected].]

1 0

Does phenix have any utility to do the 2mFo-DFc composite omit map calculations?
by zhangh1＠umbc.edu 03 Nov '10

03 Nov '10

Hi, Does phenix have any utility to do the 2mFo-DFc composite omit map calculations (similar to OMIT in CCP4)? Thanks in advance! Hailiang

4 7

phenix.model_vs_data : resolution cutoff from .pdb
by Bryan Lepore 02 Nov '10

02 Nov '10

[ phenix dev-568 ] TIL using phenix.model_vs_data, if the "REMARK 3" records in .pdb are adjusted e.g. to set a resolution limit, e.g REMARK 3 RESOLUTION RANGE HIGH (ANGSTROMS) : [change to whatever high resolution limit] REMARK 3 RESOLUTION RANGE LOW (ANGSTROMS) : 42.511 phenix.model_vs_data will pick this up and compute the R-factors for that range. Thanks Pavel -bryan

1 0

Riding H
by Jianghai Zhu 01 Nov '10

01 Nov '10

Hi all, If I use the riding H atoms come out of phenix.refine in molprobity, it gives me a decent clash score. But if I strip the H atoms and let molprobity to add H atoms, I get a much worse clash score. I saw the manual of phenix recommends the first method. So what is the cause of the difference here? I know long time ago that phenix and molprobity used different bond length to add the riding H atoms. I don't know what the situation is now. Does phenix.refine refine the riding H atom positions? Thanks. -- Jianghai

4 7

cancel that ...
by Phoebe Rice 01 Nov '10

01 Nov '10

Despite the error messages, it really did make restraints that seem to work. I just didn't realize they'd be in a hidden directory. Can I put in a request for no hidden defaults? They make me feel like a trained monkey. thanks, Phoebe

2 1

error in making links
by Phoebe Rice 01 Nov '10

01 Nov '10

Hi BB, I wanted to run phenix on a structure from the database with covalent serine-DNA links, 1zr4. The "ready-set" feature to prepare restraints from the pdb link records sounded lovely, but fails with the message "Error interpreting command line argument parament definition: "1ZR4.link.edits" RuntimeError: Unexpected end of input" This happens with either version 1.6.4-486 or a nightly build installed a few weeks ago (phenix 551; installed on 10/7/2010.) Am I asking to much to get nice Ser-DNA links without sweating blood to read the fine print? thanks, Phoebe ===================================== Phoebe A. Rice Dept. of Biochemistry & Molecular Biology The University of Chicago phone 773 834 1723 http://bmb.bsd.uchicago.edu/Faculty_and_Research/01_Faculty/01_Faculty_Alph… http://www.rsc.org/shop/books/2008/9780854042722.asp

2 1

refine : "merged" in .def
by Bryan Lepore 29 Oct '10

29 Oct '10

[ dev-564 on GNU/Linux ] why does phenix - apparently - keep putting the label "merged" in output .def's, when neither the input .def or .mtz has this label - and what generally is the significance? the .mtz has separate Bijvoets, and I's. e.g: input .def: (concise) labels = "I(+),SIGI(+),I(-),SIGI(-)" label = "R-free-flags" output .def: labels = "I(+),SIGI(+),I(-),SIGI(-),merged" label = "R-free-flags" i re-checked the .mtz lacks "merged" (not shown). -bryan

3 3

Autobuild error - bash/dash/sh
by Stephen Graham 29 Oct '10

29 Oct '10

Hi there, Apologies if this has been reported before, but the autobuild job I tried to run last night failed after a short while... The job was started from the command line using: phenix.autobuild data=../mydata_free.mtz seq_file=../protein.fasta map_file=../ccp4i/protein_phaser.1_parrot2.mtz input_map_labels='parrot.F_phi.F parrot.F_phi.phi None' And it got to the point where it was trying to submit the three jobs before it died as follows: Splitting work into 3 jobs and running with 1 processors using sh background=False in /home/stepheng/Projects/my_protein/autobuild/AutoBuild_run_3_/TEMP0 Starting job 1...Log will be: /home/stepheng/Projects/my_protein/autobuild/AutoBuild_run_3_/TEMP0/RUN_FILE_1.log *************ERROR ENDING ******************* ******************************************************************************** Failed to carry out AutoBuild_build_cycle: child process stderr output: command: 'sh /home/stepheng/Projects/my_protein/autobuild/AutoBuild_run_3_/TEMP0/RUN_FILE_1' /home/stepheng/Projects/my_protein/autobuild/AutoBuild_run_3_/TEMP0/RUN_FILE_1.sh: 5: Syntax error: Bad fd number ******************************************************************************** *************ERROR ENDING ******************* I think the problem might be that the RUN_FILE_1.sh script seems to use a bash command/redirection that is not implemented in plain old 'sh' spec, and thus not supported by dash (the lean version of sh that ships with ubuntu)? If I change the #! line of RUN_FILE_1.sh to call /bin/bash rather than /bin/sh the script seems to run to completion... Is there any easy way to change the script template to use /bin/bash or remove the offending command? I'm running phenix dev-550 on an x86 ubuntu 10.10 machine. Thanks, Stephen P.S. The reason I'm trying to run autobuild from the command line rather from the GUI is that the GUI seems to crash every time I give autobuild a sequence file in either the FASTA or PIR format. Is there another format I should be using? -- Dr Stephen Graham 1851 Research Fellow Cambridge Institute for Medical Research Wellcome Trust/MRC Building Addenbrooke's Hospital, Hills Road Cambridge, CB2 0XY, UK Phone: +44 1223 762 638

4 5

Twin Law question
by Young-Jin Cho 28 Oct '10

28 Oct '10

Dear Phenix users, Thanks in advance to answer my question. As my recent data did not give me perfect statistics, I was wondering what the problem was (one of things is ice ring though). When I ran Xtriage with my struggling current data set and looked like I got twinned data. I reran PHASER as was suggested and confirmed that current space group is right. My question is now how I can handle this problem. First, as was recommended I input twin_law, however not all the cases (I have a couple of data sets) were working. One favorite(!) error message was as following: "Map type 'Fc' not supported for twinned structures. Allowed types: Fo-Fc, Fobs-Fobs-Fmodel, 2mFo-DFc, 2mFobs-DFmodel, mFo-DFc, mFobs-DFmodel, gradient, m_gradient." I got this from both GUI and shell environment. . 1. Now from now on, whenever I want to phenix.refine, should I include twin_law command such as 'twin_law=-h,l,k'? 2. In this case, what do I expect to see such as lowering Rfree or else? If it doesn't help lowering Rfree, what to do next? 3. Do I have to report I did with twinlaw option in case of preparing a paper or presentation? I mean I also heard that if one got twinned data, they re-collected data. How much can I convince my refinement is OK with twined data? Additionally, I checked other data, and applied twin_law but the result R/Rfree were exactly same as the one I refined without twin_law option, it this normal case? With many thanks in advance, Young-Jin Cho

3 2

maximum allowable B-factor value
by Christian Roth 27 Oct '10

27 Oct '10

Dear all, I did a refinement run with a built in ligand. Within the ligand the B-factor jumps between two connected atoms from 60 up ~110, which is i think a bit too much. Is there a way to restrain the B-factors that they could not jump so much. I would accept negative difference density for the atom. Maybe the ligand is a bit degraded. Best regards Christian @Pavel Tanks for your suggestions with occupancy refinement and B-factors. -- Christian Roth Institut für Bioanalytische Chemie Biotechnologisch-Biomedizinisches Zentrum Fakultät für Chemie und Mineralogie Universität Leipzig Deutscher Platz 5 04103 Leipzig Telefon: +49 (0)341 97 31316 Fax: +49 (0)341 97 31319

2 3

Re: [phenixbb] Sigma-A vs Resolution plot
by Bryan Lepore 25 Oct '10

25 Oct '10

On Thu, Oct 21, 2010 at 12:04 PM, Pavel Afonine <pafonine(a)lbl.gov> wrote: > the plot "Sigma-A vs Resolution" will be available in the next PHENIX > nightly build: Looks nice - I greatly appreciate this - thanks Pavel, Nat, anyone else I forgot. -Bryan

1 0

nightly builds
by Bryan Lepore 22 Oct '10

22 Oct '10

is it me, or was there just a dev-562? ------------------------------------------------------------------------------- PHENIX: Python-based Hierarchical ENvironment for Integrated Xtallography Version: dev Release tag: 562 Platform: intel-linux-2.6 linux User: ------------------------------------------------------------------------------- but latest is : dev-555 2010-10-10 successful cctbx, cci_apps, phaser, phenix, misc docs; changelog

2 1

R-free problem with phenix.refine
by Vinson LIANG 22 Oct '10

22 Oct '10

Dear all, Recently I deposited a structure to PDB. However, I was told >> The R value of SF check are different from your reported value. >> >> R work in SF check = 0.2859 >> your reported R work = 0.1474 >> Also there are no R free flags in your sending files, but you reported R >> free value in your all deposition. So, I rerun my deposited file with phenix.refine, and it gave a R-work=0.1475, R-free=0.1742. However, when I tried to run the same file with refmac5, error came with "Free R label has not been set". I guess phenix.refine gives Rree-R label that other program could not recognize. Could anyone please help me to solve this problem? Thank you very much in advance. I don't know why structure refined with phenix.refine has such problem. I wish this problem will disappear in the future, since I so like phenix.refine. Best, Vinson,Liang

4 7

refine : peculiar TLS tensors
by Bryan Lepore 21 Oct '10

21 Oct '10

[ phenix 1.6.4-486 ] [ ALSO : tlsanl, refmac, etc. from ccp4 6.1.3 ] [ molscript/povscript ] the one-liner : i have some TLS tensors that sometimes do not compare at all well with refmac/anisoanl tensors. sometimes, the principal axes of a TLS tensor refines *very* far away from the center of mass(=center of reaction), or the axes from tlsanl are exceptionally small as if not scaled properly. the eigenvalues are really small too. see [*] for troubleshooting and strategy used. i suspect i am missing something. i also understand this might not be a phenix question. -bryan [*] i compared another .def that gave sensible TLS results with a different dataset/model, and my {} are all good AFAIK. anisoanl/refmac/parvati results are reasonable. the phenix refinement appears well behaved. these are all using tlsextract to get the eigenvalues/origins, and molscript to check the axes/centers of mass. input B factors are Beq unless i made an error. the tls group selections can vary - however, simple groups also can give these results. in general, using strategy (excerpt) : refine { strategy = *individual_sites individual_sites_real_space rigid_body \ *individual_adp group_adp *tls *occupancies *group_anomalous sites { individual = None torsion_angles = None rigid_body = None } adp { individual { isotropic = None anisotropic = None } group_adp_refinement_mode = *one_adp_group_per_residue \ two_adp_groups_per_residue group_selection group = None tls = "chain A and (resid 1:41 or resid 42:159 or resid 195:237 or resid 263:420)" tls = "chain A and resid 160:194" tls = "chain A and resid 245:252" } [...cut off ... ] example output of an ill-behaved S tensor from tlsanl : SCREW ROTATION AXIS ABSOLUTE POSITION (A) PITCH (A) X Y Z 1 70.667 154.197 308.627 -144.435 2 977.334 819.588******** 6222.947 3 640.7191431.084 471.341 269.382

3 5

maps : sigma-a
by Lepore, Bryan 21 Oct '10

21 Oct '10

how might one start to go about getting phenix to automatically write out the value of sigma_A vs. resolution? e.g. when making maps, to get that nice table in Randy Read's SIGMAA program output of ln(sigma_A) vs. resolution. -bryan

4 6

Bad ligand parameters
by Joe Krahn 21 Oct '10

21 Oct '10

I found that residue TRS has different dihedral target angles for the three hydroxyl oxygens, which are actually all equivalent. They have a periodicity of 1, but should have a periodicity of 3. What is the best way to handle such discrepancies? Most ligand parameters are probably auto-generated based on a reference structure, and there is probably a lot of conformation bias. In this case, it seems that parameterization should have been able to detect symmetry. It would be nice if parameters could include information about whether a human has done any validation. Actually, a ligand Wiki might be nice, where people can easily put comments, even if they are not sure about how to improve the parameters. Joe Krahn

9 13

refine : a priori TLS tensor eigenvalues
by Bryan Lepore 21 Oct '10

21 Oct '10

[ phenix 1.6.4-486 ] i gather that phenix.refine will not take a priori refined TLS tensor eigenvalues - either from input or from the pdb header. is this right? -bryan

2 3

Joint X-ray/Neutron Refinement
by Louis Lazar 21 Oct '10

21 Oct '10

Hello, Sorry for reposting, I am trying to get a better response with a more detailed subject line. I am seeing an internal error while trying to conduct joint x-ray/neutron refinement, and any help with the issue below would be greatly appreciated. Thank you very much in advance: ================================ ADP refinement=============================== Traceback (most recent call last): File "/sware/phenix-1.6.4-486/linux64/phenix-1.6.4-486/phenix/phenix/command_line/refine.py", line 11, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/sware/phenix-1.6.4-486/linux64/phenix-1.6.4-486/phenix/phenix/refinement/command_line.py", line 89, in run call_back_handler=call_back_handler) File "/sware/phenix-1.6.4-486/linux64/phenix-1.6.4-486/phenix/phenix/refinement/driver.py", line 1182, in run call_back_handler = call_back_handler) File "/sware/phenix-1.6.4-486/linux64/phenix-1.6.4-486/phenix/phenix/refinement/strategies.py", line 735, in refinement_machine h_params = h_params) File "/sware/phenix-1.6.4-486/linux64/phenix-1.6.4-486/cctbx_project/mmtbx/refinement/adp_refinement.py", line 192, in __init__ all_params = all_params) File "/sware/phenix-1.6.4-486/linux64/phenix-1.6.4-486/cctbx_project/mmtbx/refinement/adp_refinement.py", line 216, in refine_adp model.set_refine_individual_adp() File "/sware/phenix-1.6.4-486/linux64/phenix-1.6.4-486/cctbx_project/mmtbx/model.py", line 1035, in set_refine_individual_adp iselection = selection_iso.iselection()) RuntimeError: cctbx InternalError: /net/anaconda/scratch1/phenix/phenix-1.6.4-486/cctbx_project/cctbx/xray/scatterer_flags.h(417): CCTBX_ASSERT(f.use_u_iso()) failure. This error occurs during the 2nd macro cycle of refinement using this general command: phenix.refine model.pdb data_xray.hkl neutron_data.file_name=data_neutron.hkl input.xray_data.labels=FOBSx input.neutron_data.labels=FOBSn Thanks again. Regards, Louis

2 1

Stretched Maps
by Joel Weadge 21 Oct '10

21 Oct '10

Hi all,I have just finished doing molecular replacement and an initial round of refinement (with simulated annealing) on a number of different datasets. They have all be straightforward since most of them are either point mutants and/or co-crystal structures. However, two of the structures are giving me some trouble. The data appears to be of good quality and resolution (2.25 and 2.55Å). The MR proceeds as normal with good scores and then in the subsequent round of refinement the maps (2Fo-Fc and Fo-Fc) look strange. The density is stretched or flattened in one direction. Can anyone help me with what I may be overlooking?Thanks in advanceJW

5 4

occupancy refinement of a metal ion
by Christian Roth 20 Oct '10

20 Oct '10

Dear all, if I understood the manual correctly I could refine the occupancy and set the B-Factor to a constant value, or vice versa. Is there a possibility that phenix refine the B-Factor automatically to the value of the coordinating residues and adjust than the occupancy to a good value? I seems to be that the metal site is not fully occupied. I guess I have to set the B-value manually to the value of the surrounding residues and do than a refinement, right? Best regards Christian

2 1

fix rotamers question
by Peter Grey 20 Oct '10

20 Oct '10

Dear Phenix community, I would like to have some additional information regarding fix_rotamers. In particular : - How Phenix chooses which rotamers need to be fixed ? - Does Phenix really re-calculate the three maps after correcting each rotamer ? For very large structures this is very time-consuming and one can settle for calculating a new map every 10-20 residues. Perhaps this could be a user-controlled parameter ? -Can you please explain how fix_rotamers works together with ncs restraints ? Thank you in advance for shedding light on this, Peter

3 3

about low resolution refinement
by crystallogrphy 18 Oct '10

18 Oct '10

Hi, I am refining a low resolution data. According to sigma vs. resolution, my resolution is 4.7A. Because there are four identical molecular per asymetric unit, I can use positional refinement. There are about 65% alpha-helices in the structure, so I want to use secondary structure restrain. Does phenix has this restrain in the refinement? I am using Phenix-1.4 version.

4 4