- phenixbb - phenix-online.org

Auto builiding
by Jonathan Grimes 11 Sep '09

11 Sep '09

Dear Developers, I have a molecular replacement solution with a model that has 58 seq identity to my protein. There is 1 insertion in my protein of 1 residue and 1 deletion of 4 residues compared to the model. There are 2 mols in the asym unit. When I run autobuild with standard settings (but either Auto or Yes for the rebuild in place) it falls over with this error. **************************************** AutoBuild Input failed Sorry, the PDB file and sequence file could not be aligned (with no gaps and >50.0% identity) Please restart the wizard... You have several possibilities to try .... You can set rebuild_in_place=No You can set input_sequence_file=None You can edit your model to make sure it has no residues extending N-terminal of the beginning of the corresponding sequences in your sequence file You can set min_seq_identity_percent to a lower value You can set highest_resno to a higher value You can specify start_chains_list if your PDB file sequence starts with a residue number greater than the number of residues in the sequence. **************************************** It runs with rebuild_in_place=No but given the low resolution of the data (3.5A) I thought it would be best to keep as much of the model as possible. What am I doing wrong. Thanks Jon -- Dr. Jonathan M. Grimes, University Research Lecturer Division of Structural Biology Wellcome Trust Centre for Human Genetics University of Oxford Roosevelt Drive, Oxford OX3 7BN, UK Email: Jonathan(a)strubi.ox.ac.uk, Web: www.strubi.ox.ac.uk Tel: (+44) - 1865 - 287561, FAX: (+44) - 1865 - 287547

2 1

Re: [phenixbb] SA omit map
by Gino Cingolani 11 Sep '09

11 Sep '09

Dear Pascal, my suggestion to convince yourself that your ligand is there ..... 1. omit the atoms corresponding to your ligand from your pdb 2. compute electron density difference maps using data between 15-2.5A. High res maps (1.3-2.5A) will look worst than medium res maps if your ligand is partially occupied/mobile (typically bound on the surface of your proteins as opposed to a cavity). 3. Also, make sure you aren't wiping off the ligand by using an improper solvent masks during solvent flattening. If your ligand is really there, but not quite there.... you may try to pressure freeze your crystals instead of using conventional cryo-protectants. Good luck. Gino

1 0

how to generate anomalous difference map in phenix
by crystallogrphy 10 Sep '09

10 Sep '09

Hi, Does anyone know how to generate anomalous difference map by phenix? I have the data file with the scalepack format and the refined model against this data. Thanks!

3 2

Changing the resolution cutoff for water picking
by Sam Stampfer 10 Sep '09

10 Sep '09

Hi, I'm refining a 2.88 resolution structure and would like to pick waters. I assume the default low_resolution is set to 2.8 Angstroms and so my lower resolution map is not yielding any waters. How do I change the resolution cutoff for water picking to be 3.0 Angstroms? I tried adding "refinement.ordered_solvent.low_resolution=3.0" to my command line arguments and phenix is reading this, but the refinement crashes about 15 seconds in with the message: ================== Extract refinement strategy and selections ================= Sorry: Selection string 'water' results in empty selection (selects no atoms). Thanks for your help! -Sam

2 3

Ligandfit problems with ligand geometry
by Kendall Nettles 10 Sep '09

10 Sep '09

I just tried ligandfit for the first time and it did not maintain proper geometry of the ligand. There were problems with both bond lengths and planarity. Is there a way to specify a .cif file? The ligand pdb was generated with PROGRD2, and I have had good success refining a structure that contains a very close analog with phenix.refine. Kendall Nettles

2 1

Phaser and anisotropy question
by Peter Grey 10 Sep '09

10 Sep '09

Dear Phenix users, I have a very anisotropic data as phaser reports anisotropic deltaB = 60.2. I would be grateful for advice of several issues. 1.Could you please tell me if phaser map coefficients FWT,PHWT take into account the anisotropic scaling ? 2.This means that these coefficients will be different from those calculated from a partial model in sigmaa because sigmaa has no anisotropic scaling (and no bulk solvent correction) ? 3.In the case of such severe anisotropy can the scaling diminish too strongly the well measured high resolution reflections ? If so should I calculate the coefficients my self by sigmaa and not use pahser mtz output or is there a better solution ? Many thanks in advance for sharing your thoughts and experience, Peter.

4 5

Re: [phenixbb] Phenix crashes
by Katya Heldwein 10 Sep '09

10 Sep '09

Nat, 1.4-160 works just fine. Thank you. Katya Ekaterina Heldwein, Ph.D. Assistant Professor Department of Molecular Biology and Microbiology Tufts University School of Medicine 136 Harrison Ave Boston, MA 02111 phone: 617-636-0858 fax: 617-636-0337 e-mail: katya.heldwein(a)tufts.edu phenixbb-request(a)phenix-online.org wrote: > Send phenixbb mailing list submissions to > phenixbb(a)phenix-online.org > > To subscribe or unsubscribe via the World Wide Web, visit > http://www.phenix-online.org/mailman/listinfo/phenixbb > or, via email, send a message with subject or body 'help' to > phenixbb-request(a)phenix-online.org > > You can reach the person managing the list at > phenixbb-owner(a)phenix-online.org > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of phenixbb digest..." > > > Today's Topics: > > 1. Phenix crashes (Katya Heldwein) > 2. Re: Phenix crashes (sbiswas2(a)ncsu.edu) > 3. Re: Phenix crashes (Nathaniel Echols) > 4. Re: Phenix crashes (Nathaniel Echols) > 5. specifying column labels for input MTZ file > (Tirumala Kumar Chowdary) > 6. Re: specifying column labels for input MTZ file > (Schubert, Carsten [PRDUS]) > 7. Re: specifying column labels for input MTZ file (Pavel Afonine) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 09 Sep 2009 12:02:36 -0400 > From: Katya Heldwein <Katya.Heldwein(a)tufts.edu> > Subject: [phenixbb] Phenix crashes > To: phenixbb(a)phenix-online.org > Message-ID: <20090909120236.c2mlsf9lww00kwsk(a)webmail.tufts.edu> > Content-Type: text/plain; charset=ISO-8859-1; format="flowed" > > Hi, > > I am running phenix-1.4-159 on an Intel Mac. I get the following error: > > Bus error (Python call stack above) > This crash may be due to a problem in any imported > Python module, including modules which are not part > of the cctbx project. To disable the traps leading > to this message, define these environment variables > (e.g. assign the value 1): > BOOST_ADAPTBX_FPE_DEFAULT > BOOST_ADAPTBX_SIGNALS_DEFAULT > This will NOT solve the problem, just mask it, but > may allow you to proceed in case it is not critical. > > If I fix it by > > setenv BOOST_ADAPTBX_FPE_DEFAULT 1 > > I still get the same error. > > If I also do > > setenv BOOST_ADAPTBX_SIGNALS_DEFAULT 1 > > I don't get the error but phenix crashes nonetheless. > > How do I get around this problem? > > Thanks, > > Katya > > > > > ------------------------------ > > Message: 2 > Date: Wed, 9 Sep 2009 12:17:36 -0400 (EDT) > From: sbiswas2(a)ncsu.edu > Subject: Re: [phenixbb] Phenix crashes > To: "PHENIX user mailing list" <phenixbb(a)phenix-online.org> > Message-ID: <1476.152.1.132.71.1252513056.squirrel(a)webmail.ncsu.edu> > Content-Type: text/plain;charset=iso-8859-1 > > Hi, > This looks like a version problem, try installing the latest nightly builds. > Shya > > > >> Hi, >> >> I am running phenix-1.4-159 on an Intel Mac. I get the following error: >> >> Bus error (Python call stack above) >> This crash may be due to a problem in any imported >> Python module, including modules which are not part >> of the cctbx project. To disable the traps leading >> to this message, define these environment variables >> (e.g. assign the value 1): >> BOOST_ADAPTBX_FPE_DEFAULT >> BOOST_ADAPTBX_SIGNALS_DEFAULT >> This will NOT solve the problem, just mask it, but >> may allow you to proceed in case it is not critical. >> >> If I fix it by >> >> setenv BOOST_ADAPTBX_FPE_DEFAULT 1 >> >> I still get the same error. >> >> If I also do >> >> setenv BOOST_ADAPTBX_SIGNALS_DEFAULT 1 >> >> I don't get the error but phenix crashes nonetheless. >> >> How do I get around this problem? >> >> Thanks, >> >> Katya >> >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://www.phenix-online.org/mailman/listinfo/phenixbb >> >> > > > > ------------------------------ > > Message: 3 > Date: Wed, 9 Sep 2009 09:37:50 -0700 > From: Nathaniel Echols <NEchols(a)lbl.gov> > Subject: Re: [phenixbb] Phenix crashes > To: PHENIX user mailing list <phenixbb(a)phenix-online.org> > Message-ID: <2C512891-F9F0-40FB-AC0A-FA8C810F65AF(a)lbl.gov> > Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes > > On Sep 9, 2009, at 9:02 AM, Katya Heldwein wrote: > >> I am running phenix-1.4-159 on an Intel Mac. I get the following >> error: >> Bus error (Python call stack above) >> >> If I fix it by >> setenv BOOST_ADAPTBX_FPE_DEFAULT 1 >> I still get the same error. >> If I also do >> >> setenv BOOST_ADAPTBX_SIGNALS_DEFAULT 1 >> >> I don't get the error but phenix crashes nonetheless. >> >> How do I get around this problem? >> > > > Which programs are you running, and did you see a stack trace before > the error message? I just noticed the same thing running > phenix.refine from the same installer on my Mac, so it is probably an > undiagnosed bug. > > ------------------- > Nathaniel Echols > Lawrence Berkeley Lab > 510-486-5136 > NEchols(a)lbl.gov > > > > > > ------------------------------ > > Message: 4 > Date: Wed, 9 Sep 2009 09:46:24 -0700 > From: Nathaniel Echols <NEchols(a)lbl.gov> > Subject: Re: [phenixbb] Phenix crashes > To: PHENIX user mailing list <phenixbb(a)phenix-online.org> > Message-ID: <6207C287-B0BB-4A7B-AD8E-2290222334FE(a)lbl.gov> > Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes > > On Sep 9, 2009, at 9:17 AM, sbiswas2(a)ncsu.edu wrote: > >> Hi, >> This looks like a version problem, try installing the latest nightly >> builds. >> Shya >> > > > Actually, 1.4-159 was the latest build until ten minutes ago, but I > just checked 1.4-160 and it does not appear to have this problem, so > it's now ready to download. Sorry for the trouble. > > -Nat > > ------------------- > Nathaniel Echols > Lawrence Berkeley Lab > 510-486-5136 > NEchols(a)lbl.gov > > > > > > ------------------------------ > > Message: 5 > Date: Wed, 09 Sep 2009 14:10:27 -0400 > From: Tirumala Kumar Chowdary <Tirumala.Chowdary(a)tufts.edu> > Subject: [phenixbb] specifying column labels for input MTZ file > To: phenixbb(a)phenix-online.org > Message-ID: <20090909141027.belblmob7s4s440w(a)webmail.tufts.edu> > Content-Type: text/plain; charset=ISO-8859-1; format="flowed" > > Hi, > > I am trying to run phenix.refine with a sharp output MTZ file (column > labels H K L NUMORB FB PHIB FOM HX HY FH PHIH FP SIGFP HLA HLB HLC HLD > HL0 Fcent PHIcent FOMcent Fmap PHImap FHref PHIHref Fav SIGFav). > > Phenix.refine stopped at > > 'Multiple equally suitable arrays of observed xray data found. > > Possible choices: > sharp_sad_se_eden.mtz:FP,SIGFP > sharp_sad_se_eden.mtz:Fav,SIGFav > > Please use refinement.input.xray_data.labels' > > Can someone tell me on how to use this > 'refinement.input.xray_data.labels' and how to specify which column > labels to use. There is nothing written in the documentation about it. > > Thanks > Tirumal > > > > > > > ------------------------------ > > Message: 6 > Date: Wed, 9 Sep 2009 14:22:40 -0400 > From: "Schubert, Carsten [PRDUS]" <CSCHUBER(a)its.jnj.com> > Subject: Re: [phenixbb] specifying column labels for input MTZ file > To: "PHENIX user mailing list" <phenixbb(a)phenix-online.org> > Message-ID: > <204D6E1BD9ACEA40B3F76FE9C125528F035DE9FF(a)JNJUSSHGMS01.na.jnj.com> > Content-Type: text/plain; charset="us-ascii" > > Tirumal: > > refinement.input.xray_data.labels="FP,SIGFP" > > should do the trick. > > Cheers, > > Carsten > > > >> -----Original Message----- >> From: phenixbb-bounces(a)phenix-online.org [mailto:phenixbb- >> bounces(a)phenix-online.org] On Behalf Of Tirumala Kumar Chowdary >> Sent: Wednesday, September 09, 2009 2:10 PM >> To: phenixbb(a)phenix-online.org >> Subject: [phenixbb] specifying column labels for input MTZ file >> >> Hi, >> >> I am trying to run phenix.refine with a sharp output MTZ file (column >> labels H K L NUMORB FB PHIB FOM HX HY FH PHIH FP SIGFP HLA HLB HLC >> > HLD > >> HL0 Fcent PHIcent FOMcent Fmap PHImap FHref PHIHref Fav SIGFav). >> >> Phenix.refine stopped at >> >> 'Multiple equally suitable arrays of observed xray data found. >> >> Possible choices: >> sharp_sad_se_eden.mtz:FP,SIGFP >> sharp_sad_se_eden.mtz:Fav,SIGFav >> >> Please use refinement.input.xray_data.labels' >> >> Can someone tell me on how to use this >> 'refinement.input.xray_data.labels' and how to specify which column >> labels to use. There is nothing written in the documentation about it. >> >> Thanks >> Tirumal >> >> >> >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://www.phenix-online.org/mailman/listinfo/phenixbb >> > > > > > ------------------------------ > > Message: 7 > Date: Wed, 09 Sep 2009 11:24:46 -0700 > From: Pavel Afonine <PAfonine(a)lbl.gov> > Subject: Re: [phenixbb] specifying column labels for input MTZ file > To: PHENIX user mailing list <phenixbb(a)phenix-online.org> > Cc: Tirumala.Chowdary(a)tufts.edu > Message-ID: <4AA7F2EE.6080106(a)lbl.gov> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > Hi Tirumala, > > I agree, you are not the first who confused by this and an example right > in the message would help. > > First, your MTZ file contains the following data: (FP,SIGFP) and > (Fav,SIGFav) and it is your choice to select the one that you want to > use in refinement. > > Second, once decided on the data - for example it is FP, then just add > to your command line this: > > refinement.input.xray_data.labels=FP > > or, which is the same > > xray_data.labels="FP,SIGFP" > > Pavel. > > > > On 9/9/09 11:10 AM, Tirumala Kumar Chowdary wrote: > >> Hi, >> >> I am trying to run phenix.refine with a sharp output MTZ file (column >> labels H K L NUMORB FB PHIB FOM HX HY FH PHIH FP SIGFP HLA HLB HLC HLD >> HL0 Fcent PHIcent FOMcent Fmap PHImap FHref PHIHref Fav SIGFav). >> >> Phenix.refine stopped at >> >> 'Multiple equally suitable arrays of observed xray data found. >> >> Possible choices: >> sharp_sad_se_eden.mtz:FP,SIGFP >> sharp_sad_se_eden.mtz:Fav,SIGFav >> >> Please use refinement.input.xray_data.labels' >> >> Can someone tell me on how to use this >> 'refinement.input.xray_data.labels' and how to specify which column >> labels to use. There is nothing written in the documentation about it. >> >> Thanks >> Tirumal >> >> >> >> >> _______________________________________________ >> 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 > > > End of phenixbb Digest, Vol 46, Issue 14 > **************************************** >

4 4

run autosol with rough phases from shelx?
by Francis E Reyes 09 Sep '09

09 Sep '09

Hi all Is it possible to input the phs file from a shelx run into an autosol run for some phase improvement with solve/resolve?? I tried putting in the HA sites but then the resultant map is worse. Also what are acceptable skews for low resolution phasing (worse than 4A) ? I understand the docs say 0.1 but does that also apply for low res? Thanks FR --------------------------------------------- Francis Reyes M.Sc. 215 UCB University of Colorado at Boulder gpg --keyserver pgp.mit.edu --recv-keys 67BA8D5D 8AE2 F2F4 90F7 9640 28BC 686F 78FD 6669 67BA 8D5D

1 0

Phenix crashes
by Katya Heldwein 09 Sep '09

09 Sep '09

Hi, I am running phenix-1.4-159 on an Intel Mac. I get the following error: Bus error (Python call stack above) This crash may be due to a problem in any imported Python module, including modules which are not part of the cctbx project. To disable the traps leading to this message, define these environment variables (e.g. assign the value 1): BOOST_ADAPTBX_FPE_DEFAULT BOOST_ADAPTBX_SIGNALS_DEFAULT This will NOT solve the problem, just mask it, but may allow you to proceed in case it is not critical. If I fix it by setenv BOOST_ADAPTBX_FPE_DEFAULT 1 I still get the same error. If I also do setenv BOOST_ADAPTBX_SIGNALS_DEFAULT 1 I don't get the error but phenix crashes nonetheless. How do I get around this problem? Thanks, Katya

3 3

question about pick water in autobuild
by crystallogrphy 08 Sep '09

08 Sep '09

Hi, I am building a model against a 2.1A data set using phenix.autobuild. Right now Rfree is about 35%. I do not want to pick water in this stage because I see denisty at N-terminal of my model which need to build more residues. However, the program still picks water even after I say 'place_waters=False'. Does anyone know why? phenix.autobuild data=exptl_fobs_phases_freeR_flags.mtz model=AutoBuild_run_4_/overall_best-090905a.pdb seq_file=../ztb-sd.seq refinement.place_waters=False refinement.refine_with_ncs=True nproc=2

5 6

RuntimeError
by Young-Jin Cho 08 Sep '09

08 Sep '09

Hi everyone, while I was running phenix, I got the following error message. Could some one give me any tips to deal with such this problem? Thanks in advace, YoungJin Traceback (most recent call last): File "/home/yjcho2/Phenix/phenix-1.4-95/phenix/phenix/command_line/refine.py", line 19, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/home/yjcho2/Phenix/phenix-1.4-95/phenix/phenix/refinement/command_line.py", line 68, in run log=log) File "/home/yjcho2/Phenix/phenix-1.4-95/phenix/phenix/refinement/command_line.py", line 227, in __init__ self.pdb_file_names) File "/home/yjcho2/Phenix/phenix-1.4-95/cctbx_project/mmtbx/utils.py", line 876, in process_pdb_files raw_records = raw_records) File "/home/yjcho2/Phenix/phenix-1.4-95/cctbx_project/mmtbx/utils.py", line 907, in _process_pdb_file log = self.log) File "/home/yjcho2/Phenix/phenix-1.4-95/cctbx_project/mmtbx/monomer_library/pdb_interpretation.py", line 2912, in __init__ log=log) File "/home/yjcho2/Phenix/phenix-1.4-95/cctbx_project/mmtbx/monomer_library/pdb_interpretation.py", line 1962, in __init__ self.site_symmetry_table().special_position_indices() File "/home/yjcho2/Phenix/phenix-1.4-95/cctbx_project/mmtbx/monomer_library/pdb_interpretation.py", line 2217, in site_symmetry_table sites_cart=self.sites_cart) File "/home/yjcho2/Phenix/phenix-1.4-95/cctbx_project/cctbx/crystal/__init__.py", line 394, in site_symmetry_table assert_min_distance_sym_equiv=self.assert_min_distance_sym_equiv()) RuntimeError: cctbx Error: site_symmetry: min_distance_sym_equiv too large. --

2 1

Re: [phenixbb] question about pick water in autobuild
by Thomas C. Terwilliger 08 Sep '09

08 Sep '09

Hi Qun, With your command you need to add one more keyword to get rid of all the waters: place_waters_in_combine=False I'm sorry that doesn't happen automatically if you say place_waters=False ...I will make that happen in the future. Note however that these waters should not interfere with building...as waters are always deleted prior to building model. I note also that your command does not add new residues; it only adjusts the ones that are there. To build new model add: rebuild_in_place=False All the best, -Tom T >> Hi, >> >> from what you explain, the main issue I see here is that the model building program does not treat water simply as density peaks, and so removes them (waters) if necessary to build a protein chain instead. >> >> Otherwise, it's probably a good idea in general to fill in all unmodeled >> places with some placeholders like water, or probably isolated N atoms would work even better (if I recall correctly). >> >> Pavel. >> >> >> On 9/7/09 6:18 AM, crystallogrphy wrote: >>> Hi, Pavel: >>> The problem here is that autobuild placed waters into the densities where they are the N-terminal residues. Without the waters here, autobuild will build residues, not waters, to decrease R and Rfree with less model bias. Now with water picking, autobuild does not build more residues, just water. >>> I still do not know why 'place_waters=False' cannot work. Tom has not given me the answer yet. >>> On Mon, Sep 7, 2009 at 1:45 AM, Pavel Afonine <PAfonine(a)lbl.gov <mailto:[email protected]>> wrote: >>> Hi, >>> I hope Tom gives you the answer about how to stop water picking in autobuild. >>> From my side, I'm wondering why it bothers you (I'm not saying that not working option is good)? Why don't you consider those waters as a tool for improving your map? You are free to remove any of them them later on. This can be viewed as a prototype of what ARP/wARP does, or other "density growing" techniques... In your case you just temporarily place waters where the model is still absent, and when you build those missing places just consider waters as density peaks (I hope model building programs assume waters as just density peaks and remove waters if they are on the macromolecule's way). On the other hand, model >>> incompleteness may make the map ugly even for relatively well defined places. >>> Pavel. >>> On 9/6/09 7:38 AM, crystallogrphy wrote: >>>> Hi, >>>> I am building a model against a 2.1A data set using >>>> phenix.autobuild. Right now Rfree is about 35%. I do not want to pick water in this stage because I see denisty at N-terminal of my model which need to build more residues. However, the program still picks water even after I say 'place_waters=False'. Does anyone know why? >>>> phenix.autobuild data=exptl_fobs_phases_freeR_flags.mtz >>>> model=AutoBuild_run_4_/overall_best-090905a.pdb >>>> seq_file=../ztb-sd.seq refinement.place_waters=False >>>> refinement.refine_with_ncs=True nproc=2 >>>> ------------------------------------------------------------------------ _______________________________________________ >>>> phenixbb mailing list >>>> phenixbb(a)phenix-online.org <mailto:[email protected]> http://www.phenix-online.org/mailman/listinfo/phenixbb >>

1 0

about kick map generation
by crystallogrphy 04 Sep '09

04 Sep '09

Hi, I am generating a averaged kick maps with the following script. According to the script, there should be both kick map and normal 2fofc maps. However, there is only one map files created. I am not sure it is kick map or routine 2fofc map. Does anyone knows? electron_density_maps { map_format = *xplor map_coefficients_format = *mtz phs suppress = None map { mtz_label_amplitudes = "2FOFCWT" mtz_label_phases = "PH2FOFCWT" likelihood_weighted = True obs_factor = 2 calc_factor = 1 kicked = False fill_missing_f_obs_with_weighted_f_model = True } map { mtz_label_amplitudes = "FOFCWT" mtz_label_phases = "PHFOFCWT" likelihood_weighted = True obs_factor = 1 calc_factor = 1 kicked = False fill_missing_f_obs_with_weighted_f_model = True } map { mtz_label_amplitudes = "2FOFCWT_kick" mtz_label_phases = "PH2FOFCWT_kick" likelihood_weighted = True obs_factor = 2 calc_factor = 1 kicked = True fill_missing_f_obs_with_weighted_f_model = False } map { mtz_label_amplitudes = "FOFCWT_kick" mtz_label_phases = "PHFOFCWT_kick" likelihood_weighted = True obs_factor = 1 calc_factor = 1 kicked = True fill_missing_f_obs_with_weighted_f_model = False } map { mtz_label_amplitudes = "2FOFCWT_no_fill" mtz_label_phases = "PH2FOFCWT_no_fill" likelihood_weighted = True obs_factor = 2 calc_factor = 1 kicked = False fill_missing_f_obs_with_weighted_f_model = False } map { mtz_label_amplitudes = "FOFCWT_no_fill" mtz_label_phases = "PHFOFCWT_no_fill" likelihood_weighted = True obs_factor = 1 calc_factor = 1 kicked = False fill_missing_f_obs_with_weighted_f_model = False } anomalous_difference_map { mtz_label_amplitudes = "ANOM" mtz_label_phases = "PHANOM" } grid_resolution_factor = 0.25 region = *selection cell atom_selection = None atom_selection_buffer = 3 apply_sigma_scaling = True apply_volume_scaling = False }

2 1

assertion error with autobuild (phenix 1.4.3)
by Jan Stout 04 Sep '09

04 Sep '09

Dear, I am trying to run the Autobuild Wizard to make a iterative-build omit map. I get the following errors after running Autobuild using the input file autobuild.eff generated by the wizard itself. I am working with phenix 1.4.3. Any input will be greatly appreciated. Many thanks, Jan Stout Traceback (most recent call last): File "/usr/local/PHENIX/phenix-1.4.3/phenix-1.4-3/phenix/phenix/command_line/autobuild.py", line 339, in <module> run_autobuild=run_autobuild(args=sys.argv[1:]) File "/usr/local/PHENIX/phenix-1.4.3/phenix-1.4-3/phenix/phenix/autosol/wizard_command_line.py", line 598, in __init__ self.RunWizard() File "/usr/local/PHENIX/phenix-1.4.3/phenix-1.4-3/phenix/phenix/autosol/wizard_command_line.py", line 973, in RunWizard self.runWizard.run() File "/usr/local/PHENIX/phenix-1.4.3/phenix-1.4-3/phenix/phenix/autosol/wizard_command_line.py", line 230, in run self.Wizard.Run(info='script') File "/usr/local/PHENIX/phenix-1.4.3/phenix-1.4-3/phenix/phenix/autosol/AutoBaseExtend.py", line 726, in Run ok=self.CollectInputs()# raises AssertionError on failure File "/usr/local/PHENIX/phenix-1.4.3/phenix-1.4-3/phenix/phenix/autosol/AutoBaseExtend.py", line 1790, in CollectInputs getattr(self,"finish_"+str(self.application_method))() File "/usr/local/PHENIX/phenix-1.4.3/phenix-1.4-3/phenix/phenix/wizards/AutoBuild.py", line 1074, in finish_AutoBuild_start all_files_ok=self.finish_Build_start() File "/usr/local/PHENIX/phenix-1.4.3/phenix-1.4-3/phenix/phenix/autosol/BuildBase.py", line 382, in finish_Build_start all_files_ok=self.Build_finish_up_start() File "/usr/local/PHENIX/phenix-1.4.3/phenix-1.4-3/phenix/phenix/autosol/BuildBase.py", line 508, in Build_finish_up_start str(self.Facts['input_data_file'])+"'") File "/usr/local/PHENIX/phenix-1.4.3/phenix-1.4-3/phenix/phenix/autosol/UserMethods.py", line 84, in Error raise AssertionError,self.ErrorLine AssertionError: Sorry, cannot find the input_data_file ''

2 1

question about phenix.refine
by crystallogrphy 03 Sep '09

03 Sep '09

Hi, I want to use phenix to refine a model right after autobuild with the following script: phenix.refine xray_data.file_name=exptl_fobs_phases_freeR_flags.mtz overall_best-090903a.pdb strategy=rigid_body+individual_sites simulated_annealing=true simulated_annealing.start_temperature=8000 sites.rigid_body="chain A" sites.rigid_body="chain B" However, the program says " F-obs: exptl_fobs_phases_freeR_flags.mtz:FP,SIGFP If previously used R-free flags are available run this command again with the name of the file containing the original flags as an additional input. If the structure was never refined before, or if the original R-free flags are unrecoverable, run this command again with the additional definition: refinement.input.xray_data.r_free_flags.generate=True when I use phenix.mtz.dump exptl_fobs_phases_freeR_flags.mtz, it indicates that it has FreeR_flag=1. Can anyone tell me why there is FreeR_flag in my data, while the refinement program does not recognize it? label #valid %valid min max type H 35619 100.00% 0.00 22.00 H: index h,k,l K 35619 100.00% 0.00 45.00 H: index h,k,l L 35619 100.00% -33.00 33.00 H: index h,k,l FP 35619 100.00% 0.00 114.22 F: amplitude SIGFP 35619 100.00% 0.15 6.13 Q: standard deviation PHIM 35619 100.00% 0.00 0.00 P: phase angle in degrees FOMM 35619 100.00% 0.00 0.00 W: weight (of some sort) FreeR_flag 35619 100.00% 0.00 1.00 I: integer FC 35619 100.00% 0.00 0.00 F: amplitude

2 1

target omit map
by Susan Tsutakawa 03 Sep '09

03 Sep '09

I am trying to generate a simulated annealing omit map of a Mg-water cluster in the active site of my protein but when I run the SA omit map, I can only see the Mg when I use the Mgwater cluster as the target pdb. When I was building the model which did not have the metal present (phases solved by MR), I observed this Mg-water cluster in the Fo-Fc maps. Geometry and distance were ideal, etc. Now that I'm at the final stages of refinement, I would like to do a SA omit map as a figure for the paper. However, if I goto the lowest sigma, I only see the Mg. I'm using coot with the weights included to view the SA omit maps. I noticed that there is very little noise at the lowest sigma and so I'm wondering if the 0 sigma is truly 0 sigma and any suggestions on what to do. Thank you in advance. Sincerely, Susan

2 1

reflection file utils : .cif?
by Bryan W. Lepore 03 Sep '09

03 Sep '09

[ phenix 1.4r3 cctbx tag: 2008_12_07_1353 intel-linux-2.6 redhat-e4.7 ] i obtained a recent .cif from the PDB, then run : phenix.reflection_file_converter --mtz=output.mtz test.cif file_name: test.cif file_type: None No reflection data found in input file. .... not sure where to go from here (a clip of the file is below, hope not too much) - i know i can run ccp4:cif2mtz, but are .cif part of phenix.reflection_file_converter? -bryan (sample of test.cif): data_testsf # _diffrn_radiation.id 1 _diffrn_radiation.pdbx_wavelength 1.00600 # _audit.revision_id 1_0 _audit.creation_date 2009-07-07 _audit.update_record 'Initial release' # _symmetry.Int_Tables_number 1 _symmetry.space_group_name_H-M 'P 1' # _cell.length_a 79.6610 _cell.length_b 104.1490 _cell.length_c 154.0250 _cell.angle_alpha 81.9570 _cell.angle_beta 75.9250 _cell.angle_gamma 73.7280 # _reflns.d_resolution_high 2.967 _reflns.d_resolution_low 46.392 _reflns.limit_h_max 26 _reflns.limit_h_min -25 _reflns.limit_k_max 34 _reflns.limit_k_min -34 _reflns.limit_l_max 51 _reflns.limit_l_min 0 _reflns.number_all 93210 _reflns.number_obs 85990 # loop_ _symmetry_equiv.id _symmetry_equiv.pos_as_xyz 1 'X, Y, Z' # loop_ _refln.crystal_id _refln.wavelength_id _refln.scale_group_code _refln.status _refln.index_h _refln.index_k _refln.index_l _refln.F_meas_au _refln.F_meas_sigma_au 1 1 1 o -25 -16 1 218.4 77.9 1 1 1 o -25 -15 1 250.7 87.9 1 1 1 o -25 -15 2 217.9 78.8 (... etc... here is a free set hkl :) 1 1 1 f -25 -7 3 315.6 89.7 (... etc... until the end :) #END OF REFLECTIONS

2 1

Custom Planar Restraints
by cheung＠lmb.uni-muenchen.de 03 Sep '09

03 Sep '09

Hi all - is it possible to specify custom planar restraints in phenix.refine? I'm refining a 3.8A protein-dna complex and although i've set hydrogen bonding restraints between base pairs, they still become very twisted and buckled after refinement. Alan

2 1

question about autobuild
by crystallogrphy 03 Sep '09

03 Sep '09

Hi, I have used phenix.automr followed by phenix.autobuild to build a 264 residue model against my 2.1A resolution data set. The first run of autobuild gives me about 180 residues built. I have adjusted the model according to the denity. The N-terminal of the model has OK density to add more residues to the model. Can phenix.autobuild has some command to just put more residues to the N-terminal of my current model according to the density, so that the cost of CPU time will be much less than rerun the job completely? Thanks!

2 1

Chain id of water is not reading by coot
by Raja Dey 02 Sep '09

02 Sep '09

Dear Friends, I have 7 chains in my pdb file. When I am reading the pdb file in COOT, for some reason the chain id(W) of water is not coming in the list. I clicked on 'draw' and then 'go to atom'. I can see all other chains but not W. I created water in PHENIX Is there any wrong in my pdb? HETATM 7562 O HOH W 542 -9.043 -2.133 -12.633 1.00 34.55 W O HETATM 7563 O HOH W 543 -22.004 8.228 -21.695 1.00 34.55 W O HETATM 7564 O HOH W 544 5.926 38.977 -11.538 1.00 34.55 W O HETATM 7565 O HOH W 545 -31.162 5.813 -31.403 1.00 34.55 W O Raja

2 1

Improving picture quality of ED map in PYMOL
by Raja Dey 01 Sep '09

01 Sep '09

Dear Friends, Sorry, for the of topic question. Does anyone can tell how to improve the quality of a picture showing the electron density map of dna in a complex using PYMOL. I created the figure, but not happy with the quality as it is now. Might be there is a way to improve which I don't know. I want the lines in the mesh to be more sharp and fine. Any suggested program to prepare a publication quality figure to show the ED map would be of great help. Thanking you in advance... Raja

6 5

Source for Nightlys Location?
by Andrew Ring 01 Sep '09

01 Sep '09

Is it possible to get a source distribution for the nightlys? If so, where is it located? I can only finding the builds on the nightlys downloads page. For my systems and uses, the --openmp switch on version 1.4-3 has given encouraging results and I would like to investigate with something newer. -- =================================== Andrew Ring System Administrator Kuriyan Laboratory http://jkweb.qb3.berkeley.edu/ Doudna Laboratory http://rna.berkeley.edu/ University of California, Berkeley Office: 542 Stanley Hall Shipping: 176 Stanley Hall, QB3 Berkeley, CA 94720-3220 tel: (510) 643 0166 fax: (510) 643 2352 ====================================

2 2

R-work is too low
by Sam Stampfer 01 Sep '09

01 Sep '09

Hi, I'm refining a 2.88 Angstrom structure and am getting a difference of greater than 6% between the R-work and R-free. The best result (I think) is what I have from minimization alone: Start R-work = 0.3219, R-free = 0.3505 (no bulk solvent and anisotropic scale) Final R-work = 0.2424, R-free = 0.3174 (no bulk solvent and anisotropic scale) Start R-work = 0.3220, R-free = 0.3500 Final R-work = 0.2411, R-free = 0.3176 The results from simulated annealing are the following: Start R-work = 0.3219, R-free = 0.3505 (no bulk solvent and anisotropic scale) Final R-work = 0.2454, R-free = 0.3322 (no bulk solvent and anisotropic scale) Start R-work = 0.3220, R-free = 0.3500 Final R-work = 0.2403, R-free = 0.3333 I have also done the refinement with ncs but ncs is probably not appropriate for my structure. The final R-free ends up being higher (though the difference between R-free and R-work ends up being 6-7%) My PI (who is on vacation until next week) said that it was a problem if R-work was more than 0.06 lower than R-free. Any recommendations for me? Is it possible to reduce this difference? Is it worth it to rebuild the refined pdb if the R-work is this different from the R-free? My refinement script is: phenix.refine ss46_p3_cv.mtz ss46_rebuilt.pdb \ output.prefix=results/ss46_rebuilt_sa \ strategy=individual_sites+group_adp cif_link.params nag-nag.cif \ twin_law="h,-h-k,-l" \ simulated_annealing=true (the cif_link.params and nag-nag.cif files are made to handle the N-glycosylation sites, I think). Thanks for your help! -Sam

3 2

error TLS refinement
by Marjolein Thunnissen 01 Sep '09

01 Sep '09

I'm trying to finish refining a structure that has been almost refined completely about a year ago with an older version of Phenix.refine. When I start TLS refinement I get the following error message: RuntimeError : mean_weighted() argument is an empty array Traceback: Traceback (most recent call last): File "/usr/local/phenix-1.4-3/phenix/wxGUI2/Refine.py", line 259, in OnExcept raise event.data RuntimeError: mean_weighted() argument is an empty array I do not get this error when refining without TLS, only when I switch it on. The program advises me to post the error here so... Thanks Marjolein ________________________________________________________________________ Marjolein Thunnissen Phone +46-(0)46-22 24584 Associate Professor Fax +46-(0)46-22 24692 Dept of Molecular Biophysics, Lund University http://www.mbfys.lu.se PO-Box 124 S-221 00 Lund, Sweden Scientific coordinator I911 (Max-lab): MAD and fixed-wavelength stations for macromolecular crystallography ________________________________________________________________________

2 1

Specifying the number of processors to run phenix with
by Sam Stampfer 31 Aug '09

31 Aug '09

Hi- anyone know how to specify the number of processors to run phenix under in linux? Or how to tell what the default is? I'm running linux (fedora) with 4 processors and can't tell whether phenix is using all of them or not. I do the refinement with .csh files and tried doing "myfile.csh --nproc=4". The file ran but no faster- I can't tell if it is using 4 by default or if that's the incorrect command to specify the use of 4 processors. Thanks for your help! -Sam

3 3