May 2012 - phenixbb - phenix-online.org

links in phenix summary
by Christian Roth 04 Apr '13

04 Apr '13

Dear all, thanks to the phenix team which helped me a lot to resolve all issues regarding the correct link definitions. I have summarised a few things which one might consider if nonstandard links are necessary. Try phenix ready set to get your links out of the pdb. Do not forget to change the possible LINKR to LINK! Inspect the links and decide if they are useful! Peptide links are automatically incorporated if residues are close (inspect geo file) Non standard peptide links could be created using using an now available alpha release script: elbow.ligand_linking model.pdb Inspect the generated cif file and correct angles and incorporate dihedrals and so on if necessary. Sugar links which exists in the mon_lib.cif should be incorporated using the apply_cif link option like refinement.pdb_interpretation.apply_cif_link { data_link = ALPHA1-4 residue_selection_1 = chain B and resname GLC and resseq 1 residue_selection_2 = chain B and resname GLC and resseq 2 } According to the link definition the C1 oxygen of the following sugar with the higher number must be removed. The automatic incorporation of sugar links will be probably soon available in phenix. Nigel is already working on this. If one want use hydrogens in the refinement, they have to be added previously with ready_set and afterwards the OH4 sugar hydrogen has to be removed before the refinement. For links with symmetry related atoms the symmetry operator must be included in a custom bond definition efinement.geometry_restraints.edits { bond { action = *add delete change atom_selection_1 = chain A and resname HEM and resid 154 and name FE atom_selection_2 = chain A and resname MTO and resid 155 and name O symmetry_operation = -x,y,-z distance_ideal = 2.2 sigma = 0.01 } } The operator can be obtained using iotbx.show_distances your.pdb > all_distances Search for keyword "sym" or from Coot clicking on the symm option. Coot print out the operator and the translation term separately. Both has to be merged in the final operator used in phenix! Christian

3 2

Re: [phenixbb] SA omit map
by Huang, Raven H 01 Jun '12

01 Jun '12

Hi all, Was the issue of SA omit map debated last July (involving Tiaard Pining, Pavel, and Folmer) resolved in the newer version of Phenix? Specifically, We'd like to generate the same map for a ligand. It is not for the refinement, but for publication (a reviewer asked us to cover the ligand with the map in a figure). Previously, I used cns to do that (sa_omit_map.inp). But we refined our current structures solely in Phenix. If there is a choice, I'd rather not to move everything to cns and set it up just for the sake of generating a omit map. Thanks, Raven ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Raven H. Huang, Associate Professor Department of Biochemistry Center for Biophysics & Computational Biology University of Illinois at Urbana-Champaign 411 Roger Adams Lab., 600 S. Mathews Ave. Urbana, IL 61801 Phone: 217-333-3967; Fax: 217-244-5858

3 3

Occupancy refinement
by Jerome Nwachukwu 01 Jun '12

01 Jun '12

Hi: I have a 2.2Å structure containing a ligand in at least 2 alternate conformations. When the structure is refined via XYZ refinement + occupancy using phenix.refine, both ligand conformers are kicked out of the electron density. Is there a way to prevent this? Thank you, -Jerome Sent from my iPhone

3 3

TLS makes my r-factors shoot through the roof
by Morten Groftehauge 30 May '12

30 May '12

Hi Phenixbb, 3.2A data, 110 Wilson B, I41, twinned. Imperfect inverted repeat dsDNA with 4 protein monomers, NCS not used in refinement. If I don't include TLS I get a nice stable R-factor while RMS goes down. If I include TLS (whether as one group per chain or as chosen by Phenix) this happens: [image: Inline images 1] Density disappears around some of the DNA bases as an obvious difference. If I run Refmac with TLS it just never finishes. While I could just refine without TLS this does make me worried that there is something wrong with my model or my data (besides being terrible). Cheers, Morten -- Morten K Grøftehauge, PhD Pohl Group Durham University

2 5

Wilson B and REMARK 3
by Michael Thompson 28 May '12

28 May '12

Hello Everyone, I have a question regarding the Wilson B information reported in REMARK 3 of a pdb file. I am preparing a pdb file for deposition, and I am wondering which is the most appropriate value to report as the Wilson B in REMARK 3? My inclination is to report the Wilson B for the full dataset calculated by XSCALE. This value differs slightly from the Wilson B value calculated by phenix.refine and printed to the refinement log file (XSCALE=27.6 vs. phenix.refine=19.8). I assume this difference is due to phenix.refine using a low-resolution cutoff of ~20A. My concern stems from the fact that REMARK 3 "presents information on refinement program(s) used and the related statistics," so I'm not sure if it's best to include the Wilson B calculated by the refinement software, or the Wilson B for the whole data set. Thanks, Mike -- Michael C. Thompson Graduate Student Biochemistry & Molecular Biology Division Department of Chemistry & Biochemistry University of California, Los Angeles miket(a)chem.ucla.edu

2 1

Phenix crashes during refinement of a large structure
by Gino Cingolani 18 May '12

18 May '12

Hi all, PHENIX keeps crashing while refining a large structure (>18,000 residues/asym unit). Can anybody make sense of the following error message (see below)? Thanks so much, Gino ------------------------------------------------------------------------------------------------------------------------------------------------ File "/usr/local/phenix-1.7.3-928/cctbx_project/mmtbx/model_statistics.py", line 112, in __init__ self.clashscore_obj.analyze_clashes(hierarchy = pdb_hierarchy) File "/usr/local/phenix-1.7.3-928/cctbx_project/mmtbx/validation/clashscore.py", line 176, in analyze_clashes stdin_lines=tmp_r.as_pdb_string()) File "/usr/local/phenix-1.7.3-928/cctbx_project/libtbx/easy_run.py", line 147, in __init__ close_fds=not subprocess.mswindows) File "/usr/local/phenix-1.7.3-928/cctbx_project/libtbx/subprocess_with_fixes.py", line 598, in __init__ errread, errwrite) File "/usr/local/phenix-1.7.3-928/cctbx_project/libtbx/subprocess_with_fixes.py", line 1015, in _execute_child self.pid = os.fork() ------------------------------------------------------------------------------------------------------------------------------------------------ ****************************************************************************** Gino Cingolani, Ph.D. Associate Professor Thomas Jefferson University Dept. of Biochemistry & Molecular Biology 233 South 10th Street - Room 826 Philadelphia PA 19107 Office (215) 503 4573 Lab (215) 503 4595 Fax (215) 923 2117 E-mail: gino.cingolani(a)jefferson.edu Website: http://www.cingolanilab.org ****************************************************************************** "Nati non foste per viver come bruti, ma per seguir virtute e canoscenza" ("You were not born to live like brutes, but to follow virtue and knowledge") Dante, The Divine Comedy (Inferno, XXVI, vv. 119-120)

2 4

iterative-build omitmap
by ARJUN＠pitt.edu 17 May '12

17 May '12

Hi I am trying to calculate iterative-build omitmap. Initially program run all 126 omit region and created map coeffs file in corresponding TEMPO subdirectory and unable to create final combined map and this reduce program is running three days without creating the final map. I am trying to run combine_map script separately with all omit region coeffs file and I got the following error: >combine_map overall_best_denmod_map_coeffs.mtz_OMIT_REGION_126 combining map from binary file overall_best_denmod_map_coeffs.mtz_OMIT_REGION_126 >omit Full cross-validation will be done Setting n_small_cycle to 3 All done with inputs Building protein Setting up for building protein Standard fill-in residues: ALA A GLY G Standard main-chain atoms: N CA C O CB All main-chain atoms: N CA C O Max dist between adjacent N atoms: 6.000000 List of all residues: GLY ALA SER VAL ILE LEU MET CYS PHE TYR LYS ARG TRP HIS GLU ASP GLN ASN PRO THR List of 1-letter code of residues: GASVILMCFYKRWHEDQNPT List of number of atoms in side-chains: 0 1 2 3 4 4 4 2 7 8 5 7 10 6 5 4 5 4 3 3 checking license file.. solve2.access The date today is 17-may-12. Your license is good until 15-dec-96. >>>>>> CCP4 library signal library_file:Cannot open file (Error) raised in ccp4_file_open2 <<<<<< >>>>>> CCP4 library signal mtz:Cannot open file (Error) raised in MtzGet <<<<<< Error: failed to open file for read! Sorry, unable to open mtz file: exptl_fobs_phases_freeR_flags.mtz Sorry, unable to open mtz file resolve malloc count: 0 resolve exit_info: How to overcome this error problem Thank you P.Arjunan

2 1

New UPDATE-- Difficult dataset and refinement--P422? I422?Twinning?
by Zhou, Tongqing (NIH/VRC) [E] 17 May '12

17 May '12

Dear All, Kay helped me reprocess the data set with XDS, the program pick up I422 as space group. It turned out that when I originally used XDS, I used spindle axis -1, 0, 0 for the APS ID22 detector, it should be 1, 0, 0. The output from XDS: SUBSET OF INTENSITY DATA WITH SIGNAL/NOISE >= -3.0 AS FUNCTION OF RESOLUTION RESOLUTION NUMBER OF REFLECTIONS COMPLETENESS R-FACTOR R-FACTOR COMPARED I/SIGMA R-meas CC(1/2) Anomal SigAno Nano LIMIT OBSERVED UNIQUE POSSIBLE OF DATA observed expected Corr 8.87 4414 841 858 98.0% 2.7% 3.2% 4402 45.77 3.0% 99.9* -18 0.664 577 6.33 7883 1374 1385 99.2% 3.7% 3.9% 7872 35.76 4.1% 99.9* -17 0.687 1108 5.18 10205 1724 1734 99.4% 5.1% 4.9% 10195 29.38 5.6% 99.8* -15 0.748 1468 4.49 12194 2042 2055 99.4% 5.0% 4.7% 12185 30.77 5.4% 99.8* -20 0.727 1764 4.02 13830 2284 2297 99.4% 6.2% 5.9% 13822 24.55 6.7% 99.8* -15 0.736 2019 3.67 15377 2501 2530 98.9% 9.4% 9.3% 15368 16.77 10.3% 99.6* -12 0.752 2245 3.40 16709 2702 2724 99.2% 14.3% 14.5% 16700 11.31 15.6% 99.2* -5 0.761 2430 3.18 18091 2904 2934 99.0% 23.3% 24.6% 18084 7.14 25.4% 98.4* -3 0.726 2635 3.00 18105 2972 3107 95.7% 37.3% 40.4% 18071 4.42 40.7% 96.4* -6 0.665 2624 total 116808 19344 19624 98.6% 8.2% 8.4% 116699 18.79 9.0% 99.8* -10 0.724 16870 Attached please also find the xtriage log for data processed with XDS. Final model with I422 data set at 3 A resolution looked fine to me, r-free went down from ~0.32 to 0.298, the map even let me correctly fitted one place that supposed to have a Cis-Pro. The final numbers for the model is: r_work = 0.2611 r_free = 0.2983 bonds = 0.002 angles = 0.568 Thank you all for all the help! Regards, Tongqing Tongqing Zhou, Ph.D. Staff Scientist Structural Biology Section Vaccine Research Center, NIAID/NIH Building 40, Room 4609B 40 Convent Drive, MSC3027 Bethesda, MD 20892 (301) 594-8710 (Tel) (301) 793-0794 (Cell) (301) 480-2658 (Fax) ****************************************************************************** The information in this e-mail and any of its attachments is confidential and may contain sensitive information. It should not be used by anyone who is not the original intended recipient. If you have received this e-mail in error please inform the sender and delete it from your mailbox or any other storage devices. National Institute of Allergy and Infectious Diseases shall not accept liability for any statements made that are sender's own and not expressly made on behalf of the NIAID by one of its representatives. ****************************************************************** -----Original Message----- From: Pavel Afonine [mailto:[email protected]] Sent: Tuesday, May 08, 2012 2:07 PM To: phenixbb(a)phenix-online.org Subject: Re: [phenixbb] UPDATE-- Difficult dataset and refinement--P422? I422?Twinning? Hi, >> 1. Is it really twinned since including twin-law in refinment gave much >> better numbers? R-factors may not be directly comparable: G.Murshudov, Appl. Comput. Math., V.10, N.2, 2011, pp.250-261. The strategy is: - check to see if Xtriage tells it's twinned (based on reflection statistics, not just predicted twin laws based on geometry); - use phenix.model_vs_data to see if it's worth using it. Pavel _______________________________________________ phenixbb mailing list phenixbb(a)phenix-online.org<mailto:[email protected]> http://phenix-online.org/mailman/listinfo/phenixbb

2 1

annealing omit map
by r n 16 May '12

16 May '12

Hi all I have omitted some region manually from the pdb and did annealing - cartesian using phenix.refinment, the resulting map has density (1.2 sigma) in 2fofc and also seen 2.8 sigma negative density on the omitted region? any suggestions? Rab

1 1

Re: [phenixbb] mr.rosetta gui
by Terwilliger, Thomas C 16 May '12

16 May '12

Hi Alex, OK, the problem should be fixed in tomorrow's build. Thank you for pointing it out! Let us know if that does not do it! All the best, Tom T ________________________________ From: Terwilliger, Thomas C Sent: Wednesday, May 16, 2012 10:45 AM To: PHENIX user mailing list Cc: Terwilliger, Thomas C Subject: RE: [phenixbb] mr.rosetta gui Hi Alex, Yes, I can duplicate this...I'm checking on what is happening now. All the best, Tom T ________________________________ From: phenixbb-bounces(a)phenix-online.org [phenixbb-bounces(a)phenix-online.org] on behalf of Alexander Batyuk [batyuk(a)bioc.uzh.ch] Sent: Wednesday, May 16, 2012 1:54 AM To: PHENIX user mailing list Subject: [phenixbb] mr.rosetta gui Dear developers, While trying to run mr.rosetta from the GUI, I've encountered an error (see screenshot). I have un-checked the box Include solvation energy. The file rosetta_weights.wts has been created in WORK_1, fa_sol=0.0. Inspection of mr_rosetta_params.eff in GROUP_OF_PLACE_MODEL_1 revealed include_solvation_energy = False, followed by weights_file = "/Users/abatyuk/Desktop/o1map/MR_ROSETTA_74/WORK_1/rosetta_weights.wts", but it seems that the GUI wants either or. Is that the expected behavior? I'm running 1.8-1056, Rosetta 3.4 on OS X 10.6.8 Best wishes, Alex

1 0