- phenixbb - phenix-online.org

some general refine questions
by William Scott 05 Jul '07

05 Jul '07

Hi folks: With the helpful advice of the phenix crew, I've re-refined 2goz, an entry in the pdb in which I used refmac5 the first time around. There was a small error in the pdb: One nucleotide (C8 in chain A) was rotated 180 degrees about the glycosylic bond from where it should have been. The net result is one oxygen atom was in the wrong place in the structure, and this was apparent within the difference Fourier. After I rotated the base and re-refined in refmac, the Rfactor and other stats remained essentially unchanged (so this fix has nothing to do with the following): My refinement results in Refmac: REMARK 3 R VALUE (WORKING SET) : 0.191 REMARK 3 FREE R VALUE : 0.241 REMARK 3 RMS DEVIATIONS FROM IDEAL VALUES COUNT RMS WEIGHT REMARK 3 BOND LENGTHS REFINED ATOMS (A): 1503 ; 0.009 ; 0.021 REMARK 3 BOND ANGLES REFINED ATOMS (DEGREES): 2331 ; 1.676 ; 3.000 Then I refined in phenix. I get REMARK Final: r_work = 0.1737 r_free = 0.2224 bonds = 0.003 angles = 0.745 So the working and free R factor are both lower, and the geometry is substantially tighter. If I relax the weights in phenix, even to get rms bonds of around 0.006, the free-R goes up, which tells me this is not something I am justified in doing. Phil Evans once told me that anything tighter than about 0.01 (conservatively) was science-fiction, in that this is within the error of the idealized libraries for bond lengths. Other people have asked me about this too, and I didn't know how to reply. Is bonds = 0.003 angles = 0.745 realistic, or have I done something wrong? Even if I do structure idealization in Refmac, I can't get the rms bonds below about 0.005 (and my understanding is that it uses the same monomer library). Why is the weight multiplied by 0.5 by default? Since changing that to 1.0 slightly increases freeR, I assume it is the default for some reason. All the best, Bill PS: Could someone point me to a quick description for depositing phenix-refined coordinates? (The latest phenix installer seems to have omitted the docs.) William G. Scott contact info: http://chemistry.ucsc.edu/~wgscott

3 2

Re: [phenixbb] sugar links in phenix.refine
by Ralf W. Grosse-Kunstleve 04 Jul '07

04 Jul '07

Almost, but you need a separate scope (thing between braces) for each link, e.g.: refinement.pdb_interpretation.apply_cif_link { data_link = NAG-ASN residue_selection_1 = chain A and resname NAG and resid 900 residue_selection_2 = chain A and resname ASN and resid 329 } refinement.pdb_interpretation.apply_cif_link { data_link = NAG-NAG-b-D residue_selection_1 = chain A and resname NAG and resid 901 residue_selection_2 = chain A and resname ASN and resid 900 } Or alternatively: refinement.pdb_interpretation { apply_cif_link { data_link = NAG-ASN residue_selection_1 = chain A and resname NAG and resid 900 residue_selection_2 = chain A and resname ASN and resid 329 } apply_cif_link { data_link = NAG-NAG-b-D residue_selection_1 = chain A and resname NAG and resid 901 residue_selection_2 = chain A and resname ASN and resid 900 } } Ralf

1 0

Re: [phenixbb] Multivariate Z score L-test
by Claudia Scotti 29 Jun '07

29 Jun '07

Thank you so much. I'll try all the space groups and see what happens. Many thanks again, Claudia Claudia Scotti Dipartimento di Medicina Sperimentale Sezione di Patologia Generale Universita' di Pavia Piazza Botta, 10 27100 Pavia Italia Tel. 0039 0382 986335/8/1 Facs 0039 0382 303673 ---------------------------------------- > Date: Fri, 29 Jun 2007 07:35:22 -0700 > From: PHZwart(a)lbl.gov > To: phenixbb(a)phenix-online.org > Subject: Re: [phenixbb] Multivariate Z score L-test > > Hi Claudia, > > This could indicate that tyour data is actually P 2 1 1, P 1 2 1 or P > 1 1 2 or P 21 1 1, P 1 21 1, P 1 1 21. Try processing that data in > those spacegroups and see if that makes any difference. You might end > up havig to refine your structure in all spacegroups to find out which > possibility is the best. > > It could also indicate that the quality of your data is not very good, > either due to suboptimal data processing or because of bad diffraction > (or both). > > HTH > > Peter > > > 2007/6/29, Claudia Scotti : > > > > Dear List, > > > > I'm only at my third structure (a scFv antibody fragment) and I'm disturbing about a dataset that is giving problems at the refinement stage: it gets stuck at about Rcryst= 0.36 and Rfree=0.44. > > > > I've tried data processing with several softwares (Mosflm, XDS, Denzo) and the result is consistently P222, unit cell 69x79x90, 90, 90, 90. The data are anisotropic, with limits at 2.3 and 2.8 A. I cut the data at 3.0 A. > > Matthews coefficient suggests two molecules per asu. > > No NCS is detected by self-rotation and a native Patterson detects a small peak (10% of origin peak) at 0.19, 0.34, 0.5. > > MR using several softwares (AmoRe, Molrep and Phaser) and several models all agree that the best space group is P21212, though the systematic absences are not perfectly consistent. > > > > When refining, refinement gets stuck as I mentioned above. > > > > I've run Phenix x-triage and the results include this line (see below): "It could be worthwhile considering reprocessing the data". I was wondering if the Multivariate Z score L-test: 6.305 might suggest how to reprocess the data. > > > > Sorry for the probably trivial question and for probably missing important data. > > > > Thanks a lot for any help and suggestions, > > > > Claudia > > > > ------------------------------------------------------------------------------- > > Twinning and intensity statistics summary (acentric data): > > > > Statistics independent of twin laws > > - /^2 : 1.827 > > - ^2/ : 0.838 > > - : 0.649 > > - , : 0.417, 0.242 > > Multivariate Z score L-test: 6.305 > > The multivariate Z score is a quality measure of the given > > spread in intensities. Good to reasonable data is expected > > to have a Z score lower than 3.5. > > Large values can indicate twinning, but small values do not > > neccesarily exclude it. > > > > > > No (pseudo)merohedral twin laws were found. > > > > > > Patterson analyses > > - Largest peak height : 11.094 > > (correpsonding p value : 1.080e-01) > > > > > > The largest off-origin peak in the Patterson function is 11.09% of the > > height of the origin peak. No significant pseudotranslation is detected. > > > > The results of the L-test indicate that the intensity statistics > > are significantly different then is expected from good to reasonable, > > untwinned data. > > As there are no twin laws possible given the crystal symmetry, there could be > > a number of reasons for the departure of the intensity statistics from normality. > > Overmerging pseudo-symmetric or twinned data, intensity to amplitude conversion problems > > as well as bad data quality might be possible reasons. > > It could be worthwhile considering reprocessing the data. > > ----------------------------------------------------------------------------------------------- > > > > Claudia Scotti > > Dipartimento di Medicina Sperimentale > > Sezione di Patologia Generale > > Universita' di Pavia > > Piazza Botta, 10 > > 27100 Pavia > > Italia > > Tel. 0039 0382 986335/8/1 > > Facs 0039 0382 303673 > > > > _________________________________________________________________ > > Express yourself instantly with MSN Messenger! Download today it's FREE! > > http://messenger.msn.click-url.com/go/onm00200471ave/direct/01/ > > _______________________________________________ > > phenixbb mailing list > > phenixbb(a)phenix-online.org > > http://www.phenix-online.org/mailman/listinfo/phenixbb > > > _______________________________________________ > phenixbb mailing list > phenixbb(a)phenix-online.org > http://www.phenix-online.org/mailman/listinfo/phenixbb _________________________________________________________________ With Windows Live Hotmail, you can personalize your inbox with your favorite color. www.windowslive-hotmail.com/learnmore/personalize.html?locale=en-us&ocid=TX…

2 1

Multivariate Z score L-test
by Claudia Scotti 29 Jun '07

29 Jun '07

Dear List, I'm only at my third structure (a scFv antibody fragment) and I'm disturbing about a dataset that is giving problems at the refinement stage: it gets stuck at about Rcryst= 0.36 and Rfree=0.44. I've tried data processing with several softwares (Mosflm, XDS, Denzo) and the result is consistently P222, unit cell 69x79x90, 90, 90, 90. The data are anisotropic, with limits at 2.3 and 2.8 A. I cut the data at 3.0 A. Matthews coefficient suggests two molecules per asu. No NCS is detected by self-rotation and a native Patterson detects a small peak (10% of origin peak) at 0.19, 0.34, 0.5. MR using several softwares (AmoRe, Molrep and Phaser) and several models all agree that the best space group is P21212, though the systematic absences are not perfectly consistent. When refining, refinement gets stuck as I mentioned above. I've run Phenix x-triage and the results include this line (see below): "It could be worthwhile considering reprocessing the data". I was wondering if the Multivariate Z score L-test: 6.305 might suggest how to reprocess the data. Sorry for the probably trivial question and for probably missing important data. Thanks a lot for any help and suggestions, Claudia ------------------------------------------------------------------------------- Twinning and intensity statistics summary (acentric data): Statistics independent of twin laws - /^2 : 1.827 - ^2/ : 0.838 - : 0.649 - , : 0.417, 0.242 Multivariate Z score L-test: 6.305 The multivariate Z score is a quality measure of the given spread in intensities. Good to reasonable data is expected to have a Z score lower than 3.5. Large values can indicate twinning, but small values do not neccesarily exclude it. No (pseudo)merohedral twin laws were found. Patterson analyses - Largest peak height : 11.094 (correpsonding p value : 1.080e-01) The largest off-origin peak in the Patterson function is 11.09% of the height of the origin peak. No significant pseudotranslation is detected. The results of the L-test indicate that the intensity statistics are significantly different then is expected from good to reasonable, untwinned data. As there are no twin laws possible given the crystal symmetry, there could be a number of reasons for the departure of the intensity statistics from normality. Overmerging pseudo-symmetric or twinned data, intensity to amplitude conversion problems as well as bad data quality might be possible reasons. It could be worthwhile considering reprocessing the data. ----------------------------------------------------------------------------------------------- Claudia Scotti Dipartimento di Medicina Sperimentale Sezione di Patologia Generale Universita' di Pavia Piazza Botta, 10 27100 Pavia Italia Tel. 0039 0382 986335/8/1 Facs 0039 0382 303673 _________________________________________________________________ Express yourself instantly with MSN Messenger! Download today it's FREE! http://messenger.msn.click-url.com/go/onm00200471ave/direct/01/

2 1

sugar links in phenix.refine
by Jianghai Zhu 21 Jun '07

21 Jun '07

Hi, How does phenix.refine handle the sugar links in refinement? I suppose phenix.refine will automatically recognize NAG and MAN. Will it recognized NAS-ASN, beta1-4, alpha1-4, etc linkages between sugars? Thanks. Jianghai +++++++++++++++++++++++++++++++ Jianghai Zhu, Ph.D CBR Institute for Biomedical Research Department of Pathology Harvard Medical School 200 Longwood Ave., Boston, MA 02115 Ph: 617-278-3211 Fx: 618-278-3232 +++++++++++++++++++++++++++++++

2 1

'Tasks' tab missing from 1.3b GUI and HL coefficients from a model
by Stephen Graham 21 Jun '07

21 Jun '07

Hi all, I recently (well, yesterday) installed the latest beta of phenix and everything is working great from the command line (which is how I use phenix 99.9% of the time). However, I fired up the GUI and found that the 'Tasks' tab appears to be missing from the 'Strategy Interface'. Is this a bug, a manifestation of an installation problem, or a feature? And... I was trying to use the GUI for prime-and-switch phasing to improve the maps for a partly-refined MR solution. What is the recommended path for creating HL coefficients from a model (from simulated annealing in phenix.refine) so that I can use the GUI's resolve 'prime-and-switch' strategy? Or is there an even better way of removing model bias and density-improving a big, low-res, MR-solved structure with 4-fold NCS in phenix? Thanks, Stephen -- Dr Stephen Graham Nuffield Medical Fellow Division of Structural Biology Wellcome Trust Centre for Human Genetics Roosevelt Drive Oxford OX3 7BN United Kingdom Phone: +44 1865 287 549

5 5

harmonic restrains in phenix.refine
by Bjørn Pañella Pedersen 21 Jun '07

21 Jun '07

Hi all, Is it possible to use harmonic restrains in phenix.refine, in a manner similar to CNS v1.2? I have searched the phenix.refine documentation and the web, but can't find anything. I use version 2007_04_06_1210 of phenix. Thanks! -Bjørn -- Bjørn Pañella Pedersen Ph.D Student, MSc Department of Molecular Biology Office: +45 89425021 University of Aarhus Lab: +45 89425010 Gustav Wieds Vej 10c Email: bpp(a)mb.au.dk DK-8000 Aarhus C Lab WWW: http://www.bioxray.dk Denmark

2 1

SAD + MR
by Lucas Bleicher 19 Jun '07

19 Jun '07

I have a dataset from a Cs-soaked crystal, and the protein has a domain solved in the PDB. The problem is that the map obtained just with the experimental phasing is not so easily interpreted, and the domain is just about 25% of the whole protein. Inspection of the map generated by MR suggests that it's a good solution, with residues fitting reasonably well on the density. Given the new features on Phaser 2.0, which I suppose is included in Phenix 1.3b, can anyone suggest a protocol to profit from both things (anomalous scattering + good search model)? Thanks in advance, Lucas ____________________________________________________________________________________ Novo Yahoo! Cadê? - Experimente uma nova busca. http://yahoo.com.br/oqueeuganhocomisso

3 2

Death by Big-ass structure&[email protected]
by Kate Kavanagh 13 Jun '07

13 Jun '07

Hi Peter, The reindexing works fine but when I try to convert back to R3:H all atoms with ANISOU records are given B-factors of -79.86: ATOM 1 N MET A 1 16.402 60.069 44.768 1.00-78.96 N ATOM 2 CA MET A 1 16.316 58.621 44.978 1.00-78.96 C ATOM 3 CB MET A 1 17.622 57.922 44.573 1.00-78.96 C ATOM 4 C MET A 1 15.966 58.300 46.432 1.00-78.96 C ATOM 5 O MET A 1 16.196 59.115 47.326 1.00-78.96 O ATOM 6 N PHE A 2 15.408 57.114 46.661 1.00-78.96 N ATOM 7 CA PHE A 2 15.047 56.686 48.013 1.00-78.96 C ATOM 8 CB PHE A 2 13.750 55.866 48.016 1.00-78.96 C ATOM 9 CG PHE A 2 13.323 55.436 49.391 1.00-78.96 C ..... And since TLS refinement has been used, all protein atoms have ANISOU records. Any suggestions? Kate

3 2

Re: [phenixbb] [ccp4bb] DANO from PDB
by Peter Zwart 08 Jun '07

08 Jun '07

Hi Santosh, Here is an alternative to methods mentioned earlier. If you get the latest version of the cctbx or cci_apps (from www.phenix-online.org) you will find a new application named phenix.xmanip (or mmtbx.xmanip if you only download the cctbx). xmanip has some useful features, one of them is that one can provide pieces of python code to the interface that will be executed. You have full access to all cctbx methods available. Please find below an example parameter file that perform structure factor calculation (in less than 90 lines of code and input) that includes elements with f" and f' not being equal to zero (you have to set the values of f" and f' yourself rather then rely on lookup tables by wavelength) and allows you to set bulk solvent parameters such as ksol, bsol and the overall anisotropic B of the data (B_cart). Save the text below (between the hash-es in a file named params.def , change the unit cell and sg parameters, the input pdb file and f" and f' value dictionary to suit your needs: ( like fp_dict = { "S": 0, "C": 0 } fdp_dict = { "S": 0, "C": 2.8 } for instance } You can fool around with solvent density, solvent B, overall scale factor and overall anisotropic B value and resolution (d_min). Then run phenix.xmanip params.def and get a file named xmanip.mtz that will contain your error free data. If you want, you can add uniform or Gaussian errors yourself. Since xmanip is fairly new, you will find a bug in the current release that somehow appears when you do not read in any external xray data. This happens however right after the mtz file in question has been written out. This will be fixed in the forthcoming release (soonish). Let me know if anything is unclear or if you have trouble running it. HTH Peter Zwart ################################################################### xmanip { input { unit_cell = "80.127 80.127 71.582 90.00 90.00 120.00" space_group = R3:H model { file_name = tst.pdb } } parameters { action = reindex manipulate_pdb *manipulate_miller manipulate_miller { task = get_dano get_diso lsq_scale sfcalc *custom None output_label_root = "Imock" custom{ code = """ #------------------------------------------------------- # This you can change if desired # # Dictionary for fp and fpp values fp_dict = { "Zn": -0.189 } fdp_dict = { "Zn": 4.898 } # final resolution of the data d_min = 2.0 # anisotropic B value b_cart = [ 0, 0, 0, 0, 0, 0 ] # solvent parameters k_sol = 0.35 b_sol = 45.0 # overall scale k_overall = 1.0 #------------------------------------------------------- # below this point, things should be more or less okai # from mmtbx import f_model from cctbx.xray import observation_types #loop over all scatterers in the xray structure #print >> out, dir(xray_structure) scatterers = xray_structure.scatterers() for atom in scatterers.as_1d(): if fp_dict.has_key( atom.element_symbol() ): atom.fp = fp_dict[ atom.element_symbol() ] if fdp_dict.has_key( atom.element_symbol() ): atom.fdp = fdp_dict[ atom.element_symbol() ] #the changes have been made to fp and fdp #Now we have to do some trickery to get some fobs with bulk solvent #contribution ... #First make a 'fake' set of f_obs fake_f_obs = abs(xray_structure.structure_factors( d_min = d_min, anomalous_flag = True ).f_calc()) # make a 'fake' set of free flags fake_free_flags = fake_f_obs.generate_r_free_flags(fraction = 0.1, max_free = 99999999) fmodel = f_model.manager( xray_structure = xray_structure, r_free_flags = fake_free_flags, target_name = "ls_wunit_k1", f_obs = fake_f_obs, b_cart = b_cart, k_sol = k_sol, b_sol = b_sol, overall_scale = k_overall ) # now get the final values of F+ and F-. # Note the abs to get amplitudes only result = abs( fmodel.f_model() ).set_observation_type( observation_types.amplitude() ) # you want make this into intensities? result = result.f_as_f_sq() # lets leave the intensities alone for now, and just make 'mock' sigmas # to statisfy certain programs. # Don't forget to set the observation type! sigmas = result.data()/100.0 result = result.customized_copy(data = result.data(), sigmas=sigmas).set_observation_type(result) #By default, a parameter named result is passed back to the xmanip #routines and will be written to the final mtz file """ } } } output { logfile = "xmanip.log" hklout = "xmanip.mtz" xyzout = "xmanip.pdb" } } ################################################################### Santosh Panjikar wrote: > Hi all, > > Does anybody have a program which can calculate anomalous differences or > F+ and F- from a refined structure at given wavelength and resolution ? > > Thanks > Santosh > > > Santosh Panjikar, Ph.D. panjikar(a)embl-hamburg.de > Staff Scientist > EMBL Hamburg outstation http://www.embl-hamburg.de/~panjikar/ > DESY, Notkestrasse 85 Tel.:+49-40-89902-141 > 22603, Hamburg, Germany Fax: +49-40-89902-149

1 0

removing Van der Waals repulsion between two atoms
by Miguel Ortiz-Lombardía 06 Jun '07

06 Jun '07

Dear all, I've a question regarding phenix.refine. Is it possible in phenix.refine to remove repulsion between a couple of atoms to see to what distance they refine when unrestrained? I'm thinking in something like what in CNS/XPLOR one would do by setting: igroup interaction (selection1) (selection2) weights * 1 vdw 0. elec 0. end end I've been searching the web for a solution to this question, but haven't come up with anything. Sorry if I missed it! I would appreciate your help. Cheers, Miguel -- correo-e: ibdeno(a)gmail.com ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Je suis de la mauvaise herbe, Braves gens, braves gens, Je pousse en liberté Dans les jardins mal fréquentés! Georges Brassens

2 1

Re: [phenixbb] how to add phase info to the autosol wizard
by Thomas C. Terwilliger 06 Jun '07

06 Jun '07

Hi again Chuanbing, I'm sorry, I mean 1.3b-d5 (d4 is just about to be released). -Tom T Hi Chuanbing, It will be available in version 1.3b-d4 which has yet to be released (or even built)...it takes a little time between the time I have put the code in and the time that it gets to you, unfortunately! -Tom T > Dear Tom, > > Thank you very much! It's expected! > I am using the Version 1.24.1b, is version 1.3 the next version you > mentioned? > > Chuanbing > > > > On 6/5/07, Tom Terwilliger <terwilliger(a)lanl.gov> wrote: >> >> Hi Chuanbing, >> >> I have added this to the AutoSol wizard >> (input_phase_file=my-phases.mtz) and the next version of PHENIX will have >> it. >> >> -Tom T >> >> At 01:51 AM 5/31/2007, CB Bian wrote: >> >Dear all, >> > >> > I am trying to add phase info from externally calculated >> > phases to the autosol wizard with S-SAD data. But I can not find the option to do it in the running. Can someone help me? >> > >> > >> >Thanks in advance >> > >> >Chuanbing >> > >> >-- >> >ChuanBing Bian, Ph.D Student >> >Structural biology group >> >State Key laboratory of Structural Chemistry >> >Fujian Institute of Research on the Structure of Matter, >> >Chinese Academy of Sciences >> >155 Yang Qiao Xi Lu, Fuzhou, Fujian 350002 P. R. China >> >E-mail: >> ><mailto:[email protected]>cbbian(a)gmail.com >> ><mailto:[email protected]> cbbian(a)fjirsm.ac.cn >> >Tel:(office)+86-591-8370-5814 >> ><http://protein.fjirsm.ac.cn/>http://protein.fjirsm.ac.cn >> >_______________________________________________ >> >phenixbb mailing list >> >phenixbb(a)phenix-online.org >> >http://www.phenix-online.org/mailman/listinfo/phenixbb >> >> >> >> >> Thomas C. Terwilliger >> Mail Stop M888 >> Los Alamos National Laboratory >> Los Alamos, NM 87545 >> >> Tel: 505-667-0072 email: terwilliger(a)LANL.gov >> Fax: 505-665-3024 SOLVE web site: http://solve.lanl.gov PHENIX web site: http:www.phenix-online.org >> ISFI Integrated Center for Structure and Function Innovation web site: http://techcenter.mbi.ucla.edu >> TB Structural Genomics Consortium web site: >> http://www.doe-mbi.ucla.edu/TB >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://www.phenix-online.org/mailman/listinfo/phenixbb >> >> >> >> >> Thomas C. Terwilliger >> Mail Stop M888 >> Los Alamos National Laboratory >> Los Alamos, NM 87545 >> >> Tel: 505-667-0072 email: terwilliger(a)LANL.gov >> Fax: 505-665-3024 SOLVE web site: http://solve.lanl.gov PHENIX web site: http:www.phenix-online.org >> ISFI Integrated Center for Structure and Function Innovation web site: http://techcenter.mbi.ucla.edu >> TB Structural Genomics Consortium web site: >> http://www.doe-mbi.ucla.edu/TB >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://www.phenix-online.org/mailman/listinfo/phenixbb >> > > > > -- > ChuanBing Bian, Ph.D Student > Structural biology group > State Key laboratory of Structural Chemistry > Fujian Institute of Research on the Structure of Matter, > Chinese Academy of Sciences > 155 Yang Qiao Xi Lu, Fuzhou, Fujian 350002 P. R. China > E-mail: cbbian(a)gmail.com cbbian(a)fjirsm.ac.cn > Tel:(office)+86-591-8370-5814 > _______________________________________________ > phenixbb mailing list > phenixbb(a)phenix-online.org > http://www.phenix-online.org/mailman/listinfo/phenixbb >

2 1

Re: [phenixbb] how to add phase info to the autosol wizard
by Tom Terwilliger 06 Jun '07

06 Jun '07

Hi Chuanbing, I have added this to the AutoSol wizard (input_phase_file=my-phases.mtz) and the next version of PHENIX will have it. -Tom T At 01:51 AM 5/31/2007, CB Bian wrote: >Dear all, > > I am trying to add phase info from externally calculated > phases to the autosol wizard with S-SAD data. But I can not find > the option to do it in the running. Can someone help me? > > >Thanks in advance > >Chuanbing > >-- >ChuanBing Bian, Ph.D Student >Structural biology group >State Key laboratory of Structural Chemistry >Fujian Institute of Research on the Structure of Matter, >Chinese Academy of Sciences >155 Yang Qiao Xi Lu, Fuzhou, Fujian 350002 P. R. China >E-mail: ><mailto:[email protected]>cbbian(a)gmail.com ><mailto:[email protected]>cbbian(a)fjirsm.ac.cn >Tel:(office)+86-591-8370-5814 ><http://protein.fjirsm.ac.cn/>http://protein.fjirsm.ac.cn >_______________________________________________ >phenixbb mailing list >phenixbb(a)phenix-online.org >http://www.phenix-online.org/mailman/listinfo/phenixbb Thomas C. Terwilliger Mail Stop M888 Los Alamos National Laboratory Los Alamos, NM 87545 Tel: 505-667-0072 email: terwilliger(a)LANL.gov Fax: 505-665-3024 SOLVE web site: http://solve.lanl.gov PHENIX web site: http:www.phenix-online.org ISFI Integrated Center for Structure and Function Innovation web site: http://techcenter.mbi.ucla.edu TB Structural Genomics Consortium web site: http://www.doe-mbi.ucla.edu/TB _______________________________________________ phenixbb mailing list phenixbb(a)phenix-online.org http://www.phenix-online.org/mailman/listinfo/phenixbb Thomas C. Terwilliger Mail Stop M888 Los Alamos National Laboratory Los Alamos, NM 87545 Tel: 505-667-0072 email: terwilliger(a)LANL.gov Fax: 505-665-3024 SOLVE web site: http://solve.lanl.gov PHENIX web site: http:www.phenix-online.org ISFI Integrated Center for Structure and Function Innovation web site: http://techcenter.mbi.ucla.edu TB Structural Genomics Consortium web site: http://www.doe-mbi.ucla.edu/TB

3 2

secondary structure restraint
by Jianghai Zhu 05 Jun '07

05 Jun '07

Hi all, I am refining a low resolution structure (~3.8 A). For some of the domains in the structure, we already have high resolution structures. Despite of some conformational changes in some loops, the core structures should be the same. I am wondering if it is possible to restrain the secondary structures in phenix.refine. Thanks. Jianghai +++++++++++++++++++++++++++++++ Jianghai Zhu, Ph.D CBR Institute for Biomedical Research Department of Pathology Harvard Medical School 200 Longwood Ave., Boston, MA 02115 Ph: 617-278-3211 Fx: 618-278-3232 +++++++++++++++++++++++++++++++

3 7

phenix.refine crashed
by Jianghai Zhu 05 Jun '07

05 Jun '07

Hi, In my refinement phenix.refine crashed with the following errors. ================================ xyz refinement =============================== Traceback (most recent call last): File "/nfs/home/jzhu/cci_apps_sources/phenix/phenix/command_line/ refine.py", line 5, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/nfs/home/jzhu/cci_apps_sources/phenix/phenix/refinement/ command_line.py", line 78, in run call_back_after_monitor_collect=call_back_after_monitor_collect) File "/nfs/home/jzhu/cci_apps_sources/phenix/phenix/refinement/ driver.py", line 1310, in run log = log) File "/nfs/home/jzhu/cci_apps_sources/phenix/phenix/refinement/ strategies.py", line 733, in refinement_machine h_params = h_params) File "/nfs/home/jzhu/cci_apps_sources/mmtbx/mmtbx/refinement/ minimization.py", line 136, in __init__ ignore_line_search_failed_step_at_lower_bound = True)) File "/nfs/home/jzhu/cci_apps_sources/scitbx/scitbx/lbfgs.py", line 242, in run log) File "/nfs/home/jzhu/cci_apps_sources/scitbx/scitbx/lbfgs.py", line 121, in run_c_plus_plus f, g = target_evaluator.compute_functional_and_gradients() File "/nfs/home/jzhu/cci_apps_sources/mmtbx/mmtbx/refinement/ minimization.py", line 369, in compute_functional_and_gradients u_iso_refinable_params = u_iso_refinable_params) File "/nfs/home/jzhu/cci_apps_sources/mmtbx/mmtbx/refinement/ minimization.py", line 224, in compute_target sf_v3d_sel = sf_v3d.set_selected(self.hd_selection, [0,0,0]) RuntimeError: scitbx Internal Error: /net/sharptail/scratch1/rwgk/ auto_build/sources/scitbx/include/scitbx/array_family/boost_python/ flex_wrapper.h(473): SCITBX_ASSERT(a.accessor() == flags.accessor()) failure. What happened here? Thanks. Jianghai +++++++++++++++++++++++++++++++ Jianghai Zhu, Ph.D CBR Institute for Biomedical Research Department of Pathology Harvard Medical School 200 Longwood Ave., Boston, MA 02115 Ph: 617-278-3211 Fx: 618-278-3232 +++++++++++++++++++++++++++++++

2 1

Re: [phenixbb] elbow.builder
by William Scott 04 Jun '07

04 Jun '07

Problem solved: The OS X pre-compiled distribution comes with an executable gamess.Mar242007R1.x but the shell script seems to require gamess.00.x so in /usr/local/gamess I did this: sudo ln -s gamess.Mar242007R1.x gamess.00.x and benzene is as one would expect. BTW, when I run this under tcsh, I get a more informative error message -- it says it can't find gamess. Under zsh, it just says "Failed to read coordinates. Aborting optimisation." Anyway, it all works now. Thanks. Bill On Jun 3, 2007, at 3:08 PM, Peter Zwart wrote: > gamess is installed right? > > P > > > > 2007/6/3, William Scott <wgscott(a)chemistry.ucsc.edu>: >> actually, "small distortions" was too optimistic: >> >> http://chemistry.ucsc.edu/~wgscott/temp/benzene.png >> >> >> >> >> On Jun 3, 2007, at 2:33 PM, William Scott wrote: >> >> > Hi folks: >> > >> > I've been playing around a bit with elbow.builder, and found the >> > following: >> > >> > -------------- >> > >> > >> > elbow.builder --opt --smiles c1ccccc1 >> > >> > gives me a nice, flat, hexagonal benzene >> > >> > elbow.builder --gamess --smiles c1ccccc1 >> > >> > gives me a benzene that is slightly bent out of plane and has other >> > small distortions with respect to 6-fold symmetry. >> > >> > Is this a problem with elbow, a problem with gamess, or a Jahn- >> > Teller effect? >> > >> > >> > ------------ >> > >> > >> > >> > Also, I noticed elbow.doc requires a mozilla or netscape >> > executable. Since both of these are obsolete, maybe a firefox >> > executable should be an option? On OS X, "open" should open a URL >> > in the user's favorite browser, but it doesn't seem to work right >> > with Safari (which stupidly opens a Finder window of the directory >> > containing the html file). My workaround was to make a shell >> > script called "mozilla" that invokes Camino: >> > >> > #!/bin/zsh -f >> > /usr/bin/open -a Camino "$@" >> > >> > HTH someone. >> > >> > ------------------ >> > >> > >> > Bill >> > >> > >> > >> > >> > William G. Scott >> > >> > contact info: http://chemistry.ucsc.edu/~wgscott >> > >> > >> > >> > >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://www.phenix-online.org/mailman/listinfo/phenixbb >>

2 2

elbow.builder
by William Scott 03 Jun '07

03 Jun '07

Hi folks: I've been playing around a bit with elbow.builder, and found the following: -------------- elbow.builder --opt --smiles c1ccccc1 gives me a nice, flat, hexagonal benzene elbow.builder --gamess --smiles c1ccccc1 gives me a benzene that is slightly bent out of plane and has other small distortions with respect to 6-fold symmetry. Is this a problem with elbow, a problem with gamess, or a Jahn-Teller effect? ------------ Also, I noticed elbow.doc requires a mozilla or netscape executable. Since both of these are obsolete, maybe a firefox executable should be an option? On OS X, "open" should open a URL in the user's favorite browser, but it doesn't seem to work right with Safari (which stupidly opens a Finder window of the directory containing the html file). My workaround was to make a shell script called "mozilla" that invokes Camino: #!/bin/zsh -f /usr/bin/open -a Camino "$@" HTH someone. ------------------ Bill William G. Scott contact info: http://chemistry.ucsc.edu/~wgscott

1 1

Binary file not executable during installation?
by lkemp＠mrc-lmb.cam.ac.uk 01 Jun '07

01 Jun '07

Hello BB I am trying to install Phenix on Suse 9.3. I went for the 2.6 version as i got this: 2.6.11.4-21.17-smp i686 when i typed uname. However it required libtermcap.so.2 so i switched to the 2.4 version as suggested on a previous BB query. This however also goes wrong: This time the binary.log file says that python in /common is not executable I checked and this was due to the need for libstdc++.so.6. I made a softlink directing this to the .5 version but it doesnt like this. Does anyone have any ideas??? Thank you in advance for any help Lauris

2 1

AssertionError & TLS
by Kanchan Anand 01 Jun '07

01 Jun '07

Hello everybody, I am a new user of Phenix and presently trying to refine 3.6 A data of a protein-RNA complex. whenever i add tls into my .def file for the next run, i get this problem: Traceback (most recent call last): File "/data/software/phenix-1.3b/phenix/phenix/command_line/refine.py", line 5 , in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/data/software/phenix-1.3b/phenix/phenix/refinement/command_line.py", li ne 78, in run call_back_after_monitor_collect=call_back_after_monitor_collect) File "/data/software/phenix-1.3b/phenix/phenix/refinement/driver.py", line 130 8, in run log = log) File "/data/software/phenix-1.3b/phenix/phenix/refinement/strategies.py", line 887, in refinement_machine h_params = h_params) File "/data/software/phenix-1.3b/mmtbx/mmtbx/refinement/adp_refinement.py", li ne 160, in __init__ tools.split_u(fmodel.xray_structure, tls_selections, offset) File "/data/software/phenix-1.3b/mmtbx/mmtbx/tls/tools.py", line 468, in split _u assert (u_iso < 0.0).count(True) == 0 AssertionError I have downloaded the new 1.3b version and without tls the program runs fine! could you please help Thanks in advance Regards kanchan

1 0

textal problem
by Daniele de Sanctis 31 May '07

31 May '07

Dear all, I'm running textal.build to build a model in a 3.2 A map, and I'm using the option --input_model and --sequence. It seems to run ok until finishing the lookup but then i got the error message: doing sequence alignment against: ../pippo.seq sh: Syntax error: Bad fd number error: [Errno 2] No such file or directory: 'phentx-lookup-correct.pdb' indeed the file is not there, but i guess that is not there because o the previous Syntax error. any clue? thanks d -- Daniele de Sanctis, PhD Homo sum humani nil alienum a me puto ______________________________________________________________ phone ++ 351 21 4469 662 fax ++ 351 21 4433 664 e-mail de.sanctisd(a)gmail.com, sanctis(a)itqb.unl.pt Mailing address: ITQB, Av. República, Apartato 127 2781-901 Oeiras Portugal

1 0

Re: [phenixbb] how to add phase info to the autosol wizard
by Tom Terwilliger 31 May '07

31 May '07

Hi Chuanbing, I'm sorry, that option isn't there yet! It is one that is on my list of to-do's. -Tom T At 01:51 AM 5/31/2007, CB Bian wrote: >Dear all, > > I am trying to add phase info from externally calculated > phases to the autosol wizard with S-SAD data. But I can not find > the option to do it in the running. Can someone help me? > > >Thanks in advance > >Chuanbing > >-- >ChuanBing Bian, Ph.D Student >Structural biology group >State Key laboratory of Structural Chemistry >Fujian Institute of Research on the Structure of Matter, >Chinese Academy of Sciences >155 Yang Qiao Xi Lu, Fuzhou, Fujian 350002 P. R. China >E-mail: ><mailto:[email protected]>cbbian(a)gmail.com ><mailto:[email protected]>cbbian(a)fjirsm.ac.cn >Tel:(office)+86-591-8370-5814 ><http://protein.fjirsm.ac.cn/>http://protein.fjirsm.ac.cn >_______________________________________________ >phenixbb mailing list >phenixbb(a)phenix-online.org >http://www.phenix-online.org/mailman/listinfo/phenixbb Thomas C. Terwilliger Mail Stop M888 Los Alamos National Laboratory Los Alamos, NM 87545 Tel: 505-667-0072 email: terwilliger(a)LANL.gov Fax: 505-665-3024 SOLVE web site: http://solve.lanl.gov PHENIX web site: http:www.phenix-online.org ISFI Integrated Center for Structure and Function Innovation web site: http://techcenter.mbi.ucla.edu TB Structural Genomics Consortium web site: http://www.doe-mbi.ucla.edu/TB

1 0

how to add phase info to the autosol wizard
by CB Bian 31 May '07

31 May '07

Dear all, I am trying to add phase info from externally calculated phases to the autosol wizard with S-SAD data. But I can not find the option to do it in the running. Can someone help me? Thanks in advance Chuanbing -- ChuanBing Bian, Ph.D Student Structural biology group State Key laboratory of Structural Chemistry Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences 155 Yang Qiao Xi Lu, Fuzhou, Fujian 350002 P. R. China E-mail: cbbian(a)gmail.com cbbian(a)fjirsm.ac.cn Tel:(office)+86-591-8370-5814 http://protein.fjirsm.ac.cn

1 0

Re: [phenixbb] Death by Big-ass structure
by Petrus H Zwart 30 May '07

30 May '07

Hi Ulrich, Frank, Ralf has just placed a new version of cci_apps on our download sites (www.phenix-online.org). It contains a new application named phenix.xmanip (or mmtbx.xmanip if you only download the cctbx). It is basically a rewritten version of the reindex.py script I once wrote as an example jiffy. xmanip can read in various reflection files or mtz/cns data columns which subsequently can be used in reindexing or other tasks. The parameter file I send Frank a couple of days ago should still work, even though some rearrangements have been made . Typing phenix.xmanip reindex.def should do the trick. Other options are available as well. One could - generate map coefficients from a model (plain FC,PHIC, or sigmaA style maps if observed data is available) - get anomalous / isomorphous differences - scale two datasets together - generate symmetry mates of a model - the 'custom' option allows you to pass in a snippet of python code and have it executed. One has access to all datasets and molecular model read in. You could for instance add two datasets or multiply them if it makes you happy. Full documentation is under way. If you have any question, please let me know. Peter

1 0

Re: [phenixbb] phenix.refine questions: weights, docs, etc.
by Ralf W. Grosse-Kunstleve 29 May '07

29 May '07

Hi Dale, > The sets of parameters that can be refined simultaneously is a result of the > choice of optimization method, not the choice of the function being optimized. > The crucial difference is how the second derivatives are handled. In the > methods used in Phenix and CNS the seconds derivatives of all parameters are > assumed to be equal and uncorrelated. This is true only for the first iteration of a minimization in a phenix macro cycle. After that the L-BFGS method builds up an approximation of the second derivative information. > To ensure that this assumption is true > only parameters of the same category can be varied in a single cycle. This > means that the coordinates can be varied, but the B factors and scale factors > have to be held fixed. When the B factors are varied the coordinates and > scale factors must be constant. When the scale factors are varied everything > else must be fixed. > > Other refinement programs use second derivatives in a more explicit way > than Phenix and do vary more types of parameters in a single cycle. Shelxd is > the most powerful and, by default, varies all parameters of all classes each > cycle. Refmac, to the best of my knowledge, refines both coordinates and ADPs > together, but does refine TLS parameters in a separate step. As far as I know, Shelx doesn't use second derivatives at all. Shelx is shaped completely by the least-squares approach, using the Jacobian. You couldn't plug another target function into Shelx as you can do in the common macromolecular programs. I.e. there is in fact a link between the minimization method and choice of the target function. REFMAC uses a block-diagonal approach that only takes the correlations of the parameters of one atom into account. phenix approximates these with the L-BFGS method. We have not actually tried to refine coordinates and B-factors simultaneously since we've been too busy with other priorities. In theory it should be possible, but I'm sure it will take a lot of fine-tuning. > Since all the parameters of our models are correlated with one another, it > is better to refine as many of them at once as possible. Implementing > the join refinement of all these kinds of parameters is difficult, so to save > programmers' time approximations are sometimes made. We'll spend the time when we find it. :-) But I don't expect spectacular gains, no matter how long we beat on the second derivatives. Ralf

2 1

phenix.refine questions: weights, docs, etc.
by Frank von Delft 29 May '07

29 May '07

Hi, four questions on phenix.refine: 1) What's the reason for refining bss, xyz and adp in separate cycles? CNS does that too while refmac (I believe) refines them simultaneously. I suspect I missed something, but I thought the power of ML refinement is that no refinable parameters are held constant if they're not supposed to be. (At least, wasn't that the big deal in Sharp when it came out?) 2) How does one figure the right restraint weights? (wxc and wxc). In refmac, I've always calibrated the weight with rmsd(bonds), tightening it until rmsd(bonds) is ~0.018 (or below for lower reso); invariably this is where Rfree will decrease furthest and smoothly as well (i.e. without that silly rise at the end.) With phenix.refine I'm finding that I only get Rfree to drop monotonously when I tighten the weights (both wxc and wxu) so much that rmsd(bonds) goes down to ~0.002. Would there be a different way to do this? 3) There isn't currently torsional refinement, is there? For high-temperature simulated annealing, I mean. I'm having the same weight-adjustment question for SA, except is harder to fiddle with because jobs take so long, but with torsional that becomes moot. Or should one always put very tight weights on for SA? 4) Is there good documentation for all the loads and loads of keywords? The names themselves are just descriptive enough to be interesting, but obviously not enough to be informative if I don't just want the black box. Nice program, though, got to say... phx.

3 2

Death by Big-ass structure
by Frank von Delft 25 May '07

25 May '07

Hi, we get the memory error below for a structure of 30603 atoms 273217 reflections 2.1A H3 193.868 193.868 352.223 90.00 90.00 120.00 10x NCS Anything I said? It's not THAT large a problem, I'd have thought, or should we go to the rhombohedral setting instead? (Hmmmm, this rings a bell, wasn't there a posting about this 3 days ago? Nothing in the archives, was it CCP4bb?) Cheers phx Traceback (most recent call last): File "/usr/local/phenix/phenix-1.26b/phenix/phenix/command_line/refine.py", line 5, in ? command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/usr/local/phenix/phenix-1.26b/phenix/phenix/refinement/command_line.py", line 75, in run refine_object.run( File "/usr/local/phenix/phenix-1.26b/phenix/phenix/refinement/driver.py", line 1236, in run log = log) File "/usr/local/phenix/phenix-1.26b/phenix/phenix/refinement/strategies.py", line 273, in refinement_machine log = log) File "/usr/local/phenix/phenix-1.26b/phenix/phenix/refinement/weight_xray_chem.py ", line 111, in __init__ self.fmodel_copy = self.fmodel.deep_copy() File "/usr/local/phenix/phenix-1.26b/mmtbx/mmtbx/twinning/twin_f_model.py", line 1003, in deep_copy map_types = self.map_types File "/usr/local/phenix/phenix-1.26b/mmtbx/mmtbx/twinning/twin_f_model.py", line 888, in __init__ self.update_f_mask() File "/usr/local/phenix/phenix-1.26b/mmtbx/mmtbx/twinning/twin_f_model.py", line 1370, in update_f_mask self.f_mask_array = mask.structure_factors( self.miller_set ) File "/usr/local/phenix/phenix-1.26b/mmtbx/mmtbx/masks.py", line 103, in structure_factors flex.grid(fft_manager.m_real()).set_focus(fft_manager.n_real())) MemoryError

4 6

Re: [phenixbb] Question about restraints in individualB-factor re finement
by Schubert, Carsten [PRDUS] 25 May '07

25 May '07

Thanks Pavel,Nigel: so it looks as if neither approach works and I would need to split the ligands into seperate PDB files for a truly independent B-factor refinement. Thanks for the hint about the non-bonded interactions. That was useful. Cheers Carsten -----Original Message----- From: phenixbb-bounces(a)phenix-online.org [mailto:[email protected]]On Behalf Of Pavel Afonine Sent: Friday, May 25, 2007 12:51 PM To: PHENIX user mailing list Subject: Re: [phenixbb] Question about restraints in individualB-factor refinement Hi Carsten, for the refinement of b-factors there is no difference between two options. We use the sphere b-factors restraints which are calculated based on interatomic distances and the information about alternative conformations is not used at that point. However if you refine the coordinates, the difference is significant: if you refine it as two separate residues then they will "see" each other for non-bonded interactions term calculation and hence will be pushed apart. Please let me know if you still have questions! Pavel. Schubert, Carsten [PRDUS] wrote: Hi, I have a question about the handling of restraints in the individual B-factor refinement routine. What I'd like to do is to refine a ligand, which can be present in 2 or more different conformations/orientations in its binding site. I'd like to use B-factor refinement on the various instances of the ligand, which one is the most relevant one, assuming that the most relevant conformation/orientation is associated with the lowest B-factor (Validity of that assumption set aside ...) My question is now is there any difference in the restraints applied to the b-factors in the scenarios where A) the ligand is modeled as alternate conformations i.e. ATOM 2724 C01AINH I 1 27.808 26.376 23.301 0.50 27.77 I C ATOM 2733 C01BINH I 1 30.898 22.496 17.340 0.50 22.15 I C ... vs. scenario B) where the same ligand is modeled as 2 different residues ATOM 2724 C01 INH I 1 27.808 26.376 23.301 0.50 27.77 I C ATOM 2733 C01 INH I 2 30.898 22.496 17.340 0.50 22.15 I C .... Basically, what I am trying to achieve is to uncouple the B-factor refinement of each individual instance of the ligand from its other instance(s). Are there any hidden pitfalls between these 2 scenarios I should be aware of? Many thanks for any input. Cheers Carsten _______________________________________________ phenixbb mailing list phenixbb(a)phenix-online.org http://www.phenix-online.org/mailman/listinfo/phenixbb

1 0

Question about restraints in individual B-factor refinement
by Schubert, Carsten [PRDUS] 25 May '07

25 May '07

Hi, I have a question about the handling of restraints in the individual B-factor refinement routine. What I'd like to do is to refine a ligand, which can be present in 2 or more different conformations/orientations in its binding site. I'd like to use B-factor refinement on the various instances of the ligand, which one is the most relevant one, assuming that the most relevant conformation/orientation is associated with the lowest B-factor (Validity of that assumption set aside ...) My question is now is there any difference in the restraints applied to the b-factors in the scenarios where A) the ligand is modeled as alternate conformations i.e. ATOM 2724 C01AINH I 1 27.808 26.376 23.301 0.50 27.77 I C ATOM 2733 C01BINH I 1 30.898 22.496 17.340 0.50 22.15 I C ... vs. scenario B) where the same ligand is modeled as 2 different residues ATOM 2724 C01 INH I 1 27.808 26.376 23.301 0.50 27.77 I C ATOM 2733 C01 INH I 2 30.898 22.496 17.340 0.50 22.15 I C .... Basically, what I am trying to achieve is to uncouple the B-factor refinement of each individual instance of the ligand from its other instance(s). Are there any hidden pitfalls between these 2 scenarios I should be aware of? Many thanks for any input. Cheers Carsten

2 1

probable RTFM question
by William Scott 17 May '07

17 May '07

Howdie: I just used phenix.geometry_minimization to try to see the effect of deleting some disulphide bonds in a protein structure. Is there a way to do the equivalent of CNS's model_anneal (i.e., simulated annealing with just a model)? It wasn't apparent from reading the phenix.refine docs. Thanks. Bill William G. Scott Contact info: http://chemistry.ucsc.edu/~wgscott/

3 2

Re: [phenixbb] statistics output by phenix.refine
by Ralf W. Grosse-Kunstleve 17 May '07

17 May '07

> Are these maximum likelihood targets the normal maximum likelihood > function or the minus one? We are using the maximum likelihood targets described here: Lunin, V.Y., Afonine, P.V. & Urzhumtsev, A. (2002). Acta Cryst. A58, 270-282. Lunin, V.Y. & Skovoroda, T.P. (1995). Acta Cryst. A51, 880-887. > Should I expect an increase or decrease of the target? Decrease. But see below... > If the part of the restraints is not combined with > the x-ray part, the target is showing is not the actual target the > program is trying to maximize or minimize. Yes, that's right. > The values don't have a trend, at least I can't see a trend in my > log file. The x-ray target may go up since we refine x-ray target + geometry target. The total will always go down (that's how the minimizer works), but an individual contribution may go up a bit. Ralf

1 0

statistics output by phenix.refine
by Jianghai Zhu 16 May '07

16 May '07

Hi, At the end of the log file output by phenix.refine, only R and Rfree are output. Is it possible to generate a statistical table which includes R, Rfree, FOM, LLG (total), LLG (free), rmsd (bond), rmsd (angle), ... for every refinement cycle? I want to monitor the progress of my refinement. Thanks. Jianghai +++++++++++++++++++++++++++++++ Jianghai Zhu, Ph.D CBR Institute for Biomedical Research Department of Pathology Harvard Medical School 200 Longwood Ave., Boston, MA 02115 Ph: 617-278-3211 Fx: 618-278-3232 +++++++++++++++++++++++++++++++

2 3

Using Phenix to calculate Atom Maps
by James Fraser 14 May '07

14 May '07

Hello, Is there a way to use Phenix to calculate maps that are just representative of the model? This would be similar to the ATMMAP mode of SFALL or using model_fcalc.inp/model_phase.inp in CNS. Thanks, James

3 3

Reflection output
by Rife, Christopher L 09 May '07

09 May '07

Hi, In trying to get a reflection file written out after refinement with phenix.refine (version 2007_03_13_2155), we found a couple of problems. 1. The documentation gives the example of including "export_final_f_model=mtz" in the command line input, however, phenix complains that it does not know this command. It is possible to work around this by changing "write_refined_fmodel_to_file" in the param file to true. 2. When the model is refined with twinning, trying to write out the reflection file causes it to crash with the following error: > Writing refined fmodel to plain text file: > /data5/Completed/PH10062D-47717/Phenix/model_001_fmodel.hkl > Traceback (most recent call last): > File > "/home/sdcsoftware/linux/cci_apps_sources/phenix/phenix/comman > d_line/refine.py", line 5, in <module> > command_line.run(command_name="phenix.refine", args=sys.argv[1:]) > File > "/home/sdcsoftware/linux/cci_apps_sources/phenix/phenix/refine > ment/command_line.py", line 86, in run > if(output.write_refined_fmodel_to_file): > refine_object.write_refined_fmodel(f_hkl) > File > "/home/sdcsoftware/linux/cci_apps_sources/phenix/phenix/refine > ment/driver.py", line 927, in write_refined_fmodel > self.fmodel.dump(out = f) > AttributeError: 'twin_model_manager' object has no attribute 'dump' We are trying to prepare files for deposition, so any assistance would be welcome. Thanks, Chris ----------------------------------------- Structure Determination Core Joint Center for Structural Genomics Stanford Synchrotron Radiation Laboratory

2 1

Remodelling using Phenix
by Ethayathulla A.S 04 May '07

04 May '07

hi I am doing structure at 2.3A resolution. It is a dimer not NCS. I have probelm in tracing last 10 amino acids in the C-terminus in both chains. Both chains come close each other and winds up. So Iam finding difficulty in tracing the chain. I used phenix for rebuilding with all the amino acid included. But autorebuilding is not rebuilding properly, whatever input is given it rebuilds on that. Actually it should be a helix and only at lower cuff off we can observe density for main chain so I tried phenix rebuild to improve it but it's not working. Is there any option to use phenix. I used rebuilding with option automatic and input coordinates and sequence. plz suggest do improve the phases. -- A.S.Ethayathulla Department of BIophysics All India Institute of Medical Sciences Ansai Nagar New Delhi-110029 India

2 1

Re: [phenixbb] PHENIX user mailing list <[email protected]>
by Ralf W. Grosse-Kunstleve 16 Apr '07

16 Apr '07

Hi James, > Number of atoms with unknown nonbonded energy type symbols: 14 > "CD CD D 1 " > "CD CD D 2 " This problem should be fixed in the latest phenix development release (1.26 beta @ http://www.phenix-online.org/download/index.html) and the slightly more recent CCI Apps (link at same web page). Are you using phenix.refine from an older release? Could you please try the latest phenix or CCI Apps? Ralf

1 0

Interpreting Cd atoms
by James Fraser 16 Apr '07

16 Apr '07

Hello, I'm trying to refine a structure with Cadmium bound. Is there a standard way to define "energy type symbols"? Thanks, James Number of atoms with unknown nonbonded energy type symbols: 14 "CD CD D 1 " "CD CD D 2 " "CD CD D 3 " "CD CD D 4 " "CD CD D 5 " "CD CD D 6 " "CD CD D 7 " "CD CD D 8 " "CD CD D 9 " "CD CD D 10 " ... (remaining 4 not shown) Time building chain proxies: 2.36, per 1000 atoms: 1.27 Sorry: Fatal problems interpreting PDB file: Number of atoms with unknown nonbonded energy type symbols: 14 Please edit the PDB file to resolve the problems and/or supply a CIF file with matching restraint definitions. Note that elbow.builder is available to create restraint definitions.

1 0

Re: [phenixbb] Special positions
by Ralf W. Grosse-Kunstleve 11 Apr '07

11 Apr '07

Hi Morton, > I have a symmetrical ligand, malonate -OOCCCOO-, that binds right on the > crystallographic axis of a C2 structure. Is the correct approach to set > the occupancy to 0.5 and then refine or is there a better way to handle is > in Phenix? The 0.5 trick currently doesn't work in phenix.refine since we don't provide a way to turn off nonbonded interactions between symmetry copies. However, we fully support handling of atoms on special positions and bonds involving symmetry operations. For this to work, you have to give only an asymmetric unit worth of atoms, i.e. in your case half of your molecule. Is the C in the middle exactly on the two-fold axis? If so, delete these atoms from the PDB file: -OOCCCOO- ^^^ Then use the custom bond and angle feature (for the angles you'll need the 2007_04_06 CCI Apps) to define the required bonds and angles, including the bonds to the protein. I hope it is straightforward. If not, could you send me the fragments from you pdb file with the malonate and the piece of protein it is linked to? Then I could work out the bond and angle definitions for you. Ralf

1 0

Special positions
by morteng＠bioxray.dk 11 Apr '07

11 Apr '07

Hello everybody I have a symmetrical ligand, malonate -OOCCCOO-, that binds right on the crystallographic axis of a C2 structure. Is the correct approach to set the occupancy to 0.5 and then refine or is there a better way to handle is in Phenix? Best regards -- Morten K Grøftehauge PhD student Department of Molecular Biology Gustav Wieds Vej 10 C 8000 Aarhus C - Denmark Phone: +45 89 42 50 10 Fax: +45 86 12 31 78 www.bioxray.dk

1 0

phenix.refine: How to define links and bond distances
by Matthew Bowler 03 Apr '07

03 Apr '07

Dear All, I have been trying to find the utility in elbow_builder to define bond distances and link distances with no luck, the .py file has bond order information for my ligand but nothing regarding distances etc, can anyone help? Thanks in advance, yours, Matt. Matthew Bowler MRC Dunn Human Nutrition Unit Wellcome Trust / MRC Building Hills Road Cambridge CB2 2XY Tel: 0044 (0) 1223 252826 Fax: 0044 (0) 1223 252825

2 1

Re: [phenixbb] coordination of metal ions
by Ralf W. Grosse-Kunstleve 29 Mar '07

29 Mar '07

Hi Jesper, could you send me the fragments from your PDB file with the metal ion and the bound residues? I'll work this into an example for making the required cif definitions manually. Ralf

1 0

coordination of metal ions
by Jesper Lykkegaard Karlsen 29 Mar '07

29 Mar '07

Dear Phenix-devels Is it possible set up geometric restrains for metal ions bound to protein in phenix.refine? elbow.builder seems only to work on separate molecules and not on metal ions with its coordinating sidechains. I found the feature to introduce custom bond distance between two atoms, but what about custom angles? And by the way will this these custom setting ignore Van Der Waals radii? or is there a way to set these for special atoms? Thanks! Jesper

2 1

riding hydrogen and mlhl
by Jianghai Zhu 28 Mar '07

28 Mar '07

Hi All, I just found a refinement.hydrogens entry in the phenix.refine. Can we use phenix.refine to generate riding hydrogens to help the geometry in low- or mid-resolution refinement now? Using hydrogens in refmac5 really gives a better result than not using hydrogens. Another question. Is there a way to blur the experimental or density modified phases when using mlhl target? Thanks. Jianghai +++++++++++++++++++++++++++++++ Jianghai Zhu, Ph.D CBR Institute for Biomedical Research Department of Pathology Harvard Medical School 200 Longwood Ave., Boston, MA 02115 Ph: 617-278-3211 Fx: 618-278-3232 +++++++++++++++++++++++++++++++

3 4

NCS weighting
by Jianghai Zhu 28 Mar '07

28 Mar '07

Hi, I am wondering how the NCS weighting in phenix.refine works. I guess the default values (coordinate_sigma = 0.05, b_factor_weight = 10) are pretty tight. What kind values would be appropriate for medium and loose NCS? Jianghai +++++++++++++++++++++++++++++++ Jianghai Zhu, Ph.D CBR Institute for Biomedical Research Department of Pathology Harvard Medical School 200 Longwood Ave., Boston, MA 02115 Ph: 617-278-3211 Fx: 618-278-3232 +++++++++++++++++++++++++++++++

4 6

Re: [phenixbb] ls_wunit_k1
by Ralf W. Grosse-Kunstleve 26 Mar '07

26 Mar '07

> All right. The target for my TLS is also LS. Is that intentional too? Pavel, who wrote the TLS code, is current offline. I don't know if Pavel ran systematic tests to determine the default target for TLS. A couple of months ago I ran systematic tests exercising rigid body refinement using both the LS and ML targets. Overall it turned out that ML is slightly better than LS. For an individual structure and a particular starting point it can be the other way around, but averaged over about 70 structures with 100 random starting points each ML generally outperforms LS. We should have updated the default accordingly, but it fell through the cracks. If you have the time to try out both ML and LS in the refinement of your structure, we'd be interested to know the results. Ralf P.S.: What's the version of phenix.refine you are using? -- I'm still puzzled why you didn't get the "automatic adjustment" messages.

4 5

Re: [phenixbb] ls_wunit_k1
by Ralf W. Grosse-Kunstleve 26 Mar '07

26 Mar '07

> You are right. Somehow, the default target function in > refinement.bulk_solvent_and_scale is LS, not ML. The default (phenix.refine --show-defaults=all) is actually ml. However, if the main target is set to ls the bulk_solvent_and_scale target is automatically reset to ls, too. There should be corresponding messages in the log file ("automatic adjustment"). Ralf

3 5

ls_wunit_k1
by Jianghai Zhu 26 Mar '07

26 Mar '07

Hi, I happen to notice that the target function of "Bulk solvent correction and scaling" in a phenix.refine run is ls_wunit_k1. Does that mean LS function for bulk solvent correction and scaling is the default in phenix.refine? If I uncheck ls_wunit_k1 in refinement.ls_target_names in my parameter file, will phenix.refine use the ML function or something bad will happen? Jianghai +++++++++++++++++++++++++++++++ Jianghai Zhu, Ph.D CBR Institute for Biomedical Research Department of Pathology Harvard Medical School 200 Longwood Ave., Boston, MA 02115 Ph: 617-278-3211 Fx: 618-278-3232 +++++++++++++++++++++++++++++++

2 3

Re: [phenixbb] The performance of phenix.refine
by Ralf W. Grosse-Kunstleve 22 Mar '07

22 Mar '07

> Did anybody actually compare the performance among phenix.refine, CNS > and refmac5? I am trying the phenix.refine from the latest CCI > bundle for my refinement. I have a relative large structure (800 kD > unique) and unit cell (140, 160, 550, 90, 90, 90). The data is at > low resolution ( 3.8 A) with experimental phases. My impression is > that phenix.refine takes a much longer time to finish the jobs, say > rigid body refinement or SA refinement, and it also takes a lot more > memories ( in my case, phenix.refine uses more than 700 MB, but other > programs use less than 200 MB). We have been concentrating on features and a high level of automation much more than raw speed and memory efficiency. Given that memory is so cheap these days (ca. $80/GB) there isn't much incentive for us to trade developer time for a smaller memory footprint. Speed optimizations are a bit higher up on the list. > Is this because the python code in > phenix.refine? CNS and refmac5 are written in FROTRAN, I think. Just > curious. phenix.refine doesn't spend much time in the Python code. It is a long while ago I did careful timings, but I'd guess it is significantly less than 10% of the time. C++ is intrinsically almost as fast as FORTRAN if the same development time is invested in the code. I.e. what you see is not a language issue, just a tradeoff features vs. optimization. We'll work on making phenix.refine a bit faster in the future, but at the minute more features are still our focus. Ralf

1 0

The performance of phenix.refine
by Jianghai Zhu 22 Mar '07

22 Mar '07

Hi, Did anybody actually compare the performance among phenix.refine, CNS and refmac5? I am trying the phenix.refine from the latest CCI bundle for my refinement. I have a relative large structure (800 kD unique) and unit cell (140, 160, 550, 90, 90, 90). The data is at low resolution ( 3.8 A) with experimental phases. My impression is that phenix.refine takes a much longer time to finish the jobs, say rigid body refinement or SA refinement, and it also takes a lot more memories ( in my case, phenix.refine uses more than 700 MB, but other programs use less than 200 MB). Is this because the python code in phenix.refine? CNS and refmac5 are written in FROTRAN, I think. Just curious. Jianghai +++++++++++++++++++++++++++++++ Jianghai Zhu, Ph.D CBR Institute for Biomedical Research Department of Pathology Harvard Medical School 200 Longwood Ave., Boston, MA 02115 Ph: 617-278-3211 Fx: 618-278-3232 +++++++++++++++++++++++++++++++

1 0

Re: [phenixbb] freeR syntax question in phenix.refine
by Ralf W. Grosse-Kunstleve 20 Mar '07

20 Mar '07

> I was running refinement from an mtz file from sigmaa in ccp4 with an > appropriate FreeR flag designated in the MTZ file. phenix.refine > worked fine and seemed to exclude the test reflections, but I wanted > to know what these commands do. Don't! :-) > r_free_flags { > file_name = "S19_sigmaa1.mtz" > label = "FreeR_flag" > test_flag_value = 0 --> don't know what this means? > disable_suitability_test = False --> implications of this? phenix.refine is designed to figure everything out automatically. You will have to specify file_name and label only if you have multiple R-free-flag arrays in the reflection files you give to phenix.refine, which is very unusual. phenix.refine does some sanity checks to determine if arrays that "look like" R-free-flag arrays are actually suitable for refinement, and reject them if "it doesn't think so". It will then tell you that you can override the sanity checks by setting disable_suitability_test to True. Again, this is very unusual. The test_flag_value could potentially be important. phenix.refine analyses the putative R-free-flag array(s) to guess if they came from CNS, where 0 usually means work set, 1 free set, or CCP4, which mostly uses values from 0 to 19, 0 is the default free set in REFMAC. Therefore, roughly, if phenix.refine "thinks" it is a CNS R-free-flag array it assigns test_flag_value=1, or 0 if it "thinks" it is a CCP4 array. Only if you used REFMAC before for refinement, and only if you manually instructed REFMAC to use a value different from 0 for the free set, will you have to manually tell phenix.refine about it, e.g. via r_free_flags.test_flag_value=19. If you are interested in more details, inspect the script iotbx/iotbx/reflection_file_utils.py, which does all the guesswork. Look for "get_r_free_flags_scores".

1 0

freeR syntax question in phenix.refine
by Simon Newstead 20 Mar '07

20 Mar '07

Dear BB, Sorry for the obvious question, but I couldn't find a description of the syntax for the FreeR flag function in phenix.refine. Is there a paragraph on the phenix web site that explains what I should put in these fields: r_free_flags { file_name = "S19_sigmaa1.mtz" label = "FreeR_flag" test_flag_value = 0 --> don't know what this means? disable_suitability_test = False --> implications of this? I was running refinement from an mtz file from sigmaa in ccp4 with an appropriate FreeR flag designated in the MTZ file. phenix.refine worked fine and seemed to exclude the test reflections, but I wanted to know what these commands do. Cheers, Simon.

2 1