- phenixbb - phenix-online.org

Using the Same Test Set in AutoBuild and Phenix.Refine
by Dale Tronrud 02 Jan '08

02 Jan '08

Hi all, I have another problem, I'm afraid. I have built a model using phenix.autobuild and now want to run some refinement. While in the long run I'll do some manually rebuilding using Coot I just wanted to run a test of phenix.refine to ensure I have the script right and have a baseline to compare against later. My autobuild script is: phenix.autobuild model=AutoMR_run_4_/MR.1-protein.pdb data=1M50-2.mtz \ input_refinement_labels="FP SIGFP None None None None None None FreeR_flag" \ map_file=AutoMR_run_4_/MR.MAP_COEFFS.1.mtz \ seq_file=../fmo-ct.pir \ resolution=2.2 dmax=20 refinement_resolution=2.2 \ cif_def_file_list=/usr/users/dale/geometry/chromophores/bcl_tnt.cif \ input_lig_file_list=AutoMR_run_4_/MR.1-Bchl-a.pdb \ rebuild_in_place=Yes My refine script is phenix.refine AutoBuild_run_12_/overall_best.pdb \ refinement.input.xray_data.file_name=1M50-2.mtz 1M50-2.mtz \ refinement.main.high_resolution=2.2 refinement.main.low_resolution=20 \ /usr/users/dale/geometry/chromophores/bcl_tnt.cif As you can guess, my test set flags are in the same mtz file as the amplitudes. I'm feeding exactly the same file into both runs. Despite this I get in my output ******************************************************************************* ******************************************************************************* The MD5 checksum for the R-free flags array summarized above is: 785fd03f6881898dcd91bc5f8c3e5b26 The corresponding MD5 checksum in the PDB file summarized above is: 6f86ee71dbcd2dc1f5282cf18547c79b These checksums should be identical but are in fact different. This is because the R-free flags used at previous stages of refinement are different from the R-free flags summarized above. As a consequence, the values for R-free will be biased and misleading. It is best to avoid this situation by consistently using the same R-free flags throughout the refinement of a model. If the previously used R-free flags are still available run this command again with the name of the file containing the original flags as an additional input. If the original R-free flags are unrecoverable, remove the REMARK r_free_flags.md5.hexdigest 6f86ee71dbcd2dc1f5282cf18547c79b record from the input PDB file to proceed with the refinement. In this case the values for R-free will become meaningful only after many cycles of refinement. ******************************************************************************* ******************************************************************************* Sorry: Please resolve the R-free flags mismatch. While I'm glad that Phenix is checking to ensure that I haven't goofed and tried to switch test sets, I believe I'm being unjustly accused. Why does phenix.refine think I'm a bad boy? Dale Tronrud

2 5

Re: [phenixbb] Autosol Feedback and questions
by Ralf W. Grosse-Kunstleve 02 Jan '08

02 Jan '08

Hi Partha, > 3. I see that phenix.refine has been used in the internal runs, I have used > phenix.refine earlier. But after playing with the side chains in Coot, when > I tried to use the modified pdb file and resolve.mtz, it complained about > the good old Free_R flag. Is this the MD5 error message? > However, it ran directly in Refmac. Is it a known > problem or maybe I did something stupid? I could do a dry run and change > the read mtz and not cns if that would solve the problem. I'm not sure what's going wrong. To be sure, I'd need a reproducer, or at least an exact description of what happened in which order. > 4. Lastly, is it possible to use multiple processors to run the job? No, sorry. Parallelization is a big, complex project. The Cambridge group has started some work in this direction, but it is still very experimental. Ralf

1 0

Fatal Error: Incomplete libtbx environment!
by Yong-Fu Li 01 Jan '08

01 Jan '08

Hi, I just installed Phenix 1.3b, however got this error message when I tried to run phenix. ******************************************* Fatal Error: Incomplete libtbx environment! ******************************************* Please re-run the libtbx/configure.py command. My system is: 2.6.23.8-34.fc7 i686 Here is what I did: 1) downloaded phenix-installer-1.3b-rc6-intel-linux-2.6.tar 2) tar -xvf ... 3) In directory: phenix-installer-1.3b/, execute ./install 4) moved phenix-1.3b-rc6/ from /usr/local to /usr/local/progs where all of my programs dwell 5) source phenix_env in: /usr/local/progs/Phenix/phenix-1.3b-rc6 6) then run: ./phenix, and got the above message 7) as root go to /libtbx, chmod +x configure.py, and ./configure.py (I don't know about python, and my Python is 2.5), it hangs there for very long (>10 min) with another error: ./configure.py: line 3: syntax error near unexpected token `(' ./configure.py: line 3: `def run():' I checked the bb for help but didn't find any useful info (not finished reading all the messages). I wonder whether someone could help! Thanks and happy new year! Yongfu

3 2

heavy atom sites in autobuild
by Bryan W. Lepore 29 Dec '07

29 Dec '07

i'm confused about treatment of heavy atom sites in autobuild. couple questions : 1. since heavy atom sites can be treated as single atoms, should phenix.autobuild simply refine the positions within a protein model? i.e. heavy atom sites should not need to be treated as ligands, with a .cif -- right? 2. if i want autobuild to build around the heavy atom site, can i just use input_pdb_file=heavy.pdb? input_lig_file_list= is not working. -bryan

2 1

Se SAD data to include Sulphurs?
by Partha Chakrabarti 28 Dec '07

28 Dec '07

Hiya, I have a SAD dataset collected at 0.979A to very good completeness and redundancy, upto 2.1A. Two molecules in AU, each containing 2 Se and the solvent content is 40%. The Signal to noise is more than 1 upto 2.1A , however solve resolve and phaser did not work. I am not an expert but what is getting somewhere is if I use SHELEX (or hyss) to find the sites, pick the best 3 of them, use MLPHARE or so to refine and Addsolve to find the other Se. However, looks like the fourth Se is less ordered. Now, the SHELEX file also finds several others which I guess are Sulphurs. I can see helix like density in the map after Addsolve (& resolve) or SHELEXE but have not been able to place any aromatic residue yet. I guess one problem is lack of proper NCS and low solvent content. Could someone point me to a method in which I can improve the map by including the Sulphurs? For example, using the already found Se as seeds? I am trying the autosol option anyway. Cheers, Partha

2 1

Re: [phenixbb] How do I gather the pieces after AutoBuild?
by Thomas C. Terwilliger 27 Dec '07

27 Dec '07

Hi Dale, If you run AutoBuild with "rebuild_in_place=Yes" then it will keep the pieces that you have, try to rebuild them (in nearby locations), and not reassemble the molecule. Rebuilding in place is supposed to be a way to finish off your structure, and that's why it doesn't move it around. In your case, you are rebuilding in place after Phaser molecular replacement. Usually Phaser will make a compact molecule during MR, so this should not normally be a problem. I gather that in this case it didn't do so, however. With rebuild_in_place=No then AutoBuild should try and make a compact molecule. If you have NCS then it will try to make a compact molecule with NCS. All the best, Tom T > Hi all, > > I have finished my Autobuild but the resulting model has some residues > that are scattered in distant asymmetric units. Coot shows their symmetry > images where they belong, next to the rest of the protein. I could > rebuild these residues by hand in the right place, or collecte together > all the symmetry operators and write a script to gather them, but being > lazy would like a more automated solution. Surely this happens to > other people and there is a tool that solves this problem. I've read > a bunch of the web pages for AutoBuild but didn't find a solution there. > > My AutoBuild command was > > phenix.autobuild model=AutoMR_run_4_/MR.1-protein.pdb data=1M50-2.mtz \ > map_file=AutoMR_run_4_/MR.MAP_COEFFS.1.mtz \ > seq_file=../fmo-ct.pir \ > resolution=2.2 dmax=20 refinement_resolution=2.2 \ > cif_def_file_list=/usr/users/dale/geometry/chromophores/bcl_tnt.cif > \ > input_lig_file_list=AutoMR_run_4_/MR.1-Bchl-a.pdb \ > rebuild_in_place=Yes > > Shouldn't AutoBuild build a compact model by default? > > Dale Tronrud > _______________________________________________ > phenixbb mailing list > phenixbb(a)phenix-online.org > http://www.phenix-online.org/mailman/listinfo/phenixbb >

2 1

How do I gather the pieces after AutoBuild?
by Dale Tronrud 23 Dec '07

23 Dec '07

Hi all, I have finished my Autobuild but the resulting model has some residues that are scattered in distant asymmetric units. Coot shows their symmetry images where they belong, next to the rest of the protein. I could rebuild these residues by hand in the right place, or collecte together all the symmetry operators and write a script to gather them, but being lazy would like a more automated solution. Surely this happens to other people and there is a tool that solves this problem. I've read a bunch of the web pages for AutoBuild but didn't find a solution there. My AutoBuild command was phenix.autobuild model=AutoMR_run_4_/MR.1-protein.pdb data=1M50-2.mtz \ map_file=AutoMR_run_4_/MR.MAP_COEFFS.1.mtz \ seq_file=../fmo-ct.pir \ resolution=2.2 dmax=20 refinement_resolution=2.2 \ cif_def_file_list=/usr/users/dale/geometry/chromophores/bcl_tnt.cif \ input_lig_file_list=AutoMR_run_4_/MR.1-Bchl-a.pdb \ rebuild_in_place=Yes Shouldn't AutoBuild build a compact model by default? Dale Tronrud

1 0

Re: [phenixbb] autosol w/ "no merge original index"?
by Tom Terwilliger 19 Dec '07

19 Dec '07

Hi Brian, Thanks for sending me the top of your hkl file....I added some random data after that and it runs on my version just fine. Can you possibly try again with the most recent version of phenix on the web site (rc6)? Thanks! -Tom T At 01:24 PM 12/19/2007, you wrote: >phenix 1.3b of 21-10-2007 > >confusion w/ "no merge original index" data file in phenix.autosol. this >should work, right? > >caveat - the 'nmoi' hkl is actually from scala, but that usually is fine >w/ e.g. solve-2.13. 'nmoi' requires declaration of cell='(numbers)' on >command line. > >phenix.autosol seems to go through all the usual checks fine, then >reports (in AutoSol_run_2_/LAST.LOG): >".... > The error seems to be at line 0 > Here is what the first few lines of this file look like: > 1 > -987 > 86.052 86.052 132.608 90.000 90.000 90.000 p41212 > You have set the keywords readdenzo and unmerged > Please check whether you should have set readformatted or readtrek or >premerged >...." > >... 'readdenzo' and 'unmerged' are just right for 'nmoi' files. etc. > >-bryan >_______________________________________________ >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] autosol w/ "no merge original index"?
by Tom Terwilliger 19 Dec '07

19 Dec '07

Hi Bryan, Yes, "No merge original index" files should work fine. This file should be converted by the autosol wizard to premerged scalepack during the run. In your case this seems to be so as the file header that is being read during the middle of the run 1 -987 86.052 86.052 132.608 90.000 90.000 90.000 p41212 is premerged scalepack. Now what does not make sense is the error message "You have set the keywords readdenzo and unmerged". I am not sure where that has come from. I just ran a SAD dataset to test this with an unmerged original index file, and it ran fine for me. If this doesn't help, could you send me (terwilliger(a)lanl.gov) your input script and datafile ( or even the first 1000 lines of the datafile) and I can try to debug it. Thanks! -Tom T At 01:24 PM 12/19/2007, you wrote: >phenix 1.3b of 21-10-2007 > >confusion w/ "no merge original index" data file in phenix.autosol. this >should work, right? > >caveat - the 'nmoi' hkl is actually from scala, but that usually is fine >w/ e.g. solve-2.13. 'nmoi' requires declaration of cell='(numbers)' on >command line. > >phenix.autosol seems to go through all the usual checks fine, then >reports (in AutoSol_run_2_/LAST.LOG): >".... > The error seems to be at line 0 > Here is what the first few lines of this file look like: > 1 > -987 > 86.052 86.052 132.608 90.000 90.000 90.000 p41212 > You have set the keywords readdenzo and unmerged > Please check whether you should have set readformatted or readtrek or >premerged >...." > >... 'readdenzo' and 'unmerged' are just right for 'nmoi' files. etc. > >-bryan >_______________________________________________ >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

autosol w/ "no merge original index"?
by Bryan W. Lepore 19 Dec '07

19 Dec '07

phenix 1.3b of 21-10-2007 confusion w/ "no merge original index" data file in phenix.autosol. this should work, right? caveat - the 'nmoi' hkl is actually from scala, but that usually is fine w/ e.g. solve-2.13. 'nmoi' requires declaration of cell='(numbers)' on command line. phenix.autosol seems to go through all the usual checks fine, then reports (in AutoSol_run_2_/LAST.LOG): ".... The error seems to be at line 0 Here is what the first few lines of this file look like: 1 -987 86.052 86.052 132.608 90.000 90.000 90.000 p41212 You have set the keywords readdenzo and unmerged Please check whether you should have set readformatted or readtrek or premerged ...." ... 'readdenzo' and 'unmerged' are just right for 'nmoi' files. etc. -bryan

1 0