- phenixbb - phenix-online.org

elbow builder / phenix question
by Sudharsan Sridharan 20 Sep '07

20 Sep '07

Hi all, Is it ok to use the cif files generated by PRODRG server directly in phenix refine ? I fit a ligand molecule in density in Coot using this cif file. but after refinement in phenix using the same cif file the ligand coordinates does not look right. Then I created a new cif file using elbow.builder (with final geometry option) for the coordinates fit in Coot. But I see differences in the bond properties in the new cif file compared with the one from PRODRG server. e.g. please see below for parts of the cif generated in these two. could you please advice how to get the right cif file using the server or phenix and do the refinement. in a related question, is the command to use for building the cif in phenix phenix.elbow or elbow.builder. the former is given in http://www.phenix-online.org/new-site/documentation/elbow.htm but didn't work (i'm using 1.3 beta) and the latter does work (which is good !). Thanks, Sid. from prodrg: DRG OAB CAK deloc 1.250 0.020 DRG CAK OAD deloc 1.250 0.020 DRG CAK CAM single 1.530 0.025 DRG CAM CAA double 1.390 0.020 DRG CAM OAJ single 1.360 0.022 DRG OAJ CAO single 1.435 0.033 DRG CAO CAI single 1.390 0.020 DRG CAO CAP single 1.530 0.025 DRG CAI CAG double 1.390 0.020 DRG CAG CAH single 1.390 0.020 DRG CAH CAN double 1.390 0.020 DRG CAN CAL single 1.530 0.025 DRG CAN CAP single 1.530 0.025 DRG CAL OAE deloc 1.250 0.020 DRG CAL OAC deloc 1.250 0.020 DRG CAP OAF single 1.430 0.025 DRG OAF HAA single 1.000 0.027 DRG CAA 1HAA single 1.000 0.020 DRG CAA 2HAA single 1.000 0.020 DRG CAO 1HAO single 1.000 0.020 DRG CAI 1HAI single 1.000 0.020 DRG CAG 1HAG single 1.000 0.020 DRG CAH 1HAH single 1.000 0.020 DRG CAP 1HAP single 1.000 0.020 from elbow.builder: ISC OAB CAK double 1.23599 0.02 ISC CAK OAD single 1.23956 0.02 ISC CAK CAM single 1.51571 0.02 ISC CAM CAA single 1.38593 0.02 ISC CAM OAJ single 1.35985 0.02 ISC OAJ CAO double 1.45557 0.02 ISC CAO CAI aromatic 1.38350 0.02 ISC CAO CAP aromatic 1.53664 0.02 ISC CAI CAG aromatic 1.37763 0.02 ISC CAG CAH aromatic 1.37668 0.02 ISC CAH CAN aromatic 1.37587 0.02 ISC CAN CAL single 1.53109 0.02 ISC CAN CAP aromatic 1.52972 0.02 ISC CAL OAE double 1.22830 0.02 ISC CAL OAC single 1.23738 0.02 ISC CAP OAF double 1.41071 0.02 ISC 1HAD OAD single 1.08000 0.02 ISC 1HAM CAM single 1.08000 0.02 ISC 1HAA CAA single 1.08000 0.02 ISC 2HAA CAA single 1.08000 0.02 ISC 3HAA CAA single 1.08000 0.02 ISC 1HAH CAH single 1.08000 0.02 ISC 1HAC OAC single 1.08000 0.02

2 1

wxu scale in phenix.refine
by Jan Abendroth 18 Sep '07

18 Sep '07

Hi all, in the final stages of refining a structure, I would like to optimize the restraint weights for geometry and B factors. It seems to work fine for geometry by defining " refinement.target_weights.wxc_scale=1.0". However, when the respective definition for B-factors does not seem to have an impact at all. I tried both: "refinement.target_weights.wxu_scale=2.0" in the parameter file or "wxu_scale=2.0" in the command line. The main reason why I try to fiddle around with the B-weights is that coot gives me a whole bunch of red bars in the Temperature factor variance analysis... I am using phenix 1.3d (downloaded past Sunday), mac-intel-osx, osx-10.4.10 Any ideas? Cheers Jan

3 5

Re: [phenixbb] Problem in Autosol
by Eric Tee 16 Sep '07

16 Sep '07

Hi Tom, Many thanks for your reply. I was running phenix in my d_drive (a vfat partition) and got the error message. However, when I moved the mtz file to a directory under home directory (home/user/phenix-local), The program work without any problem. Let me explain how this work: I'm running phenix on fedora7 (a dual boot machine with Windows). I installed phenix under home/user/program/phenix (not /usr/local). I have a vfat partition (as d_drive) where I keep all my data and basically process everything there. I got the error message if I run phenix in the vfat partition but the program can be run successfully under the home directory /home/user/phenix-local). Following is what I got if I run phenix in the vfat partition (d_drive): WORKDIR: /mnt/d_drive/ccp4_data/SAD/phenix-localmachine Checking for KnownWizards: ['AutoSol', 'AutoBuild', 'ImportRawData'] application method: Null Getting inputs on first pass... SETTING UP INPUTS AutoSol Wizard Screen 1 of 4 AutoSol: Choose data files for dataset #1 A "dataset" is a group of data with a single set of heavy-atom sites You will be able to combine the results of all datasets later Note: Data files can be any format. If they are not all Scalepack premerged or all Scalepack unmerged then all will be converted to Scalepack premerged SUMMARY Input files dataset #1 Se_pk.mtz Creating new Data Storage /mnt/d_drive/ccp4_data/SAD/phenix-localmachine/PDS/AutoSol.1 sh: TEMP0/temporary.com: Permission denied Could be something related to permission problem but I have set permission to read/write/execute in vfat partition by adding the following line to /etc/fstab: /dev/sda6 /mnt/d_drive vfat users,uid=eric,gid=users,umask=0002,utf8=true 0 0 Any suggestion to enable running phenix in the vfat partition? Eric On 9/16/07, Tom Terwilliger <terwilliger(a)lanl.gov> wrote: > > Hi Eric, > > Can you possibly check to see that you have a recent version of phenix: > > phenix.version will do this. > > > Then if that is ok, I'm not sure what is causing the problem. This > message means that solve could not run with the data it got, but > usually an error message should come up earlier if there is a > problem. Have a look at LAST.LOG in your AutoSol directory and see > if it has any hints. Also look at AutoSol_run_1_.output and see the > end of this log file. > > If none of those work, if you can send me your input data file I can > see if there is a problem with it (all I would need is the first 100 > lines or so). > -Tom T > >

1 0

Problem in Autosol
by Eric Tee 15 Sep '07

15 Sep '07

Hi All, I was trying to run SAD phasing with phenix GUI for the first time. After conversion of mtz to sca file, at the step of autosol_scale_and_analyze_mad, I got the following message: failed to carry out Autosol_scale_and_analyze_mad None of the solve versions worked I'm pretty sure I choose the right package. Should I do something to configure solve? Eric

1 0

elbow.builder to generate charges on individual atoms?
by Francis E Reyes 13 Sep '07

13 Sep '07

Ok, this is probably an unsupported feature(TM),but is there anyway to hack elbow.builder to output an electronic charge (presumably from the AM1 calculation) for each atom of a ligand? (by say substituting in for the Bfactor)) Thanks FR --------------------------------------------- Francis Reyes M.Sc. 215 UCB University of Colorado at Boulder gpg --keyserver pgp.mit.edu --recv-keys 67BA8D5D 8AE2 F2F4 90F7 9640 28BC 686F 78FD 6669 67BA 8D5D

2 1

Position for PSI Knowledgebase development
by Paul Adams 12 Sep '07

12 Sep '07

======================================================================= Lawrence Berkeley National Laboratory Software Developer II ======================================================================= Software Developer II Job ID: 20930 Division: Physical Biosciences Date Opened: 8/7/2007 Summary This position will be part of the Computational Crystallography Initiative at Lawrence Berkeley Lab and will be responsible for developing a web site to disseminate information about the Technologies developed by the NIH Protein Structure Initiative (PSI). This program funds the majority of structural genomics efforts within the US and is responsible for the development of many new technologies that are useful to macromolecular crystallographers. This technology resource is part of the larger PSI Knowledge Base that will provide an information portal for the scientific community. Duties Essential: Under general direction of the Principal Investigator: Develop and maintain web site documenting PSI developed technologies. Interact with PSI researchers to collect information about new technologies. Create a relational database to store crystallographic data. Install and maintain SVN repositories for the PSI community. Marginal: Perform routine analysis of crystallographic data. Maintain web site computing hardware and software. Other duties as assigned. Qualifications Essential: Experience with writing and maintaining HTML web pages. Knowledge of standard programming methods, and Python in a UNIX environment. Strong analytical and organizational skills. Ability to handle multiple tasks in parallel. Excellent interpersonal skills for working in a diverse team environment, ability to interact effectively with other scientists, and to address their needs in an effective and timely manner. Marginal: B.S. Degree or higher in biophysics, computer science, chemistry, electrical engineering, physics or related field preferred or equivalent experience, with demonstrated programming experience. Experience maintaining computing hardware and software. Experience with analysis of crystallographic data. NOTE: This is a one year Term appointment with the possibility of renewal under the same terms and conditions, contingent upon continued funding and availability of work. Lawrence Berkeley National Laboratory is a world leader in science and engineering research, with 11 Nobel Prize recipients over the past 75 years, and 59 present members of the National Academy of Sciences. LBNL conducts unclassified research across a wide range of scientific disciplines and hosts four national user facilities. AA/ EEO employer committed to the development of a safe and diverse workforce. Learn more at http://www.lbl.gov. For more information about the CCI visit: http://cci.lbl.gov/ To apply: visit http://cjo.lbl.gov/ and search for the job number 20930. Please see http://cjo.lbl.gov/app_instr.html for details of how to apply. -- Paul Adams Senior Scientist, Physical Biosciences Division, Lawrence Berkeley Lab Adjunct Professor, Department of Bioengineering, U.C. Berkeley Head, Berkeley Center for Structural Biology Deputy Principal Investigator, Berkeley Structural Genomics Center Building 64, Room 248 Tel: 510-486-4225, Fax: 510-486-5909 http://cci.lbl.gov/paul Lawrence Berkeley Laboratory 1 Cyclotron Road BLDG 64R0121 Berkeley, CA 94720, USA. --

1 0

pseudo-merohedral twinning in P212121
by Schneider Sabine 12 Sep '07

12 Sep '07

Hi everyone, I have a data set to 2.8A, with most likely 2 molecules in the asu and a molecular replacement model with 40% sequence identity. (Cell 71.3 71.3 227.9 90 90 90) Initially the data appeared to be either P4322 or P43212 (Rsym = 7.4%; Cell 71.3 71.3 227.9 90 90 90), but the molecular replacement solutions had major clashes with the symmetry related molecules and a lower LLG for the solution of the second molecule. No twinning was detected by Xtriage. Processing the data in P212121 (Rsym 5.6%, cell 71.3 71.4 228.2 90 90 90) gave me a better molecular replacement solution (LLG 121/373), but I had 5 clashes with the second molecule in the asu. (I tested all alternative spacegroups) I also have a MAD dataset (very low quality!) and I combined the phases with the one obtained from molecular replacement in Sharp. After NCS averaging the maps had improved (still not great). I tried to improve things in the model and run Phaser again and the solution had better statistics (LLG 338/1223) and the clashes were reduced to 2. Xtriage detected pseudo-merohedral twinning with a twin fraction of 43% with twin operator -k, -h, -l. So I tried refining the model with phenix telling it the twinning operator and fraction, but R/Rfree are stuck around 45%. I also tried simulated annealing but that gave even higher R/Rfree. Is there anything I can do, or is it a hopeless case? Thanks for your help in advance! Sabine This message has been checked for viruses but the contents of an attachment may still contain software viruses, which could damage your computer system: you are advised to perform your own checks. Email communications with the University of Nottingham may be monitored as permitted by UK legislation.

2 1

phenix.refine error
by Jianghai Zhu 11 Sep '07

11 Sep '07

Hi, I was trying the twinning refinement by phenix.refine and I got the following error messages during TLS refinement. Traceback (most recent call last): File "/nfs/home/jzhu/cci_apps_sources/phenix/phenix/command_line/ refine.py", line 5, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/nfs/home/jzhu/cci_apps_sources/phenix/phenix/refinement/ command_line.py", line 81, in run call_back_after_monitor_collect=call_back_after_monitor_collect) File "/nfs/home/jzhu/cci_apps_sources/phenix/phenix/refinement/ driver.py", line 1244, in run log = log) File "/nfs/home/jzhu/cci_apps_sources/phenix/phenix/refinement/ strategies.py", line 866, in refinement_machine h_params = h_params) File "/nfs/home/jzhu/cci_apps_sources/mmtbx/mmtbx/refinement/ adp_refinement.py", line 175, in __init__ macro_cycle = macro_cycle) File "/nfs/home/jzhu/cci_apps_sources/mmtbx/mmtbx/tls/tools.py", line 600, in __init__ fmodel.info().show_targets(text = prefix+" start model", out = out) AttributeError: 'info' object has no attribute 'show_targets' Could somebody help me? Thanks. Jianghai PS. the version is 2007_08_18_1856. +++++++++++++++++++++++++++++++ Jianghai Zhu, Ph.D CBR Institute for Biomedical Research Department of Pathology Harvard Medical School 200 Longwood Ave., Boston, MA 02115 Ph: 617-278-3211 Fx: 618-278-3232 +++++++++++++++++++++++++++++++

2 1

What does mmtbx.twin_map_utils output?
by Francis E Reyes 11 Sep '07

11 Sep '07

Specifically I'm asking about the obs_and_calc.mtz file it produces. Thanks FR --------------------------------------------- Francis Reyes M.Sc. 215 UCB University of Colorado at Boulder gpg --keyserver pgp.mit.edu --recv-keys 67BA8D5D 8AE2 F2F4 90F7 9640 28BC 686F 78FD 6669 67BA 8D5D

2 3

So I think I have a twinning problem... ?
by Francis E Reyes 10 Sep '07

10 Sep '07

Hi all The short of the long is that phenix.xtriage on a C2 data set shows twinning with 3 twinning operators (-h,-k-l -k-h-l and k,h,-l). Pushing it through twin_utils, the soution k,h,-l gives an r/rfree of 0.30 and 0.3014 when the twin fraction is 0.352 . The question is where to go from here? I try to phenix.refine the same model in phenix.xtriage (except the mtz file is from the latest round of refinement of that model) and the R/Rfree on macro cycle =2 is 30 and 34 . Unfortunately, it doesn't seem to finish the refinement, I get a seg fault. In the def file for the refinement I set twinning = True and twin_law = k,h,-l . I tried to subject the same model and the .ref file to an autobuild passing in a .eff file for refinement specifying the twin_law, and again I get a segfault / crash when refining the input model. Please advise! --------------------------------------------- Francis Reyes M.Sc. 215 UCB University of Colorado at Boulder gpg --keyserver pgp.mit.edu --recv-keys 67BA8D5D 8AE2 F2F4 90F7 9640 28BC 686F 78FD 6669 67BA 8D5D

2 1

Twinning in phenix.autobuild?
by Francis E Reyes 10 Sep '07

10 Sep '07

Hello all So I discovered from xtriage that my crystal is twinned. I have a pretty decent model that xtriage says has an r/rfree of 30/30 to a resolution of 2.6. I wanted to try a phenix.autobuild at 2.8A. I couldn't find a parameter flag for autobuild that turns on twinning. I already autobuilt this model once, and the xtriage log did find that the crystal was twinned, but twinning in refinement of the model during rebuilding wasn't turned on. Is there a way to tell autobuild to rebuild using twinning? Thanks FR --------------------------------------------- Francis Reyes M.Sc. 215 UCB University of Colorado at Boulder gpg --keyserver pgp.mit.edu --recv-keys 67BA8D5D 8AE2 F2F4 90F7 9640 28BC 686F 78FD 6669 67BA 8D5D

3 3

crash during refinement with twinned data
by Scott White 10 Sep '07

10 Sep '07

Dear PhenixBB, I'm trying to refine a structure against potentially twinned data. I got the twin law from xtriage. phenix.refine starts, and at some point crashes with the messages below. Refinement proceeds to the end if I don't include twinning. Any suggestions? Thanks, Scott Traceback (most recent call last): File "/software/phenix-1.3b/phenix/phenix/command_line/refine.py", line 5, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/software/phenix-1.3b/phenix/phenix/refinement/command_line.py", line 81, in run call_back_after_monitor_collect=call_back_after_monitor_collect) File "/software/phenix-1.3b/phenix/phenix/refinement/driver.py", line 1325, in run log = log) File "/software/phenix-1.3b/phenix/phenix/refinement/strategies.py", line 559, in refinement_machine log = log) File "/software/phenix-1.3b/mmtbx/mmtbx/refinement/rigid_body.py", line 347, in __init__ fmodel_copy = fmodel.resolution_filter(d_min = res) File "/software/phenix-1.3b/mmtbx/mmtbx/twinning/twin_f_model.py", line 1042, in resolution_filter sf_and_grads_accuracy_params = dc.sfg_params File "/software/phenix-1.3b/mmtbx/mmtbx/twinning/twin_f_model.py", line 890, in __init__ u_part=self.scaling_parameters.u_part ) RuntimeError: cctbx Internal Error: /net/anaconda/scratch1/phenix/phenix-1.3b/cctbx/include/cctbx/xray/f_model.h(91): CCTBX_ASSERT((f_mask.size() == 0) || (hkl.size() == f_mask.size())) failure. -- _______________________________________ Dr Scott A. White Protein Crystallography Group School of Biosciences University of Birmingham B15 2TT E: s.a.white(a)bham.ac.uk T: +44-(0)121-414-7534 F: +44-(0)121-414-5925 W: http://www.bip.bham.ac.uk/cryst ______________________________________

2 1

Re: [phenixbb] Unknown atoms revisitied
by Ralf W. Grosse-Kunstleve 06 Sep '07

06 Sep '07

> Does the atom type have to be a known element or can it something like > element X as seen in some pdb files? > > An example from 1aqn: > HETATM 1951 UNK UNX 277 -4.541 50.540 38.754 0.25 13.93 X > HETATM 1952 UNK UNX 278 -4.443 49.701 37.169 1.00 15.64 X Pavel is correct, but I'll add a some more details. If you want to use these atoms in refinement, you have to decide what scattering type to use. That has to be one of the symbols published in the International Tables for Crystallography, Volume C (see cctbx/eltbx/xray_scattering/it1992.cpp under $PHENIX or cci_app_sources). Regarding the .cif file: data_comp_SIN loop_ _chem_comp_atom.comp_id _chem_comp_atom.atom_id _chem_comp_atom.type_symbol _chem_comp_atom.type_energy _chem_comp_atom.partial_charge SIN N N N 0.000 | | | | | this must be a valid "energy type" (see below) | this must be a valid scattering type symbol (Int. Tab.) this can be anything you want mon_lib/ener_lib.cif defines energy types (phenix.where_mon_lib_list_cif shows the path to mon_lib). You can define other energy types in additional .cif files you give to phenix.refine, using the same syntax as in the central ener_lib.cif. Regarding the .pdb file: HETATM 1952 UNK UNX 278 -4.443 49.701 37.169 1.00 15.64 X | | | | this must match the data_comp_UNX | | this must match _chem_comp_atom.type_symbol this must match _chem_comp_atom.atom_id You cannot use X in columns 77-78 because it must match _chem_comp_atom.type_symbol, and that in turn must be a valid scattering type symbol. However, you can leave the column blank. I think I'll change phenix.refine to treat the X like a blank, then you don't have to remove it. But in the current release you have no choice. Ralf

1 0

Re: [phenixbb] Unknown atoms revisitied
by Ralf W. Grosse-Kunstleve 06 Sep '07

06 Sep '07

Hi Jeff, > I was wondering if there were any way to have a single unknown atom > recognized by phenix for refinement. Make up any residue name you like, e.g. SIN (for single N). Create a .cif file like this: data_comp_SIN loop_ _chem_comp_atom.comp_id _chem_comp_atom.atom_id _chem_comp_atom.type_symbol _chem_comp_atom.type_energy _chem_comp_atom.partial_charge SIN N N N 0.000 In your PDB file add: ATOM 68 N SIN 15 20.000 20.000 20.000 1.00 20.00 N Add the name of the .cif file to the phenix.refine command line. You can have more than one data_comp_ in the same cif file, if you want to use different element types. Ralf

1 0

Unknown atoms in refinement
by Jeff Habel 06 Sep '07

06 Sep '07

What's the best way of dealing with unknown atoms (UNK) in phenix.refine? I have a long disordered loop that runs basically the length of the protein on the back side and there is a single "amino acid" in the middle that I'm trying to model as UNK. I can run a round of refmac and auto generate a lib file with my ligand and the unk atoms but phenix.refine balks at that with respect to the Van der Waals radii. Thanks for the help, -Jeff

4 5

wxu_scale during TLS + INDIVIDUAL_ADP refinement
by Jose M de Pereda 06 Sep '07

06 Sep '07

Dear colleagues, I am refining a structure with phenix.refine (cci_apps 2007_08_18_1856). The strategy includes refinement of tls, individual_adp, and individual_sites (during the same run.) In order to test different restraints for ADP´s refinement, I have changed the value of the wxu_scale from the default. Yet, the value of the wxu_scale parameter is always automatically reset to the same value as indicated in the log file under the section "Automatic adjustment due to TLS refinement". That is, phenix.refine uses in all cases the same wxu_scale independently of the value I set in the definitions file. Of course, when TLS refinement is not done, the wxu_scale value remains as specified. I would appreciate if someone could indicate me how to modify the ADPs restraints when doing TLS refinement. Cheers Jose -- -------------------------------------------------- Jose M de Pereda, PhD Centro de Investigación del Cáncer Spanish National Research Council - U. Salamanca Campus Unamuno s/n E-37007 Salamanca, Spain Phone: +34-923-294819 Fax: +34-923-294795 http://xtal.cicancer.org/ -------------------------------------------------

2 1

statistics for highest shell in phenix.refine
by michael hothorn 05 Sep '07

05 Sep '07

Dear all, For the crystallographic table in my manuscript, can I define a highest shell in phenix.refine and have it output the refinement statistics for this range? If not, how can I calculate the highest shell statistics using the same parameters as used for the overall refinment? Thanks Michael

2 1

Re: [phenixbb] selective nonbonded restraint removal
by Cathy Lawson 05 Sep '07

05 Sep '07

Thanks to all who replied. And thank you Ralf for your suggestion (quoted below), this worked very nicely! Somehow I missed this possibility in the online documentation. To reply to Pavel's question regarding how I figured out that tris is present at 1/2 occupancy: (1) density was first modelled with waters but they were too close together, and there was a fairly clear tetrahedral shape. (2) phosphate didn't refine well -- lead to huge negative difference peak at the P atom. (3) only other tetrahedral shape molecule in the soup was tris (TAM), and after refining just the tetrahedral core it was possible to see where the three arms extended. (4) b-factors for the tris were high relative surrounding protein atoms, and extraneous density suggested two alternate conformations for the nearby arg side-chain, one of which is too close to realistically exist next to the tris ligand. the other clue was that trial occupancy refinement of all atoms put most tris atoms in the 0.4-0.6 occupancy range. (5) tris with 0.5 occupancy, arg with two alt confs, plus a few waters at 0.5 occupancy is so far the best model with no significant different densities for this region. (6) the modelling will need to be confirmed with an SA omit ... > Hi Cathy, > > > The density suggests that when tris is not there, > > a side chain moves very close to that position > > You could model this with alternative conformations. If you assign > "altloc" A and B to two conformers of your protein, and just, say, B > to your TAM, then the ligand would only "see" your B protein. If you > don't define the B conformer (no atoms) then TAM wouldn't see that > part of the protein at all. > Actually, this approach seems better to me than working with > atom selections. > > Ralf

1 0

Re: [phenixbb] selective nonbonded restraint removal
by Cathy Lawson 05 Sep '07

05 Sep '07

Hi Ralf, That is disappointing. Lovely maps produced after a few phenix.refine runs for a 1.6 A structure revealed an unexpected tris ligand (TAM) with occupancy ~0.5. The density suggests that when tris is not there, a side chain moves very close to that position, thus my question. In other regards phenix seems a great step forward in automating refinement, and I look forward to new developments. best wishes, Cathy > Hi Cathy, > > > Is it possible to turn off nonbonded restraints for particular > > residues/monomers, or for atoms with occupancy less than 1.0? > > Nope, sorry. It is on the (long) to-do list. > BTW: each time someone asks the priority goes up. > > Ralf

4 3

selective nonbonded restraint removal
by Cathy Lawson 04 Sep '07

04 Sep '07

Hi Phenix bber's, I have been trying out phenix.refine and like it very much so far. Is it possible to turn off nonbonded restraints for particular residues/monomers, or for atoms with occupancy less than 1.0? Thanks, Cathy Lawson /////////////////////////////////////////////////////// http://rutchem.rutgers.edu/~lawson

2 1

SAD related query
by tara kashav 04 Sep '07

04 Sep '07

Dear all, I am a novice user of phenix. I am trying to obtain phases for SAD dataset collected at 2.5 angs from autosol.But when i give the phases to autobuild, the R-factor is not decreasing below 49%. Also, in warp, it is unable to built with a message "encounterd an unknown element". Can someone suggest me what could be the problem. Thanks for any suggestion in advance. Tara Kashav On 9/3/07, phenixbb-request(a)phenix-online.org < phenixbb-request(a)phenix-online.org> wrote: > > Send phenixbb mailing list submissions to > phenixbb(a)phenix-online.org > > To subscribe or unsubscribe via the World Wide Web, visit > http://www.phenix-online.org/mailman/listinfo/phenixbb > or, via email, send a message with subject or body 'help' to > phenixbb-request(a)phenix-online.org > > You can reach the person managing the list at > phenixbb-owner(a)phenix-online.org > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of phenixbb digest..." > > > Today's Topics: > > 1. command-line Patterson maps (Bryan W. Lepore) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sun, 2 Sep 2007 00:07:51 -0500 (CDT) > From: "Bryan W. Lepore" <bryanlepore(a)mail.utexas.edu> > Subject: [phenixbb] command-line Patterson maps > To: phenixbb(a)phenix-online.org > Message-ID: > < > Pine.LNX.4.64.0709020001190.16588(a)cpe-70-116-17-26.austin.res.rr.com> > Content-Type: TEXT/PLAIN; charset=US-ASCII; format=flowed > > will phenix calculate Pattersons from trial sites or a reflection file on > the command line? i.e. something besides hacking what is already there. > > e.g. cns has predict_patterson.inp. i can't seem to find a way to do it > from the command line - such as `phenix.patterson --sg=94 > reflections.mtz`. i saw phenix.maps but that looks like electron density > only. > > -bryan > > > ------------------------------ > > _______________________________________________ > phenixbb mailing list > phenixbb(a)phenix-online.org > http://www.phenix-online.org/mailman/listinfo/phenixbb > > > End of phenixbb Digest, Vol 22, Issue 1 > *************************************** > -- Tara Kashav, Dr. S. Gourinath's Lab, Lab No 430, School of Life Sciences, Jawaharlal Nehru University, New Delhi - 110067

3 2

command-line Patterson maps
by Bryan W. Lepore 03 Sep '07

03 Sep '07

will phenix calculate Pattersons from trial sites or a reflection file on the command line? i.e. something besides hacking what is already there. e.g. cns has predict_patterson.inp. i can't seem to find a way to do it from the command line - such as `phenix.patterson --sg=94 reflections.mtz`. i saw phenix.maps but that looks like electron density only. -bryan

2 1

ligand restraints
by Jeff Habel 30 Aug '07

30 Aug '07

I have a ligand with two phenyl rings where 2Fo and difference density clearly shows some torsion angle between the two rings. However, no matter how many times I tilt one ring refinement always makes them coplanar. Is there an option that I can use in phenix.refine to prevent the co-planarity or how do people typically handle this situation? Thanks for the help. All the best, -Jeff **************************************** Jeff Habel PhD Scripps Florida Drug Discovery 5353 Parkside Dr, RF-1 Rm 102 Jupiter, FL 33458 Phone: 561-799-8849 Fax: 561-799-8957 Email: jhabel(a)scripps.edu *****************************************

2 2

Turing off Auto-Modify initial isotropic ADP
by James Fraser 30 Aug '07

30 Aug '07

I'm trying to do an Fcalc with "main.number_of_macro_cycles=0 main.bulk_solvent_and_scale=false export_final_f_model=mtz ", but the b-factors get randomized before the 0th cycle of refinement (see excerpt from log below). Is there anyway to turn this off? Thanks, jf <<LOG FILE>> ----------Auto-Modify initial isotropic ADP---------- Randomizing |-ADP statistics (Wilson B = -4.823)------------------------------------------| .... -- James Fraser jamesfraser(a)berkeley.edu Alber Lab 356 Stanley Hall, QB3 UC Berkeley Berkeley, CA 94720 http://ucxray.berkeley.edu/~jfraser/

2 1

occupancy refinement question
by John Bruning 27 Aug '07

27 Aug '07

Hi, I have a ligand with dual conformations. I want to do group occupancy refinement on the conformers (each conformer is one group). I keep having trouble. The ligand is one chain named chain E and conformer 1 has residue name ADRG and conformer 2 has residue name BDRG. This is what I suppose my .def file should have in it: occupancies { individual = None group = "chain E and resname ADRG" group = "chain E and resname BDRG" one_occupancy_group_per_residue = True } When I run this it tells me I have an ambiguous selection. Any ideas why this isn't working and how I can change my .def file to make it work?

2 1

Re: [phenixbb] occupancy refinement question
by Paul Adams 27 Aug '07

27 Aug '07

Hi is it possible that the A and B in the residue name are the alternate conformation ids? If so your selection should use the DRG resname and the altloc to select the appropriate residues. Paul Adams Lawrence Berkeley Laboratory -----Original Message----- From: John Bruning <jbruning(a)gmail.com> Date: Monday, Aug 27, 2007 1:38 pm Subject: [phenixbb] occupancy refinement question To: phenixbb(a)phenix-online.orgReply-to: PHENIX user mailing list <phenixbb(a)phenix-online.org> Hi, > >I have a ligand with dual conformations. I want to do group occupancy refinement on the conformers (each conformer is one group). I keep >having trouble. The ligand is one chain named chain E and conformer 1 has residue name ADRG and conformer 2 has residue name BDRG. This is what I suppose my .def file should have in it: > >occupancies { > individual = None > group = "chain E and resname ADRG" > group = "chain E and resname BDRG" > one_occupancy_group_per_residue = True > } > >When I run this it tells me I have an ambiguous selection. > >Any ideas why this isn't working and how I can change my .def file to make it work? >_______________________________________________ >phenixbb mailing list >phenixbb(a)phenix-online.org >http://www.phenix-online.org/mailman/listinfo/phenixbb >

1 0

PX Jobs Team Leader, Postdoc, Tech (SGC-Oxford).txt
by Frank von Delft 24 Aug '07

24 Aug '07

My group has several vacancies to fill immediately: a Team Leader for Infrastructure and Methods Development; two postdoc positions; and a technician. This is the Protein Crystallography group of the Structural Genomics Consortium, Oxford. For details, please see respectively: Team Leader: http://www.sgc.ox.ac.uk/jobs/07057.html Postdocs: http://www.sgc.ox.ac.uk/jobs/07058.html Technician: http://www.sgc.ox.ac.uk/jobs/07067.html Feel free to contact me for details. The remit of the Structural Genomics Consortium (SGC) is to solve human proteins of medical relevance and place them in the public domain without restrictions; it is funded by a consortium of public and industrial funders, and has independently operating sites in Oxford and Toronto Universities and Karolinska Institutet (Stockholm). The Oxford site solved nearly 200 such structures in its first three years, and has now entered Phase II, with 4 years of funding till June 2011; it now has an increased emphasis on chemical biology and membrane proteins. The Protein Crystallography group collaborates tightly with the 5 other groups to get their purified proteins crystallized and solved -- five per month, to be precise. As importantly, with all the robotics and toys you'd care for, and access not only to stacks of real, interesting samples but also to heaps of structured data accumulated during Phase I, we're ideally positioned for development and testing of crystallographic methods, which will be the major scientific thrust of the group in this phase. Enthusiastic crystallographers, especially those that really enjoy the nuts and bolts of the technique and mucking around with tricks and toys, both crystallization and algorithms, are encouraged to apply. phx.

1 0

Help making selections
by James Fraser 23 Aug '07

23 Aug '07

I am trying to refine a structure where I want chains C and D refined anisotropically. Chain A refined isotropically and as a TLS group. Chain B refined isotropically and as a (separate) TLS group. Chain S (waters) just refined isotropically. I've tried playing around with selections on the command line and in the .def file. Yet always end up not getting the refinement to run correctly or have phenix.refine throw an "AssertionError". I'm using the latest CCIapps phenix.refine. Can someone help me out. Below are the excerpt from the DEF file and the error message thrown to the screen. Thanks, James --FROM DEF FILE-- refine { strategy = *individual_sites rigid_body *individual_adp group_adp *tls \ individual_occupancies group_occupancies group_anomalous none sites { individual = None rigid_body = None } adp { individual { isotropic = None anisotropic = "chain c or chain d" } group = None one_adp_group_per_residue = True tls = "chain a" tls = "chain b" } occupancies { individual = None group = None one_occupancy_group_per_residue = False } anomalous_s --FROM LOG FILE-- ----------Target weights (before refinement)---------- Traceback (most recent call last): File "/programs/cci_apps_2007_08_18_1856/cci_apps_sources/phenix/phenix/command_line/refine.py", line 5, in ? command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/programs/cci_apps_2007_08_18_1856/cci_apps_sources/phenix/phenix/refinement/command_line.py", line 80, in run refine_object.run( File "/programs/cci_apps_2007_08_18_1856/cci_apps_sources/phenix/phenix/refinement/driver.py", line 1244, in run log = log) File "/programs/cci_apps_2007_08_18_1856/cci_apps_sources/phenix/phenix/refinement/strategies.py", line 268, in refinement_machine log = log) File "/programs/cci_apps_2007_08_18_1856/cci_apps_sources/phenix/phenix/refinement/weight_xray_chem.py", line 137, in __init__ compute_gradients = True) File "/programs/cci_apps_2007_08_18_1856/cci_apps_sources/mmtbx/mmtbx/restraints.py", line 126, in energies_adp_aniso xray_structure = xray_structure, File "/programs/cci_apps_2007_08_18_1856/cci_apps_sources/cctbx/cctbx/adp_restraints/__init__.py", line 113, in __init__ self.check_flags(fl_i) File "/programs/cci_apps_2007_08_18_1856/cci_apps_sources/cctbx/cctbx/adp_restraints/__init__.py", line 181, in check_flags assert not fl.use_u_aniso() AssertionError -- James Fraser jamesfraser(a)berkeley.edu Alber Lab 356 Stanley Hall, QB3 UC Berkeley Berkeley, CA 94720 http://ucxray.berkeley.edu/~jfraser/

2 1

Re: [phenixbb] density modification?
by Tom Terwilliger 22 Aug '07

22 Aug '07

Hi John, Yes, you have several choices for density modification in phenix: 1. Use your experimental Fobs and phases as a starting point for the AutoBuild Wizard. After the first cycle, you will have density-modified phases. In this example below model-building is turned off with maps_only=True to make this go quickly: phenix.autobuild my_fobs_exptl_phases.mtz input_seq_file=seq.dat \ maps_only=True \ n_cycle_build_max=1 n_cycle_rebuild_max=0 Your density-modified phases will be in AutoBuild_run_1_/overall_best_denmod_map_coeffs.mtz after this run. 2. You can also use resolve directly in phenix with a command like this phenix.resolve<<EOD hklin my_fobs_exptl_phases.mtz labin FP=FP PHIB=PHIB FOM=FOM HLA=HLA HLB=HLB HLC=HLC HLD=HLD no_build solvent_content 0.45 EOD Either of these will do solvent flattening and histogram matching. If you have ncs you can specify it with a set of matrices or with a set of heavy-atom positions in a PDB file that show the NCS with: ha_file heavy_atoms_with_ncs.pdb All the best, Tom T . 03:26 PM 8/22/2007, you wrote: >Is there a way to do density modification in Phenix (solvent and or >histogram)? If so how? > >Thanks, >John Bruning >_______________________________________________ >phenixbb mailing list >phenixbb(a)phenix-online.org >http://www.phenix-online.org/mailman/listinfo/phenixbb Thomas C. Terwilliger Mail Stop M888 Los Alamos National Laboratory Los Alamos, NM 87545 Tel: 505-667-0072 email: terwilliger(a)LANL.gov Fax: 505-665-3024 SOLVE web site: http://solve.lanl.gov PHENIX web site: http:www.phenix-online.org ISFI Integrated Center for Structure and Function Innovation web site: http://techcenter.mbi.ucla.edu TB Structural Genomics Consortium web site: http://www.doe-mbi.ucla.edu/TB

1 0

density modification?
by John Bruning 22 Aug '07

22 Aug '07

Is there a way to do density modification in Phenix (solvent and or histogram)? If so how? Thanks, John Bruning

2 1

Re: [phenixbb] Question about antibumping restraints and restraints on dihedrals
by Ralf W. Grosse-Kunstleve 22 Aug '07

22 Aug '07

Hi Madhu, > 1. Is there a way of applying antibumping restraints to decrease the > number of clashes in phenix.refine, I looked at older posts but couldn't > get a satisfactory guideline regarding the same. The antibumping restraints (nonbonded interactions) are on by default. Currently there isn't a syntax to turn them off or to change the default parameterization. If you refine for a few cycles with the weight for the x-ray target set to zero (wxc_scale=0) the clashes should definitely go away. Then go back to normal refinement (wxc_scale=0.5). I'd hope that the clashes don't come back. Otherwise that's what the x-ray data guide the minimization to, even though the restraints are pushing back. > 2. The default setttings for geometry refinement looks like this: > refinement.geometry_restraints.remove { > angles = None > dihedrals = None > chiralities = None > planarities = None > } > Do I have to change "dihedrals=True" This will not work: Sorry: Invalid atom selection: refinement.geometry_restraints.remove.dihedrals="True" (RuntimeError: Syntax error.) You have to specify atom selections to *remove* default restraints. > if I do not want to put restraints on > dihedrals. I was thinking that one puts restraints on psi/phi angles only > by setting "discard_psi_phi = False" because when I set it to > "discard_psi_phi = True", the Ramachandran plot was much better. But I did > not want to do this because then the Ramachandran plot is no more a > validation tool. I just wanted to know if Unknowingly I did that by > setting up "dihedrals = None". No, if you leave the geometry_restraints.remove block untouched you'll get all restraints incl. dihedrals, except for the peptide psi and phi angles, which are controlled by the dedicated parameter you quote. Cheers, Ralf

1 0

Question about antibumping restraints and restraints on dihedrals
by Madhumati Sevvana; Biotechnik 22 Aug '07

22 Aug '07

Hi! I could successfully refine a merohedral twin in phenix.refine but I have a few unanswered questions: 1. Is there a way of applying antibumping restraints to decrease the number of clashes in phenix.refine, I looked at older posts but couldn't get a satisfactory guideline regarding the same. 2. The default setttings for geometry refinement looks like this: refinement.geometry_restraints.remove { angles = None dihedrals = None chiralities = None planarities = None } Do I have to change "dihedrals=True" if I do not want to put restraints on dihedrals. I was thinking that one puts restraints on psi/phi angles only by setting "discard_psi_phi = False" because when I set it to "discard_psi_phi = True", the Ramachandran plot was much better. But I did not want to do this because then the Ramachandran plot is no more a validation tool. I just wanted to know if Unknowingly I did that by setting up "dihedrals = None". Thank you in advance. Regards, Madhu. ------------------------------------------ Madhumati Sevvana, PhD Lehrstuhl fuer Biotechnik Friedrich-Alexander-Universitaet Henkestr. 91, D-91052 Erlangen, Germany Tel: 0049-9131-8523083 Email: msevvana(a)biologie.uni-erlangen.de -------------------------------------------

1 0

Re: [phenixbb] low res structures and Phenix
by Tom Terwilliger 21 Aug '07

21 Aug '07

Hi John, I'll answer the first of your questions... -Tom T >I have a couple of questions regarding low res structures in Phenix >(about 3.5 angstroms res). > >1. Can Phenix do B-factor sharpening for improved maps? If not >what program do people usually use for this? If you are using the AutoSol wizard, you can use the keywords: correct_aniso=True b_overall=20.0 to force correction of any anisotropy in your data and to set the overall B to 20.0. Note that this will affect all the maps produced by the AutoSol wizard. I have not yet found this to improve the quality of the resulting model; if it works for you please let me know! You can also use phenix.xtriage to do this: phenix.xtriage perfect.mtz hklout=output.mtz aniso.action=True aniso.final_b=user_b_iso aniso.b_iso=77 will set the b_iso of perfect.mtz to 77 >2. How can one optimize Ksol and Rprobe (bulk solvent model) in Phenix? > >3. Are there any other suggestions for dealing with low resolution >structures in Phenix? > >Thanks, >John Bruning >_______________________________________________ >phenixbb mailing list >phenixbb(a)phenix-online.org >http://www.phenix-online.org/mailman/listinfo/phenixbb Thomas C. Terwilliger Mail Stop M888 Los Alamos National Laboratory Los Alamos, NM 87545 Tel: 505-667-0072 email: terwilliger(a)LANL.gov Fax: 505-665-3024 SOLVE web site: http://solve.lanl.gov PHENIX web site: http:www.phenix-online.org ISFI Integrated Center for Structure and Function Innovation web site: http://techcenter.mbi.ucla.edu TB Structural Genomics Consortium web site: http://www.doe-mbi.ucla.edu/TB

3 2

cis peptide and pdb head
by Jianghai Zhu 21 Aug '07

21 Aug '07

Hi, How does phenix.refine treat cis peptide? Will it automatically detect it and restrain it like refmac5? or do I have to input the information somewhere? or generally how does phenix.refine treat the pdb head? I remember that some earlier version manual stated that all the head information are discarded, is that correct? Jianghai +++++++++++++++++++++++++++++++ Jianghai Zhu, Ph.D CBR Institute for Biomedical Research Department of Pathology Harvard Medical School 200 Longwood Ave., Boston, MA 02115 Ph: 617-278-3211 Fx: 618-278-3232 +++++++++++++++++++++++++++++++

3 3

low res structures and Phenix
by John Bruning 20 Aug '07

20 Aug '07

Hi, I have a couple of questions regarding low res structures in Phenix (about 3.5 angstroms res). 1. Can Phenix do B-factor sharpening for improved maps? If not what program do people usually use for this? 2. How can one optimize Ksol and Rprobe (bulk solvent model) in Phenix? 3. Are there any other suggestions for dealing with low resolution structures in Phenix? Thanks, John Bruning

2 1

Re: [phenixbb] autobuild wizard question
by Tom Terwilliger 19 Aug '07

19 Aug '07

Hi Jeff, I'm sorry for the trouble with the AutoBuild Wizard! I was not able to duplicate the problem with my version and putting in a model with 2 chains and trying to do rebuild_in_place of just a few residues, unfortunately. Everything seemed to work as it is supposed to. Can you check a couple things: 1. Make sure you have the latest version of PHENIX (development version) from the website. 2. Make sure that the range of residues you want to rebuild is at least 2 in from each end of the chain (so if the chain is 1-350, you can rebuild 3-6 or 347-348 but not 340-350). 3. If you accidentally give the Wizard a range that is too big, it should give you an error message and stop, but if you fix the range then it should go on. 4. Try running the Wizard from the command line, and specify "verbose=True debug=True" which will help because it gives you a lot of output that may help figure out what is going on. You will need to add also " |tee debug.log" to catch the output from this run. If none of those things help, then if you can send me the AutoBuild_run_xx_.output file from your GUI-version run and the output from running from the command line; that may help me figure out what is going wrong. Thanks! -Tom T Jeff Habel wrote: >Im a new PHENIX user and having trouble with the Autobuild wizard. > >The structure has 2 molecules in the asymmetric >unit and I am trying to build a small portion >(286-297) of both chains. I have tried building >both at the same time and trying just chain A >but always receive the error about needing 2 >amino acids on either end for anchoring the >build. Every time I alter the build range (284 >299 282 301 etc) the error persists. > >Any suggestions would be appreciated. > >Thanks, > >-Jeff Thomas C. Terwilliger Mail Stop M888 Los Alamos National Laboratory Los Alamos, NM 87545 Tel: 505-667-0072 email: terwilliger(a)LANL.gov Fax: 505-665-3024 SOLVE web site: http://solve.lanl.gov PHENIX web site: http:www.phenix-online.org ISFI Integrated Center for Structure and Function Innovation web site: http://techcenter.mbi.ucla.edu TB Structural Genomics Consortium web site: http://www.doe-mbi.ucla.edu/TB

1 0

New CCI Apps
by Ralf W. Grosse-Kunstleve 19 Aug '07

19 Aug '07

On behalf of the Phenix development team I'm pleased to announce the release of a new CCI Apps version: http://phenix-online.org/download/cci_apps/ Send bug reports to: bugs(a)phenix-online.org For help write to: help(a)phenix-online.org Version 2007_08_18_1856 phenix.reduce updated to produce PDB-V3 hydrogen names (Duke group) phenix.refine updated to handle the new atom names from the remediated PDB-V3. Currently protein and RNA/DNA only. For ligands with new names, use elbow.builder to generate CIF files with matching restraints. phenix.refine manual significantly re-worked phenix.refine bug fixes: - Proper handling of chain breaks, to avoid incorrect angle, dihedral, chirality and planarity restraints. This bug was likely to lead to bad restraints for structures with alternative conformations. - Avoid Wilson B crash for no data if low resolution omitted - Proper handling of NCS + remove_h=true - Proper handling of negative residue numbers - Avoid segfault in twin refinement - Cumulative fixes and enhancements in picking and refinement of interatomic scatterers (IAS) phenix.refine vdw_1_4_factor changed from 2/3 to 0.8 (to improve geometry and MolProbity clash scores in particular) phenix.refine obsolete remove_hydrogens option removed New phenix.pdbtools (use --help) Command aliases removed (i.e. mmtbx.xtriage is now only available as phenix.xtriage) Most major commands renamed to phenix.*, e.g. iotbx.emma -> phenix.emma iotbx.pdb.hierarchy -> phenix.pdb.hierarchy Use phenix<TAB> (in most shells) to get a full list of available phenix commands. eLBOW - Perform simple eLBOW jobs in COOT - Automatically perform the appropriate calculations on all the unknown ligands in a PDB file and combine the CIF results into one file - Covalent bonding of ligands to macromolecules (phase I implementation) - Print the sequence of PDB file (elbow.print_sequence) - User control of addition and writing of Hydrogens - Output of TRIPOS mol2 files - User control over naming of output files - elbow.join_cif_files has re-ordered arguments (target as last argument) - elbow.link_edits will generate "edits" from PDB LINK records for input to phenix.refine - elbow.metal_coordination will generate phenix.refine "edits" for metal coordination spheres including angles

1 0

autobuild wizard question
by Jeff Habel 17 Aug '07

17 Aug '07

I'm a new PHENIX user and having trouble with the Autobuild wizard. The structure has 2 molecules in the asymmetric unit and I am trying to build a small portion (286-297) of both chains. I have tried building both at the same time and trying just chain A but always receive the error about needing 2 amino acids on either end for anchoring the build. Every time I alter the build range (284... 299... 282... 301... etc) the error persists. Any suggestions would be appreciated. Thanks, -Jeff **************************************** Jeff Habel PhD Scripps Florida Drug Discovery 5353 Parkside Dr, RF-1 Rm 102 Jupiter, FL 33458 Phone: 561-799-8849 Fax: 561-799-8957 Email: jhabel(a)scripps.edu *****************************************

1 0

Re: [phenixbb] TLS refinement
by Pavel Afonine 16 Aug '07

16 Aug '07

Hi John, > I was wondering if there is an easy way to take the output from TLSmd > and input the definitions into phenix.refine. Unfortunately, I do not know since it never happened to me to use it. It is our long term plan to implement our own automatic TLS domains identification option. > If any one has an example file of a phenix TLS file with multiple TLS > definitions per chain I would love to have it (The manual just gives > an example of TLS definitions of one per chain, but I am getting > confused on how to define mutlitple). I'm loading two files at the web address below: one with TLS selections for multiple groups, another one is the PDB file with refined model and TLS records in its header. This is actually the most complex selection for TLS groups I ever used -:) http://cci.lbl.gov/~afonine/tls_selections/ Please let me know if you have any more questions. Pavel.

1 0

TLS refinement
by John Bruning 16 Aug '07

16 Aug '07

Hi, I was wondering if there is an easy way to take the output from TLSmd and input the definitions into phenix.refine. If any one has an example file of a phenix TLS file with multiple TLS definitions per chain I would love to have it (The manual just gives an example of TLS definitions of one per chain, but I am getting confused on how to define mutlitple). Thanks, John Bruning

2 1

Re: [phenixbb] Question about precedence of command line parameters vs. .def file
by Ralf W. Grosse-Kunstleve 14 Aug '07

14 Aug '07

Hi Carsten, > I see how the loading of both pdb files can be a desired feature, but I > also recall reading somewhere that command line params take precedence > over .def files, which in turn takes precedence over defaults. There > seems to be some deviation from that philosophy. Each parameter has a bunch of associated attributes to control the behavior. The attribute that causes trouble here is .multiple=True. It means that you can specify multiple values. Values accumulate in the order given. .multiple=True is used in many places and in general this feature works well, but in the case of the pdb file name it sometimes causes confusion. Unfortunately I don't see an obvious way around this. It would be possible to support something like pdb.file_name==your.pdb, where == means discard all accumulated values and start over with the new value, but that wouldn't be very intuitive. Or we'd have to stop supporting reading of multiple pdb files. I wonder if you could achieve your goals with a different approach. The idea is not to use the .def files produced by phenix.refine. I'm guessing you're doing this because they contain all your accumulated customizations defined via the command line. Could you organize your workflow to accumulate the customizations in a separate file, then pass that file as another argument? E.g.: Content of "params" file: refinement.main.macro_cycles=5 refinement.main.high_resolution=3 phenix.refine my.pdb params > This also would help future automation we are implementing, because > it would eliminate the need to read and parse an existing .def. file. Parsing and rewriting a .def file is very easy from Python. If you have phenix, you also have both Python and a parameter file parser and writer. The "phil" parser/writer is documented in this newsletter article: http://cci.lbl.gov/publications/download/iucrcompcomm_jan2005.pdf I could help you getting started if that's an interesting alternative. If you can think of a more intuitive way of handling the interaction of parameters with multiple values and command-line parameters, please let me know! Ralf

1 0

Question about precedence of command line parameters vs. .def file
by Schubert, Carsten [PRDUS] 14 Aug '07

14 Aug '07

Hi, I have a question / request about parameters defined on the command line vs. defined in the .def file are handled. Specifically the input.pdb.file_name parameter. I noticed that if I have already a input.pdb.file_name defined in the .def file and call phenix.refine with another pdb file, both files are loaded. If both pdb files represent the same protein in different stages of refinement that obviously causes the job to fail due to close non-bonding interactions. I see how the loading of both pdb files can be a desired feature, but I also recall reading somewhere that command line params take precedence over .def files, which in turn takes precedence over defaults. There seems to be some deviation from that philosophy. For me this is more a work-flow issue (which becomes important if one does to many structures per day ...) than a bug. I'd really like to minimize the time spend on editing the def file and just define new job parameters from the command line and have the def file fill in the missing bits. This also would help future automation we are implementing, because it would eliminate the need to read and parse an existing .def. file. There are ways around this but not as clean, but that is another matter .... Would it be possible to implement a switch which would check for duplicate parameters on the command line vs. .def file, give precedence to the command line and run the job like that? That way we would not break the existing scheme, which other people may have already grown accustomed too. Other suggestions are welcome too. Cheers, Carsten

2 1

conversion from h3 to r3 with phenix apps
by Green, Todd 14 Aug '07

14 Aug '07

Hello all, I have what i believe to be a rhombohedral crystal that has an insanely large cell with the hexagonal setting: 533.026 533.026 652.887 90.000 90.000 120.000 and a modestly large cell comparatively with the rhombohedral setting: 377 377 377 90 90 90 I should be able to easily reindex to the smaller cell in scalepack but for some reason i'm not getting it to work correctly. Rather than struggle further on a friday afternoon, i figured that i'd give a phenix app a try. i assume xtriage can do this, can someone point me in the correct direction? thanks in advance- Todd

3 7

refinement questions
by Sudharsan Sridharan 10 Aug '07

10 Aug '07

Hi all, I'm starting to refine a structure using phenix.refine. Almost all my model residues have been built starting from a Sharp map. I've been trying phenix.refine with and without SA and getting bond angle rmsds of 1.0 -1.1 (2.5 A data, using mlhl targets; example below) (yet to add water and also trying water picking in phenix.refine). I'm wondering if this sounds low (I have read discussions on bond and angle rmsds). More importantly I'd like to know how to work with different parameters (weights etc) during refinement in phenix. Could you please point to what params are safe to manipulate. Thanks, Sid. REMARK ******************** REFINEMENT SUMMARY: QUICK FACTS ******************* REMARK Start: r_work = 0.3637 r_free = 0.4003 bonds = 0.010 angles = 1.175 REMARK Final: r_work = 0.2607 r_free = 0.3236 bonds = 0.007 angles = 1.017 REMARK ************************************************************************

2 1

Position for PSI Knowledgebase development
by Paul Adams 07 Aug '07

07 Aug '07

======================================================================= Lawrence Berkeley National Laboratory Software Developer II ======================================================================= Software Developer II Job ID: 20930 Division: Physical Biosciences Date Opened: 8/7/2007 Summary This position will be part of the Computational Crystallography Initiative at Lawrence Berkeley Lab and will be responsible for developing a web site to disseminate information about the Technologies developed by the NIH Protein Structure Initiative (PSI). This program funds the majority of structural genomics efforts within the US and is responsible for the development of many new technologies that are useful to macromolecular crystallographers. This technology resource is part of the larger PSI Knowledge Base that will provide an information portal for the scientific community. Duties Essential: Under general direction of the Principal Investigator: Develop and maintain web site documenting PSI developed technologies. Interact with PSI researchers to collect information about new technologies. Create a relational database to store crystallographic data. Install and maintain SVN repositories for the PSI community. Marginal: Perform routine analysis of crystallographic data. Maintain web site computing hardware and software. Other duties as assigned. Qualifications Essential: Experience with writing and maintaining HTML web pages. Knowledge of standard programming methods, and Python in a UNIX environment. Strong analytical and organizational skills. Ability to handle multiple tasks in parallel. Excellent interpersonal skills for working in a diverse team environment, ability to interact effectively with other scientists, and to address their needs in an effective and timely manner. Marginal: B.S. Degree or higher in biophysics, computer science, chemistry, electrical engineering, physics or related field preferred or equivalent experience, with demonstrated programming experience. Experience maintaining computing hardware and software. Experience with analysis of crystallographic data. NOTE: This is a one year Term appointment with the possibility of renewal under the same terms and conditions, contingent upon continued funding and availability of work. Lawrence Berkeley National Laboratory is a world leader in science and engineering research, with 11 Nobel Prize recipients over the past 75 years, and 59 present members of the National Academy of Sciences. LBNL conducts unclassified research across a wide range of scientific disciplines and hosts four national user facilities. AA/ EEO employer committed to the development of a safe and diverse workforce. Learn more at http://www.lbl.gov. For more information about the CCI visit: http://cci.lbl.gov/ To apply: visit http://cjo.lbl.gov/ and search for the job number 20930. Please see http://cjo.lbl.gov/app_instr.html for details of how to apply. -- Paul Adams Senior Scientist, Physical Biosciences Division Head, Berkeley Center for Structural Biology Deputy Principal Investigator, Berkeley Structural Genomics Center Building 64, Room 248 Tel: 510-486-4225, Fax: 510-486-5909 http://cci.lbl.gov/ Lawrence Berkeley Laboratory 1 Cyclotron Road BLDG 64R0121 Berkeley, CA 94720, USA. --

1 0

ncs restraints weight
by Daniele de Sanctis 06 Aug '07

06 Aug '07

Dear all, i'm refining a low resolution structure (3.2A) with four molecules in the asymmetric unit. I saw that phenix.refine doesn't allow constraints, so i'd like to change the restraints weight but i can't find the right keyword... thanks for your help d -- Daniele de Sanctis, PhD Homo sum humani nil alienum a me puto ______________________________________________________________ phone ++ 351 21 4469 662 fax ++ 351 21 4433 664 e-mail de.sanctisd(a)gmail.com, sanctis(a)itqb.unl.pt Mailing address: ITQB, Av. República, Apartato 127 2781-901 Oeiras Portugal

2 1

Re: [phenixbb] set_biso_to_wilson_b gone?
by Pavel Afonine 05 Aug '07

05 Aug '07

Hi Luca, > I've just noticed that in the last version of PHENIX (1.3b-d7) this is > not recognised anymore: > > refinement.modify_start_model.adp.set_biso_to_wilson_b = true > > is this intentional? I'm afraid, yes. There are some technical reasons for this. You still have at least two ways of doing this: 1) When running phenix.refine : % phenix.refine data.mtz model.pdb modify_start_model.adp.set_b_iso=xxx where xxx is Wilson B values (that you may get from previous phenix.refine runs; phenix.refine reports it). 2) We recently added a tool for PDB file(s) manipulations. Currently it is very simple but will grow based on needs and user's requests. For example, it can read in a model and for selected or all atoms: - modify ADP (randomize, shift, set to a value, scale, convert to/from iso(aniso)tropic); - modify sites (shake, apply rotation / translation); - modify occupancies; - remove selected atoms; - write out a PDB file with modified model. To get started type: phenix.pdbtools -h I'm attaching a copy of a quick manual for phenix.pdbtools (it will be a part of new updated PHENIX documentation). Pavel.

3 3

Re: [phenixbb] Problematic dataset
by Peter Zwart 05 Aug '07

05 Aug '07

Hi Lucas, Your P21 just looks twinned. Have a look at the last section of the xtriage output: ------------------------------------------------------------------------------- Twinning and intensity statistics summary (acentric data): Statistics independent of twin laws - <I^2>/<I>^2 : 1.915 - <F>^2/<F^2> : 0.821 - <|E^2-1|> : 0.677 - <|L|>, <L^2>: 0.442, 0.267 Multivariate Z score L-test: 3.623 The multivariate Z score is a quality measure of the given spread in intensities. Good to reasonable data is expected to have a Z score lower than 3.5. Large values can indicate twinning, but small values do not neccesarily exclude it. Statistics depending on twin laws ------------------------------------------------------------------ | Operator | type | R obs. | Britton alpha | H alpha | ML alpha | ------------------------------------------------------------------ | h,-k,-h-l | PM | 0.090 | 0.406 | 0.417 | 0.370 | ------------------------------------------------------------------ Patterson analyses - Largest peak height : 6.379 (correpsonding p value : 6.273e-01) The largest off-origin peak in the Patterson function is 6.38% of the height of the origin peak. No significant pseudotranslation is detected. The results of the L-test indicate that the intensity statistics are significantly different then is expected from good to reasonable, untwinned data. As there are twin laws possible given the crystal symmetry, twinning could be the reason for the departure of the intensity statistics from normality. It might be worthwhile carrying out refinement with a twin specific target function. ------------------------------------------------------------------------------- You might have NCS parallel to the twin axis. If you run mmtbx.xtwin_map_utils on your P21 data, you get a mtzfile called obs_and_calc.mtz push that through xtriage (make sure to specify labels for observed and calculated data) to get the RvsR statistic. When you do that, please also provide the keyword 'perform=True' on the xtriage command line for my pleasure and send me the resulting log file. I guess the twin refinement will push the r values down a bit. Note that for futurte reference, xtriage produces a logfile as well, usually named logfile.log unless you ask it to be a different name. logfile.log contains some ccp4 style graphs you can view with loggraph (such as the L, NZ, Britton and H curves). If the above is too cryptic and you need more info, please let me know. HTH Peter Lucas Bleicher wrote: > Hi, Peter! > > I am sending the xtriage output for both space groups. > > So, the correct approach would be to use the MTZ > processed in P21 with the twin law from the xtriage > output from this MTZ file? > > Thanks in advance, > Lucas > > > Alertas do Yahoo! Mail em seu celular. Saiba mais em http://br.mobile.yahoo.com/mailalertas/ > ------------------------------------------------------------------------ > > ############################################################# > ## mmtbx.xtriage ## > ## ## > ## P.H. Zwart, R.W. Grosse-Kunstleve & P.D. Adams ## > ## ## > ############################################################# > Date 2007-07-31 Time 16:43:49 BRT -0300 > > > ##-------------------------------------------## > ## WARNING: ## > ## Number of residues unspecified ## > ##-------------------------------------------## > > > Effective parameters: > scaling.input { > parameters { > asu_contents { > n_residues = None > n_bases = None > n_copies_per_asu = None > } > misc_twin_parameters { > missing_symmetry { > tanh_location = 0.08 > tanh_slope = 50 > } > twinning_with_ncs { > perform_analyses = False > n_bins = 7 > } > twin_test_cuts { > low_resolution = 10 > high_resolution = None > isigi_cut = 3 > completeness_cut = 0.85 > } > } > reporting { > verbose = 1 > log = "logfile.log" > ccp4_style_graphs = True > } > } > xray_data { > file_name = "LamNatP1High_scala1.mtz" > obs_labels = "MEAN_lnls,SIGIMEAN_lnls" > calc_labels = None > unit_cell = 52.16379929 64.59190369 108.2256012 89.99259949 89.99089813 \ > 66.10500336 > space_group = "P 1" > high_resolution = None > low_resolution = None > } > } > <scitbx_array_family_flex_ext.double object at 0xb5e96dac> > > > Symmetry, cell and reflection file content summary > > Miller array info: LamNatP1High_scala1.mtz:IMEAN_lnls,SIGIMEAN_lnls > Observation type: xray.amplitude > Type of data: double, size=94082 > Type of sigmas: double, size=94082 > Number of Miller indices: 94082 > Anomalous flag: False > Unit cell: (52.1638, 64.5919, 108.226, 89.9926, 89.9909, 66.105) > Space group: P 1 (No. 1) > Systematic absences: 0 > Centric reflections: 0 > Resolution range: 26.0819 1.85 > Completeness in resolution range: 0.854134 > Completeness with d_max=infinity: 0.853839 > > ##----------------------------------------------------## > ## Basic statistics ## > ##----------------------------------------------------## > > > Matthews coefficient and Solvent content statistics > > Number of residues unknown, assuming 50% solvent content > > ---------------------------------------------------------------- > | Best guess : 1220 residues in the asu | > ---------------------------------------------------------------- > > > Completeness and data strength analyses > > The following table lists the completeness in various resolution > ranges, after applying a I/sigI cut. Miller indices for which > individual I/sigI values are larger than the value specified in > the top row of the table, are retained, while other intensities > are discarded. The resulting completeness profiles are an indication > of the strength of the data. > > ---------------------------------------------------------------------------------------- > | Res. Range | I/sigI>1 | I/sigI>2 | I/sigI>3 | I/sigI>5 | I/sigI>10 | I/sigI>15 | > ---------------------------------------------------------------------------------------- > | 26.08 - 4.55 | 97.0% | 96.3% | 95.5% | 93.5% | 87.3% | 80.3% | > | 4.55 - 3.62 | 96.5% | 96.1% | 95.4% | 93.8% | 87.8% | 81.2% | > | 3.62 - 3.16 | 95.7% | 94.5% | 93.2% | 89.4% | 78.1% | 66.3% | > | 3.16 - 2.87 | 93.1% | 90.7% | 87.1% | 79.1% | 60.8% | 44.5% | > | 2.87 - 2.67 | 88.4% | 84.0% | 78.5% | 66.7% | 44.0% | 27.9% | > | 2.67 - 2.51 | 85.5% | 78.6% | 71.1% | 57.7% | 32.6% | 17.9% | > | 2.51 - 2.38 | 82.4% | 74.1% | 65.3% | 49.3% | 24.8% | 12.2% | > | 2.38 - 2.28 | 81.6% | 72.3% | 62.4% | 46.0% | 20.4% | 9.0% | > | 2.28 - 2.19 | 78.5% | 67.6% | 55.1% | 37.2% | 14.8% | 5.7% | > | 2.19 - 2.12 | 75.5% | 61.9% | 49.0% | 30.7% | 9.9% | 3.5% | > | 2.12 - 2.05 | 72.1% | 55.6% | 41.8% | 23.7% | 6.3% | 1.9% | > | 2.05 - 1.99 | 68.6% | 50.2% | 35.8% | 17.9% | 4.1% | 1.0% | > | 1.99 - 1.94 | 63.2% | 42.9% | 28.5% | 13.7% | 3.1% | 0.5% | > | 1.94 - 1.89 | 56.6% | 36.7% | 23.2% | 10.4% | 1.7% | 0.3% | > ---------------------------------------------------------------------------------------- > > The completeness of data for which I/sig(I)>3.00, exceeds 85% for > for resolution ranges lower than 2.87A. > The data is cut at this resolution for the potential twin tests > and intensity statistics. > > > Maximum likelihood isotropic Wilson scaling > ML estimate of overall B value of LamNatP1High_scala1.mtz:IMEAN_lnls,SIGIMEAN_lnls: > 15.59 A**(-2) > Estimated -log of scale factor of LamNatP1High_scala1.mtz:IMEAN_lnls,SIGIMEAN_lnls: > -2.45 > > > Maximum likelihood anisotropic Wilson scaling > ML estimate of overall B_cart value of LamNatP1High_scala1.mtz:IMEAN_lnls,SIGIMEAN_lnls: > 15.61, 1.39, -1.33 > 11.40, 0.27 > 21.10 > Equivalent representation as U_cif: > 0.18, -0.04, -0.02 > 0.14, 0.00 > 0.27 > > ML estimate of -log of scale factor of LamNatP1High_scala1.mtz:IMEAN_lnls,SIGIMEAN_lnls: > -2.45 > Correcting for anisotropy in the data > > Some basic intensity statistics follow. > > > > Low resolution completeness analyses > > The following table shows the completeness > of the data to 5 Angstrom. > unused: - 26.0822 [ 0/38 ] 0.000 > bin 1: 26.0822 - 10.5555 [528/541] 0.976 > bin 2: 10.5555 - 8.4726 [551/564] 0.977 > bin 3: 8.4726 - 7.4299 [530/547] 0.969 > bin 4: 7.4299 - 6.7636 [568/582] 0.976 > bin 5: 6.7636 - 6.2860 [525/540] 0.972 > bin 6: 6.2860 - 5.9200 [547/559] 0.979 > bin 7: 5.9200 - 5.6266 [530/546] 0.971 > bin 8: 5.6266 - 5.3839 [519/532] 0.976 > bin 9: 5.3839 - 5.1783 [568/584] 0.973 > bin 10: 5.1783 - 5.0009 [516/532] 0.970 > unused: 5.0009 - [ 0/0 ] > > > > Mean intensity analyses > Analyses of the mean intensity. > Inspired by: Morris et al. (2004). J. Synch. Rad.11, 56-59. > The following resolution shells are worrisome: > ------------------------------------------------ > | d_spacing | z_score | compl. | <Iobs>/<Iexp> | > ------------------------------------------------ > | 9.993 | 7.15 | 0.97 | 0.426 | > | 8.448 | 7.29 | 0.98 | 0.559 | > | 7.451 | 6.09 | 0.97 | 0.676 | > | 4.303 | 4.64 | 0.97 | 1.223 | > | 3.410 | 4.52 | 0.96 | 1.176 | > | 3.333 | 6.63 | 0.96 | 1.262 | > | 3.261 | 5.36 | 0.96 | 1.205 | > | 3.194 | 5.48 | 0.96 | 1.203 | > | 3.071 | 4.87 | 0.94 | 1.183 | > | 2.962 | 6.43 | 0.94 | 1.245 | > | 2.911 | 4.61 | 0.93 | 1.160 | > | 2.454 | 4.96 | 0.87 | 0.872 | > | 2.425 | 6.89 | 0.87 | 0.825 | > | 2.397 | 4.68 | 0.87 | 0.873 | > | 2.182 | 6.79 | 0.84 | 0.830 | > | 2.142 | 5.54 | 0.84 | 0.856 | > | 1.870 | 5.46 | 0.60 | 1.205 | > | 1.857 | 11.65 | 0.54 | 1.570 | > ------------------------------------------------ > > Possible reasons for the presence of the reported > unexpected low or elevated mean intensity in > a given resolution bin are : > - missing overloaded or weak reflections > - suboptimal data processing > - satelite (ice) crystals > - NCS > - translational pseudo symmetry (detected elsewhere) > - outliers (detected elsewhere) > - ice rings (detected elsewhere) > - other problems > Note that the presence of abnormalities > in a certain region of reciprocal space might > confuse the data validation algorithm throughout > a large region of reciprocal space, even though > the data is acceptable in those areas. > > > Possible outliers > Inspired by: Read, Acta Cryst. (1999). D55, 1759-1764 > > Acentric reflections: > > ----------------------------------------------------------------- > | d_space | H K L | |E| | p(wilson) | p(extreme) | > ----------------------------------------------------------------- > | 2.114 | -18, -9, 35 | 3.79 | 5.57e-07 | 5.08e-02 | > | 3.066 | -16, -8, 12 | 3.93 | 1.90e-07 | 1.77e-02 | > | 3.066 | 16, 8, 12 | 3.95 | 1.69e-07 | 1.57e-02 | > | 2.115 | 18, 9, 35 | 3.79 | 5.65e-07 | 5.15e-02 | > ----------------------------------------------------------------- > > p(wilson) : 1-(1-exp[-|E|^2]) > p(extreme) : 1-(1-exp[-|E|^2])^(n_acentrics) > p(wilson) is the probability that an E-value of the specified > value would be observed when it would selected at random from > the given data set. > p(extreme) is the probability that the largest |E| value is > larger or equal then the observed largest |E| value. > > Both measures can be used for outlier detection. p(extreme) > takes into account the size of the dataset. > > > Centric reflections: > > None > > > Ice ring related problems > > The following statistics were obtained from ice-ring > insensitive resolution ranges > mean bin z_score : 3.46 > ( rms deviation : 1.89 ) > mean bin completeness : 0.89 > ( rms deviation : 0.09 ) > > The following table shows the z-scores > and completeness in ice-ring sensitive areas. > Large z-scores and high completeness in these > resolution ranges might be a reason to re-assess > your data processsing if ice rings were present. > > ------------------------------------------------ > | d_spacing | z_score | compl. | Rel. Ice int. | > ------------------------------------------------ > | 3.897 | 0.51 | 0.97 | 1.000 | > | 3.669 | 1.87 | 0.97 | 0.750 | > | 3.441 | 4.52 | 0.96 | 0.530 | > | 2.671 | 3.84 | 0.89 | 0.170 | > | 2.249 | 3.61 | 0.87 | 0.390 | > | 2.072 | 3.23 | 0.82 | 0.300 | > | 1.948 | 0.62 | 0.76 | 0.040 | > | 1.918 | 0.44 | 0.73 | 0.180 | > | 1.883 | 5.46 | 0.60 | 0.030 | > ------------------------------------------------ > > Abnormalities in mean intensity or completeness at > resolution ranges with a relative ice ring intensity > lower then 0.10 will be ignored. > > No ice ring related problems detected. > If ice rings were present, the data does not look > worse at ice ring related d_spacings as compared > to the rest of the data set > > > > > Basic analyses completed > > ##----------------------------------------------------## > ## Twinning Analyses ## > ##----------------------------------------------------## > > > > Using data between 10.00 to 2.87 Angstrom. > > Determining possible twin laws. > > The following twin laws have been found: > > --------------------------------------------------------------------------------- > | Type | Axis | R metric (%) | delta (le Page) | delta (Lebedev) | Twin law | > --------------------------------------------------------------------------------- > | PM | 2-fold | 0.009 | 0.010 | 0.000 | -h,-k,l | > | PM | 2-fold | 0.090 | 0.079 | 0.001 | -h,-h+k,-l | > | PM | 2-fold | 0.099 | 0.080 | 0.001 | h,h-k,-l | > --------------------------------------------------------------------------------- > M: Merohedral twin law > PM: Pseudomerohedral twin law > > 0 merohedral twin operators found > 3 pseudo-merohedral twin operators found > In total, 3 twin operator were found > > > Number of centrics : 0 > Number of acentrics : 27599 > > Largest patterson peak with length larger then 15 Angstrom > > Frac. coord. : 0.119 -0.276 -0.000 > Distance to origin : 16.337 > Height (origin=100) : 5.893 > p_value(height) : 7.365e-01 > > The reported p_value has the following meaning: > The probability that a peak of the specified height > or larger is found in a Patterson function of a > macro molecule that does not have any translational > pseudo symmetry is equal to 7.365e-01 > p_values smaller then 0.05 might indicate > weak translation pseudo symmetry, or the self vector of > a large anomalous scatterer such as Hg, whereas values > smaller then 1e-3 are a very strong indication for > the presence of translational pseudo symmetry. > > > > Wilson ratio and moments > > Acentric reflections > <I^2>/<I>^2 :1.932 (untwinned: 2.000; perfect twin 1.500) > <F>^2/<F^2> :0.820 (untwinned: 0.785; perfect twin 0.885) > <|E^2 - 1|> :0.678 (untwinned: 0.736; perfect twin 0.541) > > > > NZ test (0<=z<1) to detect twinning and possible translational NCS > > > ----------------------------------------------- > | Z | Nac_obs | Nac_theo | Nc_obs | Nc_theo | > ----------------------------------------------- > | 0.0 | 0.000 | 0.000 | 0.000 | 0.000 | > | 0.1 | 0.046 | 0.095 | 0.000 | 0.248 | > | 0.2 | 0.127 | 0.181 | 0.000 | 0.345 | > | 0.3 | 0.210 | 0.259 | 0.000 | 0.419 | > | 0.4 | 0.287 | 0.330 | 0.000 | 0.474 | > | 0.5 | 0.361 | 0.394 | 0.000 | 0.520 | > | 0.6 | 0.427 | 0.451 | 0.000 | 0.561 | > | 0.7 | 0.486 | 0.503 | 0.000 | 0.597 | > | 0.8 | 0.544 | 0.551 | 0.000 | 0.629 | > | 0.9 | 0.597 | 0.593 | 0.000 | 0.657 | > | 1.0 | 0.643 | 0.632 | 0.000 | 0.683 | > ----------------------------------------------- > | Maximum deviation acentric : 0.055 | > | Maximum deviation centric : 0.683 | > | | > | <NZ(obs)-NZ(twinned)>_acentric : -0.024 | > | <NZ(obs)-NZ(twinned)>_centric : -0.467 | > ----------------------------------------------- > > > L test for acentric data > > using difference vectors (dh,dk,dl) of the form: > (2hp,2kp,2lp) > where hp, kp, and lp are random signed integers such that > 2 <= |dh| + |dk| + |dl| <= 8 > > Mean |L| :0.441 (untwinned: 0.500; perfect twin: 0.375) > Mean L^2 :0.267 (untwinned: 0.333; perfect twin: 0.200) > > The distribution of |L| values indicates a twin fraction of > 0.08. Note that this estimate is not as reliable as obtained > via a Britton plot or H-test if twin laws are available. > > > > --------------------------------------------- > Analysing possible twin law : -h,-k,l > --------------------------------------------- > > > Results of the H-test on a-centric data: > > (Only 50.0% of the strongest twin pairs were used) > > mean |H| : 0.053 (0.50: untwinned; 0.0: 50% twinned) > mean H^2 : 0.005 (0.33: untwinned; 0.0: 50% twinned) > Estimation of twin fraction via mean |H|: 0.447 > Estimation of twin fraction via cum. dist. of H: 0.449 > > > > Britton analyses > > Extrapolation performed on 0.47 < alpha < 0.495 > Estimated twin fraction: 0.442 > Correlation: 0.9953 > > R vs R statistic: > R_abs_twin = <|I1-I2|>/<|I1+I2|> > Lebedev, Vagin, Murshudov. Acta Cryst. (2006). D62, 83-95 > > R_abs_twin observed data : 0.054 > > R_sq_twin = <(I1-I2)^2>/<(I1+I2)^2> > R_sq_twin observed data : 0.004 > No calculated data available. > R_twin for calculated data not determined. > > > Maximum Likelihood twin fraction determination > Zwart, Read, Grosse-Kunstleve & Adams, to be published. > > > The estimated twin fraction is equal to 0.415 > > > --------------------------------------------- > Analysing possible twin law : -h,-h+k,-l > --------------------------------------------- > > > Results of the H-test on a-centric data: > > (Only 50.0% of the strongest twin pairs were used) > > mean |H| : 0.093 (0.50: untwinned; 0.0: 50% twinned) > mean H^2 : 0.015 (0.33: untwinned; 0.0: 50% twinned) > Estimation of twin fraction via mean |H|: 0.407 > Estimation of twin fraction via cum. dist. of H: 0.407 > > > > Britton analyses > > Extrapolation performed on 0.44 < alpha < 0.495 > Estimated twin fraction: 0.399 > Correlation: 0.9959 > > R vs R statistic: > R_abs_twin = <|I1-I2|>/<|I1+I2|> > Lebedev, Vagin, Murshudov. Acta Cryst. (2006). D62, 83-95 > > R_abs_twin observed data : 0.093 > > R_sq_twin = <(I1-I2)^2>/<(I1+I2)^2> > R_sq_twin observed data : 0.011 > No calculated data available. > R_twin for calculated data not determined. > > > Maximum Likelihood twin fraction determination > Zwart, Read, Grosse-Kunstleve & Adams, to be published. > > > The estimated twin fraction is equal to 0.361 > > > --------------------------------------------- > Analysing possible twin law : h,h-k,-l > --------------------------------------------- > > > Results of the H-test on a-centric data: > > (Only 50.0% of the strongest twin pairs were used) > > mean |H| : 0.085 (0.50: untwinned; 0.0: 50% twinned) > mean H^2 : 0.013 (0.33: untwinned; 0.0: 50% twinned) > Estimation of twin fraction via mean |H|: 0.415 > Estimation of twin fraction via cum. dist. of H: 0.417 > > > > Britton analyses > > Extrapolation performed on 0.46 < alpha < 0.495 > Estimated twin fraction: 0.411 > Correlation: 0.9961 > > R vs R statistic: > R_abs_twin = <|I1-I2|>/<|I1+I2|> > Lebedev, Vagin, Murshudov. Acta Cryst. (2006). D62, 83-95 > > R_abs_twin observed data : 0.086 > > R_sq_twin = <(I1-I2)^2>/<(I1+I2)^2> > R_sq_twin observed data : 0.009 > No calculated data available. > R_twin for calculated data not determined. > > > Maximum Likelihood twin fraction determination > Zwart, Read, Grosse-Kunstleve & Adams, to be published. > > > The estimated twin fraction is equal to 0.350 > > > > Exploring higher metric symmetry > > The point group of data as dictated by the space group is P 1 > the point group in the niggli setting is P 1 > The point group of the lattice is Hall: C 2 2 (x+y,2*y,z) > A summary of R values for various possible point groups follow. > > ----------------------------------------------------------------------------------------------- > | Point group | mean R_used | max R_used | mean R_unused | min R_unused | choice | > ----------------------------------------------------------------------------------------------- > | Hall: C 2y (x-y,2*x,z) | 0.086 | 0.086 | 0.054 | 0.054 | | > | P 1 | None | None | 0.070 | 0.054 | | > | Hall: C 2y (x+y,2*y,z) | 0.093 | 0.093 | 0.054 | 0.054 | | > | Hall: C 2 2 (x+y,2*y,z) | 0.070 | 0.086 | None | None | | > | P 1 1 2 | 0.054 | 0.054 | 0.093 | 0.093 | <--- | > ----------------------------------------------------------------------------------------------- > > R_used: mean and maximum R value for symmetry operators *used* in this point group > R_unused: mean and minimum R value for symmetry operators *not used* in this point group > The likely point group of the data is: P 1 1 2 > > Possible space groups in this point groups are: > Unit cell: (52.1638, 108.226, 64.5258, 90, 113.762, 90) > Space group: P 1 2 1 (No. 3) > > Unit cell: (52.1638, 108.226, 64.5258, 90, 113.762, 90) > Space group: P 1 21 1 (No. 4) > > Note that this analyses does not take into account the effects of twinning. > If the data is (allmost) perfectly twinned, the symmetry will appear to be > higher than it actually is. > > > > > > ------------------------------------------------------------------------------- > Twinning and intensity statistics summary (acentric data): > > Statistics independent of twin laws > - <I^2>/<I>^2 : 1.932 > - <F>^2/<F^2> : 0.820 > - <|E^2-1|> : 0.678 > - <|L|>, <L^2>: 0.441, 0.267 > Multivariate Z score L-test: 3.807 > The multivariate Z score is a quality measure of the given > spread in intensities. Good to reasonable data is expected > to have a Z score lower than 3.5. > Large values can indicate twinning, but small values do not > neccesarily exclude it. > > > Statistics depending on twin laws > ------------------------------------------------------------------- > | Operator | type | R obs. | Britton alpha | H alpha | ML alpha | > ------------------------------------------------------------------- > | -h,-k,l | PM | 0.054 | 0.442 | 0.449 | 0.415 | > | -h,-h+k,-l | PM | 0.093 | 0.399 | 0.407 | 0.361 | > | h,h-k,-l | PM | 0.086 | 0.411 | 0.417 | 0.350 | > ------------------------------------------------------------------- > > Patterson analyses > - Largest peak height : 5.893 > (correpsonding p value : 7.365e-01) > > > The largest off-origin peak in the Patterson function is 5.89% of the > height of the origin peak. No significant pseudotranslation is detected. > > The results of the L-test indicate that the intensity statistics > are significantly different then is expected from good to reasonable, > untwinned data. > As there are twin laws possible given the crystal symmetry, twinning could > be the reason for the departure of the intensity statistics from normality. > It might be worthwhile carrying refinement with a twin specific target function. > > Note that the symmetry of the intensities suggest that the assumed space group > is too low. As twinning is however suspected, it is not immediuatly clear if this > is the case. Carefull reprocessing and (twin)refinement for all cases might resolve > this question. > > -------------------------------------------------------------------------------

1 0

R, Rfree, and wxu_scale
by Jianghai Zhu 05 Aug '07

05 Aug '07

Hi, Just finished a round of refinement. Rfree went up about 1% and R went down a bit. I am wondering what happened here. Also Rfree always went down after TLS refinement but always went up a lot more after the individual ADP refinement. Is there anything I can do here? Should I use weaker B restraints? I tried to assign a wxu_scale value. But phenix.refine automatically adjust it to 1.81 every time. ------------------------------------------------------------------------ R-factors, x-ray target values and norm of gradient of x-ray target stage r-work r-free xray_target_w xray_target_t 0 : 0.3982 0.5701 6.156859e+00 6.238509e+00 1_bss: 0.2008 0.2494 5.949117e+00 6.088855e+00 1_xyz: 0.2030 0.2497 5.953829e+00 6.091696e+00 1_adp: 0.1998 0.2559 5.952833e+00 6.099340e+00 2_bss: 0.1996 0.2557 5.952517e+00 6.099196e+00 2_xyz: 0.1984 0.2576 5.953246e+00 6.104085e+00 2_adp: 0.1974 0.2578 5.950434e+00 6.102336e+00 3_bss: 0.1972 0.2574 5.950333e+00 6.102222e+00 3_xyz: 0.1997 0.2576 5.956524e+00 6.103647e+00 3_adp: 0.1987 0.2590 5.955554e+00 6.104434e+00 3_bss: 0.1986 0.2588 5.955445e+00 6.104227e+00 ------------------------------------------------------------------------ Weights for target T = Exray * wxc * wxc_scale + Echem * wc and angles between gradient vectors, eg. (d_Exray/d_sites, d_Echem/d_sites) stage wxc wxu wxc_sc wxu_sc /_gxc,gc /_gxu,gu 0 : 5.5611e-01 7.2249e-01 0.050 1.810 96.334 101.736 1_bss: 5.5611e-01 7.2249e-01 0.050 1.810 96.334 101.736 1_xyz: 4.5001e-01 6.1023e-01 0.050 1.810 96.591 102.016 1_adp: 4.5001e-01 6.1023e-01 0.050 1.810 96.591 102.016 2_bss: 4.5001e-01 6.1023e-01 0.050 1.810 96.591 102.016 2_xyz: 5.1326e-01 6.5356e-01 0.050 1.810 112.802 141.393 2_adp: 5.1326e-01 6.5356e-01 0.050 1.810 112.802 141.393 3_bss: 5.1326e-01 6.5356e-01 0.050 1.810 112.802 141.393 3_xyz: 3.3568e-01 7.0695e-01 0.050 1.810 111.975 135.760 3_adp: 3.3568e-01 7.0695e-01 0.050 1.810 111.975 135.760 3_bss: 3.3568e-01 7.0695e-01 0.050 1.810 111.975 135.760 ------------------------------------------------------------------------ Jianghai +++++++++++++++++++++++++++++++ Jianghai Zhu, Ph.D CBR Institute for Biomedical Research Department of Pathology Harvard Medical School 200 Longwood Ave., Boston, MA 02115 Ph: 617-278-3211 Fx: 618-278-3232 +++++++++++++++++++++++++++++++

3 9

Postdoctoral Position for PHENIX development
by Paul D. Adams 02 Aug '07

02 Aug '07

======================================================================= Lawrence Berkeley National Laboratory Computational Postdoctoral Fellow ======================================================================= Computational Postdoctoral Fellow Job ID: 20907 Division: Physical Biosciences Date Opened: 8/2/2007 Summary We are seeking an intelligent, motivated scientist to develop graphical tools for the PHENIX software system. The successful candidate will work under the direction of the Principal Investigator for this project. Primary duties and responsibilities will include development and implementation of graphical interfaces, user workflow tools, and mechanisms for data tracking and storage. Duties Essential - Primary duties and responsibilities will include development and implementation of computational algorithms for structure determination using X-ray crystallography. Develop graphic interfaces to components of the PHENIX software system. Create interfaces between PHENIX and other software packages such as Coot. Implement software code for capturing and tracking data and other information as it is created by the PHENIX software. This position requires interaction with P.I.s, postdocs, Scientists and external collaborators involved in the PHENIX project. Publication of research results in scientific journals and presentation of results at scientific meetings. Marginal - Actively participate in group meetings. Assume shared responsibility for maintenance of computational equipment with other group members. Qualifications Essential - Ph.D. within the last three years in computational crystallography or a related field. Demonstrated experience with graphical interface development. First-author publications in macromolecular crystallography or closely related field. Familiarity with software development tools such as C++, C or Fortran programming languages and scripting languages such as Python or Perl. Strong programming skills. Solid interpersonal skills and the ability to work in a team environment are critical. Marginal - Experience with the practical aspects of protein crystallography: protein production, crystallization and structure solution. NOTE: This is a one year Term appointment with the possibility of renewal under the same terms and conditions, contingent upon continued funding and availability of work. Lawrence Berkeley National Laboratory is a world leader in science and engineering research, with 11 Nobel Prize recipients over the past 75 years, and 59 present members of the National Academy of Sciences. LBNL conducts unclassified research across a wide range of scientific disciplines and hosts four national user facilities. AA/ EEO employer committed to the development of a safe and diverse workforce. Learn more at http://www.lbl.gov. For more information about PHENIX visit: http://www.phenix-online.org/ For more information about the CCI visit: http://cci.lbl.gov/ To apply: visit http://cjo.lbl.gov/ and search for the job number 20907. Please see http://cjo.lbl.gov/app_instr.html for details of how to apply. -- Paul Adams Senior Scientist, Physical Biosciences Division Head, Berkeley Center for Structural Biology Deputy Principal Investigator, Berkeley Structural Genomics Center Building 64, Room 248 Tel: 510-486-4225, Fax: 510-486-5909 http://cci.lbl.gov/ Lawrence Berkeley Laboratory 1 Cyclotron Road BLDG 64R0121 Berkeley, CA 94720, USA. --

1 0