- phenixbb - phenix-online.org

Refinement with new mtz.file
by Aleksandar Bijelic 14 Jul '15

14 Jul '15

Dear CCP4 user and experts, I refined (with PHENIX) a 3.0 A dataset obtaining Rfree of about 0.24 (with good geometry according to Ramachandran, Beta Outliers, etc.) .... Everything seems to be ok (especially in relation to the resolution) .... because the mtz.file I used was quite old and I cannot find my xscale.hkl file, I processed the data set again (this time with optimizng and polishing) and received a "better" file according to almost everything (resolution limit, I/sigma, CC(1/2), Rmeas), thus I used this new mtz.file and put it in my last refinement step (the refinement which led to the above mentioned Rfree = 0.24). Suprisingly, the refinement starts at Rwork = 0.18 and Rfree = 0.19 but ending up with 0.20 and 0.22, respectively. So I wanted to know if this is usual? I was expecting my data to become slightlly better but what is irritating me is the starting R-values of the refinement and that it get worse during refinement. Maybe I did something wrong? Is it reasonable to replace the mtz.file with a new one in the last refinement step or should I start the refinement from the scratch? In the sake of completeness, I deleted the header of the used pdb.file because of the R-flag error which occurs, since Phenix reconizes that the pdb.file was already used with other Rflags. Sorry, but I am still a beginner in this field, so I would be very grateful if somebody could explain me this situation and my mistake and if I need to start refinement from the beginning. Thank you in advance! Best Regards. Aleksandar -- ------------------------------------------- Aleksandar Bijelic, MSc. Institut für Biophysikalische Chemie Universität Wien Althanstrasse 14 A-1090 Wien Tel: +43 1 4277 52536 e-Mail: aleksandar.bijelic(a)univie.ac.at --------------------------------------------

7 10

Set up Simulated Annealing starting tempreture
by Xiao Lei 13 Jul '15

13 Jul '15

Dear Phenixbb members, I'm learning phenix this days. I use phenix refinement and found there is a simulated annealing option (SA Cartesian or SA torsion option). However, I found the SA cartesian starts at 5000K and SA torsion starts at 2500K by default, and I'd like to change the starting temperature by myself. I search phenix manual and found: To change starting SA temperature: % phenix.refine data.hkl model.pdb simulated_annealing=true \ strategy=individual_sites simulated_annealing.start_temperature=10000 I do not understand this, does this mean I have to type the above command in somewhere in phenix? is there a place I can simply put SA start temperature numbers in? Xiao

2 2

NCS and R-Values Hung up at 0.42/0.46 for 2.1Å resolution
by John Bruning 11 Jul '15

11 Jul '15

I collected two datasets of the same ligand bound structure at 2.4 and 2.1 Å resolution. The apo structure has already been solved to 2.2Å resolution in P31 (and P3121 for SeMet). After processing each dataset as both P31 and P3121 in MOSFLM, I proceeded with MR with Phaser using sequential partial search models to place either 3 copies in the ASU for P3121, or 6 copies for P31. (The previously solved apo structure had 6 copies of the same protein in the ASU). *** in summary two datasets, processed as both P31 or P3121 After running xtriage for either .mtz, I get this warning about NCS: Distance to origin : 56.728 Height relative to origin : 26.388 % Translational pseudo-symmetry is very likely present in these data. Be aware that this will change the intensity statistics and may impact subsequent analyses, and in practice may lead to higher R-factors in refinement. Regardless, I can proceed with the MR using sequential partial search models to place either 3 copies (P3121) or 6 copies (P31) with decent solutions (i.e. packing looks good upon inspection, scoring looks appropriate -> SOLU SET RFZ=2.8 TFZ=8.3 PAK=13*LLG*=-37 *LLG*=2199 as an example). However, I did get the following somewhat expected warnings: Number of known components (1) not divisible by 2: Translational *NCS* correction not applied Large non-origin Patterson peak indicates that translational *NCS* is present Upon any attempt at initial refinement (with a number of different strategies), the R-factors are hung up > 0.42 (for example in the P3121, 2.1Å case): Start R-work = 0.4966, R-free = 0.5133; Final R-work = 0.4231, R-free = 0.4477 My concern is that there is some pathology that is not being recognized (merohedral twinning) or that the NCS/pseudosymmetry is not being addressed appropriately. -- Cheers, John

1 0

Re: [phenixbb] Fwd: [ccp4bb] Phenix; converting a map to structure factors
by Pavel Afonine 11 Jul '15

11 Jul '15

Hi Mark, Tom copied me your email. > As I'm sure you are aware there is tremendous interest in building > into and refining against low resolution (by X-ray standards) EM maps. > Could you comment to the phenix BB about whether its appropriate to > convert maps to reciprocal space and use ML refinement? While it is a possible way of doing this (and this is how it's been done in the past) I don't think it is a good route to follow. In case of cryo-EM your data is map, so why not just refine against the map? Why take a detour through reciprocal-space by converting map into structure factors (=manipulating your original data) and doing reciprocal space using refinement programs that were designed for crystallography? With this in mind we have developed a new refinement tool in Phenix (it is still in active development!) called phenix.real_space_refine that is designed to refine atomic models directly into maps (any map: X-ray, neutron, cryo-EM, low or high resolution), for quick summary see http://phenix-online.org/presentations/real_space_refine.pdf as well as documentation page (at the bottom): http://phenix-online.org/documentation/ and newsletter article (page 43): http://phenix-online.org/newsletter/CCN_2013_07.pdf Let me know if you have any questions or comments! Pavel

2 4

atom vdW dictionary for Python
by David Shin 11 Jul '15

11 Jul '15

Hi - I'm writing a Python script to measure the volume of pockets within protein structures - plus extra bells and whistles... thus a custom job. I was wondering just for speed, if someone could point me to values for van der Waals radii that are somewhat agreed upon for protein atoms, or if there is some dictionary already made that someone could point to that resides in phenix or some other structural biology package that would be ok for me to copy/use. Thanks, Dave

4 3

Re: [phenixbb] [ccp4bb] Phenix; converting a map to structure factors
by Terwilliger, Thomas Charles 10 Jul '15

10 Jul '15

Hi Reza, Just use: phenix.map_to_structure_factors map.ccp4 d_min=2.5 or something like that, All the best, Tom T On Jul 10, 2015, at 11:27 AM, Reza Khayat wrote: Hi, Is it possible to convert a map to structure factors using Phenix? Thanks. Best wishes, Reza Reza Khayat, PhD Assistant Professor City College of New York Department of Chemistry New York, NY 10031

3 4

NCS and R-Values Hung up at ~0.42/0.46
by John Bruning 10 Jul '15

10 Jul '15

Hi Everyone, I collected two datasets of the same ligand bound structure at 2.4 and 2.1 Å resolution. The apo structure has already been solved to 2.2Å resolution in P31 (and P3121 for SeMet). After processing each dataset as both P31 and P3121 in MOSFLM, I proceeded with MR with Phaser using sequential partial search models to place either 3 copies in the ASU for P3121, or 6 copies for P31. (The previously solved apo structure had 6 copies of the same protein in the ASU). *** in summary two datasets, processed as both P31 or P3121 After running xtraige for either .mtz, I get this warning about NCS: Largest Patterson peak with length larger than 15 Angstrom: Frac. coord. : 0.333 -0.333 -0.153 Distance to origin : 56.728 Height relative to origin : 26.388 % p_value(height) : 2.659e-03symmetry. Translational pseudo-symmetry is very likely present in these data. Be aware that this will change the intensity statistics and may impact subsequent analyses, and in practice may lead to higher R-factors in refinement. Regardless, I can proceed with the MR using sequential partial search models to place either 3 copies (P3121) or 6 copies (P31) with decent solutions (i.e. packing looks good upon inspection, scoring looks appropriate -> SOLU SET RFZ=2.8 TFZ=8.3 PAK=13*LLG*=-37 *LLG*=2199 as an example). However, I did get the following somewhat expected warnings: Number of known components (1) not divisible by 2: Translational *NCS* correction not applied Large non-origin Patterson peak indicates that translational *NCS* is present Upon any attempt at initial refinement (with a number of different strategies), the R-factors are hung up > 0.42 (for example in the P3121, 2.1Å case): Start R-work = 0.4966, R-free = 0.5133 Final R-work = 0.4231, R-free = 0.4477 *** best R-factors from ~ 15 different strategies My concern is that there is some pathology that is not being recognized (merohedral twinning) or that the NCS/pseudosymmetry is not being addressed appropriately. Any suggestions in troubleshooting this would be much appreciated, as I have pretty much exhausted everything I could think of. -- Cheers, John

2 1

(no subject)
by Ashok Nayak 10 Jul '15

10 Jul '15

Dear Phenixites, I have a helical bundle structure, where in I wish to understand by what rotation and about which axis the dimers are inter-related. What I get from the superposition with Superpose is the transformation matrix and not the axis, and the angles are also elusive. Ashok Nayak CSIR-CDRI Lucknow, India --

2 1

Axis of Rotation
by Ashok Nayak 10 Jul '15

10 Jul '15

Dear Phenixites, I have a helical bundle structure, where in I wish to understand by what rotation and about which axis the dimers are inter-related. What I get from the superposition with Superpose is the transformation matrix and not the axis, and the angles are also elusive. Thanks in advance -- Ashok Senior Research Fellow - Dr JV Pratap, Lab No-LSN 008 Molecular and Structural Biology Division Central Drug Research Institute, Janakipuram Extension Lucknow-226031 India

2 1

rosetta.build_phenix_interface problem
by Appu kumar 09 Jul '15

09 Jul '15

Dear All, I have installed Rosetta on Linux machine. When i am issuing the command " rosetta.build_phenix_interface" it is giving me following error scons: Building targets ... scons: *** [build/src/release/linux/2.6/64/x86/gcc/4.4/python/libObjexxFCL.so] Implicit dependency `/usr/local/lib/libdl.so' not found, needed by target `build/src/release/linux/2.6/64/x86/gcc/4.4/python/libObjexxFCL.so'. scons: building terminated because of errors. Would you please give me your valuable suggestions to successfully carry out the Rosetta refinement. Thank you Appu

1 0

Rosetta Refinement problem
by Appu kumar 09 Jul '15

09 Jul '15

Dear All, I have installed Rosetta successfully on Linux machine and after installation it say building done. While refining a structure with "phenix.rosetta_refine" it is giving following problem 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. In my .cshrc file i have setup the path as follow setenv PHENIX_ROSETTA_PATH /home/appu/2014_09_29_NSLS_X29/rosetta_bin_linux_2015.25.57927_bundle setenv ROSETTA_BIN /home/appu/2014_09_29_NSLS_X29/rosetta_bin_linux_2015.25.57927_bundle/main/source/bin setenv phenix /home/appu/software/phenix/phenix-1.9-1692 setenv ROSETTA3_DB /home/appu/2014_09_29_NSLS_X29/rosetta_bin_linux_2015.25.57927_bundle/main/database I would appreciate any help and suggestions in executing the Rosetta refinement Thank you Appu

1 0

Problem with python while running Autobuild
by Tam＠em.uni-frankfurt.de 09 Jul '15

09 Jul '15

Dear all, I am currently facing a problem of python while running the Autobuild in Phenix. phenix.autobuild data=data1.mtz model=model.pdb seq_file=protein.dat rebuild_in_place=False \ composite_omit_type=iterative_build_omit then it stopped with error Error in `/home/mauricetam/software/phenix-1.9-1692/build/intel-linux-2.6-x86_64/base/bin/python2.7': free(): invalid next size (normal): 0x0000000003e100a0 *** I wonder what's wrong with it. Regards Heng-Keat

2 1

problems with rosetta
by Benny Belviso 09 Jul '15

09 Jul '15

Dear all, I am trying to use rosetta with the Phenix Gui (Phenix is installed on Debian). I have downloaded the tar.gz file named "rosetta_bin_linux_2015.25.57927_bundle" that contains source and binary files. Becouse the tar.gz contains pre-compiled files, I have not compiled rosetta again (however, I have the same problems also in the case of tar.gz containing source after the compilation). According to the instructions showed on http://strucbio.biologie.uni-konstanz.de/ccp4wiki/index.php/Phenix#Installa… I have setup my environmental variables in .bashrc as: export PHENIX_ROSETTA_PATH=/home/bbdanilo/program/rosetta_bin_linux_2015.25.57927_bundle export ROSETTA3_DB=$PHENIX_ROSETTA_PATH/main/database However, by opening MR-Rosetta on Phenix GUI, I get a first warning message that says me that rosetta directory seems to be not correct. By running rosetta, the program says that it is not able to find PHENIX_ROSETTA_PATH. I have tried to change the path also in the "Directory setup" by typing /home/bbdanilo/program/rosetta_bin_linux_2015.25.57927_bundle as directory for rosetta path /main/source/bin as Rosetta binary directory name mr_protocols.linuxgccrelease as binary of rosetta In this case, the program says that does not find the binary for rosetta (it requires mr_rosetta.default, instead I should use mr_protocols.linuxgccrelease). What am I doing wrong? Any help is welcome. Thanks. Danilo This message and any attachment are intended solely for the addressee and may contain confidential information. If you have received this message in error, please send it back to me, and immediately delete it. Please do not use, copy or disclose the information contained in this message or in any attachment. Any views or opinions expressed by the author of this email do not necessarily reflect the views of the University of Nottingham. This message has been checked for viruses but the contents of an attachment may still contain software viruses which could damage your computer system, you are advised to perform your own checks. Email communications with the University of Nottingham may be monitored as permitted by UK legislation.

2 1

Re: [phenixbb] on determination whether a protein is a monomer or complex by phenix
by Valerie Pye 07 Jul '15

07 Jul '15

Hi Smith, Your solvent content is most likely the 66% one (if you have diffraction to about 1A there's a chance it could be 33%.....in which case try to solve it with Balbes (ccp4)). You can not tell what the crystal packing is like until you solve the structure! Even if you know how many molecules you have in your a.s.u. you can not tell if the contacts are formed because of crystal packing or oligo formation! Its much better to look at your protein in solution for this (as i said before!) Cheers, val On 7 July 2015 at 03:17, Smith Liu <smith_liu123(a)163.com> wrote: > Dear Valerie and all, > > My space group is I422 (attached image for I422). There are 2 solutions > for Mathews coefficient, one is solvent content 66.76%, the other is > 33.52%. Based on this message and the attached images, how can we know > whether my protein is a monomer (for example, monomer protein only exist in > the "+" position), a tetramer (or crystal packing as a "tetramer", each > chain exist in the "+" position, and each chain so close to each other so > that they contact each other and a tetramer forms), a 8-subunit protein > complex (or crystal packing as a " 8-subunit protein complex ", each 2 > chains exist in the "+" position, and each 2 chians so close to each other > so that they contact each other and a 8-subunit protein complex form)? > > I am looking forward to getting a reply from you. > > Best regards. > > Smith > > > > > > At 2015-07-07 00:34:32, "Valerie Pye" <valpye(a)googlemail.com> wrote: > > Hi Smith, > > I'm fairly sure you can not tell the oligomeric state of your protein from > an mtz and your sequence. > What you can do is calculate the Matthews co-efficient with your given > unit cell parameters and size of monomer to see how many copies of your > monomer you may have in the a.s.u. > The difficulty lies in the fact that if say you have 4 copies of your > protein in the a.s.u - this doesn't necessarily mean you have a tetramer - > they could be 4 monomers related by n.c.s - or 2 dimers.... > Is it possible to solve your structure? Then you confirm your space group > and number of molecules in the a.s.u. - then look at the area of contact > between monomers to see if they are likely to be because of multimerisation > or crystal packing. > > Totally agree with Tim that a better way (which you'd probably want to do > to back-up a multimeric crystal structure) is by SEC or native PAGE. You > could also try SAXS, SEC-MALS, EM, or DLS depending on the size of protein > you're looking at. You'd also want to know that you have a clean oligomer > and not an aggregate ;) > > Good luck, > > val > > On 6 July 2015 at 12:21, Smith Liu <smith_liu123(a)163.com> wrote: > >> Dear Tim, >> >> I mean based on the mtz file, can I know whether I will get a PDB file of >> a monomer or dimer, or trimer, or tetramer or oligomer? >> >> Best regards. >> >> Smith >> >> >> >> >> >> At 2015-07-06 18:14:03, "Tim Gruene" <tg(a)shelx.uni-ac.gwdg.de> wrote: >> >Dear Smith, >> > >> >I would say that size exclusion chromatography or native PAGE are better >> >experimental techniques than crystallography for this type of question. >> > >> >Best regards, >> >Tim >> > >> >On 07/06/2015 11:20 AM, Smith Liu wrote: >> >> Dear All, >> >> >> >> >> >> Suppose I have got a mtz file, and I know the aa sequence of the protein from which I got the crystal and the corresponding mtz file, will you please tell me which phenix function (program) can tell me whether my protein crystal was formed by protein monomer or by protein complex, and if formed by protein complex, how many subunits existing in the complex? >> >> >> >> >> >> Best regards. >> >> >> >> >> >> Smith >> >> >> >> >> >> >> >> _______________________________________________ >> >> phenixbb mailing list >> >> phenixbb(a)phenix-online.org >> >> http://phenix-online.org/mailman/listinfo/phenixbb >> >> Unsubscribe: phenixbb-leave(a)phenix-online.org >> >> >> > >> >-- >> >-- >> >Dr Tim Gruene >> >Institut fuer anorganische Chemie >> >Tammannstr. 4 >> >D-37077 Goettingen >> >phone: +49 (0)551 39 22149 >> > >> >GPG Key ID = A46BEE1A >> > >> > >> >> >> >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://phenix-online.org/mailman/listinfo/phenixbb >> Unsubscribe: phenixbb-leave(a)phenix-online.org >> > > > >

1 0

Weird problem of Coot fetch PDB&Map using PDB_redo
by Xiao Lei 07 Jul '15

07 Jul '15

Hi All, Today I found one weird problem of coot. I have coot installed in different platforms, I have coot 0.8.1 version in Ubuntu 15.04 and coot 0.8.1 version in Mac OS X 10.10.3, and I also have wincoot 0.8.1 in Windows 8.1. When I tried to fetch pdb&map using EDS, coot in all three platforms works very well. However, when I tried to fetch pdb&map using PDB_redo, only the windows wincoot works, which shows "successfully read coordinates...". The coot in mac and Ubuntu failed to fetch and shows "failed to read coordinates file...". I do not think the problem comes from PDB_redo as I tried to fetch the same pdb to compare coot in three platforms, only wincoot works well. I do not know if anyone else has this problem and how to fix it. Thanks ahead. Xiao

2 1

on determination whether a protein is a monomer or complex by phenix
by Smith Liu 06 Jul '15

06 Jul '15

Dear All, Suppose I have got a mtz file, and I know the aa sequence of the protein from which I got the crystal and the corresponding mtz file, will you please tell me which phenix function (program) can tell me whether my protein crystal was formed by protein monomer or by protein complex, and if formed by protein complex, how many subunits existing in the complex? Best regards. Smith

5 6

Re: [phenixbb] phenix rosetta refine
by kaiser 06 Jul '15

06 Jul '15

Thanks, Nigel! <div>-------- Original message --------</div><div>From: Nigel Moriarty <nwmoriarty(a)lbl.gov> </div><div>Date:2015/07/06 14:11 (GMT-08:00) </div><div>To: Jens Kaiser <kaiser(a)caltech.edu> </div><div>Cc: phenixbb(a)phenix-online.org </div><div>Subject: Re: [phenixbb] phenix rosetta refine </div><div> </div>Jens I believe that is all true and correct. 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 On Mon, Jul 6, 2015 at 1:55 PM, Jens Kaiser <kaiser(a)caltech.edu> wrote: Nigel, Sorry for the late reply, I was away from the office last week. Yes, I did follow these instructions. At the moment I'm using rosetta_bin_linux_2015.22.57859_bundle and phenix-1.10pre-2084. I requested 20 models (which I get) but there are rosetta_##.sh files created that are empty and subdirectories _ros_temp?????? that contain files with _adp_refine in their names. I'm just wondering if it is really running as intended. Cheers, Jens On Mon, 2015-06-29 at 10:25 -0700, Nigel Moriarty wrote: > Jens > > > Have you followed the instructions? > > > http://www.phenix-online.org/documentation/reference/rosetta_refine.html > > > 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 > > On Sat, Jun 27, 2015 at 2:38 AM, Jens Kaiser <kaiser(a)caltech.edu> > wrote: > Hey there, > who would be the person to talk about setting up > phenix.rosetta_refine? I tried our SBgrid installation with > overrides to > "nightly" for phenix and rosetta; I tried a clean install of > rosetta and > phenix - production version and nightly's in all combinations > - and all > combinations of compiling (rosetta first, then > rosetta.build_phenix_interface, phenix first, then > rosetta.build_phenix_interface). The problem is that > "phenix.rosetta_refine" seems to run, and produce a log file > that looks > reasonable, but geometry params are worse than with an > ordinary refine. > I think I tracked it down to the fact that the "rosetta_XX.sh" > files are > empty, and rosetta actually does not run. > > Any help is welcome - and any support will be extended. > > Cheers, > > Jens > > _______________________________________________ > 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 rosetta refine
by Jens Kaiser 06 Jul '15

06 Jul '15

Hey there, who would be the person to talk about setting up phenix.rosetta_refine? I tried our SBgrid installation with overrides to "nightly" for phenix and rosetta; I tried a clean install of rosetta and phenix - production version and nightly's in all combinations - and all combinations of compiling (rosetta first, then rosetta.build_phenix_interface, phenix first, then rosetta.build_phenix_interface). The problem is that "phenix.rosetta_refine" seems to run, and produce a log file that looks reasonable, but geometry params are worse than with an ordinary refine. I think I tracked it down to the fact that the "rosetta_XX.sh" files are empty, and rosetta actually does not run. Any help is welcome - and any support will be extended. Cheers, Jens

2 3

Re: [phenixbb] Optimal settings for finding NCS from density
by Terwilliger, Thomas Charles 05 Jul '15

05 Jul '15

Hi William, For find_ncs_from_density the shape of the molecule itself is not so important. What matters is if there is some local part of the molecule that is recognizable in both copies. So I would try default parameters first. I assume you did that and it didn't give you a useful answer...so then I would do these things: 1. try many locations : n_center_use=20 This is the best option...it will try 20 times using density at 20 different locations to search. If you have an older version of phenix you may need to download a nightly build to do this. 2. Try varying the resolution 3. Try varying the angular spacing of the search (delta_phi) All the best, Tom T ________________________________ From: CCP4 bulletin board [CCP4BB(a)JISCMAIL.AC.UK] on behalf of William Chao [William.Chao(a)CRICK.AC.UK] Sent: Friday, July 03, 2015 10:24 AM To: CCP4BB(a)JISCMAIL.AC.UK Subject: [ccp4bb] Optimal settings for finding NCS from density Dear all, What would be the best settings to look for NCS using phenix.find_ncs_from_density on an elongated (~150A) molecule? No coordinate is available unfortunately and there should be 2 copies in a P21 space group. Many thanks in advance! William --- The Francis Crick Institute Limited is a registered charity in England and Wales no. 1140062 and a company registered in England and Wales no. 06885462, with its registered office at 215 Euston Road, London NW1 2BE.

2 2

labelit access in Phenix?
by Keller, Jacob 03 Jul '15

03 Jul '15

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

3 2

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