February 2010 - phenixbb - phenix-online.org

riding model and deposition
by Alex Theodossis 19 Feb '10

19 Feb '10

Dear phenixbb, I have refined a 1.3A structure using individual anisotropic B-factors (Rfactor, 15.5%; Rfree, 17.7%). Additional use of the riding model during refinement has improvement the R values (Rf, 14.5%; Rfree, 16.6%), although bond length and angle outliers did increase (anyone know why that would be?). I have decided to deposit the model obtained using the riding model, but am unsure whether to include the hydrogen coordinates in the deposition. The Phenix.refine documentation advises that they are left in, as they are needed to reproduce the refinement statistics. However, reviewing the high resolution structures in the PDB I found that: 1. In the case of structures refined using REFMAC or SHELX, the riding model does not seem to be included, although its use during refinement is mentioned in the pdb header. 2. In the case of structures refined using PHENIX, use of the riding model is sometimes indicated in the header, but hydrogen coordinates only appear to be included in the pdb in the case of structures refined to a resolution better than 1.3A. As far as I can understand the purpose of the riding model is to provide additional restraints for refinement, so as long as its use is made evident in the pdb header (and maybe the methods section of the publication) that should be enough. Inclusion of hydrogens in the model, when you can't see them seems inappropriate. Nevertheless, I would appreciate a more authoritative perspective on this issue. Thank you in advance, Alex Theodossis

5 7

refinement error with ligand
by Christian Roth 18 Feb '10

18 Feb '10

Dear all, I started a refinement of one of my structures using the GUI of phenix (Version 1.6). I used the options of adding hydrogens and dis a suimulated annealing. I used the output file and corrected it manually using Coot. I incorporated a few Ethylenglycol Monomers which might have bound. I merged everything wrote the file out and loaded this pdb in the GUI. I got an error which might be du to the fact that I didnt specify a CIF file for the ligand. I started ready set from the GUI to check the pdb and used REEL to create a CIF File. I added the CIF-File to the File list and tried to start again phenix.refine using the hydrogen option but no sim. annealing. However I got always the following error log: # Date 2010-02-17 Time 12:00:35 CET +0100 (1266404435.31 s) #phil __OFF__ Command line arguments: "/homes/roth/clcf/049b32mlsyncref/phenix/phenix/project_data/refine_15.eff" "-- overwrite" HOST = bbzws339 HOSTTYPE = i486-linux VENDOR = intel USER = roth PID = 29983 ------------------------------------------------------------------------------- PHENIX: Python-based Hierarchical ENvironment for Integrated Xtallography Version: 1.6 Release tag: 289 Platform: intel-linux-2.6 linux User: roth ------------------------------------------------------------------------------- phenix.refine: Macromolecular Structure Refinement Monomer Library directory: "/software/linuxsoft/phenix-1.6-289/chem_data/mon_lib" Total number of atoms: 16881 Number of models: 1 Model: "" Number of chains: 19 Chain: "A" Number of atoms: 1575 Number of conformers: 1 Conformer: "" Number of residues, atoms: 94, 1575 Classifications: {'peptide': 94} Modifications used: {'COO': 1, 'NH3': 1} Link IDs: {'PTRANS': 6, 'TRANS': 87} and so on for other chains Chain: "R" Number of atoms: 7 Number of conformers: 1 Conformer: "" Number of residues, atoms: 1, 7 Unusual residues: {'EDO': 1} Unexpected atoms: {'EDO,HO1': 1, 'EDO,HO2': 1} Classifications: {'undetermined': 1} Number of atoms with unknown nonbonded energy type symbols: 16 "HETATM16821 HO2 EDO K 1 .*. H " "HETATM16826 HO1 EDO K 1 .*. H " "HETATM16828 HO2 EDO L 1 .*. H " "HETATM16834 HO1 EDO L 1 .*. H " "HETATM16836 HO2 EDO M 1 .*. H " "HETATM16840 HO1 EDO M 1 .*. H " "HETATM16842 HO2 EDO N 1 .*. H " "HETATM16848 HO1 EDO N 1 .*. H " "HETATM16850 HO2 EDO O 1 .*. H " "HETATM16857 HO1 EDO O 1 .*. H " ... (remaining 6 not shown) I looked in the pdb for this atoms, but I cannot find them. When I use a cif file created directly by phenix elbow using command line and when I did not use the update hydrogen button the refinement works. But the new residues as well as the EDO molecules have no hydrogens in riding positions. Do I have to remove all hydrogens and attach them again for every cycle? One Refinement runs for more then 5 hours regardless I use sim annealing or not. This seems to be very long for me. Maybe I want too many optimization steps? Could anyone explain me what goes wrong here. I am no so much experienced with phenix and now I run out of ideas. Thanks in advance Christian

4 4

phenix.elbow and compound refinement
by YS. Chung 17 Feb '10

17 Feb '10

Dear all, I am having some difficulties in refining my compound within my structure. I was able to create the compound and the cif file using PRODRG and was able to fit it in coot perfectly fine. However, as pointed out by many users of phenix.refine, the cif file from PRODRG does not work well with phenix.refine. I have created a cif file using phenix.elbow and the refinement itself works fine. However, when I try to use the cif file from phenix elbow in coot, it pulls apart my compound and it's extremely difficult to work with. This still happens when I create the cif file using the --final_geometry option. Is there a way around this or is there something I can do? Any help will be much appreciated. Thank you, Yu Seon Chung

2 1

refinement at 4.4 A => side chains or not
by jonathan elegheert 13 Feb '10

13 Feb '10

Dear bb, could anyone please share his/her opinion on following matter; I am refining a 4.4 A protein complex structure which I could solve using MR with the high resolution models of the separate components. I also have another 3.5 A crystal form of the complex. My question is whether or not to rebuild/leave "positionally known" side chains in the 4.4 A structure, because of course 2/3 of the time there is of course no density for it. However, a review of the PDB learnes that almost all models deposited between 4 and 5 A contain full side chains! Personally my opinion is not to include them because a model should be a representation of the information content of the experiment and thus density, but maybe that deviates from the consensus. On a more technical note: until now the trimming procedure has only sporadically led to Fo-Fc 'blobs' in the hydrophobic protein interior; are there proven modifications one has to apply to the bulk solvent model mask to avoid these even more? Thanks a lot for your insights, Jonathan -- Jonathan Elegheert Ph.D. Student Unit for Structural Biology& Biophysics http://www.lprobe.ugent.be/xray.html Lab for Protein Biochemistry and Biomolecular Engineering Department of Biochemistry& Microbiology Ghent University, Belgium e-mail: jonathan.elegheert(a)ugent.be

3 2

multiple model refinement
by chene＠umbi.umd.edu 13 Feb '10

13 Feb '10

Hi, PHENIX users, Does anyone know if multiple model refinement is implemented in the latest version? I am using version 1.3. It was not implemented yet in that version. Thanks Chen Chen Chen, PhD, Research Associate Center for Advanced Research in Biotechnology (CARB), University of Maryland Biotechnology Institute (UMBI), 9600 Gudelsky Drive, Rockville, MD 20850 240-314-6122

2 1

iMosflm .mtz
by anon 12 Feb '10

12 Feb '10

Hi, i was under the impression Phenix (AutoSol) could handle .mtz files from iMosflm, however when i input the .mtz file into AutoSol, no labels appears in the pull-down menu, and pressing Run results in a dialogue box appearing stating you must select a column label for the reflection file, if no labels are present it may not be a valid format. If i input the iMosflm .mtz file into phenix.reflection_file_editor, and copy over the I, SIGI, M_ISYM columns, i get the following error; A Python error was detected. This is probably a bug; click "OK" to send a bug report to the PHENIX developers. TypeError : not all arguments converted during string formatting Traceback (most recent call last): File "/usr/local/phenix-1.6-289/phenix/wxGUI2/MainWindow.py", line 880, in OnRun self.run() File "/usr/local/phenix-1.6-289/phenix/wxGUI2/Programs/ HKLTools.py", line 81, in run symmetry_callback=None) # TODO File "/usr/local/phenix-1.6-289/cctbx_project/iotbx/ reflection_file_editor.py", line 358, in __init__ "related reflections are not identical.)" % array_name) TypeError: not all arguments converted during string formatting When i use the scaled .mtz file produced by scala, AutoSol works fine. Just wondering what's going on! Also, how do you delete individual Jobs in the Job History window (ie jobs that have failed or you no longer want to keep etc) Thanks for the help

3 2

phenix.refine frozen
by Roberto Steiner 12 Feb '10

12 Feb '10

Hi there, I would appreciate some feedback on the following: I am playing around with phenix.refine which comes with Phenix distro 1.6-289 for Mac-Intel. When running a normal refinement job (ind site + grouped ADPs + ncs+TorsionSimulatedAnnealing) phenix.refine gets stuck after bulk solvent modelling and scaling. The last things I see in the log page are: ====================== bulk solvent modeling and scaling ====================== |--(resolution: 3.50 - 58.55 A; n_refl. = 55140)------------------------------| | | | r_work= 0.2687 r_free= 0.2800 ksol= 0.38 Bsol= 10.09 scale= 1.148 | | | | overall anisotropic scale matrix (Cartesian basis; B11,B22,B33,B12,B13,B23):| | (-2.67,-2.67,5.34,0.00,0.00,-0.00); trace/3= 0.00 | | | | maximum likelihood estimate for coordinate error: 0.41 A | | x-ray target function (ml) for work reflections: 7.216427 | | -----------------------------------------------------------------------------| ========================== Target weights: x-ray data ========================= and that's it. For completeness, I must report that I changed the defaults of wxc_scale and wxu_scale (but surely that cannot be the problem...). Switching off the TorsionSimAnnealing does not help. Additionally, if I click on the abort button on the bottom-left of the page nothing happens. On the other hand if I click on the abort icon at the top of the page I get the message "You are not currently running any calculations. Do you want to close this window?" If I say Ok the following message appears: "There is at least one process running; are you sure you want to exit before it is finished" Other tasks (phenix.xtriage for example) seem to run fine. Any suggestion? Best Roberto

3 2

phenix.model_vs_data or something else for 'most disagreeable reflections' ?
by Francis E Reyes 11 Feb '10

11 Feb '10

Hi phenixer's I'm looking for a utility that'll calculate Fc given a refined model and compare it to Fo not on a global scale, but on an individual reflection scale. ShelX does this a reflection list of 'most disagreeable reflections'. Is this possible in phenix? Or another utility (besides shelx)? 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

6 8

Re: [phenixbb] Omit map Problem
by Carruthers, Carl W. 10 Feb '10

10 Feb '10

Coot works fine. When I click the GUI button after the run, it loads the pdb file that I used as a start model and nothing else. When I try opening the composite map manually (using coot: "open mtz", "auto open mtz" or even "open map" ) it says at the bottom "read map /file path/ resolve_composite_map.mtz failed." Thanks again! Carl Methodist. Leading Medicine. Recognized by U.S.News & World Report as an "Honor Roll" hospital Named to FORTUNE magazine's "100 Best Companies to Work For" list five years in a row Designated as a Magnet hospital for excellence in nursing Visit us at methodisthealth.com Follow us at twitter.com/MethodistHosp and www.facebook.com/methodisthospital ***CONFIDENTIALITY NOTICE*** This e-mail is the property of The Methodist Hospital and/or its relevant affiliates and may contain restricted and privileged material for the sole use of the intended recipient(s). Any review, use, distribution or disclosure by others is strictly prohibited. If you are not the intended recipient (or authorized to receive for the recipient), please contact the sender and delete all copies of the message. Thank you.

2 1

Omit map Problem
by Carruthers, Carl W. 10 Feb '10

10 Feb '10

When making an composite omit map in Phenix 1.6 GUI, the final "resolve_composite_map.mtz" file won't load in coot or pymol. When I click the little magnifying glass on the Summary page, Phenix states that it can't find the file, even though I can see it when using Finder (Mac OS). The final lines of the run status are below: Putting all omit results in AutoBuild_run_61_/OMIT OMIT results in AutoBuild_run_61_/OMIT : ['overall_best.pdb_OMIT_REGION_1', 'overall_best.pdb_OMIT_REGION_10', 'overall_best.pdb_OMIT_REGION_11', 'overall_best.pdb_OMIT_REGION_12', 'overall_best.pdb_OMIT_REGION_2', 'overall_best.pdb_OMIT_REGION_3', 'overall_best.pdb_OMIT_REGION_4', 'overall_best.pdb_OMIT_REGION_5', 'overall_best.pdb_OMIT_REGION_6', 'overall_best.pdb_OMIT_REGION_7', 'overall_best.pdb_OMIT_REGION_8', 'overall_best.pdb_OMIT_REGION_9', 'overall_best_denmod_map_coeffs.mtz_OMIT_REGION_1', 'overall_best_denmod_map_coeffs.mtz_OMIT_REGION_10', 'overall_best_denmod_map_coeffs.mtz_OMIT_REGION_11', 'overall_best_denmod_map_coeffs.mtz_OMIT_REGION_12', 'overall_best_denmod_map_coeffs.mtz_OMIT_REGION_2', 'overall_best_denmod_map_coeffs.mtz_OMIT_REGION_3', 'overall_best_denmod_map_coeffs.mtz_OMIT_REGION_4', 'overall_best_denmod_map_coeffs.mtz_OMIT_REGION_5', 'overall_best_denmod_map_coeffs.mtz_OMIT_REGION_6', 'overall_best_denmod_map_coeffs.mtz_OMIT_REGION_7', 'overall_best_denmod_map_coeffs.mtz_OMIT_REGION_8', 'overall_best_denmod_map_coeffs.mtz_OMIT_REGION_9', 'resolve_composite_map.mtz'] Combining all 12 omit regions in AutoBuild_run_61_/OMIT No background map in composite map Composite omit result is in /Users/carlcarruthers/Crystallography/AutoBuild_run_61_/OMIT/resolve_composite_map.mtz Possible operation: finished Required antecedents: [] Expected results: [] Application method: finished Yes, have antecedents... Resulting expectation:10000.0 Am I doing something wrong? Is it a bug? Thanks for helping a crystallogaophy newb (this is my first model!)... Carl Methodist. Leading Medicine. Recognized by U.S.News & World Report as an "Honor Roll" hospital Named to FORTUNE magazine's "100 Best Companies to Work For" list five years in a row Designated as a Magnet hospital for excellence in nursing Visit us at methodisthealth.com Follow us at twitter.com/MethodistHosp and www.facebook.com/methodisthospital ***CONFIDENTIALITY NOTICE*** This e-mail is the property of The Methodist Hospital and/or its relevant affiliates and may contain restricted and privileged material for the sole use of the intended recipient(s). Any review, use, distribution or disclosure by others is strictly prohibited. If you are not the intended recipient (or authorized to receive for the recipient), please contact the sender and delete all copies of the message. Thank you.

2 1