- phenixbb - phenix-online.org

Re: [phenixbb] OS X installation problem
by Ralf W. Grosse-Kunstleve 17 Mar '07

17 Mar '07

I have not seen this error before. Could you send me the output of c++ --version Note that we have binary bundles for PowerPC Macs. Did you try that one? Ralf

2 1

OS X installation problem
by Arthur Glasfeld 17 Mar '07

17 Mar '07

Hi, I'm following Bill Scott's instructions for installing cci_apps on a G5 mac. I've attempted it from the self-extracting file and from the gzipped file. I type the following commands: /System/Library/Frameworks/Python.framework/Versions/Current/bin/ python ../cci_apps_sources/libtbx/configure.py phenix source setpaths.sh num_cpu=$( sysctl hw.logicalcpu | awk '{print $NF}' ) libtbx.scons -j $num_cpu . and things begin OK, but then after a bit of compiling, I get: /usr/local/cci_apps/cci_apps_sources/boost/libs/python/src/object/ inheritance.cpp:120: instantiated from here /usr/local/cci_apps/cci_apps_sources/boost/boost/graph/ graph_concepts.hpp:74: error: base 'boost::concepts::IncidenceGraph<boost::reverse_graph<boost::<unnamed>:: cast_graph, const boost::<unnamed>::cast_graph&> >' with only non- default constructor in class without a constructor scons: *** [boost/libs/python/src/object/inheritance.o] Error 1 scons: building terminated because of errors. usr+sys time: 5.92 seconds wall clock time: 37.99 seconds I had no problems in an intel mac. Anything obvious I should have done? Thanks in advance for any assistance, Arthur Arthur Glasfeld Department of Chemistry Reed College 3203 SE Woodstock Blvd. Portland, OR 97202 USA

1 0

Re: [phenixbb] phenix.refine crash -- other error
by Peter Zwart 16 Mar '07

16 Mar '07

Good to know. Let me know if the refinement is any faster. Cheers Peter Ulrich Baumann wrote: > Yop, that did it --- I sourced an old phenix setpaths and not the cci_apps. > > Thanks a lot, > > Ulrich > > Quoting Peter Zwart <PHZwart(a)lbl.gov>: > >> Hi Ulrich, >> >> Please make sure that the phenix.python used is actually the >> phenix.python you installed with cci_apps. >> >> The behavoir seen might be due to using an other phenix.python >> installation (on that loads in old python libraries). This can be >> avoided by making sure that the last crystallographic package sourced in >> your .chsrc file, is the latest phenix/cci_apps installation. >> >> A which phenix.python might give a clue what is run, and hitting the tab >> key or CRTl-D after typing phenix.refine might auytocomplete or give a >> list of possibilities. >> >> Let me know what happened. >> >> >> Peter >> >> >> >> >> >> >> >> >> >> Ulrich Baumann wrote: >>> Hi Peter, >>> >>> I hate to bother you guys with this but .... >>> >>> dcb-macbeth:[lip3] > phenix.python >>> /usr/local/cci_apps/cci_apps_sources/mmtbx/mmtbx/examples/reindex.py >>> data.file=lip_r3_1.8A.mtz model.file=in.pdb action=reindex >> standard_laws=niggli >>> #phil __OFF__ >>> ================= >>> REINDEX >>> A reindexing tool >>> ================= >>> >>> Traceback (most recent call last): >>> File >> "/usr/local/cci_apps/cci_apps_sources/mmtbx/mmtbx/examples/reindex.py", >>> line 468, in ? >>> reindex_utils(sys.argv[1:]) >>> File >> "/usr/local/cci_apps/cci_apps_sources/mmtbx/mmtbx/examples/reindex.py", >>> line 260, in reindex_utils >>> combined_xs = crystal.select_crystal_symmetry( >>> AttributeError: 'module' object has no attribute 'select_crystal_symmetry' >>> >>> >>> I seem to hae to set a spacegroup? Sorry, I could not find a documentation >> for >>> the reindex utility. >>> >>> Many thanks, >>> >>> Ulrich >>> >>> Quoting Peter Zwart <PHZwart(a)lbl.gov>: >>> >>>> Hi Ulrich, >>>> >>>> You could try refining your structure in a primitive setting, this might >>>> alleviate memory usage a bit (not sure though). >>>> >>>> try >>>> >>>> phenix.python $MMTBX_DIST/mmtbx/examples/reindex.py data.file=mydata.mtz >>>> model.file=mymodel.pdb action=reindex standard_laws=niggli >>>> >>>> You will get a new mtz file and a new pdb file. Unfortunately, the free >>>> flags are not copied over, so you would have to reindex those separately >>>> using iotbx.reflection_file_converter and read that file in separately >>>> in phenix.refine. >>>> >>>> This might be a bit criptic, if so, let me know and I will post more >>>> details. >>>> >>>> I am not sure this does the trick though, I hope Ralf or Pavel can give >>>> some suggestions. >>>> >>>> Cheers >>>> >>>> Peter >>>> >>>> >>>> >>>> >>>> Ulrich Baumann wrote: >>>>> Hi there, >>>>> >>>>> if you remove the ANISOU card from the input pdb file, phenix will run >>>> fine >>>>> ....sometimes not. >>>>> >>>>> phenix.refine ( Version: 2007_01_24_2154) crashed for me as well: we >> have >>>> a >>>>> rather large unit cell and good resolution. The program crashes >>>> reproducibly >>>>> (see below) when using resolution limits better than 3 A, although >>>> sometimes >>>>> this only occurs in the 2nd macrocycle. >>>>> >>>>> Working crystal symmetry after inspecting all inputs: >>>>> Unit cell: (202.401, 202.401, 317.731, 90, 90, 120) >>>>> Space group: R 3 :H (No. 146) >>>>> >>>>> >>>>> We have 2 GB memory and could stuff in more but maybe we can limit the >>>> memeory >>>>> use in a different way? >>>>> >>>>> >>>>> Thanks a lot for your help, >>>>> >>>>> Ulrich >>>>> >>>>> >>>>> ============================= ml refinement start >>>> ============================= >>>>> ----------structure factors based statistics (before >>>> refinement)---------- >>>>> Traceback (most recent call last): >>>>> File >>>> "/usr/local/cci_apps/cci_apps_sources/phenix/phenix/command_line/refine.p >>>>> y", line 5, in <module> >>>>> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) >>>>> File >>>> "/usr/local/cci_apps/cci_apps_sources/phenix/phenix/refinement/command_li >>>>> ne.py", line 75, in run >>>>> call_back_after_monitor_collect=call_back_after_monitor_collect) >>>>> File >>>> "/usr/local/cci_apps/cci_apps_sources/phenix/phenix/refinement/driver.py" >>>>> , line 1108, in run >>>>> log = log) >>>>> File >>>> "/usr/local/cci_apps/cci_apps_sources/phenix/phenix/refinement/strategies >>>>> .py", line 174, in refinement_machine >>>>> abcd = abcd) >>>>> File "/usr/local/cci_apps/cci_apps_sources/mmtbx/mmtbx/f_model.py", >> line >>>> 185, >>>>> in __init__ >>>>> b_cart = b_cart) >>>>> File "/usr/local/cci_apps/cci_apps_sources/mmtbx/mmtbx/f_model.py", >> line >>>> 466, >>>>> in update_xray_structure >>>>> f_mask = bulk_solvent_mask_obj.structure_factors(miller_set= >>>> self.f_obs) >>>>> File "/usr/local/cci_apps/cci_apps_sources/mmtbx/mmtbx/masks.py", line >>>> 97, in >>>>> structure_factors >>>>> flex.grid(fft_manager.m_real()).set_focus(fft_manager.n_real())) >>>>> MemoryError >>>>> >>>>> ---------------------------------------------------------- >>>>> ------------------------------------------------------- >>>>> Quoting Pavel Afonine <pafonine(a)lbl.gov>: >>>>> >>>>>> Hi, >>>>>> >>>>>> this is known problem and it will be fixed ASAP; sorry for this. >>>>>> This presumably happens because an atom selected to be refined as >>>>>> individual anisotropic is included into a TLS group. >>>>>> I will let you know once the fixed version of CCI Apps is available. >>>>>> >>>>>> Pavel. >>>>>> >>>>>> >>>>>> >>>>>> Jianghai Zhu wrote: >>>>>>> Hi, >>>>>>> >>>>>>> I was using phenix.refine from the latest cci apps bundle to do >>>>>>> tls+individual_sites+individual_adp refinement. The program crashed >>>>>>> with the following error message. >>>>>>> >>>>>>> ----------Target weights (before 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 75, 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 1108, in run >>>>>>> log = log) >>>>>>> File >>>>>>> >> "/nfs/home/jzhu/cci_apps_sources/phenix/phenix/refinement/strategies.py", >>>>>>> line 245, in refinement_machine >>>>>>> log = log) >>>>>>> File >>>>>>> >> "/nfs/home/jzhu/cci_apps_sources/phenix/phenix/refinement/weight_xray_chem.py", >>>>>>> line 222, in __init__ >>>>>>> compute_gradients = True) >>>>>>> File "/nfs/home/jzhu/cci_apps_sources/mmtbx/mmtbx/restraints.py", >>>>>>> line 144, in energies_adp_aniso >>>>>>> xray_structure = xray_structure, >>>>>>> File >>>>>>> >> "/nfs/home/jzhu/cci_apps_sources/cctbx/cctbx/adp_restraints/__init__.py", >>>>>>> line 109, in __init__ >>>>>>> self.check_flags(fl_i) >>>>>>> File >>>>>>> >> "/nfs/home/jzhu/cci_apps_sources/cctbx/cctbx/adp_restraints/__init__.py", >>>>>>> line 177, in check_flags >>>>>>> assert not fl.use_u_aniso() >>>>>>> AssertionError >>>>>>> >>>>>>> Anyone has any idea what is happening? 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 >>>>>>> +++++++++++++++++++++++++++++++ >>>>>>> >>>>>>> >>>>>>> >>>>>>> >> ------------------------------------------------------------------------ >>>>>>> _______________________________________________ >>>>>>> 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 >>>> >>> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://www.phenix-online.org/mailman/listinfo/phenixbb >> > >

1 0

Re: [phenixbb] pdb waters
by Pavel Afonine 16 Mar '07

16 Mar '07

Hi, yes, the waters in your file do not match the Monomer Library definition. If you simply replace O1 with O that will fix the problem. For example: GOOD record: ATOM 22792 O HOH X 1 37.309 34.504 -14.174 1.00 38.28 O BAD record: ATOM 22792 O1 HOH X 1 37.309 34.504 -14.174 1.00 38.28 O > I have tried running elbow on my pdb already with no luck. using your > command I get the error message: elbow.builder: error: no such option: > --do-all That makes me thinking that you are using an old version of CCI Apps or PHENIX. Are you using the latest one? Is the command you run looks like: % elbow.builder x.pdb --do-all Nigel: do you have any comments? Pavel.

5 5

Re: [phenixbb] Few questions -- FreeR flag specification
by Tom Terwilliger 15 Mar '07

15 Mar '07

Hi John, I'm sorry for the slow reply; I was on a site visit and out of email contact. The AutoBuild wizard in the current version (which I am afraid you probably do not have!) does try to take the FreeR flags from your input file. You can also tell it specifically to do so with the GUI and in the command-line with refinement_file=ref.mtz input_refinement_labels="FP SIGFP None None None None None None FreeR" If you have version 1.25a , then you still should be able to specify the column for FreeR flags in the GUI window that asks which column labels to use (I could be wrong about that, however, as I don't have this version on my computer any longer). You should also be able to specify them with a script: phenix.runWizard AutoBuild AutoBuild.start_from_model terwill@rosie> cat $PHENIX/phenix/phenix/wizard_scripts/AutoBuild.start_from_model # AutoBuild.start_from_model # script to run a AutoBuild starting from a model # Version 1.26 of 26-Jan-2007 # run with: # phenix.runWizard AutoBuild AutoBuild.start_from_model [restart] # if you specify restart it will forget everything, # otherwise it will continue on from where it is now. input_ncs_file NONE # heavy-atom or PDB file for finding NCS input_data_file $PHENIX/examples/p9-sad/solve_1.mtz # file with experimental data. # Can be any format if only FP SIGFP; must be mtz if phases or HL or FreeR input input_labels FP SIGFP None None None None None None FreeR_flag # optional # line defining the 9 columns in input_mtz_file corresponding to # FP SIGFP PHIB FOM HLA HLB HLC HLD FreeR (use "None" for missing) input_map_file $PHENIX/examples/p9-sad/resolve_1.mtz # Optional map-coefficients file for first rebuild input_map_labels FP PHIM FOMM # optional line defining the 3 columns # in input_map_file corresponding to FP PHIB FOM ("None" for missing) input_seq_file $PHENIX/examples/p9-sad/seq.dat # file with sequence in 1-letter form # separate chains with >>> resolution 2.5 # high-resolution limit solvent_fraction 0.60 # you can put any other parameters here input_pdb_file $PHENIX/examples/p9-sad/p9.pdb # input PDB file for rebuilding rebuild_in_place True # rebuild my_model.pdb in place # (as opposed to building from scratch) start_chains_list 4 # chain 1 starts at residue 4 # rebuild residues 21-27 of chain "" and 31-42 of chain "" # rebuild_res_start_list 21 31 # rebuild_res_end_list 27 42 # rebuild_chain_list "" "" # n_cycle_rebuild_min 5 # n_cycle_rebuild_max 10 # AutoBuild will run automatically from here on with no more inputs I hope that helps! The new version of PHENIX that will be out very soon fixes many other problems as well! All the best, Tom T At 03:26 AM 3/13/2007, john kryst wrote: >Hi phenixbb !! > > I have few questions. > >1. Is there any utility in phenix to write out pdb formatted (which >is ready to deposit) output file.. like xtal_pdbsubmission utility in CNS. >2. In Autobuild wizard how to ask the program to use the Free R >flags present in the file. It always generates new Free R flags >during refinement. >One general question. > >Does the rmsd estimation (for eg. rmsd_bonds ) depends on the >program we use ?? > >I have been refining in CNS. then i used Autobuild wizard to build >my model. when i went back to cns the rmsd of bond increased ( >0.0218 A). But phenix.refine says it it 0.007A. > >Thanks in advance for your kind inputs. > >regards >john >_______________________________________________ >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

Re: [phenixbb] phenix.refine problems with my pdb
by Nicholas Noinaj 15 Mar '07

15 Mar '07

matt, can't offer much help, but have you tried to remove the ligand and try refinement? it would let you know if the pdb file or the presence of the ligand is the problem. at least let you would know where you should focus your attention. if all else fails, have you tried getting parameters for your ligand from the PRODRG Server and using that for inputing the ligand parameters into phenix.refine directly when you setup the refinement job? (I think i have done this before???) here is the website in case you aren't familiar... http://davapc1.bioch.dundee.ac.uk/programs/prodrg/ Again, Best of luck! cheers, nick -----Original Message----- From: Matthew Bowler <mwb(a)mrc-dunn.cam.ac.uk> To: phenixbb(a)phenix-online.org Date: Thu, 15 Mar 2007 13:00:42 +0000 Subject: [phenixbb] phenix.refine problems with my pdb Dear All, I am trying to run phenix.refine but I keep getting the message: Sorry: Fatal problems interpreting PDB file: Number of atoms with unknown nonbonded energy type symbols: 308 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. I have run my ligands through elbow and used the output .cif but I still get the same message. I can find no info on this problem in the manual or the BB, 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 _______________________________________________ phenixbb mailing list phenixbb(a)phenix-online.org http://www.phenix-online.org/mailman/listinfo/phenixbb

6 6

Re: [phenixbb] phenix.refine problems with my pdb
by Nicholas Noinaj 15 Mar '07

15 Mar '07

matt, >>> Number of atoms with unknown nonbonded energy type symbols: 308 just curious, did you have 308 waters? as peter said, hard to determine exact problem without some files or examples to look at. also phenix does a good job with water picking, in my little experience. just one more question, when you do get the refinement to work, is your ligand(s) included with the protein model? just no waters, right? it wasn't completely clear from your message. if the ligands aren't included, i recommend adding it(them) back and see if refinement works. then if that works, try water picking within phenix. again, best of luck. cheers, nick -----Original Message----- From: Matthew Bowler <mwb(a)mrc-dunn.cam.ac.uk> To: noinaj(a)uky.edu, PHENIX user mailing list <phenixbb(a)phenix-online.org> Date: Thu, 15 Mar 2007 14:15:45 +0000 Subject: Re: [phenixbb] phenix.refine problems with my pdb Thanks for replies so far.... I removed all ligands and then waters - the waters seem to be the problem as the protein model is now refining, do they need to be treated/marked in a different way for phenix? Thanks again, Matt. On 15 Mar 2007, at 13:32, Nicholas Noinaj wrote: > matt, > > can't offer much help, but have you tried to remove the ligand and > try refinement? it would let you know if the pdb file or the > presence of the ligand is the problem. at least let you would know > where you should focus your attention. > > if all else fails, have you tried getting parameters for your > ligand from the PRODRG Server and using that for inputing the > ligand parameters into phenix.refine directly when you setup the > refinement job? (I think i have done this before???) > > here is the website in case you aren't familiar... > http://davapc1.bioch.dundee.ac.uk/programs/prodrg/ > > Again, Best of luck! > > > > cheers, > nick > > > -----Original Message----- > From: Matthew Bowler <mwb(a)mrc-dunn.cam.ac.uk> > To: phenixbb(a)phenix-online.org > Date: Thu, 15 Mar 2007 13:00:42 +0000 > Subject: [phenixbb] phenix.refine problems with my pdb > > Dear All, > I am trying to run phenix.refine but I keep getting the message: > > > Sorry: Fatal problems interpreting PDB file: > Number of atoms with unknown nonbonded energy type symbols: 308 > 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. > > I have run my ligands through elbow and used the output .cif but I > still get the same message. I can find no info on this problem in > the manual or the BB, 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 > > _______________________________________________ > 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 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

1 0

Few questions
by john kryst 15 Mar '07

15 Mar '07

Hi phenixbb !! I have few questions. 1. Is there any utility in phenix to write out pdb formatted (which is ready to deposit) output file.. like xtal_pdbsubmission utility in CNS. 2. In Autobuild wizard how to ask the program to use the Free R flags present in the file. It always generates new Free R flags during refinement. One general question. Does the rmsd estimation (for eg. rmsd_bonds ) depends on the program we use ?? I have been refining in CNS. then i used Autobuild wizard to build my model. when i went back to cns the rmsd of bond increased (0.0218 A). But phenix.refine says it it 0.007A. Thanks in advance for your kind inputs. regards john

4 4

phenix.refine problems with my pdb
by Matthew Bowler 15 Mar '07

15 Mar '07

Dear All, I am trying to run phenix.refine but I keep getting the message: Sorry: Fatal problems interpreting PDB file: Number of atoms with unknown nonbonded energy type symbols: 308 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. I have run my ligands through elbow and used the output .cif but I still get the same message. I can find no info on this problem in the manual or the BB, 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

Two groups per residue and multiple trajectories for simulated annealing
by Stephen Graham 02 Mar '07

02 Mar '07

Hi all, I am taking my first tentative steps into the world of phenix.refine with a somewhat problematic 3 A data set and I had a few questions that Google couldn't answer. I am running phenix.refine version 1.24.1b (intel-linux-2.6) from the command line. Firstly, is it possible to do 'grouped' adp refinement with 2 groups per residue (like the good old bgroup.inp of CNS)? Secondly, is it possible to specify multiple starting trajectories for simulated annealing? I often find it handy to do 10 runs of simulated annealing and overlay the output co-ordinates to see where the 'wobbly' bits are. Lastly, can the latest version of phenix.refine do SA omit maps? Thanks, Stephen P.S. My .def file is the default (from phenix.refine --dry-run) w/ the following changes: 86,87c86,87 < refine_adp_individual = True < refine_adp_group = False --- > refine_adp_individual = False > refine_adp_group = True 92c92 < simulated_annealing = False --- > simulated_annealing = True -- 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

2 1

Phased anomalous difference map
by Luca Jovine 28 Feb '07

28 Feb '07

Hello, Is it possible to calculate a phased anomalous difference map with phenix.refine? Thanks, Luca

1 0

Re: [phenixbb] phenix.refine crash
by Ralf W. Grosse-Kunstleve 24 Feb '07

24 Feb '07

> Anyone has any idea what is happening? Thanks. It is a bug that we've fixed already. We will release a new cci apps in the next couple of days that will include the fix. In the meantime, could you try removing the ANISOU cards (grep -v ANISOU xxx.pdb)? I think that should work around the problem. Ralf

1 0

Re: [phenixbb] input HA
by Tom Terwilliger 13 Feb '07

13 Feb '07

Hi Morton, I'm sorry to make it difficult! Here are 2 formats that should work: 1. PDB format (I'm not sure why the one from cns failed...can you send me a line or two from that so I can make it work?): HETATM 1 SE SE 3 73.920 0.324 3.857 0.92 54.60 Se HETATM 2 SE SE 4 39.591 8.262 2.158 0.55 60.00 Se 2. Plain fractional xyz, with or without the word "xyz". For example see: $PHENIX/examples/gene-5-mad/sites.xyz : xyz 0.0271 0.8074 0.0657 xyz 0.0239 0.0116 -0.0965 Please let me know if this doesn't do it! -Tom T At 03:52 AM 2/13/2007, Morten Grøftehauge wrote: >Ok, I give up. > >What is the correct format for a heavy atom file >for AutoSOL? I've tried both the outputs from >FFT, CNS pdbs and sdbs and I've tried just >putting them on a line like you do for SHARP: >ATOM Br - 0.0120 0.5731 0.2740 >but none of them are accepted by Phenix. >Any help please? > >-- >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 ><http://www.bioxray.dk>www.bioxray.dk > >_______________________________________________ >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

input HA
by Morten Grøftehauge 13 Feb '07

13 Feb '07

Ok, I give up. What is the correct format for a heavy atom file for AutoSOL? I've tried both the outputs from FFT, CNS pdbs and sdbs and I've tried just putting them on a line like you do for SHARP: ATOM Br -0.0120 0.5731 0.2740 but none of them are accepted by Phenix. Any help please? -- 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

CCI Apps
by Noinaj 26 Jan '07

26 Jan '07

Hi, Quick question. I am trying to refine occupancy of a ligand and was attempting to do this using the method outlined under the CCI Apps documentation for phenix.refine. a.. Refinement of one occupancy factor per selected group of atoms (group occupancy refinement): % phenix.refine data.hkl model.pdb strategy=group_occupancies \ refine.occupancies.group="chain A and resid 1"Yet, when I type the commands, I keep getting an error saying that something is wrong with my commands: strategy, group, etc. I have not (at least I don't recall) loading the updated CCI Apps and only have what was installed with Phenix 1.24b installation. Am I doing something wrong here or do I need to update the CCI Apps? Thanks in advance for all feedback! Cheers, Nick ________________________________________ Nicholas Noinaj University of Kentucky School of Medicine Department of Molecular and Cellular Biochemistry Biomedical Biological Sciences Research Building, Rm 236 741 S. Limestone Lexington, Ky 40536 Lab: 859-323-8183 Home: 859-228-0978 noinaj(a)uky.edu www.noinaj.com

2 1

New CCI Apps
by Ralf W. Grosse-Kunstleve 25 Jan '07

25 Jan '07

On behalf of the Phenix development team I'm pleased to announce the release of a new CCI Apps version: http://phenix-online.org/download/cci_apps/ Version 2007_01_24_2154 Simultaneous use of NCS restraints and automatic water picking. Individual anisotropic ADP restraints (similar to SHELXL SIMU). Bug fix in processing of multiple apply_cif_link blocks applied to the same residue. Send bug reports to: bugs(a)phenix-online.org For help write to: help(a)phenix-online.org

1 0

questions
by Jianghai Zhu 17 Jan '07

17 Jan '07

Hi, Got some questions in using phenix.refine. 1) Can phenix.refine use the ligand library generated by refmac5? 2) How do I use elbow to generate the ligand library when I have more than one kind of ligands? When I used --residue=xxx more than once, elbow will only generated library for the last ligand. If I use elbow multiple times to generate multiple libraries, how do I read the libraries in? Should I combine them to one file? 3) When I ran phenix.refine, I found the following message in the log. Number of residues, atoms: 32, 418 Unusual residues: {'NAG': 22, 'MAN': 10} Classifications: {'undetermined': 32} Link IDs: {None: 31} Unresolved non-hydrogen bonds: 32 Unresolved non-hydrogen angles: 64 Unresolved non-hydrogen dihedrals: 22 Unresolved non-hydrogen chiralities: 32 Number of residues, atoms: 18, 28 Unusual residues: {'GOL': 2, ' MG': 2, ' CA': 14} Classifications: {'undetermined': 18} Link IDs: {None: 17} I already gave the cif library for NAG. Looks like nothing happened. How come phenix.refine doesn't recognize Glycerol, Mg, and Ca? 4) The refinement (TLS + ML + B individual) went through, I got reasonable R, Rfree, rmsdBOND, rmsdANGLE. But the B factors are pretty low. The B factor of the backbone is much lower than the side chain, some have numbers like 4. Some metal atoms also have B factors around 4. What did I do wrong? Thanks. Jianghai

6 13

Re: [phenixbb] phenix.hyss
by Ralf W. Grosse-Kunstleve 17 Jan '07

17 Jan '07

Hi Jianghai, I think the large memory consumption is due to the large number of triplet phase relations. Using the --real_space_squaring option may help. > The MAD data only goes to 4 A. I tried SOLVE and based on the CC > table output from SOLVE, the anomalous signal goes to 4.3 A. I > believe the cut of CC in SOLVE is 0.3. Unfortunately, the native > data is not isomorphous. Finding 50 Se given only 4.3 A resolution is difficult. I'm sorry to say I'm not very hopeful, but I'd try all tools out there. I'd also play with the resolution cutoff. hyss used 4 A in your case. I'd try 4.5 and maybe even all cutoffs 4.4, 4.3, 4.2, 4.1 if you can afford the CPU time, with --real_space_squaring first, without it if nothing worked still. Another question is the space group. You cannot distinguish I222 and I212121 from the systematic absences. Are you certain it is I212121? To find a cluster, lower the resolution significantly. I've had success with 6 A cutoff (--resolution=6). The atom type doesn't really matter, Ta should be fine. If it works, hyss will give you atoms around the centers of the clusters. You have to fit the whole cluster yourself. I'm not sure about the best approach for this (inspecting the map with coot and place the cluster manually?). Maybe someone else can help? How many clusters do you expect? In general, the more you expect the more data you need. I.e. a 6 A cutoff may work for finding a couple of sites, but for say 10 sites you'd need 5 A. -- These numbers are just to convey the idea. What works exactly depends highly on the available data. Hope this helps. Cheers, Ralf

1 0

Re: [phenixbb] phenix.hyss
by Ralf W. Grosse-Kunstleve 17 Jan '07

17 Jan '07

Hi Jianghai, phenix.hyss isn't normally using that much memory. To help me understand what is happening, could you please email me the complete output (I will keep it confidential if you email to rwgk(a)cci.lbl.gov directly). Thanks! Ralf

1 0

mmtbx.twin_map_utils utility
by Jorge Iulek 03 Jan '07

03 Jan '07

Hallo, This must be a naive question, but following Dr. Zwart's previous post and as I found in the web (http://www.phenix-online.org/phenix_wiki/mmtbx.twin_map_utils) this should be a utility to calculate a detwinned map. I just installed phenix and everything seemed to be fine during installation. But, where is this utility ? I cannot find it. Thanks, *********************************************************************** Prof. Jorge Iulek, Ph. D. Protein Purification and 3D Structure Determination Group State University of Ponta Grossa Telephone: ++ 55 - 42 - 3220-3062 Fax: ++ 55 - 42 - 3220-3262 iulek(a)uepg.br iulek(a)interponta.com.br *********************************************************************** !!!!!!!!!!!!!!!!!!!!!!!!!!Hi Brian, There are map utilities for twinned data indeed. mmtbx.twin_map_utils carries out refinement of the overall, anisotropic scale, bulk solvent parameters and twin fraction. It spits out the maps afterwards. Unfortunately, refniement against twinned data is not yet fully implemented in phenix.refine but we are working on that. HTH Peter > i see there are map utilities for twinned data. > > but before i go try refinement of a twin - will phenix.refine refine > against twinned data? > > if so, is alpha determined each time, thus, alpha is not specified > in the > .def? > > -bryan > _______________________________________________ > phenixbb mailing list > phenixbb at phenix-online.org > http://www.phenix-online.org/mailman/listinfo/phenixbb > -- Esta mensagem foi verificada pelo sistema de antivírus e acredita-se estar livre de perigo.

1 0

full omit maps
by Miguel Ortiz Lombardia 29 Dec '06

29 Dec '06

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hello everyone, I'm trying to obtain a full omit map as it was described by Tom Terwilliger in a previous post: http://www.phenix-online.org/pipermail/phenixbb/2006-October/000108.html However, when I run the AutoBuild strategy asking for no rebuilding / no refinement all I get in the FULL_OMIT directory is a collection of pdb files named 'overall_best.pdb_OMIT_REGION_x' where 'x' correspond to each of the 'omit regions'. But I get no 'resolve_composite_map.mtz' whatsoever. Surely I'm missing something. I would appreciate your guidance. Have all of you an excellent year 2007! All the best, Miguel - -- Miguel Ortiz Lombardía Centro de Investigaciones Oncológicas C/ Melchor Fernández Almagro, 3 28029 Madrid, Spain Tel. +34 912 246 900 Fax. +34 912 246 976 email: molatwork(a)yahoo.es www: http://www.ysbl.york.ac.uk/~mol/ ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Le travail est ce que l'homme a trouvé de mieux pour ne rien faire de sa vie. (Raoul Vaneigem) -----BEGIN PGP SIGNATURE----- Version: GnuPG v1.4.2.2 (GNU/Linux) iD8DBQFFlNB7F6oOrDvhbQIRAuBpAJwK4/+sj+5PowBBIvNVE0Tvyin7wwCdHFOS QHrtAM+9afN+2RTGshFB8KU= =IXga -----END PGP SIGNATURE-----

4 4

Re: [phenixbb] CCI apps bundles upgrading
by Ralf W. Grosse-Kunstleve 20 Dec '06

20 Dec '06

Hi Miguel, Thanks for asking the question! I just added this information to the CCI Apps download page: Please do NOT install CCI Apps inside an existing Phenix installation. It is best to install the CCI Apps bundle in a new, separate, empty directory. To override the Phenix commands with more recent equivalents in CCI Apps (e.g. phenix.refine), source cci_apps_build/setpaths.csh AFTER source phenix_env. Sorry for not explaining this earlier. It is just one of the things we don't usually think about as developers. Please let us know if there are other issues like this. Cheers, Ralf

2 1

CCI apps bundles upgrading
by Miguel Ortiz Lombardia 20 Dec '06

20 Dec '06

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Dear all, This is a pretty newbie question. I have phenix installed at: /usr/local/phenix-1.24.1b and have downloaded the new version of the CCI Apps. I would like to install/upgrade them in my current installation (I have a backup, of course ;-) Should I simply run: python /PATH-TO/cctbx_sources/libtbx/configure.py phenix in the /usr/local/phenix-1.24.1b directory and then do the: source setpaths.csh libtbx.scons or it should be done differently to keep the phenix installation working? Thank you! Cheers, Miguel - -- Miguel Ortiz Lombardía Centro de Investigaciones Oncológicas C/ Melchor Fernández Almagro, 3 28029 Madrid, Spain Tel. +34 912 246 900 Fax. +34 912 246 976 email: molatwork(a)yahoo.es www: http://www.ysbl.york.ac.uk/~mol/ ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Le travail est ce que l'homme a trouvé de mieux pour ne rien faire de sa vie. (Raoul Vaneigem) -----BEGIN PGP SIGNATURE----- Version: GnuPG v1.4.2.2 (GNU/Linux) iD8DBQFFiVtPF6oOrDvhbQIRAkcKAJ414Qrsa+8DSj67uFQ6gv56Sc7TRQCfY8xG 7jXKhu+xv6nCc90xYjXt0T8= =2vZq -----END PGP SIGNATURE-----

2 1

Anisotropic ADP
by James Fraser 16 Dec '06

16 Dec '06

Hello, We've just started using phenix.refine and I'm a bit confused about whether anisotropic ADP's are actually being used or not. On the wiki (http://www.phenix-online.org/phenix_wiki/phenix.refine) it says: << phenix.refine is the macromolecular refinement engine of the Phenix suite. It features * Individual as well as grouped B factor refinement * Anisotropic or isotropic+anisotropic B factor refinement * TLS refinement >> However, the documentation says (http://www.phenix-online.org/download/documentation/phenix/phenix/structure…): << Current limitations * No restrained individual anisotropic ADP refinement >> In my .def file I can set: adp { individual { isotropic = None anisotropic = chain a force_all_to_be_refined_isotropically = True } and get ANISO lines in my output pdb. But I'm confused as to whether this is actually refining with individual anisotropic b's for chain a (analagous to refining with anisotropic temperature factors in refmac). Should I be resetting the force_all... to false? Can someone with more knowledge clarify this? Thanks, the program seems to be working great otherwise! James Alber Lab UC Berkeley

2 1

Re: [phenixbb] questions / TLS+NCS bug
by Pavel Afonine 14 Dec '06

14 Dec '06

Hi Jianghai, Thanks for the .log file ! I looked at it and now I think I now what's the problem: this is the bug in phenix.refine that we will fix for the next release. The bug arises when you try to refine using TLS and NCS at once. Sorry for this. The possible solution is to split your refinement into 2 parts: 1) First, refine coordinates and isotropic B-factors. Use NCS restraints. DO NOT refine TLS. 2) At the end, as the final tune up, refine ONLY isotropic B-factors + TLS and do not refine coordinates. For this run, please REMOVE all NCS information (NCS selections) from the command file. This refinement run should be the final one. Once again, sorry for the inconvenience. This will be fixed in the next release of CCI APPS and PHENIX. Once this fixed, you will be able to run everything in "one go" and in any combination. Please let me know if what I suggested helped you. Any further questions are welcome! Cheers, Pavel. Jianghai Zhu wrote: > Here is the log file. Thanks. > > Jianghai > > = > ------------------------------------------------------------------------ > > > > On Dec 14, 2006, at 12:50 PM, Pavel Afonine wrote: > >> Hi Jianghai, >> >> could you please send us .log file from your refinement run, so we >> can analyze what's going on. >> >> In general, "bad" B-factors can be: >> - misplaced model; >> - inadequate TLS model (= domains chosen for TLS do not correspond to >> the reality). >> >> If you are using like 2 months old or older version of phenix.refine, >> you may want to get the latest CCI APPS since we made lots of >> improvements. Just goto http://www.phenix-online.org/download/cci_apps/ >> >> Pavel. >> >> >> Jianghai Zhu wrote: >>> The resolution is 2.5 A. The wilson B is about 50. I know B factor >>> of the backbone is lower than that of the sidechian. But a B factor >>> like 4 is definitely wrong. >>> >>> Jianghai >>> >>> >>> On Dec 14, 2006, at 12:11 PM, Peter Zwart wrote: >>> >>>> >>>>>> 4) The refinement (TLS + ML + B individual) went through, I got >>>>>> reasonable R, Rfree, rmsdBOND, rmsdANGLE. But the B factors are >>>>>> pretty low. The B factor of the backbone is much lower than the >>>>>> side >>>>>> chain, some have numbers like 4. Some metal atoms also have B >>>>>> factors around 4. What did I do wrong? >>>> >>>> What is the resolution of your data? >>>> >>>> Backbone B-s usually are lower than the main chain. >>>> >>>> What is the Wilson B value reported by phenix.refine? >>>> You could re-refine and randomize all B-values and see what happens (I >>>> have to get back to you to to get the exact command for this). >>>> Maybe it >>>> is useful to obtain a copy of the latest verison of phenix.refine by >>>> downloading cci_apps from our server http://www.phenix-online.org. >>>> >>>> >>>> If your B-values still come out lowish, try growing crystals that >>>> do not >>>> diffract very well, that usually does the trick. >>>> >>>> HTH >>>> >>>> Peter >>>> >>>> >>>> >>>> >>>> >>>> >>>> >>>> >>>> _______________________________________________ >>>> phenixbb mailing list >>>> phenixbb(a)phenix-online.org <mailto:[email protected]> >>>> http://www.phenix-online.org/mailman/listinfo/phenixbb >>> >>> ------------------------------------------------------------------------ >>> _______________________________________________ >>> phenixbb mailing list >>> phenixbb(a)phenix-online.org >>> http://www.phenix-online.org/mailman/listinfo/phenixbb >>> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org <mailto:[email protected]> >> http://www.phenix-online.org/mailman/listinfo/phenixbb > > =

1 0