- phenixbb - phenix-online.org

metal refinement: adps and restraints
by Alexandra Marques 02 Jul '15

02 Jul '15

Hi, I am refining a protein with a Mo in the active site and which is bound to a cofactor and to a water molecule. After each refinement step there is always a lot of positive and negative density around the Mo and the cofactor molecules, although their occupancy refined to 1. I read this message http://www.phenix-online.org/pipermail/phenixbb/2015-May/022031.html, but I still have doubts about how to choose the refinement parameters. I know I have to select the Mo atoms to be refined anisotropically. And the other atoms, do I have to specifically select them to be refined isotropically?. For instance: Isotropic atoms: not element Mo Anisotropic atoms: element Mo I am asking that because I have been using TLS for the protein atoms during refinement, so I am not sure if it is correct to select the other atoms as isotropic or if I leave it in blank… My other problem is with the metal restraints. I am not sure about the ideal distance between the Mo and O atom but based on homologous protein structures I choose the following: refinement.geometry_restraints.edits { bond { action = *add atom_selection_1 = name MO and chain A and resname 4MO and resseq 1 atom_selection_2 = name O and chain A and resname O and resseq 2 distance_ideal = 1.900000 sigma = 0.050 } } However after the refinement this ideal distance interval (1.95-1.85) is not respected; I have 4 active sites and in some of them the Mo-O distance refines to 0.83 A although there is a lot of density for the water molecule to be at 1.9 A. How to make the restrainst work? Am I missing something? Thank you very much, Alexandra

3 2

accept MR solution even if the number of clashes is above the threshold
by Almudena Ponce Salvatierra 02 Jul '15

02 Jul '15

Hi everyone, I am running Phaser (from Phenix) and while checking the .log file (it is still running) I realize that it found some solutions with a TFZ score over 7, but it won't take them I guess because the number of clashes is higher than allowed (>12, I guess they're not so many either). My models are poor and for so I would like to check these solutions that for now Phaser will reject. However in the input I don't want to increase the allowed number of clashes because this will then affect all the other solutions and this will take forever... So my question is: can I tell phaser somehow to write down/save the solutions for which the TFZ score will be higher than 7? Even though the clashes Thanks a lot in advance All the best, Almu -- Almudena Ponce-Salvatierra Macromolecular crystallography and Nucleic acid chemistry Max Planck Institute for Biophysical Chemistry Am Fassberg 11 37077 Göttingen Germany

3 2

C-beta RMSD
by Kaushik Hatti 02 Jul '15

02 Jul '15

Hi, Is there a tool which can calculate C-beta RMSD for equivalent amino acids of homologous structures, post C-alpha superposition? Sorry if its off topic, Thanks, Kaushik -- Stupidity is everyone’s birthright. However, only the learned exercise it! --Kaushik (28Oct2014)

3 2

Simulated annealing refinement for just one domain
by Appu kumar 30 Jun '15

30 Jun '15

Hello All, I was wondering if there is a simulated annealing refinement for the just one domain and keeping the other domain fix. I have a structure where the electron density is poorly defined which could be because of the some conformation of domain is different. Although the resolution of the data is 4A but channel has well defined density in the TM domain. Thanks in advance for your advice and suggestions. Appu

2 1

while doing tar xvzf of nightly build 1.10pre-2089
by Murpholino Peligro 29 Jun '15

29 Jun '15

I get an error >>>> gzip: stdin: unexpected end of file phenix-installer-1.10pre-2089-intel-linux-2.6-x86_64-centos5/modules/phenix_regression/wizards/command_line_rosetta_tests/test_mr_rosetta_ncs.tgz tar: Unexpected EOF in archive tar: Unexpected EOF in archive tar: Error is not recoverable: exiting now <<<< Should I try another version?

1 0

Re: [phenixbb] labelit access in Phenix?
by R.D. Oeffner 28 Jun '15

28 Jun '15

This should have been forward to the mailing list. Rob -------- Original Message -------- Subject: Re: [phenixbb] labelit access in Phenix? Date: 2015-06-26 21:46 From: "R.D. Oeffner" <rdo20(a)cam.ac.uk> To: "Keller, Jacob" <kellerj(a)janelia.hhmi.org> Copy: Aaron Brewster <asbrewster(a)lbl.gov>, Nigel W Moriarty <NWMoriarty(a)lbl.gov> Hi Jacob, Labelit is not present on Windows in the last official release of Phenix from May last year, version 1.9-1692. But it is present in recent pre-releases of Phenix as command line tools. When you install Phenix on Windows a startup group is created that includes a command line prompt which recognises Phenix commands. From there you can execute all labelit command line tools. I think Nigel is about to put the recent Windows pre-release online on the Phenix website. Feel free to get in touch with further questions. Regards, Rob On 2015-06-26 18:44, Aaron Brewster wrote: > Hi Jacob, replying off of the list for now. There is no labelit gui > support in the phenix gui. There may be a batch script you have to > run to get the phenix command line binaries in the path for the cmd > program in win7, but I don't know it. If you can run phenix.refine on > the command line but not labelit, that would be instructive. > > cc'ing Rob who has been working on the phenix windows build. > > -Aaron > > On Thu, Jun 25, 2015 at 7:58 PM, Keller, Jacob > <kellerj(a)janelia.hhmi.org> wrote: > >> Dear Phenix list, >> >> Does anyone know how to access labelit through Phenix to determine >> whether the highest possible symmetry has been found? I exclusively >> use the GUI, and it's on win7. >> >> In the labelit documentation/homepage it mentions the command: >> >> labelit.check_pdb_symmetry [pdb coordinate file] [data=mtz file] >> >> But I did not find any executables called labelit, nor can I figure >> out how to run Phenix from command line in Win7, short of installing >> Cygwin or something. >> >> JPK >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://phenix-online.org/mailman/listinfo/phenixbb [1] >> Unsubscribe: phenixbb-leave(a)phenix-online.org > > > > Links: > ------ > [1] http://phenix-online.org/mailman/listinfo/phenixbb

1 0

phenix-1.9-1692 pymol problem
by ChenWeiFei 27 Jun '15

27 Jun '15

Dear all, I just install phenix and I found phenix-1.9-1692/pymol/install and it says: INSTALLATION of the PyMOL Executable Build for Unix 1) Extract the appropriate archive. tar -zxvf pymol-0_99-bin-linux-x86-glibc23.tgz 2) Change directory into the root level of the archive. cd pymol 3) Run "./setup.sh" which will create "./pymol" with a record of the current working directory. ./setup.sh 4) Link or copy ./pymol into a directory that is in your path. ln -s `pwd`/pymol /usr/local/bin/pymol (or perhaps) cp pymol /usr/local/bin/pymol 5) Now try it! ./pymol I follow it and when I type pymol there comes an error: $ pymol /usr/local/bin/pymol: /usr/local/phenix-1.9-1692/pymol/pymol.exe: /lib/ld-linux.so.2: bad ELF interpreter: No such file or directory Is any one konw how to deal with it? Best, Weifei

1 0

phenix.composite_omit_map
by Dyda 26 Jun '15

26 Jun '15

Dear Phenix list, I have generated an mtz file with phenix.composite_omit_map. It has labels assigned as 2mFo-DFC for F and PH2mFo-DFC for P. Now, I would like to generate a map file. I was trying to convince phenix.maps to use this mtz file to do that, however it didn't seem to like the labels. I tried both the GUI interface of phenix.maps and also tried to force feed it through the command line with a maps.param file. Neither of these worked, so I'm obviously missing something. Btw., Coot can open the mtz file and display the map, and it does look reasonable. This is phenix-1.9-1692 running under RHEL6. Thank you for your help. Fred [32m******************************************************************************* Fred Dyda, Ph.D. Phone:301-402-4496 Laboratory of Molecular Biology Fax: 301-496-0201 DHHS/NIH/NIDDK e-mail:[email protected] Bldg. 5. Room 303 Bethesda, MD 20892-0560 URGENT message e-mail: 2022476710(a)mms.att.net Google maps coords: 39.000597, -77.102102 http://www2.niddk.nih.gov/NIDDKLabs/IntramuralFaculty/DydaFred *******************************************************************************[m

2 1

Re: [phenixbb] [ccp4bb] phenix installtion problem
by Hattne, Johan 26 Jun '15

26 Jun '15

> On Jun 26, 2015, at 09:05, ChenWeiFei <weifei_c(a)OUTLOOK.COM> wrote: > > Dear all, > Sorry to disturb you. > I have install Phenix to my computer Centos 7.0. > When I want to use it I need to > cd /usr/local/phenix-1.9-1692/ > and > . /usr/local/phenix-1.9-1692/phenix_env.sh > then type phenix in the command line every time. > I wish to use phenix in every terminate window and just type phenix. > Some told me use export PATH to do it but I don't know how to do it. > I am new to linux(centos) > Thank you so much! > Best, > Weifei You could edit ~/.bashrc and add a line . /usr/local/phenix-1.9-1692/phenix_env.sh This will not necessarily have an affect on newly spawned login shells. To address that, you could make sure that your ~/.bash_profile contains a line like test -r "${HOME}/.bashrc" && . "${HOME}/.bashrc" near the top. // Johan Research Specialist @ Gonen Lab ____________________________________________________ Janelia Research Campus * 19700 Helix Drive Ashburn, VA 20147 * +1 (571) 209-4000 extension 3376

1 0

Re: [phenixbb] [ccp4bb] phenix installtion problem
by Alexandre OURJOUMTSEV 26 Jun '15

26 Jun '15

Dear Weifei, I am neither in the Phenix nor in the CCP4 teams so I let me remark that it does not look reasonable, just a few minutes after Tom Terwilliger's mail : "I'll answer you on the Phenix BB as you've asked specific questions about Phenix tools." sending your mail (see below) asking explicitly phenix questions again through the CCP4bb. I do not think that you "disturb", you should not sorry for this, the BB are for this, but we need to follow the rules. I ask everybody to respect both ccp4 and phenix teams, who assures a lot of service and help for our community, and send your requests properly (for your own benefit), through the BB concerned. My another request (luckily not applicable to the Weifei mail) is that people always sign properly their mails and give their name when asking their questions : ccp4, phenix teams and other experts answer you openly and this looks quite weird seeing anonymous requests. Again, this is for your own benefit, for your better integration into the research community. With best regards, Sacha Urzhumtsev ________________________________ De : CCP4 bulletin board <CCP4BB(a)JISCMAIL.AC.UK> de la part de ChenWeiFei <weifei_c(a)OUTLOOK.COM> Envoyé : vendredi 26 juin 2015 15:05 À : CCP4BB(a)JISCMAIL.AC.UK Objet : [ccp4bb] phenix installtion problem Dear all, Sorry to disturb you. I have install Phenix to my computer Centos 7.0. When I want to use it I need to cd /usr/local/phenix-1.9-1692/ and . /usr/local/phenix-1.9-1692/phenix_env.sh then type phenix in the command line every time. I wish to use phenix in every terminate window and just type phenix. Some told me use export PATH to do it but I don't know how to do it. I am new to linux(centos) Thank you so much! Best, Weifei

1 0

anomalous phasing with HySS - phaser - Autobuild
by Almudena Ponce Salvatierra 26 Jun '15

26 Jun '15

Hi all, I have SeMet data for which HySs locates 127 Selenium atoms with a CC of 0.31, which I think is decent. Then I run Phaser to generate the first maps, and it gives a score of 25 and LLG of 190 or so... The next step would be running Autobuild, however, the first models I realize it is building have R work and R free between 50 and 55%. So I guess it is just some random stuff there. If the result from HySS and Phaser looks more or less good, why does AutoBuild show these numbers?? Any idea of what I could do to improve so? Any ideas are welcome. Best wishes, Almudena -- Almudena Ponce-Salvatierra Macromolecular crystallography and Nucleic acid chemistry Max Planck Institute for Biophysical Chemistry Am Fassberg 11 37077 Göttingen Germany

2 1

Optimal Linux configuration for Phenix.refine
by Colin Levy 24 Jun '15

24 Jun '15

Apologies for the slightly off topic question. I am configuring a Linux workstation and would like your thoughts on an optimal setup for running Phenix. A lot of my computing time is spent in refinement and I am keen to put something together that will see a considerable increase in speed from my creaking Mac Pro (2 2.4GHz Quad core Intel Xeon with 24Gb RAM). Many thanks, Colin Manchester Protein Structure Facility Dr. Colin W. Levy MIB G034 Tel. 0161 275 5090 Mob.07786 197 554 c.levy(a)manchester.ac.uk<mailto:[email protected]>

5 5

Fw:Re: Autobuild Density
by Smith Liu 19 Jun '15

19 Jun '15

Dear All, Is any method to make the density modified map (from AutoBuild) comparable to the 2mFo-DFc map (from phenix.refine)? Smith -------- Forwarding messages -------- From: "Pavel Afonine" <pafonine(a)lbl.gov> Date: 2015-06-19 00:16:45 To: phenixbb(a)phenix-online.org Subject: Re: [phenixbb] Autobuild Density Hi Appu, a guess: perhaps you are comparing density modified map (from AutoBuild) with 2mFo-DFc map (from phenix.refine). If that's the case the maps are going to look differently indeed, so it's not too unexpected. Pavel On 6/18/15 08:47, Appu kumar wrote: > Dear All, > I have query about the density map from autobuild. The density map > from the autobuild look exceptionally good at 3.8A resolution. But > when i refined the autobuild pdb against the original mtz file, the > density become worse. Is it i am doing something wrong or i have to > use some parameterization? Structure is solved with the MR. > Appu _______________________________________________ phenixbb mailing list phenixbb(a)phenix-online.org http://phenix-online.org/mailman/listinfo/phenixbb Unsubscribe: phenixbb-leave(a)phenix-online.org

2 1

Autobuild Density
by Appu kumar 18 Jun '15

18 Jun '15

Dear All, I have query about the density map from autobuild. The density map from the autobuild look exceptionally good at 3.8A resolution. But when i refined the autobuild pdb against the original mtz file, the density become worse. Is it i am doing something wrong or i have to use some parameterization? Structure is solved with the MR. Appu

2 1

Get R values without refining in phenix.refine
by Murpholino Peligro 17 Jun '15

17 Jun '15

How can I get R values without refining anything? Just to check if my newest model is better than the old one?

3 4

Re: [phenixbb] anomalous difference map calculation
by Alexandra Marques 17 Jun '15

17 Jun '15

Thank you very much for your replies. Here is what I did: 1)I rerun xdsconv again with Friedls law=false in order to get the mtz file with the F(+) and F(-) columns (new.mtz). 2) I used my pdb file from the last refinement together with the new.mtz file as input for “Calculate maps” in Phenix and selected the option for anomalous map. Now I obtained the ANOM PHANOM map. I hope this is correct… In the anomalous map the largest density is in the Mo and Fe atoms. There is no density in the S of SO3 or SO4 groups, so, this will not help me to decide whether there are SO3 groups in the structure (because for the SO4 I have no doubt). However, while searching the mailling list I foud this post http://www.phenix-online.org/pipermail/phenixbb/2009-September/013799.html which worried me. Should I have used from the beginning this “new” mtz file with F+/F- in order to obtain anomalous differences map at each refinement step? Do you think that I should go back, transfer the R-free flags to my new mtz file and repeat all the refinements? Best regards, Alexandra 2015-06-17 16:59 GMT+01:00 Pavel Afonine <pafonine(a)lbl.gov>: > Hi Alexandra, > > to get anomalous difference map you do not need to do anything special: > phenix.refine calculates this map by default as long as your input data are > Fobs(+) and Fobs(-) or Iobs(+) and Iobs(-). > > Refining occupancy of Mo may not be a bad idea given it is likely to be > partially occupied. > > Pavel > > > On 6/17/15 06:04, Alexandra Marques wrote: > > Hi, > > > I am in the last refinement steps of a MR model and I want to calculate > an anomalous difference map essentially to confirm the presence of a > sulfite molecule and to locate vanadium (present in soaking solution). I > read that it is necessary to have a mtz file with anomalous data (i.e. > F+,F- or I+,I-). However, my data was collected at “normal” wavelenght > (0.97) and it was processed with XDS considering Friedls law= true and my > mtz file contain the following columns: H K L FP SIGFP. So, can I still > calculate a anomalous difference map based on my data? > > > Since I also have a Mo atom in the active site can I try to refine its > occupancy by using the option “anomalous groups” in the refinement > strategy? > > > Thank you very much, > Alexandra > > > _______________________________________________ > phenixbb mailing [email protected]://phenix-online.org/mailman/listinfo/phenixbb > Unsubscribe: phenixbb-leave(a)phenix-online.org > > >

3 3

anomalous difference map calculation
by Alexandra Marques 17 Jun '15

17 Jun '15

Hi, I am in the last refinement steps of a MR model and I want to calculate an anomalous difference map essentially to confirm the presence of a sulfite molecule and to locate vanadium (present in soaking solution). I read that it is necessary to have a mtz file with anomalous data (i.e. F+,F- or I+,I-). However, my data was collected at “normal” wavelenght (0.97) and it was processed with XDS considering Friedls law= true and my mtz file contain the following columns: H K L FP SIGFP. So, can I still calculate a anomalous difference map based on my data? Since I also have a Mo atom in the active site can I try to refine its occupancy by using the option “anomalous groups” in the refinement strategy? Thank you very much, Alexandra

4 4

keep torsion restraint on ligand
by CPMAS Chen 17 Jun '15

17 Jun '15

Hi, All Phenix Users, I am fitting 2-bromoethanol to the density map. I want to keep its gauche conformation, so I edited the cif for this ligand to restrain the torsion angle, loop_ _chem_comp_tor.comp_id _chem_comp_tor.id _chem_comp_tor.atom_id_1 _chem_comp_tor.atom_id_2 _chem_comp_tor.atom_id_3 _chem_comp_tor.atom_id_4 _chem_comp_tor.value_angle _chem_comp_tor.value_angle_esd _chem_comp_tor.period BRJ Var_01 BR1 CB CB1 OB1 -60.00 0.0 3 But, after refinement, this angle can be close to 90 or even 120. Even worse is that the fitting at such angles are not better than at -60/60 angle. How could I fix/restrain the angle during refinement? Thanks! Charles *************************************************** Charles Chen Research Associate University of Pittsburgh School of Medicine Department of Anesthesiology ******************************************************

3 3

Re: [phenixbb] anomalous difference map calculation
by Ashok Nayak 17 Jun '15

17 Jun '15

Hello Alexandra, I believe you have to reprocess the data keeping friedel pairs (f+ and f- ) separate to make use of the little most anomalous signal (if any) at your collected wavelength. The best thing would be to collect a high redundant set of peak data (if not all the three data sets including mad remote and inflection data) If you have access to a tunable source and xtals left. Once you process differently, you can do all sorts of substructure search and later refine anomalous scatterers using wavelength, atom type and f ' and f“ values. Sulphur sad is also a good possibility if you have a decent nmbr of sulphurs in your protein. Hope this suffices your query. Regards Ashok Nayak CSIR CDRI, Lucknow, India Alexandra Marques <at.marques(a)fct.unl.pt> wrote: >_______________________________________________ >phenixbb mailing list >phenixbb(a)phenix-online.org >http://phenix-online.org/mailman/listinfo/phenixbb >Unsubscribe: phenixbb-leave(a)phenix-online.org

1 0

phenix refine of my EM data
by Smith Liu 17 Jun '15

17 Jun '15

Dear All, By phenix.real_space_refine, I have fitted the a PDB file to an EM map, with Ramachandran favored around 90% and acceptable level of Ramachandran outlier (1.5%) and rotamer outlier (1.2%). Besides visual checking indicating the fitting of the main chain and sid chains to the dnsity map is acceptable. I hope the refine results can be improved by phenix reciprocal space refine. By phenix.map_to_structure_factors I convert my mrc file to mtz file, with the above mentioned phenix.real_space_refine refined pdb amd mtz files as the input files for the graphical phenix refine, in which I never select real space refine, and I always check xyz coordinates and Group B-factors (sometimes I check xyz coordinates and individual B-factors ). Fot this kind of phenix reciprocal space refine, I find the molprobity anslysis results have improved a little, but I find both the contouring level of 2FOFCWT map and the FOFCWT map is more than thousands of thousands high (and it seems the contouring level of FOFCWT is much higher), and visual checking indicated both the main chain and side chains hardly fit with the density map. But I really find there were publications which did phenix.real_space_refine first, and then process the reciprocal space refine, which get good results as introduced in the publications. Will you please tell me why my phenix reciprocal refine after phenix.real_space_refine did not work? Best regards. Smith

1 0

on phenix.map_to_structure_factors
by Smith Liu 13 Jun '15

13 Jun '15

Dear All For phenix.map_to_structure_factors, it needs to input the resolution d_min. But for an EM map, usually the resolution is the general resolution of the whole map (or to some extent the average resolution of the whole map). Is d_min for phenix.map_to_structure_factors the general EM map resolution, or the best resolution of the whole EM map, which is only for part of the EM map with the best quality? Smith

2 1

on real_space_refine
by Smith Liu 12 Jun '15

12 Jun '15

Dear All, After fitting a PDB into a 3.5-5 A resolution EM map by phenix.real_space_refine, should I still need to do a reciprocal space refine? For crystallography, the reciprocal space is in fact not purely imaginary, it corresponds to the diffaction spot. But for the EM data refinement, what is the reciprocal space refine? And what reciporce spce refine works on for the data refinement? Smith

2 1

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

Your abstract can still be considered for an oral presentation at the American Crystallographic Assn Annual meeting
by Terwilliger, Thomas Charles 21 May '15

21 May '15

Hi Phenix users, A message from the American Crystallographic Association! YES! We are still accepting abstracts. No late fee. ACA Annual Meeting • Philadelphia, PA • July 25 - 29, 2015 Submit now to be considered for an oral presentation. Once all the lectures are scheduled, abstracts will continue to be accepted for poster presentations. see: http://www.AmerCrystalAssn.org All the best, Tom T Vice-President, ACA

1 1

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

21 May '15

Dear All, 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， 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. 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? Smith

2 1

Map quality
by mohamed noor 19 May '15

19 May '15

Dear all Is there a single (or a few) metrics that can be used to quantitatively assess map quality instead of looking at each one in Coot? For example, I want to compare the effect of having low completeness in the low resolution shells. Thanks.

6 8