- phenixbb - phenix-online.org

Re: [phenixbb] clash guarde failure phenix coot
by Pavel Afonine 19 Jan '10

19 Jan '10

Hi Georg, I think the refinement .log file should contain the lines like these below: Warning: very small nonbonded interaction distances. Nonbonded interactions: 12604 Sorted by model distance: nonbonded pdb=" CG ARG A 3 " pdb=" CG ARG A 4 " model vdw 0.000 3.840 nonbonded pdb=" NE ARG A 3 " pdb=" NE ARG A 4 " model vdw 0.000 3.200 which pinpoints you the offending atoms. Does this answer your question? Pavel. On 1/19/10 9:57 AM, Georg Mlynek wrote: > > Dear All, > > I´m refining a 2.1A with phenix 1.5-2 and get when I click on > automatically add hydrogens to model following error message: > refinement.pdb_interpretation.clash_guard.failure: Number of nonbonded > interaction distances < 0.5:391 > > I suppose that this comes, because coot mixes up the hydrogens. >

1 0

rigid body refinement problem
by Jianghai Zhu 19 Jan '10

19 Jan '10

Hi, I am using phenix.refine_dev_282 for a rigid body refinement. The output pdb file is exactly the same as the input model even though there are obviously some large movements. Here are some information from the log file. Information about total rigid body shift of selected groups: rotation (deg) translation (A) xyz total xyz total group 1: 0.051 -5.493 -1.073 5.60 0.55 -3.15 -1.35 3.47 group 2: 0.254 0.562 0.252 0.67 0.36 -3.99 0.87 4.10 group 3: 3.807 4.001 -2.262 5.92 -0.81 -3.17 1.06 3.44 group 4: 2.275 -4.215 1.909 5.12 -0.18 -2.95 -1.43 3.28 group 5: -1.902 -5.750 -3.206 6.81 -1.99 -4.44 0.27 4.88 group 6: -4.932 -3.939 -5.340 8.16 -0.18 -2.08 -0.25 2.11 group 7: -1.140 -0.769 -7.715 7.83 0.73 -2.91 0.04 3.00 group 8: 2.008 -1.858 -1.236 3.02 -1.10 -5.44 -0.82 5.61 But then the model didn't change at all. stage Deviation of refined model from start model max min mean 0 : 0.000 0.000 0.000 1_bss: 0.000 0.000 0.000 1_rbr: 0.000 0.000 0.000 1_bss: 0.000 0.000 0.000 1_adp: 0.000 0.000 0.000 1_bss: 0.000 0.000 0.000 Is this a bug? -- Jianghai

2 2

difference density map
by 丁玮 19 Jan '10

19 Jan '10

4 3

Where is phenix.refine finding starting r values?
by Bart, Aaron G 18 Jan '10

18 Jan '10

Greetings, I am fairly new to crystallography, and have just started to use Phenix for structural refinement. I am currently working on a 2.8A structure that was run through multiple rounds of rigid body refinement via Refmac and manual revision using Coot. I eventually got stuck where large portions of my molecule had very poor density and the rwork and rfree stayed around 0.3. I then performed a TLS refinement which resulted in better density, and tried running it through phenix.refine with default ADP settings, as well as simulated annealing. It worked beautifully, however when I looked at the log file for the refinement, phenix was reporting a starting rwork and rfree that represented my first stage of refinement in refmac, 0.38/0.44. My question is, is phenix somehow neglecting the refinement that was done in refmac? Regards -Aaron

2 1

Re: [phenixbb] merohedral twinned refinement
by JXQI 18 Jan '10

18 Jan '10

Zwart, Cherney and all, Thanks for the suggestions.The problem is fixed when I tried P31. Best regards J.X. Qi > > 3. If I try P3,the Matthews coeffs tell 2 or 3 molecules in the ASU,the > molecules will crash.Besides,the overall FOM ( lower than 0.4) and the > R/Rfree values (about 0.45/0.56) suggest that the P3 is the wrong sg. try P31 or otherwise P1. You should get a solution in those groups. It could of course mean that there are OD type issues, but I would not jump there yet. Can you send the intensity statistics from xtriage?

1 0

right symmetry?
by Gabriela Guédez 18 Jan '10

18 Jan '10

Hi Phenixbb, I have a data that indexes in P6122, also P61 and P31. Data analysis in phenix.xtriage says no twinning is suspected according to the L-test. However, after MR and twin RB refinement, the model improves from P6122 or P61 to P31 5% Rfacs and 0.5 the FOM. I thought the data might have RPS and twinning not detected with the L-test. So, i ran xtriage to compute the RvsR statistics using the obs_and_calc.mtz file and the results were: in P61 R_abs_twin observed data : 0.112 R_abs_twin calculated data : 0.047 in P31 R_abs_twin observed data : 0.117 R_abs_twin calculated data : 0.139 "misespecified crystal symmetry", according to Lebedev, Vagin, Murshudov et al 2006, Acta Cryst. D62,83-95. I would like to know a test that tells me the right symmetry. I appreciate your help. Thanks, Gabriela

2 4

Re: [phenixbb] merohedral twinned refinement
by JXQI 18 Jan '10

18 Jan '10

Dr.Zwart, Thanks for the reply. 2. The systematic absences and the Rsym support the P3121.I already mentioned the model fit very with the map if I use the P3121. 3. If I try P3,the Matthews coeffs tell 2 or 3 molecules in the ASU,the molecules will crash.Besides,the overall FOM ( lower than 0.4) and the R/Rfree values (about 0.45/0.56) suggest that the P3 is the wrong sg. 4. TLS refinement also did not imrove the results. 5. There are totally 88 aa in my protein, I could build 86 aa clearly. The intensity ststistics is attached in the message. J.X. QI ************************************************************************ Shell I/Sigma in resolution shells: Lower Upper % of of reflections with I / Sigma less than limit limit 0 1 2 3 5 10 20 >20 total 50.00 4.95 0.9 1.5 1.9 2.4 2.8 4.0 9.5 90.5 100.0 4.95 3.93 0.0 0.2 0.3 0.5 0.8 2.0 6.6 93.4 100.0 3.93 3.44 0.2 0.3 0.6 1.0 2.5 5.4 15.4 84.6 100.0 3.44 3.12 0.0 0.2 0.5 0.7 2.5 7.6 22.7 77.3 100.0 3.12 2.90 0.3 1.0 2.1 3.2 7.2 17.6 41.4 58.6 100.0 2.90 2.73 0.0 0.7 2.0 3.2 9.9 20.8 52.9 47.1 100.0 2.73 2.59 0.0 1.4 2.9 5.7 12.4 30.0 60.8 39.2 100.0 2.59 2.48 0.7 1.6 4.7 7.7 16.3 36.0 71.5 28.5 100.0 2.48 2.38 0.6 2.1 5.0 10.2 21.9 42.0 76.2 23.8 100.0 2.38 2.30 0.8 2.3 7.4 12.7 25.3 49.4 83.4 16.6 100.0 All hkl 0.4 1.1 2.7 4.7 10.0 21.1 43.4 56.6 100.0 Summary of reflections intensities and R-factors by shells R linear = SUM ( ABS(I - <I>)) / SUM (I) R square = SUM ( (I - <I>) ** 2) / SUM (I ** 2) Chi**2 = SUM ( (I - <I>) ** 2) / (Error ** 2 * N / (N-1) ) ) In all sums single measurements are excluded Shell Lower Upper Average Average Norm. Linear Square limit Angstrom I error stat. Chi**2 R-fac R-fac 50.00 4.95 619.8 13.8 3.5 1.039 0.046 0.050 4.95 3.93 834.6 17.5 5.0 1.156 0.048 0.052 3.93 3.44 521.9 14.6 4.5 1.026 0.059 0.062 3.44 3.12 333.7 11.1 4.1 1.073 0.071 0.068 3.12 2.90 172.5 7.1 3.6 1.153 0.098 0.096 2.90 2.73 131.4 6.0 3.4 1.256 0.116 0.108 2.73 2.59 102.8 5.2 3.4 1.306 0.140 0.130 2.59 2.48 77.7 4.6 3.4 1.381 0.173 0.158 2.48 2.38 64.1 4.3 3.3 1.355 0.195 0.177 2.38 2.30 54.7 4.1 3.3 1.345 0.221 0.197 All reflections 297.8 8.9 3.7 1.209 0.071 0.059 Intensities of systematic absences h k l Intensity Sigma I/Sigma 0 0 2 -0.2 0.4 -0.5 0 0 4 0.1 1.0 0.1 0 0 5 0.4 1.4 0.3 0 0 7 1.1 2.5 0.4 0 0 8 0.4 2.1 0.2 0 0 10 7.3 3.4 2.2 0 0 11 6.0 3.5 1.7 0 0 13 -2.3 4.9 -0.5 0 0 14 6.1 4.3 1.4 0 0 16 8.1 5.1 1.6 0 0 17 3.9 4.9 0.8 0 0 19 0.8 4.5 0.2 0 0 20 4.0 4.6 0.9 ********************************************************************* 1. How do the intensity statistic look like? 2. How did you choose P3121? 3. Did you try solving and refining in P3? 4. TLS? 5. How complete is your model 2010/1/17 JXQI <jxqi at mail.im.ac.cn>: > Dear all, > > I have a 2.3A data set processed with HKL2000. The space group is P3121 with > 63% of solvent content and 1 molecule in the ASU. I used a 40% homology > model to do molecular replacement. Both the Phaser and Molrep gave similar > solution with good packing. After further refinement with phenix, the built > model fit well with the 2mFo-DFc map. However the R/Rfree were stuck to > 0.34/0.38. > > I rechedked the data set with xtriage and found the data set is merohedral > twinned (with twin law"-h,-k,l", the estimated twin fraction is about 0.05). > When I used the twin law in further refinement, the R/Rfree were only > decreased to 0.32/0.36. Any suggestions are appreciated, thanks! > > > J.X. QI > _______________________________________________ > phenixbb mailing list > phenixbb at phenix-online.org > http://phenix-online.org/mailman/listinfo/phenixbb

3 2

merohedral twinned refinement
by JXQI 18 Jan '10

18 Jan '10

Dear all, I have a 2.3A data set processed with HKL2000. The space group is P3121 with 63% of solvent content and 1 molecule in the ASU. I used a 40% homology model to do molecular replacement. Both the Phaser and Molrep gave similar solution with good packing. After further refinement with phenix, the built model fit well with the 2mFo-DFc map. However the R/Rfree were stuck to 0.34/0.38. I rechedked the data set with xtriage and found the data set is merohedral twinned (with twin law"-h,-k,l", the estimated twin fraction is about 0.05). When I used the twin law in further refinement, the R/Rfree were only decreased to 0.32/0.36. Any suggestions are appreciated, thanks! J.X. QI

3 2

Twin fraction in phenix refine
by Pietro Roversi 17 Jan '10

17 Jan '10

Dear all, phenix.xtriage estimates a twin fraction of between 24 and 30% but then when I start the refinement I get 47-49%. In CNS if I fix the twin fraction in a series of refinements the Rfactors jump from 44% (twin fraction 0) to 37% (from 25% to 49%, only marginal improvement upon going from alpha=25% to 49%). Questions: 1. Is the refinement of alpha dependent on the model? 2. I am tempted to work with the 30% that xtriage estimates but do not know how: is there a way to fix alpha? Thanks Pietro -- Sent from my Desktop Pietro Roversi EP Abraham Fellow in Biochemistry - Lincoln College - Oxford Sir William Dunn School of Pathology, Oxford University South Parks Road, Oxford OX1 3RE, England UK Tel. 0044-1865-275385

8 13

Re: [phenixbb] [ccp4bb] Adding H in refinement
by Schubert, Carsten [PRDUS] 13 Jan '10

13 Jan '10

Hi Pavel, sorry I hijacked the thread a bit... A question about the riding Hydrogens and their B-factors. I noticed that at least with earlier versions of PHENIX the B-factor of the Hydrogens was not fixed to the B-factor of the bonded atoms. If I would generate the Hydrogens in reduce they would inherit the appropriate B-factors, but during the course of the refinement and rebuilding they would start to diverge, sometimes drastically. Are you saying the current implementation restrains the B-factors to the bonded atoms and forces them to be equal or are they just inherited and are allowed to float afterwards? Cheers, Carsten > -----Original Message----- > From: CCP4 bulletin board [mailto:[email protected]] On Behalf Of > Pavel Afonine > Sent: Tuesday, January 12, 2010 12:24 PM > To: CCP4BB(a)JISCMAIL.AC.UK > Subject: Re: [ccp4bb] Adding H in refinement > > Hi Sara, > > - what you observe should not happen since phenix.refine uses riding > model for H atoms. The hydrogen's B-factors are automatically inherited > from the atoms these hydrogens are bonded to. For example, in X-H bond > the B-factor of X should be equal to B-factor of H. > > - make sure you are using the latest version of PHENIX: > http://www.phenix-online.org/download/ > > - if this does not help, please contact me directly, and send the data > and model so I can tell you what exactly is not right (all the files > will be handled confidentially). > > - there is PHENIX bulletin board: > http://www.phenix-online.org/mailman/listinfo/phenixbb > > Pavel. > > > I did refinement with phenix of my 1.9-2.0 Angstroem structures and > included the hydrogens (riding). However, when I checked on the > statistics (refinement close to the end), the average B-factor was > extremely high (in phenix.polygon it was higher than with any other > structure in similar resolution range). It makes sense though that this > happens, if you have a residue which has a high B-factor and carries a > lot of hydrogens the average B-factor will raise quite a lot (since the > B-factor of hydrogens is calculated 1-1.5x of B-factor from the atom > it sits on), right? > > (When I removed the hydrogens again, the average B-factor was > fine...) > > > > My question is now did I do something wrong in my refinement (-> do I > have to change something that this does not happen), or is this > something everybody sees? > > > > If this is common, what would happen (during evaluation) if you want > to publish a structure and the statistics show such a high average B- > factor? Is it better to have hydrogens on, but a bad average B, or no > hydrogens on and a good average B... Obviously, I never published a > structure ;) > >

2 1

Beamtime at the ALS
by Peter Zwart 12 Jan '10

12 Jan '10

Dear All, January 15, 2010 is the deadline for the March/April 2010 Rapid Access Proposal cycle. All Berkeley Center for Structural Biology(BCSB) beamlines are equipped with ADSC Q315/Q315R detectors, automated sample changers and data collection software enabling high-throughput crystal screening and data collection. Remote data collection is available on all BCSB beamlines, providing the user with the full complement of sample visualization, sample manipulation, beamline control, data acquisition and data analysis tools exactly as they would see them if they were stationed at the beamline. This enhanced remote operation capability is coupled with 22hr onsite support by BCSB staff who are able to assist immediately with loading additional samples for remote users or troubleshoot any issues that might arise. Remote users can furthermore be kept up-to-date on changes in ring status via an SMS service (http://bcsb.als.lbl.gov/wiki/index.php/Ring_Status_Notifications_to_Cell_Ph…) Specific features are summarized below. Beamlines 8.2.1 and 8.2.2: To facilitate studies on small crystals, a microdiffractometer was installed in the beamline 8.2.1 endstation. The new equipment allows precise sample positioning to within 2 microns, excellent sample viewing of very small crystals, and an off-axis crystal positioning stage. Both beamlines 8.2.1 and 8.2.2 feature a Rigaku sample changer (Actor), allowing remote operations to now be a routine mode of access for these beamlines. Beamlines 5.0.1, 5.0.2, 5.0.3: The Berkeley Automounter sample handling system has a routine capacity of 96 samples (6 pucks). In a typical high-throughput screening mode, the mount-to-mount time is around 2.5 minutes per sample, allowing users to screen a full puck within 45 minutes. The sector 5 beamline user stations are equipped with fully high-adjustable, ergonomically friendly work stations. Data analyses in the BCSB is facilitated by software maintained by sbgrid (http://www.sbgrid.org). A 16 core linux machine is available for our users to process their data and solve/refine their structure. An additional mode of access to the BCSB beamlines is through the Collaborative Crystallography (CC) Program. Users apply for beamtime via the general user program, and collaborate with an expert crystallographer who will conduct the experiments and data reduction on behalf of the researchers. Depending on the users, structure solution, model building and refinement can be carried out as well. Please contact BCSBBeamtime(a)lbl.gov for more information. Please visit http://bcsb.lbl.gov/ for more details about the Center and its beamlines. To find out more, click on: http://www.als.lbl.gov/als/quickguide/independinvest.html We invite you to submit a proposal at: http://alsusweb.lbl.gov/ Scroll down to "Structural Biology beamines (includes protein SAXS)." Click on "New Proposal." If you'd like to apply for May/June 2010 beamtime at the Advanced Light Source, please submit a General User proposal by March 15, 2010. If you have any questions or would like to request open beamtime, please e-mail BCSBBeamtime(a)lbl.gov. Please note that executed user agreements must be received by LBNL prior to beamtime. Proprietary fees, if applicable, must be received by LBNL at least five working days prior to scheduled beamtime. -- ----------------------------------------------------------------- P.H. Zwart Beamline Scientist Berkeley Center for Structural Biology Lawrence Berkeley National Laboratories 1 Cyclotron Road, Berkeley, CA-94703, USA Cell: 510 289 9246 BCSB: http://bcsb.als.lbl.gov PHENIX: http://www.phenix-online.org CCTBX: http://cctbx.sf.net -----------------------------------------------------------------

1 0

January 15, 2009 deadline- User proposal submission for Collaborative Crystallography at BCSB
by Banumathi Sankaran 12 Jan '10

12 Jan '10

Dear Users, The deadline for Mar/April 2010 Collaborative Crystallography proposals will be *January 15, 2009. * ------------------------------------------------------------------------ Through the Collaborative Crystallography Program (CC) at the Advanced Light Source (ALS), scientists can send protein crystals to Berkeley Center for Structural Biology (BCSB) staff researchers for data collection and analysis. The CC Program can provide a number of benefits to researchers: * Obtain high quality data and analysis through collaborating with expert beamline researchers; * Rapid turn around on projects; and * Reduced travel costs. To apply, please submit a proposal through the ALS General User proposal review process for beamtime allocation. Proposals are reviewed and ranked by the Proposal Study Panel, and beamtime is allocated accordingly. BCSB staff schedule the CC projects on Beamlines 5.0.1 and 5.0.2 to fit into the available resources. Only non-proprietary projects will be accepted. As a condition of participation, BCSB staff researchers who participate in data collection and/or analysis must be appropriately acknowledged - typically being included as authors on publications and in PDB depositions. Please consult the website for additional information at: http://bcsb.als.lbl.gov/wiki/index.php/Collaborative_Crystallography ----------------------------------------------------------------------------------------- How To Apply: To learn more, go to: http://www-als.lbl.gov/als/quickguide/independinvest.html To submit a proposal, go to: http://alsusweb.lbl.gov/. Scroll down to "*Structural Biology beamlines (includes protein SAXS)*" and click on "New Proposal." Enter your proposal information, with attention to the following details: * For question 3/First choice, select "5.0.1 (Monochromatic);" for question 3/Second choice, select "5.0.2 (MAD)." * Check box (yes) in response to question (7) "Do you want collaborative crystallography (beamline 5.0.1 or 5.0.2 only) * In question 4, please describe a specific research project with a clear end point. In order to request CC time for Mar/April 2010 allocation period, proposals must be submitted by *January 15, 2009.* The deadline for CC proposals for the time period May/June 2010 will be March 15, 2010. Regards, Banumathi Sankaran

1 0

reduce and ready_set problem
by Jianghai Zhu 11 Jan '10

11 Jan '10

Hi, I believe I have reported this bug before, but it is not fixed in the phenix.reduce and phenix.ready_set yet. The problem is the hydrogens added by reduce when there is a glycan. Here is the error message. Number of atoms with unknown nonbonded energy type symbols: 12 "ATOM 8332 HD22 ASN B 99 .*. H " "ATOM 11820 HD22 ASN B 320 .*. H " "ATOM 12615 HD22 ASN B 371 .*. H " "ATOM 14058 HO4 NAG B3320 .*. H " "ATOM 14086 HO4 NAG B3321 .*. H " "ATOM 14136 HO4 NAG B3371 .*. H " "ATOM 22458 HD22 ASN D 99 .*. H " "ATOM 25955 HD22 ASN D 320 .*. H " "ATOM 26750 HD22 ASN D 371 .*. H " "ATOM 28256 HO4 NAG D3320 .*. H " ... (remaining 2 not shown) ... Sorry: Fatal problems interpreting PDB file: Number of atoms with unknown nonbonded energy type symbols: 12 Please edit the PDB file to resolve the problems and/or supply a CIF file with matching restraint definitions, along with apply_cif_modification and apply_cif_link parameter definitions if necessary. Also note that phenix.ready_set and phenix.elbow are available for creating restraint definitions (CIF files). The hydrogen left on the N atom of the ASN is correct, but should be named as HD21 as defined in the ccp4 cif library. And the HO4 for the NAG glycan is just wrong. I have to delete all these hydrogens to get phenix.refine going. -- Jianghai

2 1

Resolve after Phaser SAD module - question
by Peter Grey 10 Jan '10

10 Jan '10

Dear Phenix experts, It is recommended to use FWT/PHWT after Phaser SAD module for initial map calculation. Could you please explain how to tell Resolve to pick these fields for the initial map ? With MR it is clear since you need FWT and PHIC but with SAD you really need FWT and PHWT . I thank you in advance, Peter

2 1

Phenix.refine RuntimeError
by JXQI 09 Jan '10

09 Jan '10

Dear PHENIXBBs, I tried to refine a N-linked glycan using phenix.refine and used the data_link=NAG-ASN for ASN-NAG connection. But I got a RuntimeError of unknown atom_id: HD22. From the mon_lib_list.cif file of Phenix-1.5-2, the atom_id HD22 is already defined. Any suggestions are appreciated. Thanks J.X. QI ********************************************************* Monomer Library directory: "/home/jxqi/152Phenix/phenix-1.5-2/chem_data/mon_lib" Total number of atoms: 11844 apply_cif_link: data_link: NAG-ASN mod_id_1: DEL-O1 mod_id_2: DEL-HD22 residue_selection_1: chain A and resname ASN and resid 170 residue_selection_2: chain A and resname NAG and resid 1 apply_cif_link: data_link: NAG-ASN mod_id_1: DEL-O1 mod_id_2: DEL-HD22 residue_selection_1: chain C and resname ASN and resid 170 residue_selection_2: chain C and resname NAG and resid 1 apply_cif_link: data_link: NAG-ASN mod_id_1: DEL-O1 mod_id_2: DEL-HD22 residue_selection_1: chain E and resname ASN and resid 170 residue_selection_2: chain E and resname NAG and resid 1 Number of models: 1 Model: "" Number of chains: 7 Chain: "A" Number of atoms: 2559 Number of conformers: 1 Conformer: "" Traceback (most recent call last): File "/home/jxqi/152Phenix/phenix-1.5-2/phenix/phenix/command_line/refine.py", line 11, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/home/jxqi/152Phenix/phenix-1.5-2/phenix/phenix/refinement/command_line.py", line 68, in run log=log) File "/home/jxqi/152Phenix/phenix-1.5-2/phenix/phenix/refinement/command_line.py", line 227, in __init__ self.pdb_file_names) File "/home/jxqi/152Phenix/phenix-1.5-2/cctbx_project/mmtbx/utils.py", line 907, in process_pdb_files stop_if_duplicate_labels = stop_if_duplicate_labels) File "/home/jxqi/152Phenix/phenix-1.5-2/cctbx_project/mmtbx/utils.py", line 940, in _process_pdb_file log = self.log) File "/home/jxqi/152Phenix/phenix-1.5-2/cctbx_project/mmtbx/monomer_library/pdb_interpretation.py", line 2995, in __init__ log=log) File "/home/jxqi/152Phenix/phenix-1.5-2/cctbx_project/mmtbx/monomer_library/pdb_interpretation.py", line 2105, in __init__ log=log) File "/home/jxqi/152Phenix/phenix-1.5-2/cctbx_project/mmtbx/monomer_library/pdb_interpretation.py", line 1575, in __init__ prev_mm=prev_mm) File "/home/jxqi/152Phenix/phenix-1.5-2/cctbx_project/mmtbx/monomer_library/pdb_interpretation.py", line 543, in __init__ mod_mod_id=self.mon_lib_srv.mod_mod_id_dict[apply_data_mod]) File "/home/jxqi/152Phenix/phenix-1.5-2/cctbx_project/mmtbx/monomer_library/pdb_interpretation.py", line 819, in apply_mod mod_mon = self.monomer.apply_mod(mod_mod_id) File "/home/jxqi/152Phenix/phenix-1.5-2/cctbx_project/mmtbx/monomer_library/cif_types.py", line 359, in apply_mod result.delete_atom_in_place(mod_atom.atom_id) File "/home/jxqi/152Phenix/phenix-1.5-2/cctbx_project/mmtbx/monomer_library/cif_types.py", line 264, in delete_atom_in_place raise RuntimeError("delete_atom_in_place: unknown atom_id: %s" % atom_id) RuntimeError: delete_atom_in_place: unknown atom_id: HD22 **********************************************************

2 1

Hello
by Pankaj Chauhan 06 Jan '10

06 Jan '10

Hello, While during NCS based restrained refinement, how can we take two chains as single translation unit as the pdb model is a dimer-dimer. Thanks -- Pankaj Chauhan Senior Research Fellow, PhD Scholar Crystallography Laboratory Division of Protein Science & Engg. IMTECH , Chandigarh, INDIA

1 0

coot - display map
by r n 05 Jan '10

05 Jan '10

Hi I am using Phenix GUI, after refinement , tried to load the coot interface in Phenix, it opens up the coot, but did not load the map or model file. I am using coot0.6 version, Is there any place in phenix enviroment that I could mention the path for coot? Thanks

2 4

refinement
by r n 05 Jan '10

05 Jan '10

Hi I am doing refinement using phenix. I did simulated annealing, Rfree goes up to 10% more than before. (40%) and R stays at 35%. If use previous phenix refinment file for same pdb and data , R and Rfree is 25/27%. I thought even different r-free set after simulated annealing and 7 cycles of refinement in phenix should bring Rfree down, but not ? But same data set and model , in CNS, R/Rfree 26/28%. Any suggestions will be appreciated. I am using phenix-dev249 - linux. thanks ram

2 1

text output for model quality and validation report at phenix
by Ravindra D. Makde 04 Jan '10

04 Jan '10

Hi all, I will greatly appreciate the answers to following query: How can I get the text or excel output of the model quality and validation report of a model after phenix refinement? I am specially interested in real space CC of the various residues. I am using phenix.refine_gui at phenix-1.4-162. I have searched all output files, but couldn't find. Thanks in advance. Ravi Ravindra D. Makde, PhD Postdoctoral fellow, Tan lab, The Department of Biochemistry and Molecular Biology, The Pennsylvania State University, University Park, PA 16802 USA

3 2

optimize phenix
by r n 04 Jan '10

04 Jan '10

hi all I am trying to optimize the model, in which secondary structures are distorted and I did pdbtools , geometry regulatization , after optimization still the helices are distorted. Is it any tools in the phenix that could optimize the secondary structures alone in the pdb? Thanks ram

3 4

some sort of occupancy refinement
by Bryan W. Lepore 01 Jan '10

01 Jan '10

question is after the setup here : i have a ligand with multiple conformers. i gave each conformer a different chain label. refining occupancy for the overlapped parts kicks that part out of the density. i have remove_selection for the parts that are 1.0 and they seem OK. a single conformer stays in the density. i tried starting with and without site or occupancy refinement, no help. i do not want group occupancy refinement b/c some parts do not overlap - i.e. they "should" be occupancy 1.0 -- i tried it anyways with constrained_group =... to check. is the chain label an improper way to do this? -bryan p.s: Happy New Year with Blue Moon Eclipse (~10degW->180degE long. eastern HS)!

4 3

high resolution Anomalous map
by vennila Natesan 31 Dec '09

31 Dec '09

Hello all... How do we create anomalous map using high resolution reflections? For example, if we have 2.2A resolution data (i.e 20-2.2A), i want to know how to create the anomalous map for the same dataset but from 2.4-2.2A resolution range. Vennila.

4 8

phenix.refine change: support for residues on symmetry elements
by Ralf W. Grosse-Kunstleve 27 Dec '09

27 Dec '09

The current development version of phenix.refine (dev-268 or later) supports residues located on symmetry elements, e.g. SO4 on a three-fold axis. phenix.refine now computes structure factors slightly differently compared to older versions. The new rules are: - Atoms with occupancy == 1 are subject to automatic handling of the symmetry. Example: An atom near a three-fold axis (distance to symmetry copy <= 0.5 A) is moved onto the axis and constrained in refinement to stay on the axis. The "special position factor" (1/3 is this case) is taken into account internally. The factor does not appear in the occupancy. - Atoms with occupancy != 1 are assumed to take the special position factor into account in the occupancy. Automatic handling of the symmetry is turned off. Example: A disordered SO4 molecule near a three-fold axis should have occupancies of 0.33 for all atoms. - If all atoms in a residue have occupancies != 1, all nonbonded symmetry interactions between the atoms in the residue are disabled. These three rules are meant - to bring phenix.refine in line with other programs and the PDB. - "to do the right thing" in most cases without manual intervention. Please let me know if you see problems with the new way phenix.refine uses occupancies to determine special position constraints and exclusions for nonbonded symmetry interactions. Does anyone need atoms with partial occupancies but constrained to stay on a symmetry axis during refinement? Ralf

6 7

Re: [phenixbb] Phaser SAD question
by Kitty Pattridge 24 Dec '09

24 Dec '09

"Randy J. Read" <rjr27(a)cam.ac.uk> wrote: >Dear Peter, > >I'm interested in your observation that you get better signal from an >anomalous difference map computed with density-modified phases than you get >with an LLG map using the MR model. Presumably the density modification >just started from the model phases as well? We've never tried that, so >we'll have to look into that. Unfortunately, you can't at the moment start >the LLG completion with a "model" specified through structure factors, >though this is on our to-do list. > >With the MR+SAD option, because it is iterative, there are sites found in >subsequent cycles (where the information from some of the anomalous >scatterers has become available) and you end up with a set of reliable >sites that can be significantly larger than the set obtained from the first >LLG map. Does your anomalous difference Fourier with density modified >phases give you a larger apparently reliable substructure than the >iterative completion? You could use phenix.emma to compare the two >substructures and see how much they have in common. > >You can indeed give Phaser the set of sites that you got from the anomalous >difference map. One way to choose the cutoff would be to see how deep the >deepest hole is (i.e. how many standard deviations the map goes in the >negative direction), and choose a larger number in the positive direction. >Or you could just use the peaks above 6-7 sigma. > >I'd be interested in hearing how you get on with this approach, off-line if >you prefer. > >Good luck and Happy Christmas! > >Randy Read > >On Dec 24 2009, Peter Grey wrote: > >>Dear Phenix Experts, >> >>I use the SAD+MR option in Phenix. The ASU is huge and the number of >>scatterres is in the few hundreds range. I try to minimize the number of >>runs of Phaser SAD+MR since they take long time on my computer and would >>appreciate your insight regarding the following : >> >> - I noticed that I get much higher signal and much larger number of sites >> if I use density modified phases to calculate the anomalous difference >> map. Is it possible to input these phases instead of the model to Phaser >> ? >> >> - If not, do you recommend providing Phaser with the list of sites found >> in this difference anomalous map. Should such a list include only the >> sites you can be sure of (say above 7 sigma) or all possible sites (say >> above 5 sigma) ? >> >> - Should I use the (peak height)/(sigma of map) as the first estimate of >> occupancy ? if so do I have to normalize it (strongest site has occupancy >> of 1) or can I leave the height/sigma values as they are ? >> >> >>Many thanks, >> >>Peter. >> > >_______________________________________________ >phenixbb mailing list >phenixbb(a)phenix-online.org >http://phenix-online.org/mailman/listinfo/phenixbb > >

1 0

Re: [phenixbb] Phaser SAD question
by Kitty Pattridge 24 Dec '09

24 Dec '09

"Randy J. Read" <rjr27(a)cam.ac.uk> wrote: >Dear Peter, > >I'm interested in your observation that you get better signal from an >anomalous difference map computed with density-modified phases than you get >with an LLG map using the MR model. Presumably the density modification >just started from the model phases as well? We've never tried that, so >we'll have to look into that. Unfortunately, you can't at the moment start >the LLG completion with a "model" specified through structure factors, >though this is on our to-do list. > >With the MR+SAD option, because it is iterative, there are sites found in >subsequent cycles (where the information from some of the anomalous >scatterers has become available) and you end up with a set of reliable >sites that can be significantly larger than the set obtained from the first >LLG map. Does your anomalous difference Fourier with density modified >phases give you a larger apparently reliable substructure than the >iterative completion? You could use phenix.emma to compare the two >substructures and see how much they have in common. > >You can indeed give Phaser the set of sites that you got from the anomalous >difference map. One way to choose the cutoff would be to see how deep the >deepest hole is (i.e. how many standard deviations the map goes in the >negative direction), and choose a larger number in the positive direction. >Or you could just use the peaks above 6-7 sigma. > >I'd be interested in hearing how you get on with this approach, off-line if >you prefer. > >Good luck and Happy Christmas! > >Randy Read > >On Dec 24 2009, Peter Grey wrote: > >>Dear Phenix Experts, >> >>I use the SAD+MR option in Phenix. The ASU is huge and the number of >>scatterres is in the few hundreds range. I try to minimize the number of >>runs of Phaser SAD+MR since they take long time on my computer and would >>appreciate your insight regarding the following : >> >> - I noticed that I get much higher signal and much larger number of sites >> if I use density modified phases to calculate the anomalous difference >> map. Is it possible to input these phases instead of the model to Phaser >> ? >> >> - If not, do you recommend providing Phaser with the list of sites found >> in this difference anomalous map. Should such a list include only the >> sites you can be sure of (say above 7 sigma) or all possible sites (say >> above 5 sigma) ? >> >> - Should I use the (peak height)/(sigma of map) as the first estimate of >> occupancy ? if so do I have to normalize it (strongest site has occupancy >> of 1) or can I leave the height/sigma values as they are ? >> >> >>Many thanks, >> >>Peter. >> > >_______________________________________________ >phenixbb mailing list >phenixbb(a)phenix-online.org >http://phenix-online.org/mailman/listinfo/phenixbb > >

1 0

Phaser SAD question
by Peter Grey 24 Dec '09

24 Dec '09

Dear Phenix Experts, I use the SAD+MR option in Phenix. The ASU is huge and the number of scatterres is in the few hundreds range. I try to minimize the number of runs of Phaser SAD+MR since they take long time on my computer and would appreciate your insight regarding the following : - I noticed that I get much higher signal and much larger number of sites if I use density modified phases to calculate the anomalous difference map. Is it possible to input these phases instead of the model to Phaser ? - If not, do you recommend providing Phaser with the list of sites found in this difference anomalous map. Should such a list include only the sites you can be sure of (say above 7 sigma) or all possible sites (say above 5 sigma) ? - Should I use the (peak height)/(sigma of map) as the first estimate of occupancy ? if so do I have to normalize it (strongest site has occupancy of 1) or can I leave the height/sigma values as they are ? Many thanks, Peter.

2 1

Re: [phenixbb] phenix GUI freezing during refinement and file naming issues
by John Berrisford 23 Dec '09

23 Dec '09

Dear all. When the GUI freezes on me its a complete freeze with no pop up dialog boxes. This means the job history is not recorded and no results are recorded in the GUI. I'm not sure if this is screensaver or desktop switching related - it only seems to happen randomly - mostly when running long jobs where nothing much happens on the screen, in my experience. I'm also using an nVidia card and an upto date version of the X-server (on my ubuntu 9.10 box). I agree this is a minor issue compared to having a great program which works well. But it does mean that if its easier to use a command line all your hard work with the GUI (and its enhancements) will go to waste... John -- John Berrisford Medical Research Council Mitochondrial Biology Unit Wellcome Trust / MRC Building Hills Road Cambridge CB2 0XY WEB: www.mrc-mbu.cam.ac.uk TEL: +44 (01223) 252918

2 1

rotational pseudosymmetry in P61?
by Gabriela Guédez R. 23 Dec '09

23 Dec '09

Hi all, how could i determine the r.m.s.d. and matrix of the NCS operator that relates two molecules in the asymmetric unit in P61? i would like to know whether there is rotational pseudosymmetry in the crystal because i see peaks for two molecules in the kappa=180º SR map. Thanks a lot!

3 2

Autosol problem: Failed to carry out Autosol_scale_and_analyze_mad
by Zheng Zhou 22 Dec '09

22 Dec '09

Dear all I installed phenix nightly builds, dev-256. I have solved structures using molecular replacement. Phenix.refine runs smoothly. Thanks for your great work. When I run some test data using autosol, I met the following error message: ******************************************************************************** Failed to carry out AutoSol_scale_and_analyze_mad: None of the solve versions worked See overall log file and also solve log file (in AutoSol_run_1_/LAST.LOG) ******************************************************************************** The LAST.LOG appears to be empty. I also run the phenix.solve command. Here is the output: ********************************************************************************************* $ phenix.solve binary=/usr/local/phenix-dev-254/build/intel-linux-2.6/solve_resolve/exe/solve SOLVEDIR=/usr/local/phenix-dev-254/solve_resolve/ext_ref_files SYMOP=/usr/local/phenix-dev-254/solve_resolve/ext_ref_files/symop.lib SYMINFO=/usr/local/phenix-dev-254/solve_resolve/ext_ref_files/syminfo.lib CCP4_OPEN=UNKNOWN SOLVE_RESOLVE_MALLOC_MAX=7.48024e+09 /usr/local/phenix-dev-254/build/intel-linux-2.6/solve_resolve/exe/solve: error while loading shared libraries: /usr/local/phenix-dev-254/build/intel-linux-2.6/lib/libimf.so: cannot restore segment prot after reloc: Permission denied ***************************************************************************************************** The libimf.so file is there. I run phenix.solve as superuser, but I still have the same problem. Thanks for your help. Best, Joe

2 3

phenix GUI freezing during refinement and file naming issues
by John Berrisford 22 Dec '09

22 Dec '09

Dear all I've recently been having a few problems with the new phenix GUI freezing during refinement runs on several different systems, most often when I try and do something different on the computer which puts the GUI in the background or on a different desktop. I've seen this problem when running a local copy of phenix on a ubuntu 9.10 system or when running the GUI remotely on a mac connected to our Centos linux cluster. When I switch to a different desktop to get on with doing something else (refine jobs can take a long time) the GUI freezes - although the underlying task runs correctly. The only way out of this problem is to force quit the GUI and then restart it. However, this means that the task isn't recorded in my history and all settings used to run the job are lost from the GUI - the parameter file to run the job does exist though. This also means I lose all the nice molprobibity connections to coot which are really useful. Also if I try to run phenix on the Centos linux cluster when running a local windows system the GUI is so sllloooowwww whenever the GUI tries to open a new window that its unsuable. Has anyone else seen these problems and know a cure? If I run phenix locally on a mac I don't seem to have this problem, but it is nice to run the refine job on our linux cluster so that I can get on with doing something else without it slowing down my mac. On another note, when running phenix refine from the GUI it nicely puts all the files in a folder called "Refine_1" etc... with the number incrementing with each run. So, why are the pdb and mtz file always called "refine_1..." not "refine_2 etc... to match the folder name? This would help with keeping track of files during itterative refinement and model building. It would also stop the problem of having lots of files called "refine_1" loaded into coot when I try to compare models from different refine runs. It would also be useful to be able to be able to add a prefix to the file name / folder name / job name in the GUI so that I could distinguish between different refinement experiments. This is easy to do on the command line.... Thanks and merry christmas John -- John Berrisford Medical Research Council Mitochondrial Biology Unit Wellcome Trust / MRC Building Hills Road Cambridge CB2 0XY WEB: www.mrc-mbu.cam.ac.uk TEL: +44 (01223) 252918

3 4

Re: [phenixbb] ncs resolve
by Thomas C. Terwilliger 18 Dec '09

18 Dec '09

Hi Georg, I think that if phenix.find_ncs says No NCS found...you might try ncs_in_type=sites then it is interpreting the sites as parts of chains and is trying to find NCS in those chains, while what you want is for it to interpret these as sites and to find ncs in the sites. So I wonder if you did try adding the keyword ncs_in_type as in: phenix.find_ncs my_sites.pdb my_map_coeffs.mtz ncs_in_type=sites If that failed, then I would suggest trying fewer sites that are more separated: perhaps 3 or 4 sites per monomer, as far apart as possible. If that fails (or maybe as your first try actually), take the entire helix, and make 2 copies A and B, placing them in their appropriate places, and then feed these to phenix with: phenix.simple_ncs_from_pdb my_two_copies.pdb and that may do it for you. All the best, Tom T >> Dear phenix users, >> >> luckily autobuild could solve my SAD problem, but didn´t find any NCS. >> Although the protein is a dimer. I analysed the structure in coot now and >> SSM superposed a homolog protein (just one chain) which (fits quite good) >> on >> chain A and B. I then could find a helix which was built in chain A and B, >> I >> took the 8 corresponding aminoacids from A and B and run: phenix.find_ncs >> overall_best_denmod_map_coeffs.mtz overall_best198206ca155163b.pdb with >> and >> without ncs_parameters.force_ncs=True and didn´t get any solution. >> >> The program gives me an output saying: No NCS found...you might try >> ncs_in_type=sites >> >> Any suggestions whats wrong and a smarter way to do this. >> >> Thanks in advance Georg. >> >> -----Ursprüngliche Nachricht----- >> Von: phenixbb-bounces(a)phenix-online.org >> [mailto:[email protected]] Im Auftrag von Thomas Stout >> Gesendet: Freitag, 18. Dezember 2009 23:56 >> An: PHENIX user mailing list >> Betreff: Re: [phenixbb] queueing systems >> >> >> >> We don't use PHENIX, but we do use several queueing systems on three >> different clusters. >> >> SGE is bulletproof and highly recommended. >> >> Condor has some degree of "evilness" -- especially under Windows. I have >> found it prone to grabbing resources indiscriminately and also being very, >> very resistant to uninstallation. Your mileage may vary. >> >> Certain implementations of "torque" are broken - stay away, or at least be >> very cautious and enter with your eyes open. >> >> Cheers, >> Tom >> >> >> -----Original Message----- >> From: phenixbb-bounces(a)phenix-online.org on behalf of Phil Evans >> Sent: Fri 12/18/2009 2:47 PM >> To: PHENIX user mailing list >> Subject: Re: [phenixbb] queueing systems >> >> >> Dear Nat >> >> We use Sun grid engine >> >> Phil >> >> On 18 Dec 2009, at 18:47, Nathaniel Echols wrote: >> >>> Hi all-- >>> >>> I'd like to hear what kinds of queueing systems people are using for >>> running Phenix on managed clusters. The eventual goal - if feasible - >>> is to add support for submitting jobs directly from the Phenix GUI. >>> Anyone who is interested, please send me (not the list) a brief email >>> describing a) the OS, b) the queueing software, and c) the submission >>> method (i.e. shell script, configuration file, simple command). >>> (Also, anyone who is willing to help test this if and when it becomes >>> available, let me know - I have access to Sun Grid Engine and Condor >>> already, but it's probably going to be difficult to test anything else >>> personally.) >>> >>> thanks, >>> Nat >>> >>> -------------------- >>> Nathaniel Echols >>> Lawrence Berkeley Lab >>> 510-486-5136 >>> NEchols(a)lbl.gov >>> >>> >>> >>> >>> >>> _______________________________________________ >>> phenixbb mailing list >>> phenixbb(a)phenix-online.org >>> http://phenix-online.org/mailman/listinfo/phenixbb >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://phenix-online.org/mailman/listinfo/phenixbb >> >> This email (including any attachments) may contain material that is >> confidential and privileged and is for the sole use of the intended >> recipient. Any review, reliance or distribution by others or forwarding >> without express permission is strictly prohibited. If you are not the >> intended recipient, please contact the sender and delete all copies. >> Exelixis, Inc. reserves the right, to the extent and under circumstances >> permitted by applicable law, to retain, monitor and intercept e-mail >> messages to and from its systems. >> >> >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://phenix-online.org/mailman/listinfo/phenixbb >>

1 0

queueing systems
by Nathaniel Echols 18 Dec '09

18 Dec '09

Hi all-- I'd like to hear what kinds of queueing systems people are using for running Phenix on managed clusters. The eventual goal - if feasible - is to add support for submitting jobs directly from the Phenix GUI. Anyone who is interested, please send me (not the list) a brief email describing a) the OS, b) the queueing software, and c) the submission method (i.e. shell script, configuration file, simple command). (Also, anyone who is willing to help test this if and when it becomes available, let me know - I have access to Sun Grid Engine and Condor already, but it's probably going to be difficult to test anything else personally.) thanks, Nat -------------------- Nathaniel Echols Lawrence Berkeley Lab 510-486-5136 NEchols(a)lbl.gov

4 3

refinement in low resolution
by r n 17 Dec '09

17 Dec '09

Hi I am trying to refine the model in low resolution and there are multimer in the asu. I did rigid body refinement after MR , molecule did not move much and R/Rfree stays 45 % and If I do the overlap between the monomers, rms is 0.06. If I do simulated annealing or any type of refinement (group_adp), Rfree is high and R goes down, but the if do the superposition of the monomer, rms is 1.5 to 1.8. So like to know what kind of refinement is can be done in very low resolution between 5 to 6 ang data? Is it any strict ncs could be possible in phenix that should not move the monomer too much? suggestions are appreciated rn

2 3

B(sol) too high
by Matthew Merski 17 Dec '09

17 Dec '09

I am in the process of solving a series of protein/ligand complex structures and ran the Polygon tool on the series. I noticed that, in general, even when all the other statistics were good, the B(sol) term was a bit too high according to Polygon. (The structures with poorer statistics also had high B(sol) values as well.) Since this is a series of multiple structures of the same protein, should this be a concern for me? And if so, what are typical causes of a high B(sol) value? Matt Merski Shoichet Group UCSF

3 6

LigPlot
by Zhou, Tongqing (NIH/VRC) [E] 16 Dec '09

16 Dec '09

Dear Phenix Users, I am trying to get Roman Laskowski's Ligplot, but could not find it on the ftp site he provided in the website below: http://www.biochem.ucl.ac.uk/bsm/ligplot/ligplot.html Does anyone have any information on how and where to get it? Thanks, Tongqing Tongqing Zhou, Ph.D. Staff Scientist Structural Biology Section Vaccine Research Center, NIAID/NIH Building 40, Room 4607B 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. ******************************************************************

2 1

tls matrices etc.
by Bryan W. Lepore 16 Dec '09

16 Dec '09

does phenix compute the center of reaction pretty much like refmac does? i figure the matrices/center of reaction between refmac and phenix are identical - is that true? -bryan

2 1

Rfree refmac5 issue
by anon 16 Dec '09

16 Dec '09

Hi, has anyone had any experience using phenix with refmac5? I seem to be having issues with the Rfree flag. Phenix sets aside 8.1% of the reflection for Rfree. However when i use refmac for refinement, it reports that its using 91.4% of the reflections for Rfree........? thanks

3 2

polygon output (dev-256)
by Bryan W. Lepore 16 Dec '09

16 Dec '09

[ phenix dev 256 ] where - or should we - be looking for the polygon output .png (the one with the polygon)? its obvious in 1.5-2, so i think its dev-256 only... -bryan

2 3

.cif use/restraints
by Bryan W. Lepore 15 Dec '09

15 Dec '09

a ligand is distorting - as if it is ignoring the .cif and i already have the angles tighted up in the .def. what else should i check? what i did : generated some elbow output. the .cif was provided to phenix.refine on cmd line. the .pdb had the H's and CONNECT's removed, but the .cif retains them. the .pdb was incorporated into a single .pdb input file. -bryan

2 1

Firefox problem
by Phil Evans 15 Dec '09

15 Dec '09

Dear Gurus If I click the Help button on the Phenix GUI I get the following fail messages. What am I doing wrong? I'm running 64-bit Phenix on a 64-bit machine, but our Firefox (as I understand it) is 32-bit, since there is no 64-bit Java Best wishes Phil hex_pre> phenix Using browser /usr/bin/firefox (Gecko:3603): Gtk-WARNING **: /lmb/home/public/xtal/Phenix/Phenix64/phenix-dev-2 49/build/intel-linux-2.6-x86_64/base/lib/gtk-2.0/2.10.0/engines/libclearlooks.so : wrong ELF class: ELFCLASS64 (Gecko:3603): Pango-WARNING **: /lmb/home/public/xtal/Phenix/Phenix64/phenix-dev -249/build/intel-linux-2.6-x86_64/base/lib/pango/1.6.0/modules/pango-basic-fc.so : wrong ELF class: ELFCLASS64 Failed to load Pango module for id: 'BasicScriptEngineFc' (Gecko:3603): Pango-WARNING **: /lmb/home/public/xtal/Phenix/Phenix64/phenix-dev -249/build/intel-linux-2.6-x86_64/base/lib/pango/1.6.0/modules/pango-basic-fc.so : wrong ELF class: ELFCLASS64 Failed to load Pango module for id: 'BasicScriptEngineFc' (Gecko:3603): Pango-WARNING **: /lmb/home/public/xtal/Phenix/Phenix64/phenix-dev -249/build/intel-linux-2.6-x86_64/base/lib/pango/1.6.0/modules/pango-basic-fc.so : wrong ELF class: ELFCLASS64 Failed to load Pango module for id: 'BasicScriptEngineFc' (Gecko:3603): Pango-WARNING **: /lmb/home/public/xtal/Phenix/Phenix64/phenix-dev -249/build/intel-linux-2.6-x86_64/base/lib/pango/1.6.0/modules/pango-basic-fc.so : wrong ELF class: ELFCLASS64 .. snip ... Failed to load Pango module for id: 'BasicScriptEngineFc' (Gecko:3603): Pango-WARNING **: /lmb/home/public/xtal/Phenix/Phenix64/phenix-dev -249/build/intel-linux-2.6-x86_64/base/lib/pango/1.6.0/modules/pango-basic-fc.so : wrong ELF class: ELFCLASS64 Failed to load Pango module for id: 'BasicScriptEngineFc' (Gecko:3603): GLib-GObject-CRITICAL **: g_object_get_qdata: assertion `G_IS_OBJE CT (object)' failed (Gecko:3603): Pango-CRITICAL **: pango_font_describe: assertion `font != NULL' f ailed /usr/lib/firefox-3.0.15/run-mozilla.sh: line 131: 3603 Segmentation fault "$prog" ${1+"$@"}

2 1

phenix and refmac5
by anon 15 Dec '09

15 Dec '09

Hi, recently i used phenix (AutoSol and AutoBuild) on a SAD dataset containing Iodine. I was wanting to do my refinement using Refmac instead of phenix.refine. I was wondering where would i be able to locate the correct .mtz file to input into Refmac5 for rounds of refinement? I was under the impression the .mtz file used for refinement was the one produced after scaling, but im not sure where phenix puts this. Also I was wondering if anyone had any issues with phenix.xtriage carrying out its Matthews coefficient calculation? When using Phenix, it gave me a solution with 3 mols/AU, despite there being 4 mols/AU according to the Matthews coeffiecient calculation in CCP4. Looking back at the xtriage log in autosol it has n_residues= 200 (despite me not specifying anywhere the number of residues, my protein has 147 residues). Realising this I re-ran Autosol, inputting my .mtz file along with my sequence this time, and specifying the number of NCS copies as 4. Looking at the phenix.xtriage log this time results in it recognising the correct number of residues, however the matthews coefficient calculation is still incorrect, claiming 1 mol/Au. Any ideas what may be the problem? ##-------------------------------------------## ## Unit cell defined manually, will ignore ## specification in reflection file: ## From file : None ## From input: (140.124, 77.954, 54.519, 90, 98.549, 90) ##-------------------------------------------## ##-------------------------------------------## ## Space group defined manually, will ignore ## specification in reflection file: ## From file : C 1 2 1 ## From input: C 1 2 1 ##-------------------------------------------## Effective parameters: #phil __ON__ scaling { input { asu_contents { n_residues = 588 n_bases = None n_copies_per_asu = 4 } Matthews coefficient and Solvent content statistics ---------------------------------------------------------------- | Copies | Solvent content | Matthews Coef. | P(solvent cont.) | |--------|-----------------|----------------|------------------| | 1 | 0.448 | 2.226 | 1.000 | ---------------------------------------------------------------- | Best guess : 1 copies in the asu | ----------------------------------------------------------------

2 1

input/output error
by Geoffrey Feld 14 Dec '09

14 Dec '09

Greetings fellow crystallographers, This past weekend, I collected some new data on my new project, indexed, scaled it, and found a phaser solution. The resolution is 3.6 A and I need to find one of the protein components, so I'm slowly refining piece-by-piece so I don't overdo it and flatten out the missing protein. To do this, I simply do rigid body and tls refinements for decreasingly smaller pieces of the protein I've already found. This method worked for me in the past for a similar structure. However, this time, as the refinement is finishing up, I get an input/output error and the program shuts down. Sometimes I get similar errors, and I can fix them by simply running the job again, overwriting the failed job, but this time that doesn't seem to be working. Here's the output from my .log file: Traceback (most recent call last): File "/programs/phenix-1.5-2/phenix/phenix/command_line/refine.py", line 11, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/programs/phenix-1.5-2/phenix/phenix/refinement/command_line.py", line 89, in run call_back_handler=call_back_handler) File "/programs/phenix-1.5-2/phenix/phenix/refinement/driver.py", line 1121, in run call_back_handler = call_back_handler) File "/programs/phenix-1.5-2/phenix/phenix/refinement/strategies.py", line 428, in refinement_machine monitors = monitors) File "/programs/phenix-1.5-2/cctbx_project/mmtbx/refinement/rigid_body.py", line 559, in __init__ out = log) File "/programs/phenix-1.5-2/cctbx_project/mmtbx/refinement/rigid_body.py", line 619, in show translations=t_vec) File "/programs/phenix-1.5-2/cctbx_project/mmtbx/refinement/print_statistics.py", line 161, in show_rigid_body_rotations_and_translations % euler_angle_convention, suffix=frame) File "/programs/phenix-1.5-2/cctbx_project/libtbx/utils.py", line 652, in write file_object.write(str) IOError: [Errno 5] Input/output error Again, the only thing I've done is made a few rigid body and tls groups. Please advise! Thanks, Geoff -- Geoffrey K. Feld Department of Chemistry 492 Stanley Hall University of California, Berkeley "Vigilia pretium libertatis"

2 1

TLS-derived anisoU restraint
by Bryan W. Lepore 14 Dec '09

14 Dec '09

is wxu the easiest - perhaps indirect - way to restrain anisoU's derived from TLS? under the target_weights scope i see : wxu_scale fix_wxu optimize_wxu <-already using this .... b/c i want to pull back some anisoU deviations across C-N. -bryan

2 1

Re: [phenixbb] How to extend R-free set to higher resolution by new phenix-1.5.2
by Thomas C. Terwilliger 14 Dec '09

14 Dec '09

Hi PC, >> >> Hi, >> >> I know this topic was discussed before, but I am still confused and still >> unable to refine my map merged from hires data because phenix.refine >> complains "No array of R-free flags found" >> >> The command was used to generate a hires map. >> phenix.autobuild data=resolve_1.mtz hires_file=hires.sca >> \seq_file=sequence.dat maps_only=True >> >> So, now I have "exptl_fobs_phases_freeR_flags.mtz", "cycle_best_1.mtz", >> "overall_best_denmod_map_coeffs.mtz"...several files in my autobuild >> folder, What this command has done is to take the phases in resolve_1.mtz and combined them with amplitudes from hires.sca and added free R flags to generate exptl_fobs_phases_freeR_flags.mtz . Additionally the file overall_best_denmod_map_coeffs.mtz now contains the above amplitudes along with some kind of density-modified phases. I would caution that this is probably not what you really wanted to do. Instead: Almost always you should use for "data=" a file that contains experimental phases only (not density-modified phases). So it would be better to use the command: phenix.autobuild data=phaser_1.mtz hires_file=hires.sca \ eq_file=sequence.dat maps_only=True which would instead density modify using hl coeffs from modify phaser_1.mtz, extending phase information to the data in hires.sca . Now you will have exptl_fobs_phases_freeR_flags.mtz with (a) experimental amplitudes from hires.sca (b) experimental phases and hl coefficients (where available) from phaser_1.mtz (c) free R flags to the full resolution of hires.sca and overall_best_denmod_map_coeffs.mtz with coefficients for a density-modified map at the resolution of hires.sca. >> >> How do I merge low and high resolution data by >> "phenix.reflection_file_editor"? What files (input and output arrays) >> should I choose? I think reflection_file_editor will only copy arrays (or columns) of data, and not merge several columns of data into one. Note that merging high and low resolution data requires (1) a scale factor between the datasets and (2) possibly a b-factor or anisotropic b-factor relating the two datasets. It cannot be properly done simply by merging reflection files. The reflection_flle_editor will add free R flags if they are missing or extend them for you if they are present only to low resolution, however. This allows you to take a low-res dataset with free R flags and a high-res dataset with no free R and to transfer the existing low-res free R flags to the high-res dataset and to extend them to full resolution. If you want, you can use phenix.solve to scale low and high resolution data by pretending that they are part of a mad dataset: make a file "scale.csh" like this, after editing to match your case: #!/bin/csh -f cp $PHENIX/solve_resolve/ext_ref_files/c2.sym . phenix.solve<<EOD resolution 20 3.0 cell 76 28 42 90 103 90 symfile c2.sym readdenzo ! readformatted/readdenzo/readtrek/readccp4_unmerged premerged ! premerged/ unmerged read_intensities ! read_intensities/read_amplitudes lambda 1 mad_atom c fprprv_mad 1. fprimv_mad 1. rawmadfile low.sca rawmadfile high.sca skip_solve sad EOD Then run with chmod +x scale.csh ./scale.csh Now solve_1.sca should have your scaled, merged data. I hope that helps! -Tom T

3 2

How to extend R-free set to higher resolution by new phenix-1.5.2
by Pei-Chun Lin 12 Dec '09

12 Dec '09

Hi, I know this topic was discussed before, but I am still confused and still unable to refine my map merged from hires data because phenix.refine complains "No array of R-free flags found" The command was used to generate a hires map. phenix.autobuild data=resolve_1.mtz hires_file=hires.sca \seq_file=sequence.dat maps_only=True So, now I have "exptl_fobs_phases_freeR_flags.mtz", "cycle_best_1.mtz", "overall_best_denmod_map_coeffs.mtz"...several files in my autobuild folder, How do I merge low and high resolution data by "phenix.reflection_file_editor"? What files (input and output arrays) should I choose? Detail description will be very appreciate. Thanks for answer my question. PC _________________________________________________________________ Windows Live: Make it easier for your friends to see what you’re up to on Facebook. http://www.microsoft.com/middleeast/windows/windowslive/see-it-in-action/so…

1 0

Side chain at special position
by Nate Nicely 11 Dec '09

11 Dec '09

I have a structure wherein the data show alternate conformations for a particular Trp side chain. The problem is that each of these conformations also occupies a special position: rotamer A clashes with rotamer A of the symmetry object, and similar for rotamer B. I'm attaching a screenshot wherein I tried to capture what's going on. The prerefined structure has a yellow carbon backbone, and you can see how refinement pushes rotamer A out of its rightful place in the map. I readjusted it and added the alternate conformation as seen in the green carbon backbone, but refinement after that just pushes apart both rotamers (not shown). I'd appreciate advice how to handle refining this structure in phenix. Thanks, Nate N

4 3

Re: [phenixbb] structural validation and deposit format
by Pavel Afonine 11 Dec '09

11 Dec '09

On 12/11/09 12:55 PM, rui wrote: > SO in the output pdb file, there are another set of lines start with > ANISOU, what are those, should I keep those? These are ADPs - Atomic Displacement Parameters, or B-factors. They are as important as atomic coordinates (almost) stored in ATOM records -:) > I turned on the optimze_wxc and optimize_wxu flags in the GUI, it does > give better R/Rfree. Do I need to turn on every time or just use the > optimized weights, where I can find those optimized weight? In log file? Optimization of weights is usually time consuming and the default weights are just good most of the time. So I usually recommend running weights optimization as the last step to tune up your model. Of course there might be exceptions. Since optimal weight may be different from run to run, they are not stored in .def file for the next runs. Pavel.

1 0

structural validation and deposit format
by rui 11 Dec '09

11 Dec '09

Hi, All, Sorry for this stupid question. When I tried to use phenix to validate the structure, it generate geometry outliers, model validation, real-space correlations and a phenix_map_coeffs.mtz file which I assume could be open in coot. How should I use this generated mtz file? When I open in coot, it tellls me core dump. I'm wondering if phenix will automatically do the corrections or should I look one by one in the log file and modify? It also generate a test_reduced.pdb, should I use this reduced pdb ( H added) for next round of refinement? Thanks. Is the format from phenix.refine ready for deposit?

3 2

refinement
by r n 11 Dec '09

11 Dec '09

Hi all I do have very low resolution data (around 4 ang), what are the efficient way of doing refinement, either rigid body alone or rigidbody and tls or individual with group_adp. I did both, but not much significant changes in Rfree stays around 48/50%. Also I do have other questions 1. While doing rigid body refinement, phenix complaint about the special position and could not perform rigid body refinement, whereas individual site refinement is working fine. I do have to delete the atom in special position for rigid body refinement? Thanks

8 20

How to reload a project into the gui
by Dirk 11 Dec '09

11 Dec '09

Hi, during a phenix refinement the gui froze and I had to kill the process. After restarting the gui the project I was working with was disappeared - all other projects are still there. Is there a way to reload the project into the gui? All files including the project folder still exists but are somehow not recognized by the gui. I'm using dev-249 64bit. Greetings Dirk

2 1