- phenixbb - phenix-online.org

refinement of ligands and waters
by Alexandra Marques 10 Jun '15

10 Jun '15

Hi, I still have little experience in the refinement of protein structures. So, I have some (basic) doubts that I would like you to help me please. 1) I have a metal coordinated to a cofactor in the active center of my protein structure. The initial solution provided by phaser/autobuilding already contained the metal+cofactor modelled in the density but with zero occupancy. Then I created all the necessary restraints and cif files with phenix_ready_set and did two cycles of model building/phenix refinement. The refined structures always contained zero occupancy for those atoms. My question is: should I manually change the occupancy for the metal+cofactor to 1? I thought this would be made automatically by phenix because I selected the option “occupancies” and “individual B-factors” in the refinement strategy… 2) When manually placing waters in the density with coot, what is the maximum value of B-factor that is acceptable for a water? I am asking that because I have some water molecules with weak density and consequently B-factors between 45 and 55 but that are involved in interactions with other molecules. Is it better to keep those waters or remove them? Best regards, Alexandra

2 1

does the lastest Phenix only support windows 7
by Xiao Lei 09 Jun '15

09 Jun '15

Hi All, I am a beginner of crystallography and I'd like to install Phenix. I tried to download phenix for my Windows 8.1 system (64bit) but in the download page it says "64-bit (Windows 7 only) ". I do not know if I really have to use Windows 7. Does anyone have experience of installation Phenix on Windows 8.1? Thanks ahead. Xiao

2 2

Berkeley Server Downtime
by Paul Adams 08 Jun '15

08 Jun '15

Dear Colleagues, the CCI group in Berkeley is moving to a new building this week. As a result several structural biology web services will be offline for a short while. These include: - The Phenix web site (http://phenix-online.org/) - The DIALS West web site (http://dials.lbl.gov/) - The CCI web site (http://cci.lbl.gov/) - The Structural Biology Technology Portal (http://technology.sbkb.org/) The servers will be shutdown at 9am Pacific time on Wednesday 10th. Services should resume by the morning of Friday 12th at the latest. During the downtime bug reports and emails requesting help should be queued and received by us once the servers are back online. We greatly appreciate your understanding, and hope that this doesn’t present too much of an inconvenience. Cheers, Paul -- Paul Adams Deputy Division Director, Physical Biosciences Division, Lawrence Berkeley Lab Division Deputy for Biosciences, Advanced Light Source, Lawrence Berkeley Lab Adjunct Professor, Department of Bioengineering, U.C. Berkeley Vice President for Technology, the Joint BioEnergy Institute Laboratory Research Manager, ENIGMA Science Focus Area Building 64, Room 248 Building 80, Room 247 Building 978, Room 4126 Tel: 1-510-486-4225, Fax: 1-510-486-5909 http://cci.lbl.gov/paul Lawrence Berkeley Laboratory 1 Cyclotron Road BLDG 64R0121 Berkeley, CA 94720, USA. Executive Assistant: Louise Benvenue [ LBenvenue(a)lbl.gov ][ 1-510-495-2506 ] --

1 0

no merge original index
by Tarek DawoD 07 Jun '15

07 Jun '15

obtaining unmerged data from HKL2000 when i use "no merge original index" it either does nothing or gives this error ((could not read file /,,,/,,,,//,,///,,,,/.sca header: format is not supported by HKL-2000, File could be corrupted or it may be a result of a "no merge" command)) How could this be fixed?? did somebody else managed to get the unmerged data from HKL2000 or it is problematic? Thank you for your explanation and providing solutionTKDawoD

2 1

problem with B-factor plot within the Structure Comparison tool
by Emilia C. Arturo (Emily) 06 Jun '15

06 Jun '15

Hello All. I am trying to zoom in on a section of the B-factors v. residue plot that is generated using the Structure Comparison tool of the phenix GUI. I zoom in by selecting a narrow residue range. When I do so, the plot does not show all parts of each trace; they are broken up and/or residues are missing that are certainly in the pdb files that I supplied, and that showed up when I displayed a wider residue range. I am using phenix version 1.9-1692 on a Mac (OS X Yosemite). I'd input 21 pdb files. I've done some typical refresher-type things like closed and reopened the project, and restarted phenix (with no rationale, really), but the display problem persists. Does anyone else recognize this problem? Regards, Emily.

2 1

Model bias in MR-SAD
by mohamed noor 05 Jun '15

05 Jun '15

Dear all In the case of MR-SAD, how bad can model bias be? Does it also depend on the data resolution and resolution at which there is still enough anomalous signal? Or can the anomalous signal help in overcoming the bias inherent in the MR model? Thanks.

2 1

R free assignment in thin shells
by mohamed noor 04 Jun '15

04 Jun '15

Dear all I found a not-so-old thread about this topic: http://phenix-online.org/pipermail/phenixbb/2012-January/018259.html. Maybe there has been significant changes in phenix.refine since then. However, I couldn't find a conclusion. In my case, my data extends to 1.4 A with the R/R free being around 23/26 %. There is also pseudo symmetry, so there are 16 NCS copies in a lower symmetry group than what Xtriage suggested. The map looks fine with density for some H atoms even being visible. On what basis should I decide whether to use (or not) thin shell assignment? And what is the cut-off for number of NCS copies? Thanks.

3 2

Problems running phenix.rosetta_refine
by Klose, Thomas 04 Jun '15

04 Jun '15

Dear all, I was trying to set up phenix.rosetta_refine to test it for the refinement of a 4.5Å structure. I installed the latest phenix release and compiled the latest rosetta version (2015.19) using rosetta.build_phenix_interface. However, when I try to run phenix.rosetta_refine I get the following error message: …. Generating 5 models on 4 processors... DEBUG dumping error message ERROR: Rosetta exited with status 1 stderr output: Sorry: Wrong map type requested: %%map_type%% Allowed format is: [p][m]Fo+[q][D]Fc[kick][filled] where [p] and [q] are any numbers (optional), [m] and [D] indicate if the requested map is sigmaa (optional), Fo and Fc are Fobs and Fcalc, [kick] is for Average Kick Map (optional), [filled] is for missing Fobs filled map. Examples: 2mFo-DFc, 3.2Fo-2.3Fc, mFobs-DFcalc_kick, 2mFobs-DFcalc_filled, Fc, 2mFobs-DFcalc_kick_fill, anom, anom_diff, anomalous_difference, Fo ERROR: PhenixInterface::calculateDensityMap() : error writing map ERROR:: Exit from: src/core/scoring/cryst/PhenixInterface.cc line: 1122 This happens with our actual data as well as test data that I downloaded (PDB: 1yip). Any help in solving this problem would be greatly appreciated. Thanks, Thomas

3 4

Help with using FFT map coefficients tool
by Emilia C. Arturo (Emily) 04 Jun '15

04 Jun '15

Hello. I'd like to generate a map appropriate for viewing in PyMOL, using a model and structure factors downloaded from the PDB. But I'm unable to do so, and perhaps you can help me figure out what I might be doing wrong: To generate the map, I did the following using the phenix GUI (version 1.9-1692): 1. I used Reflection tools > Import CIF structure factors to convert the .cif format structure factors for PDB id 1J8U to an mtz file. 2. Then, to configure the job to generate the CCP4 map, I loaded this mtz file, but when I did that I got the following error: "This reflections file does not appear to contain any suitable map coefficients" I went through the same steps again for two other PDB ids, this time generating the mtz file by either letting phenix do it from a PDB id, or by supplying the tool with both a .cif file with structure factors and a pdb model file. I get the same error. Has anyone any advice? Regards, Emily.

3 5

Re: [phenixbb] Changing an amino acid to an unusual amino acid
by Schubert, Carsten [JRDUS] 03 Jun '15

03 Jun '15

Hi Nat, bonds in COOT do not mean much, as they do not translate into bonded interactions in the refinement programs. Phenix for instance tries to figure out the bonding by itself through some distance based heuristics. However in your case it looks as if the CRO residue may not be defined as a polymer, i.e. part of a protein chain in the monomer library you are using. In Windows the monomer library is located in C:\WinCoot\share\coot\lib\data\monomers. Look at the difference in the headers for CYS.cif and CRO.cif or PTR.cif (phosphotyrosine). The aminoacids or modified aminoacids considered part of the peptide chain have all 'peptide' in the header, CRO does not. That is probably the reason why coot does not make any bonds. This would probably cause problems later on with Phenix as well. You could pull out the CRO.cif into your directory, modify it to make it a peptide and read it into coot and phenix during refinement. I would also run it through REEL to make sure the definitions are kosher. HTH Carsten From: Natalia Ketaren [mailto:[email protected]] Sent: Tuesday, June 02, 2015 5:42 PM To: Schubert, Carsten [JRDUS] Subject: Re: Changing an amino acid to an unusual amino acid Hi Carsten, Th residue is a part of the correct chain. It just doesn't make an bonds with the residues on either side of it. Is this common for unusual amino acids when adding to the main chain? Cheers, Nat _____________________________ Natalia E. Ketaren, PhD Postdoctoral Associate Rout Laboratory Box 213 The Rockefeller University 1230 York Avenue New York, NY 10065 Tel: (212) 327-8136 Fax: (212) 327-7193 Natalia.Ketaren(a)rockefeller.edu<mailto:[email protected]> website: http://lab.rockefeller.edu/rout/ NCDIR: http://www.ncdir.org/ On Jun 2, 2015, at 5:26 PM, Schubert, Carsten [JRDUS] <CSCHUBER(a)its.jnj.com<mailto:[email protected]>> wrote: Hi Natalia, if the missing/wrong chain IDs are really the problem you can change them in COOT under Calculate -> Change Chain IDs. The dialog is a bit tricky, since you probably have to use the correct residue range to change an individual residue into an existing chain. In your case I would write out the PDB file and use a text editor to change the chain ID and move the residue into its correct place in the file. HTH Carsten From: phenixbb-bounces(a)phenix-online.org<mailto:[email protected]> [mailto:[email protected]] On Behalf Of Natalia Ketaren Sent: Tuesday, June 02, 2015 5:17 PM To: phenixbb(a)phenix-online.org<mailto:[email protected]> Subject: [phenixbb] Changing an amino acid to an unusual amino acid Hello All, I'm working on a GFP structure and I'm trying to incorporate the CRO (chromophore) residue into the structure using coot. My process is to use Extension->Modelling->Replace residues, where it's suppose to replace a Tyr residue. I type the 3 letter code, CRO, the residue appears, I'm able to real space refine it, but it does not incorporate CRO into the desired chain. Would someone be able to help me figure out what I'm doing wrong? And possibly walk me through a solution? Cheers, Nat _____________________________ Natalia E. Ketaren, PhD Postdoctoral Associate Rout Laboratory Box 213 The Rockefeller University 1230 York Avenue New York, NY 10065 Tel: (212) 327-8136 Fax: (212) 327-7193 Natalia.Ketaren(a)rockefeller.edu<mailto:[email protected]> website: http://lab.rockefeller.edu/rout/ NCDIR: http://www.ncdir.org/

4 3

Negative density
by mohamed noor 03 Jun '15

03 Jun '15

Dear all After refinement, there is a single negative peak in my map where nothing is modeled. What does it mean? Thanks.

3 2

Changing an amino acid to an unusual amino acid
by Natalia Ketaren 02 Jun '15

02 Jun '15

Hello All, I’m working on a GFP structure and I’m trying to incorporate the CRO (chromophore) residue into the structure using coot. My process is to use Extension->Modelling—>Replace residues, where it’s suppose to replace a Tyr residue. I type the 3 letter code, CRO, the residue appears, I’m able to real space refine it, but it does not incorporate CRO into the desired chain. Would someone be able to help me figure out what I’m doing wrong? And possibly walk me through a solution? Cheers, Nat _____________________________ Natalia E. Ketaren, PhD Postdoctoral Associate Rout Laboratory Box 213 The Rockefeller University 1230 York Avenue New York, NY 10065 Tel: (212) 327-8136 Fax: (212) 327-7193 Natalia.Ketaren(a)rockefeller.edu<mailto:[email protected]> website: http://lab.rockefeller.edu/rout/ NCDIR: http://www.ncdir.org/

2 1

Merging reflections
by mohamed noor 02 Jun '15

02 Jun '15

Dear all I have a few native datasets, one anomalous dataset with signal up to 10 A and two other datasets (also collected at the phasing wavelength but no anomalous signal seems to be present) and a sulfur SAD dataset with signal to about 7 A. All the datasets have a resolution of about 2-3 A and processed with xia2/XDS. I couldn't do fluorescence scan as the detector couldn't detect at such a high energy (20 keV). A few questions: 1. Why is there at least some signal from one dataset but not others? The signal should come from a covalently bound ligand, molybdopterin to be precise. 2. Since I(+) should not be equal to I(-) when there is anomalous signal, should I just merge the native dataset together with those collected at phasing wavelength that have no anomalous signal? Or do I merge everything together? 3. Will the anomalous signal to 10 A be useful to do MR-SAD? I could solve the structure with MR but I wanted to be sure my ligand is there. Thanks.

3 4

Allignment of multiple structures
by jens j birktoft 01 Jun '15

01 Jun '15

I apologize if this question this question has been asked before but I still need help finding an answer to the following. I am looking for a program/web-server that will calculate the superposition of multiple structures (non-protein!) Thanks -- +++++++++++++++++++++++++++++++++++++++++++++++++++ Jens J. Birktoft Structural DNA Nanotechnology Department of Chemistry New York University e-mail: jens.knold(a)gmail.com; Phone: 212-749-5057 very slow-mail: 350 Central Park West, Suite 9F, New York, NY 10025 +++++++++++++++++++++++++++++++++++++++++++++++++++

2 1

a question on phenix.real_space_refine
by Smith Liu 01 Jun '15

01 Jun '15

Dear Pavel, By phenix.real_space_refine, I fit a PDB with a 4.0 map. Most of the fitting were nice, but for the Ala residues, it indicates the o of the C=O occupies the density position which should be occupied by the CH3 side chain of the Ala, which also means the CH3 side chain of Ala cannot occupy the density position which clearly should be occupied by the CH3 side chain of Ala. Is any method or restraint to have the CH3 of Ala to ccupy its density position? Smith

1 0

alternate side chains conjugated with water molecules
by Jorge Iulek 31 May '15

31 May '15

Dear all, I have a 1.70 A resolution dataset in which I observe, with certain confidence, that some (few) side chains are alternated with waters. Though managing to tie occupancies of each alternate side chain with a water molecule to give an overall occupancy of 1.00, phenix refine still pushes them (water from side chain) apart as if there were bumping restraints (but their occupancies sum up to no more than 1.00). The closest to this I could find is: http://www.phenix-online.org/pipermail/phenixbb/2009-May/013354.html But I could not find how to turn off specific bumping restraints. I would be thankful if someone could point me to the recipe to do this. Jorge

2 2

Low resolution refinement - rotamers and Ramachandran
by mohamed noor 29 May '15

29 May '15

Dear developers I am refining a low resolution structure (3.8 A). After I fix most of the outliers (red bars) in Coot, phenix.refine is causing my structure to be worse than the start with 13 % Ramachandran outlier and 13 % rotamer outlier. Secondly, when I used phase_and_build to build my model, my R/Rfree is around 30/35 with atoms from some residues missing. When I fill them in with Coot, fix the outliers above and re-refine, R/Rfree jumps to 36/42 and the rmsd angles are too high (2.9). I have already used weight optimization and NCS. Am I doing something wrong? Thanks.

3 5

reference model restraints by segid
by jp d 28 May '15

28 May '15

Hi,I am having trouble getting phenix to used segid's in reference model restraints when my reference and target pdbs both don't use chain ID's.phenix quits with this message: *** Adding Reference Model Restraints *** Sorry: chain ' ' in /home/path/removed/agnref2.pdb contains one or more errant TER cards. Please remove and try again. my definition file looks like this:refinement { main { reference_model_restraints = True } reference_model { file = two.pdb } } refinement.reference_model.reference_group { reference = segid A16S selection = segid A16S } (above about 50 more times) this was on versions 1.9-1692 , dev-2016 and dev-2056 any suggestions ?thanksjpd

4 4

riding hydrogen model
by Filip Qbukarski 27 May '15

27 May '15

Dear phenix community I am currently refining a 1.75A structure and wanted to usethe riding hydrogen model in pehnix.refine. From what I found in the manual andprevious discussions on the mailing list and if I understand well, using ridingHs is advised for every resolution and the way of doing this in phenix.refine is notRefmac-like but instead I should 1) put Hs using ReadySet, 2) select the‘riding’ option for Hs and 3) use the PDB containing the H coordinates fordeposition, even though they were not individually refined. In this case theircontribution to scattering is also taken into account (which makes sense),which is a different way of treating riding Hs compared to Refmac, if Iunderstand well. Would it be possible to comment on this and in case I got thingswrong, correct?In fact for my structure, after adding the Hs prior to refinement using ReadySet and using theriding option, I actually got Rfree dropping by 1% (from 19.5 to 18.5) and better clash score (~0.5 according to Molprobity). Thanks in advancebestFilip

2 1

Re: [phenixbb] Refinement
by Pavel Afonine 26 May '15

26 May '15

Hi Mohamed, > How do tell phenix.refine to not move my HEM ligand around? yes, you can exclude selected atoms from refinement of coordinates in Phenix refinement GUI: Refinement settings -> Modify selection for: XYZ coordinates -> Edit -> type atom selection in "Edit selected" field (under XYZ refinement) like "chain A and resseq 1" or "chain A and resname ATP" However, I would rather find why correct geometry is not preserved. > I have a cif file describing the ligand but in my case, I can see > clear density of it not following the ideal geometry. If you send me files (off mailing list of course!) I will see if I can explain it.. You can copy to Nigel too as chances are that if you choose to send files I will be talking to him anyway. > As an example, I have attached a screenshot from Coot. It's hard to infer what's going on from a static image.. > In fact, the C atom is covalently linked to the nearest Cys, so the > HEM side chain should be 'moved' into the density and the green blob > shouldn't be there. Is this link defined anywhere in restraints used in refinement? Pavel

5 7

Fw: on the convergence of phenix refine
by Smith Liu 26 May '15

26 May '15

Dear All, For the meaning of non-convergence here, I mean in the whole 30 macro-cycles, both the Rwork and Rfree curves were in zigzag format, which I never met before. I am looking forward to getting an explanation from you. Smith -------- Forwarding messages -------- From: "Smith Liu" <smith_liu123(a)163.com> Date: 2015-05-26 09:26:23 To: phenixbb(a)phenix-online.org Subject: [phenixbb] on the convergence of phenix refine Dear All, Recently I have done a phenix refine, I find even I have 30 macro-cycles, the refine did not converge well. For example, for the last 5 cycles refine, the Rwork can still fluctuate between o.21 and 0.22, although both the final Rwork and Rfree decreased in comparison with the start Rwork and Rfree. What is wrong for my map and refine strategy? Smith

2 1

on the convergence of phenix refine
by Smith Liu 26 May '15

26 May '15

Dear All, Recently I have done a phenix refine, I find even I have 30 macro-cycles, the refine did not converge well. For example, for the last 5 cycles refine, the Rwork can still fluctuate between o.21 and 0.22, although both the final Rwork and Rfree decreased in comparison with the start Rwork and Rfree. What is wrong for my map and refine strategy? Smith

1 0

on B-factor refine of the phenix.real_space_refine
by Smith Liu 23 May '15

23 May '15

Dear Pavel, For the input, the last line of PDB is as followng, ATOM 21953 HH22 ARG H 295 53.862 115.032 84.116 1.00100.56 H For the run=ADP of the phenix.real_space_refine, the last line of the first model of all states pdb is, ATOM 21953 HH22 ARG H 295 53.862 115.032 84.116 1.00194.44 H and the last line of the second model of all states pdb is ATOM 21953 HH22 ARG H 295 53.862 115.032 84.116 1.00164.47 H. Clearly in each PDB of each model of the all states pdb produced, the b-factor is different, and they are different with the b-factors of the input pdb. Thus will you please explain why there are 2 states pdbs which seems both are refined pdbs? In addition, is it normal that the single state refined pdb cannot be opended? Next, will you please explain based on what of the mrc file, the phenix.real_space_refine run=ADP function got the b-factors? Further, I have noticed that the phenix.real_space_refine run=ADP gave a resolution which is different from the mrc general resolution given by relion or some other software. Will you please explain base on what the phenix.real_space_refine run=adp gives the resolution? Best regards. Smith -------- Forwarding messages -------- From: "Pavel Afonine" <pafonine(a)lbl.gov> Date: 2015-05-23 00:48:00 To: "Smith Liu" <smith_liu123(a)163.com>,"phenixbb(a)phenix-online.org" <phenixbb(a)phenix-online.org> Subject: Re: Fw:Re: [phenixbb] on B-factor refine of the phenix.real_space_refine For the all states PDB, it contains the input pdb for the phenix.real_space_refine and the real spaced refined PDB Yes I think so.. Is it not? Pavel

2 1

Fw:Re: on B-factor refine of the phenix.real_space_refine
by Smith Liu 22 May '15

22 May '15

Dear Pavel, It looks work. But for the refined PDB, it seems it was in read-only format, which I cannot copy. For the all states PDB, it contains the input pdb for the phenix.real_space_refine and the real spaced refined PDB, am I right? Smith -------- Forwarding messages -------- From: "Pavel Afonine" <pafonine(a)lbl.gov> Date: 2015-05-21 14:44:30 To: "Smith Liu" <smith_liu123(a)163.com>,phenixbb(a)phenix-online.org Subject: Re: [phenixbb] on B-factor refine of the phenix.real_space_refine I have a EM map (.mrc), with fitting and based on a crystal structure I got the fitted PDB file. I convert the mrc map to the CCP4 map， no need to convert .mrc into .ccp4 (isn't it the same thing?). Anyways, phenix.real_space_refine can use both. and I run phenix.real_space_refine model.pdb map.mtz run=adp, with phenix-dev-2034, in order to get my own b-factors for my fitted PDB based on my EM map. However after the above phenix.real_space_refine, I checked the refined PDB, and I find phenix.real_space_refine neither change the coordinates, nor change the b-factors for each atom. I need to see the files to be able to comment. If you share the files (off list and not sending directly via email but using things like Dropbox instead) then I will be happy to have a look. Thus will you please introduce to me how can I get my own b-factors for my fitted PDB in order to deposit my map and fitte pdb to the EM database? phenix.real_space_refine model.pdb map.mtz run=adp should just do what you want. If it didn't - see above (I need to see the files). Pavel

2 1

metal-ligand distances, and Optimize X-ray/stereochemistry weights
by Emilia C. Arturo (Emily) 21 May '15

21 May '15

I'm refining a model at 2.9 A, where four chains are in the ASU, and each active site has an iron with a few ligands; the configuration of this active site is the topic here. Currently the iron-ligand distances are not ideal (according to CheckMyMetal and the ideal distances generated by phenix.metal_coordination), but they are within a range observed previously for structures of truncated versions of this enzyme; furthermore, the more the metal-ligand distances change towards ideal distances, the more positive and negative density I observe in the Fo-Fc, so it appears that this deviation from the 'ideal' is supported by the data. However, when I use the 'Optimize X-ray/stereochemistry weights' option during refinement, the metal-ligand distances change (the metal appears to do most of the moving), getting closer to ideal, and this generates more positive and negative density in the Fo-Fc map). However, clearly it's nice that optimizing the X-ray/stereochemistry weights reduces R-free (and keeps the R-f/R-work ratio fine as well). So what is the most appropriate approach here, in order to keep the metal-ligand distances that best support the data, but still reap the benefits of the optimized stereochemistry weights? Some details of the structure and refinement are these: The iron in each of the four chains in the ASU is coordinated with atoms from two histidines, a glutamate, and a water (so four ligands per Fe) and there remains un-modeled density in the iterative build comp OMIT map used for modeling (because resolution is 2.9A and adding more water (which is my best guess, given what truncated structures at higher resolution show coordinating the Fe) does nothing to the Rf). I am working with an .edits file generated by phenix.metal_coordination. After some experimentation, I've left the sigma values (for Fe-O and Fe-N) at default values, and added a restraint for Fe-O (water), using an ideal distance of 2.3 and a sigma value of 0.1. Emily.

3 6