- phenixbb - phenix-online.org

norleucine (NLE) restraints
by Daniel Keedy 01 Dec '14

01 Dec '14

Hi everybody, First-time poster! I'm trying to re-refine 2f4k.pdb, but phenix.refine chokes on the non-standard norleucine (NLE) amino acid. I've tried phenix.elbow and phenix.ready_set using different version of PHENIX, with the NLE cards as either HETATM or ATOM, all to no avail. SHELX was used to refine the structure originally, but I'd like to use PHENIX to re-refine it. Any ideas on how to define these restraints? Daniel Keedy James Fraser Lab UCSF

3 3

来自"张士军" <[email protected]>的邮件
by 张士军 01 Dec '14

01 Dec '14

hello everyone how could use phenix phasing to phase a se substructure by search any possible spacegroup

1 0

Generate table 1-error message
by Lau Sze Yi (SIgN) 30 Nov '14

30 Nov '14

Hi, I got an error message when I tried to run the generate table 1 utility: Crystal symmetry mismatch between different files. My input files for pdb, mtz and cif were generated from my ccp4 refmac run. The unmerged data is the mtz output of ccp4 aimless. Did I upload the wrong file? Regards, Sze Yi Research Fellow Singapore Immunology Network (SigN) 8A Biomedical Grove, #03-06 Immunos, Singapore 138648 DID: +65 6407 0640 Fax: +65 6407 0640

2 1

Not a suitable array of R-free flags:
by Yamei Yu 26 Nov '14

26 Nov '14

Dear all, I use phenix.refine to refine a structure. Firstly, I used --dryrun to generate an .eff file. The command is "phenix.refine --dry_run 7_44.mtz 7_44.pdb refinement.input.xray_data.r_free_flags.generate=True refinement.input.xray_data.r_free_flags.max_free=1000” Then I use run the program by typing “phenix.refine 7_44_refine_001.eff” It runs correct! After that I use coot to modulate some parts of the structure, then I want to refine the structure again. I type “ phenix.refine 7_44_refine_002.def” This time I got the following error message: Not a suitable array of R-free flags: refinement.input.xray_data.r_free_flags.label=R-free-flags Sorry_Not_a_suitable_array: Not a suitable array of R-free flags: refinement.input.xray_data.r_free_flags.label=R-free-flags To override the suitability test define: refinement.input.xray_data.r_free_flags.disable_suitability_test=True I checked the generated file by phenix “7_44_refine_data.mtz” I does contain a column labeled as R-free-flags. I paste part if it below. * Space group = 'P1' (number 1) OVERALL FILE STATISTICS for resolution range 0.000 - 0.277 ======================= Col Sort Min Max Num % Mean Mean Resolution Type Column num order Missing complete abs. Low High label 1 ASC -32 31 0 100.00 -2.2 11.9 49.44 1.90 H H 2 NONE -46 45 0 100.00 -2.5 17.8 49.44 1.90 H K 3 NONE 0 49 0 100.00 18.8 18.8 49.44 1.90 H L 4 NONE 0.9 271.3 0 100.00 26.52 26.52 49.44 1.90 F F-obs 5 NONE 0.2 5.1 0 100.00 1.25 1.25 49.44 1.90 Q SIGF-obs 6 NONE 0.0 1.0 0 100.00 0.01 0.01 49.44 1.90 I R-free-flags No. of reflections used in FILE STATISTICS 143215 Then I regenerate a new mtz file with FreeRflag by XDSCONT, then I revise the .def file accordingly and rerun the program. I got same error message:Not a suitable array of R-free flags: refinement.input.xray_data.r_free_flags.label=FreeRflag Sorry_Not_a_suitable_array: Not a suitable array of R-free flags: refinement.input.xray_data.r_free_flags.label=FreeRflag To override the suitability test define: refinement.input.xray_data.r_free_flags.disable_suitability_test=True I checked my new mtz file, it dose contain a column labeled as FreeRflag Why there is an error message? How to solve it? Thank you so much!! Best wishes! yamei ************************************************ Yamei Yu State Key Laboratory of Biotherapy／Collaborative Innovation Center of Biotherapy, West China Hospital, Sichuan University,Chengdu,610041, P.R.China Tel: 15882013485 Email: ymyuxtal(a)gmail.com ymyuxtal(a)163.com yamei_yu(a)scu.edu.cn

3 2

error messages running phenix.mr_rosetta
by Ursula Schulze-Gahmen 26 Nov '14

26 Nov '14

I am running mr_rosetta on a 64bit Linux computer. The model optimization worked without a problem, but now in the molecular replacement step I am receiving error messages, and I am not sure what the problem is. Could anyone provide a suggestion? Thanks Ursula /usr/share/themes/Slider/gtk-2.0/gtkrc:80: error: unexpected identifier `colorize_scrollbar', expected character `}' Fontconfig warning: no <cachedir> elements found. Check configuration. Fontconfig warning: adding <cachedir>/var/cache/fontconfig</cachedir> Fontconfig warning: adding <cachedir>~/.fontconfig</cachedir> (gedit:20211): GLib-GObject-WARNING **: Two different plugins tried to register 'BasicEngineFc'. (gedit:20211): GLib-GObject-CRITICAL **: g_object_new: assertion `G_TYPE_IS_OBJECT (object_type)' failed (gedit:20211): Pango-WARNING **: Failed to load Pango module '/programs/x86_64-linux/phenix/dev-1819/phenix-dev-1819/build/intel-linux-2.6-x86_64/base/lib/pango/1.6.0/modules/pango-basic-fc.so' for id 'BasicScriptEngineFc' (gedit:20211): GLib-GObject-WARNING **: Two different plugins tried to register 'BasicEngineFc'. (gedit:20211): GLib-GObject-CRITICAL **: g_object_new: assertion `G_TYPE_IS_OBJECT (object_type)' failed (gedit:20211): GLib-GObject-WARNING **: Two different plugins tried to register 'BasicEngineFc'. (gedit:20211): GLib-GObject-CRITICAL **: g_object_new: assertion `G_TYPE_IS_OBJECT (object_type)' failed (gedit:20211): GLib-GObject-WARNING **: Two different plugins tried to register 'BasicEngineFc'. (gedit:20211): GLib-GObject-CRITICAL **: g_object_new: assertion `G_TYPE_IS_OBJECT (object_type)' failed (gedit:20211): GLib-GObject-WARNING **: Two different plugins tried to register 'BasicEngineFc'. (gedit:20211): GLib-GObject-CRITICAL **: g_object_new: assertion `G_TYPE_IS_OBJECT (object_type)' failed (gedit:20211): GLib-GObject-WARNING **: Two different plugins tried to register 'BasicEngineFc'. (gedit:20211): GLib-GObject-CRITICAL **: g_object_new: assertion `G_TYPE_IS_OBJECT (object_type)' failed (gedit:20211): GLib-GObject-WARNING **: Two different plugins tried to register 'BasicEngineFc'. (gedit:20211): GLib-GObject-CRITICAL **: g_object_new: assertion `G_TYPE_IS_OBJECT (object_type)' failed (gedit:20211): GLib-GObject-WARNING **: Two different plugins tried to register 'BasicEngineFc'. (gedit:20211): GLib-GObject-CRITICAL **: g_object_new: assertion `G_TYPE_IS_OBJECT (object_type)' failed (gedit:20211): GLib-GObject-WARNING **: Two different plugins tried to register 'BasicEngineFc'. (gedit:20211): GLib-GObject-CRITICAL **: g_object_new: assertion `G_TYPE_IS_OBJECT (object_type)' failed (gedit:20211): Gtk-WARNING **: Could not find the icon 'system-file-manager'. The 'hicolor' theme was not found either, perhaps you need to install it. You can get a copy from: http://icon-theme.freedesktop.org/releases -- Ursula Schulze-Gahmen, Ph.D. Project Scientist UC Berkeley, QB3 360 Stanley Hall #3220 Berkeley, CA 94720-3220 (510) 643 9491

1 0

residue selection for rsrl?
by Kevin Jude 26 Nov '14

26 Nov '14

I am refining a structure at 2.9 Å resolution and have two side chains that keep moving out of density during reciprocal space refinement; both side chains should fit the resulting Fo-Fc density in favored rotamers. Turning on real-space refinement fixes this problem but leads to higher Rfree and poorer geometry overall. Can I select a subset of residues for real space refinement while still performing reciprocal space refinement on the whole asymmetric unit? Thanks for any help Kevin

2 2

phenix rosetta refine
by Andreas Förster 26 Nov '14

26 Nov '14

Dear all, I'm trying to get phenix.rosetta_refine to run. I think it's all set up properly, with Rosetta 3.5 and Phenix 1810 on RHEL 6.3 64bit. However phenix.rosetta_refine model.pdb data.mtz nproc=48 protocol=full > log& fails with the following error in log: End of input processing Sorry: The RosettaScripts executable could not be located. Please add the appropriate directory to your PATH environment variable, or set the variable ROSETTA_BIN to point to it. To me it seems things are defined correctly: andreas@csbstereo in phenix: echo $ROSETTA_BIN /csb/soft/Linux64/share/rosetta3.5/rosetta_source/bin andreas@csbstereo in phenix: ls $ROSETTA_BIN | grep scripts rosetta_scripts.linuxgccrelease@ rosetta_scripts.default.linuxgccrelease@ andreas@csbstereo in phenix: ll $ROSETTA_BIN/rosetta_scripts.linuxgccrelease lrwxrwxrwx 1 csbadmin csb 93 Oct 9 2013 /csb/soft/Linux64/share/rosetta3.5/rosetta_source/bin/rosetta_scripts.linuxgccrelease -> ../build/src/release/linux/2.6/64/x86/gcc/4.4/default/rosetta_scripts.default.linuxgccrelease Any idea of what's wrong would be greatly appreciated. Andreas -- Andreas Förster X-ray Crystallography Facility Manager Centre for Structural Biology Imperial College London

3 2

Alternative conformation of ligand with near zero occupancy
by Singh, Bishal 26 Nov '14

26 Nov '14

Hello everybody, I thank you for clarifying my previous doubts. Your suggestions matter a lot to me. Here I have another doubt. I have structure of a protein in complex with triglycan (fucose-galactose-galactose). The terminal galactose of this triglycan possesses two conformations with occupancies 0.4 and 0.6. In order to satisfy the geometric parameters of this triglycan, I derived the alternative conformation of fucose and the central galactose as well. Then I performed occupancy refinement. The occupancy of alternative central galactose is 0.06 and that of alternative fucose is 0.00. Now the question is, should I delete this fucose of zero occupancy or keep it to maintain the geometric parameters of triglycan? I am looking forward for your opinion. If you like I can provide a snapshot of ligand with electron density. With regards, Bishal Singh PhD student DKFZ, University of Heidelberg, Germany

2 1

input for specifying labin in phenix.find_ncs
by Arka Chakraborty 26 Nov '14

26 Nov '14

Dear Phenix BBers, I have a seemingly trivial problem in specifying column labels in phenix.find_ncs. Following the manual, on entering labels as : phenix.find_ncs xx.pdb yy.mtz params.ncs.labin = "labin FP=F PHIB=PHI FOM=FOM" , phenix returns the error : 'Sorry, cannot interpret the keyword: params.ncs.labin=labin FP=F PHIB=PHI FOM=FOM'. Any clues on what mistake I might be making will be greatly appreciated. Thanks a lot, Regards, Arko -- *Arka Chakraborty* *ibmb (Institut de Biologia Molecular de Barcelona)* *BARCELONA, SPAIN*

3 4

Different space group of coordinate and structure factors
by Singh, Bishal 25 Nov '14

25 Nov '14

Hello everybody, For a data-set, the space group after data-processing was P4222. When I performed molecular replacement through PHASER in phenix keeping the search option 'All possible in same point group", the space group changed to P41212. Now the question is that shall I process the data again in this new space group or should I move on for refinement? Please suggest me your opinion. With regards, Bishal Singh PhD Student DKFZ, University of Heidelberg, Germany

2 2

explicit solvent in low resolution refinement
by Guenter Fritz 25 Nov '14

25 Nov '14

Dear all, we have determined the structure of a protein complex at 3.5 A resolution as well as the structures of three subunits of this complex independently at resolutions better than 2 A. I had used the models of the subunits as reference models in refinement, perfect. Now I was wondering whether one could use explicit solvent positions (e.g. inner shell water molecules) from the high resolution structures in the low resolution refinement. This should be superior to a simple bulk solvent model. I placed the inner shell water molecules from the high res structures into the low res model and added harmonic restraints to the water positions in low res refinement. I see a decrease of about 0.5% in Rwork and Rfree versus the model without water molecules. However, adding harmonic restraints excludes the option to use the reference models in refinement. Is it possible to use protein and water atoms from the reference models to generate restraints for the low resolution refinement? Thank you, Guenter

5 9

ready set and H/D addition on waters
by Maxime Cuypers 24 Nov '14

24 Nov '14

hello, after not using phenix for a while, i just realised the new versions do not include the option in ready_set gui to add H or D atoms to waters in the pdb models (and only to the waters without altering the rest of the model would be nice). of course it can be done on the command line but why has it been removed? it is a very neutron crystallographer friendly option don t you think? cheerios, -- Dr. Maxime Cuypers

3 4

Rwork > Rfree
by Singh, Bishal 21 Nov '14

21 Nov '14

Hello everybody, While depositing a structure in PDB, I observed that Rwork was more than Rfree in the highest resolution shell and in a few more shells as well. I shall be delighted if somebody can suggest me the reasons and also what should I do now. With regards, Bishal PhD student, DKFZ, University of Heidelberg, Germany

4 3

MRage against local "select" PDB database on commandline
by Kaushik Hatti 21 Nov '14

21 Nov '14

Hi, I am new to Phenix. Sorry if my question is naive. For a given mtz file, I want to pick templates from a local PDB database (its actually iTasser template library) and perform Molecular Replacement using MRage. I have already installed Blast, downloaded iTasser PDB library (http://zhanglab.ccmb.med.umich.edu/library/PDB.tar.bz2) and extracted sequences from both the downloaded PDB template library set and target protein for which mtz is available. How should I go about building pipeline to perform MR using templates from locally downloaded select PDB database? Are there any more parameters available for parameter name: services under scop: search other than "local" and "ncbi"? Thanks in advance, Kaushik -- People living deeply have no fear of death - Anais Nin Caution: I am still the dumbest person I have ever known :-)

2 6

Submissions for Comput. Cryst. Newsletter
by Nigel Moriarty 20 Nov '14

20 Nov '14

Folks Calling for articles and short communications of interest to structural biologists. The deadline for submission in the January 2015 issue is 17 December 2014 (with some leeway if you are still working on the text, especially as it can be polished over the break). The Computational Crystallography Newsletter (CCN) is an electronic newsletter for structural biologists, and is published online every 6 months. Feature articles, meeting announcements and reports can be submitted to me at any time for consideration. Submission of text by email or word-processing files using the CCN templates is requested. Past newsletters and a Microsoft Word template are available at www.phenix-online.org/newsletter. Cheers Nigel --- Nigel W. Moriarty Building 64R0246B, Physical Biosciences Division Lawrence Berkeley National Laboratory Berkeley, CA 94720-8235 Phone : 510-486-5709 Email : NWMoriarty(a)LBL.gov Fax : 510-486-5909 Web : CCI.LBL.gov

1 0

Re: [phenixbb] Rwork > Rfree
by Ed Pozharski 20 Nov '14

20 Nov '14

Generate a different test set and repeat refinement. Do it twenty times. If R turns out to be systematically higher than Rfree even in few resolution shells, you have an interesting case worth looking at. If not, it's just noise. Sent on a Sprint Samsung Galaxy S® III <div>-------- Original message --------</div><div>From: "Singh, Bishal" <b.singh(a)dkfz-heidelberg.de> </div><div>Date:11/20/2014 9:52 AM (GMT-05:00) </div><div>To: phenixbb(a)phenix-online.org </div><div>Subject: [phenixbb] Rwork > Rfree </div><div> </div>Hello everybody, While depositing a structure in PDB, I observed that Rwork was more than Rfree in the highest resolution shell and in a few more shells as well. I shall be delighted if somebody can suggest me the reasons and also what should I do now. With regards, Bishal PhD student, DKFZ, University of Heidelberg, Germany _______________________________________________ phenixbb mailing list phenixbb(a)phenix-online.org http://phenix-online.org/mailman/listinfo/phenixbb

1 0

GeoStd: was definition of glycosylation in phenix.refine
by Schubert, Carsten [JRDUS] 18 Nov '14

18 Nov '14

Hi Nigel, could you please comment on any plans to harmonize the use of these restraints for other refinement programs as well. Phenix originally appears to have implemented the restraint dictionaries from CCP4. Right now there seems to be some kind of parallel development going on with the dictionaries being in sync to some degree. Is there a plan to move to a unified restraints distribution at some point? Thanks Carsten From: phenixbb-bounces(a)phenix-online.org [mailto:[email protected]] On Behalf Of Nigel Moriarty Sent: Tuesday, November 18, 2014 12:38 PM Cc: PHENIX user mailing list Subject: Re: [phenixbb] definition of glycosylation in phenix.refine @Laurent Engin advice is sound but you are having trouble free feel to send me the inputs directly. NB. Any files sent to me will be held in strictest confidence. @Engin I think you will find what you are looking for in the GeoStd. As restraints are improved, they are moved from $PHENIX/chem_data/mon_lib to $PHENIX/chem_data/geostd (which is hosted on SourceForge for anybody's use). The directory structure is the same as the Monomer Library with a geostd_list.cif file corresponding to mon_lib_list.cif which specifies the links. Cheers Nigel --- Nigel W. Moriarty Building 64R0246B, Physical Biosciences Division Lawrence Berkeley National Laboratory Berkeley, CA 94720-8235 Phone : 510-486-5709 Email : NWMoriarty(a)LBL.gov<mailto:[email protected]> Fax : 510-486-5909 Web : CCI.LBL.gov<http://CCI.LBL.gov> On Tue, Nov 18, 2014 at 8:25 AM, Engin Özkan <eozkan(a)uchicago.edu<mailto:[email protected]>> wrote: Hi Laurent, If I am reading the error message right, you did not delete the hydrogens that get removed upon bond formation. That's what it is complaining about. I always thought Phenix probably removed these atoms (the hydroxyl and the hydrogen) automatically when you defined the link, but never tested if it really does it. I also usually remove all hydrogens right away when sugar residues are inserted. Finally, I will recommend that you use automatic linking in the phenix.refine GUI. It works great (I am on the latest nightly, dev-1839; you may need to upgrade). No need to define pesky links anymore. Thank you developers! At the risk of hijacking the thread, I could not find in my installation the cif files for N-linked sugars. Have they been moved elsewhere? Hope this helps, Engin On 11/18/14 10:02 AM, Laurent Maveyraud wrote: Hi, I am currently refining a structure of a glycosylated protein, but I fail to define the appropriate links... The glycosylatin is Asn638-NAG1061-(beta1-4)-NAG1062-(beta1-4)-BMA1063-(ALPHA1-3)-MAN1064 This is what I add in the eff file : apply_cif_link { data_link = NAG-ASN residue_selection_1 = chain D and resname NAG and resid 1061 residue_selection_2 = chain A and resname ASN and resid 638 } apply_cif_link { data_link = BETA1-4 residue_selection_1 = chain D and resname NAG and resid 1062 residue_selection_2 = chain D and resname NAG and resid 1061 } apply_cif_link { data_link = BETA1-4 residue_selection_1 = chain D and resname BMA and resid 1063 residue_selection_2 = chain D and resname NAG and resid 1062 } apply_cif_link { data_link = ALPHA1-3 residue_selection_1 = chain D and resname MAN and resid 1064 residue_selection_2 = chain D and resname BMA and resid 1063 } this is correctly read in refine: apply_cif_link: data_link: NAG-ASN mod_id_1: DEL-O1 mod_id_2: DEL-HD22 residue_selection_1: chain D and resname NAG and resid 1061 residue_selection_2: chain A and resname ASN and resid 638 apply_cif_link: data_link: BETA1-4 mod_id_1: DEL-HO4 mod_id_2: DEL-O1 residue_selection_1: chain D and resname NAG and resid 1062 residue_selection_2: chain D and resname NAG and resid 1061 apply_cif_link: data_link: BETA1-4 mod_id_1: DEL-HO4 mod_id_2: DEL-O1 residue_selection_1: chain D and resname BMA and resid 1063 residue_selection_2: chain D and resname NAG and resid 1062 apply_cif_link: data_link: ALPHA1-3 mod_id_1: DEL-HO3 mod_id_2: DEL-O1 residue_selection_1: chain D and resname MAN and resid 1064 residue_selection_2: chain D and resname BMA and resid 1063 but generates an error : Number of atoms with unknown nonbonded energy type symbols: 5 "HETATM10951 HO4 NAG D1052 .*. H " "HETATM10973 HO4 BMA D1053 .*. H " "HETATM11028 HO4 NAG D1062 .*. H " "HETATM11050 HO4 BMA D1063 .*. H " "HETATM11070 HO3 MAN D1064 .*. H " and the job stops. Any ideas on what if wrong ? thanks for your help Laurent _______________________________________________ phenixbb mailing list phenixbb(a)phenix-online.org<mailto:[email protected]> http://phenix-online.org/mailman/listinfo/phenixbb

1 0

definition of glycosylation in phenix.refine
by Laurent Maveyraud 18 Nov '14

18 Nov '14

Hi, I am currently refining a structure of a glycosylated protein, but I fail to define the appropriate links... The glycosylatin is Asn638-NAG1061-(beta1-4)-NAG1062-(beta1-4)-BMA1063-(ALPHA1-3)-MAN1064 This is what I add in the eff file : apply_cif_link { data_link = NAG-ASN residue_selection_1 = chain D and resname NAG and resid 1061 residue_selection_2 = chain A and resname ASN and resid 638 } apply_cif_link { data_link = BETA1-4 residue_selection_1 = chain D and resname NAG and resid 1062 residue_selection_2 = chain D and resname NAG and resid 1061 } apply_cif_link { data_link = BETA1-4 residue_selection_1 = chain D and resname BMA and resid 1063 residue_selection_2 = chain D and resname NAG and resid 1062 } apply_cif_link { data_link = ALPHA1-3 residue_selection_1 = chain D and resname MAN and resid 1064 residue_selection_2 = chain D and resname BMA and resid 1063 } this is correctly read in refine: apply_cif_link: data_link: NAG-ASN mod_id_1: DEL-O1 mod_id_2: DEL-HD22 residue_selection_1: chain D and resname NAG and resid 1061 residue_selection_2: chain A and resname ASN and resid 638 apply_cif_link: data_link: BETA1-4 mod_id_1: DEL-HO4 mod_id_2: DEL-O1 residue_selection_1: chain D and resname NAG and resid 1062 residue_selection_2: chain D and resname NAG and resid 1061 apply_cif_link: data_link: BETA1-4 mod_id_1: DEL-HO4 mod_id_2: DEL-O1 residue_selection_1: chain D and resname BMA and resid 1063 residue_selection_2: chain D and resname NAG and resid 1062 apply_cif_link: data_link: ALPHA1-3 mod_id_1: DEL-HO3 mod_id_2: DEL-O1 residue_selection_1: chain D and resname MAN and resid 1064 residue_selection_2: chain D and resname BMA and resid 1063 but generates an error : Number of atoms with unknown nonbonded energy type symbols: 5 "HETATM10951 HO4 NAG D1052 .*. H " "HETATM10973 HO4 BMA D1053 .*. H " "HETATM11028 HO4 NAG D1062 .*. H " "HETATM11050 HO4 BMA D1063 .*. H " "HETATM11070 HO3 MAN D1064 .*. H " and the job stops. Any ideas on what if wrong ? thanks for your help Laurent -- ---------------------------------------------------------- Laurent Maveyraud laurent.maveyraud AT ipbs DOT fr P I C T --- Plateforme Intégrée de Criblage de Toulouse Université Paul Sabatier / CNRS / I.P.B.S. UMR 5089 Département Biologie Structurale et Biophysique http://cribligand.ipbs.fr http://www.ipbs.fr 205 route de Narbonne 31077 TOULOUSE Cedex FRANCE Tél: +33 (0)561 175 435 Fax : +33 (0)561 175 994 ----------------------------------------------------------

3 3

mrage with nucleic acids?
by Paul Paukstelis 17 Nov '14

17 Nov '14

I'm trying to run mrage with nucleic acid models. With phaser, it is possible to specify the molecule type, but this option doesn't exist in mrage (from what I can see). Providing a sequence file automatically assumes protein one letter codes. Any suggestion on how to clearly define a sequence file as a nucleic acid? Thanks, --paul

2 1

Filtering CCP4 maps to lower resolution
by Barad, Ben 16 Nov '14

16 Nov '14

Hi All, I am currently trying to filter a CCP4-formatted real space map from the EMDB, which is at 3.2 Å resolution, to different resolutions from 3.2 to 5 Å, and I've been having a bit of a difficult time getting it to work. My first thought was to do an fft to get structure factors, then regenerate the map from those structure factors with a resolution cutoff. To do this, I used the tools "phenix.map_to_structure_factors" and "phenix.mtz2map" to generate a structure factor file, then convert it back into a map. When I do this, I end up with a very noisy map that is not usable for what I would like to do. Is there a better way to accomplish this in phenix (or another tool)? I've also tried some low-pass filtering with proc3d, without much success... Thanks, Ben Barad Fraser Lab @ UCSF

3 3

Re: [phenixbb] simulated annealing omit map
by Pavel Afonine 15 Nov '14

15 Nov '14

Hello Bishal, as discussed before, computing ligand-omit map is tricky. Quoting my previous email: """ There are two commonly used options for omitting the ligand in order to demonstrate its presence/absence in Fo-Fc OMIT map: 1) Physically remove the ligand from PDB file. Then do some refinement and calculate Fo-Fc map. 2) Keep ligand in the file, set its occupancy to zero. Then, again, do some refinement and calculate Fo-Fc map. In this case you may want to ask refinement program to not move the ligand or move it only a little. Now, here is why these two options are poor and will not give you what you want. In the first case the bulk-solvent mask will be set in the ligand region and therefore it will mask ligand density (bulk-solvent will be filled into the ligand region). Depending on the strength of ligand density it may be masked completely or deteriorated. If you follow the second option you will always get positive density in ligand area. This density may correspond to bulk-solvent, ligand or mixture of both. That is there will be no simple way to differentiate whether this density arises from the ligand or bulk-solvent. """ Dale Tronrud offered a great alternative option: """ An alternative you might want to consider is what I call the "discovery map". At some point in the refinement process there was a map that convinced YOU that this ligand was present. You should be the hardest person to be convinced so that map will be both an omit map (because the model had been refined without the ligand prior to this) and clear enough to satisfy the reader. """ Now, recently I improved composite OMIT map calculation in phenix.composite_omit_map tool, see picture that illustrates how it works: https://www.dropbox.com/s/dbp8e0348v5p78h/fig_6.png?dl=0 I think at present that's the best option to follow if for some reason you cannot follow Dale's suggestion. Pavel On 11/13/14 6:58 AM, Singh, Bishal wrote: > Hello everybody, > I am generating figures showing the electron density map around the > ligands. I deleted the ligands from final co-ordinate file and then > performed refinement with identical parameters as before while > keeping simulated annealing=true. I shall be thankful if someone could > suggest me whether I consider mFo-DFc map or 2mFo-DFc map or both. > Kindly also tell me the minimum acceptable contour level required for > showing difference map. All structures are at 1.4 - 2.0 angstrom > resolutions. > > Regards, > Bishal

5 5

CB outliers
by George Devaniranjan 13 Nov '14

13 Nov '14

Hi, I refined a structure using phenix.refine. The R/R free decreased slightly from 13.46/15.87 to 13.0/15.6. But, molprobity said there were 3 CB outliers (starting from 0). I looked at those residues in COOT. I could NOT really see a different between the refined and initial model. In addition I looked at COOT validate-->density fit analysis I cannot see ANY difference in the peak heights (same numerical value too). I looked at RSCC (using phenix validate tools as well). I still can't see a difference in the plots. I am not sure how to get rid of the 3 CB outliers that have crept in since I can't pick them out as worse off in the refined model. I ask since in another case > 3 CB outliers, with the same results in picking them out in COOT/RSCC. I am wondering if there is another way I can see the deviation and fix it? Thank you, George

3 6

Rosetta compile problem
by miyatake 13 Nov '14

13 Nov '14

Dear all, I tried to compile Rosetta for rosetta_refinement, on Mac OX10.10. I run " /rosetta_2014.22.56873_bundle/main/source] user% sudo scons.py But it was in failure, showing the following message: g++-4.2 -o build/external/debug/macos/14.0/64/x86/gcc/4.2/default/dbio/cppdb/atomic_counter.os -c -isystem external/boost_1_55_0/ -isystem external/include/ -isystem external/dbio/ -pipe -ffor-scope -Wall -Wextra -pedantic -Werror -Wno-long-long -Wno-strict-aliasing -m64 -march=nocona -mtune=generic -O0 -g -ggdb -ffloat-store -fPIC -DBOOST_ERROR_CODE_HEADER_ONLY -DBOOST_SYSTEM_NO_DEPRECATED -DCPPDB_EXPORTS -DCPPDB_LIBRARY_SUFFIX=\".dylib\" -DCPPDB_LIBRARY_PREFIX=\"lib\" -DCPPDB_DISABLE_SHARED_OBJECT_LOADING -DCPPDB_DISABLE_THREAD_SAFETY -DCPPDB_SOVERSION=\"0\" -DCPPDB_MAJOR=0 -DCPPDB_MINOR=3 -DCPPDB_PATCH=0 -DCPPDB_VERSION=\"0.3.0\" -DCPPDB_WITH_SQLITE3 -Ibuild/external/debug/macos/14.0/64/x86/gcc/4.2/default/dbio/sqlite3 -Iexternal/dbio/sqlite3 -Iexternal/boost_1_55_0 -Iexternal/dbio -I/usr/include -I/usr/local/include external/dbio/cppdb/atomic_counter.cpp sh: g++-4.2: command not found scons: *** [build/external/debug/macos/14.0/64/x86/gcc/4.2/default/dbio/cppdb/atomic_counter.os] Error 127 scons: building terminated because of errors. How can I resolve that? Best regards, H. Miyatake

3 2

simulated annealing omit map
by Singh, Bishal 13 Nov '14

13 Nov '14

Hello everybody, I am generating figures showing the electron density map around the ligands. I deleted the ligands from final co-ordinate file and then performed refinement with identical parameters as before while keeping simulated annealing=true. I shall be thankful if someone could suggest me whether I consider mFo-DFc map or 2mFo-DFc map or both. Kindly also tell me the minimum acceptable contour level required for showing difference map. All structures are at 1.4 - 2.0 angstrom resolutions. Regards, Bishal

1 0

Settings-Directories of MR-Rosetta wizard is altered when phenix-GUI is restarted?
by Martin Moche 13 Nov '14

13 Nov '14

Dear colleagues, When restarting phenix followed by opening the "MR-Rosetta(beta)" the Settings-Directories need to be reset each time? How can I set Rosetta binary directory name: Rosetta binary name: Rosetta script directory name: Rosetta database directory name: in Settings-Directories of "MR-Rosetta(beta)" wizard (or elsewhere) so that they do not change each time I restart the phenix-GUI? These settings are not changed: Base output folder: Rosetta path: Likely "Rosetta path" does not change between phenix-GUI restarts since PHENIX_ROSETTA_PATH is set in the shell as: echo $PHENIX_ROSETTA_PATH /software/apps/rosetta/2014.35/gcc01/rosetta_2014.35.57232_bundle The "Path to Rosetta" in "Phenix preferences - Wizards" is currently left blank since $PHENIX_ROSETTA_PATH is set in the shell. I have by accident started a few jobs with mixed (new and old) settings and phenix-GUI starts without error messages and all processors are running 100% CPU, however likely in some infinite loop since no work is performed beyond "place model"... Best regards, Martin Martin Moche, Ph.D. Head of Protein Crystallography Karolinska Institutet MBB/PSF Scheeles väg 2 171 77 Stockholm Sweden phone: +46-8-524 868 43 mobile: +46-73-322 93 27 fax: +46-8-524 868 68 email: Martin.Moche(a)ki.se

1 0