- phenixbb - phenix-online.org

Re: [phenixbb] Experimental Data Question
by Yuri 16 Jul '11

16 Jul '11

I am not positive, because I had to convert my reflections.sca (I SIGI) into reflections.mtz (F SIGF) -- Yuri Pompeu

2 1

Experimental Data Question
by Yuri 15 Jul '11

15 Jul '11

Hello All, Heres my situation: I collected data on a new crystal, indexed, scaled (2.75 A), ran molecular replacement (40% homologue) and found a solution. I did some building got a new model, did an initial refinement (using my .sca file, not data.mtz from mol rep). I then ran phase_and_build with density modification. The new maps look a lot better and I am able to see some residues I was not able to place before. I did some more building and I am ready to refine. What should I use as experimental data? phase and_build gave me a .mtz file with amplitudes (no I) and "experimental phase" coeeficients (HLAM. HLBM,...) Cheers -- Yuri Pompeu

2 1

Re: [phenixbb] ADP refinement selection
by Pavel Afonine 15 Jul '11

15 Jul '11

Hi Christian, the problem is that you did not ask phenix.refine to refine the rest of the atoms isotropically (with your command below you only refine anisotropic ADPs of atoms in chain Q, and the B-factors of all the other atoms are not refined): > adp { > individual { > anisotropic = "chain Q" > } I would do it like this: adp { individual { anisotropic = chain Q and element CL isotropic = not (chain Q and element CL) } } Pavel.

2 2

off-topic-Question
by Haytham Wahba 14 Jul '11

14 Jul '11

Dear all i have protein soaked in Dimethyltindibromide in the active site 2 cysteine and 1 aspartic acid by help of anomalous map i can see 2 Sn atoms in each subunit in AU there is difference in both subunits, Sub-A shows acetate near Sn and Sub-B show phosphate ion, even the position of one of the 2 cystein is different from unit to other 1-is it normal to see such difference between 2 subunit? also the distances between Sn and surrounding atoms does not suggest covalent binding 2-the important question is how to confirm if the 2 methyl groups or even one methyl group of dimethlytin are still there or cleaved i mean one can not differentiate between density of water and methygroup. How to define coordinate number of Heavy atom and identity of the attached atoms as you will see in the link below (with water) and (with Methyl) the map are the same and as well r-factors 3-what can i do to confirm presence or cleavage of methyl: collect more data or using other technique like exafs, Mossbauer spectroscopy could help (with water) http://dl.dropbox.com/u/12001878/Subunit-A.png http://dl.dropbox.com/u/12001878/subunit-A-anomalous.png http://dl.dropbox.com/u/12001878/Subunit-B.png http://dl.dropbox.com/u/12001878/subunit-B-anomalous.png (with methyl) http://dl.dropbox.com/u/12001878/DimethylTin-subunit-A-alt-conf.png http://dl.dropbox.com/u/12001878/DimethylTin-subunit-B.png N.B: resolution 1.7 A r factors and molprop are perfect. thank you for your time haytham biochimi canada

1 0

Parsing error
by Yuri 14 Jul '11

14 Jul '11

Nat, Thank you. I am still using dev-814. You probably have fixed it already. ill update when I am done running these jobs. -- Yuri Pompeu

1 0

Parsing error
by Yuri 14 Jul '11

14 Jul '11

Hello everyone, I cannot figure out why I keep getting this error. I dont see anything wrong with my PDB. I realize I realize I have breaks in my chain, as there are regions of unclear density as of now I am building in what I can see... Heres what the log says: ======== ================= Separating polymer chains into rigid bodies ================= (chain 'A' and segid ' ' and resseq 393:6) (chain 'B' and segid 'B ' and resseq 401:401) ================== Extract refinement strategy and selections ================= The error is initial RESSEQ index > final RESSEQ index. I do not understand why it is picking them this way It then halts at an error... Any Ideas?

3 4

restraints_edite.params
by Haytham Wahba 13 Jul '11

13 Jul '11

Dear All How to get rid of restraints_edite.params in the following run of phenix.refine, i mean after using it in previous run.and what is .def files. and bond.slack Haytham Udem Biochemi Canada

4 3

R-merge and MOLREP
by Yuri 13 Jul '11

13 Jul '11

Ed and Nathaniel, Thank you for the input. I agree with Ed. The I/sigma values for this data set are well within accepted ranges... I believe the R-merge values are due to some inherent crystal property. They grew in very strange manner. -- Yuri Pompeu

1 0

R-merge and MOLREP
by Yuri 13 Jul '11

13 Jul '11

Hello everyone, I just collected some data on a new protein. It indexes in P3 and looks like I can scale it all the way to P 63 2 2 . My questions are: 1-I see a reasonable chi-square value 1.69 but r-merge 0.15 (0.49 at 2.75 A). It was 0.20 before I excluded some bad frames. I thought this R value is a bit high... Is it? Mosaicty was around 1.00. 2- My search model has an identity of around 40%. I tried MOLREP. Seems I got a solution LLG=486 PAK=10 TFZ ~30 RFZ ~20. My R-free before any refinement is 0.54. Two questions here, a) I think I only fully understand the meaning of LLG, what about PAK and both Z-values? b)is this R-free acceptable or is it also a little high? Thank you for any input

3 2

Summer X-rays at the NSLS
by Robert Sweet 12 Jul '11

12 Jul '11

Summer X-rays at the NSLS The PXRR operates five beamlines at the NSLS: http://www.px.nsls.bnl.gov/. Two of these are the undulator beamlines X25 and X29. The dipole beamlines X12-B and X12-C are available, and are very effective for conventional crystallography. There are efficient ALS-Style Automounters available -- 8400 crystals mounted in the last 12 months. Also, at beamline X26-C, we have a unique-in-the-US station for combined x-ray diffraction and single-crystal optical spectroscopies: http://www.px.nsls.bnl.gov/x26c/x26c_info.html Notice also our world-leading Mail-In crystallography program: http://www.px.nsls.bnl.gov/Mailin.html You may request access to any of these here: http://www.px.nsls.bnl.gov/beamtime_requests.html, and can look at the schedules for all five here: http://www.px.nsls.bnl.gov/ Stay tuned for the announcement in December of RapiData, our annual course in rapid data collection and structure solving: http://www.bnl.gov/rapidata/ Bob Sweet ========================================================================= Robert M. Sweet E-Dress: sweet(a)bnl.gov Group Leader, PXRR: Macromolecular ^ (that's L Crystallography Research Resource at NSLS not 1) http://px.nsls.bnl.gov/ Biology Dept Brookhaven Nat'l Lab. Phones: Upton, NY 11973 631 344 3401 (Office) U.S.A. 631 344 2741 (Facsimile) =========================================================================

11 19

significant figures in reporting Rfree
by Jennifer Weinreich 12 Jul '11

12 Jul '11

Hi. I am wondering what is the common practice with respect to how many significant figures should be used for reporting Rfree at different resolution? Should we always use all 4 figures that are reported by the program? Thank you, Jenn

4 4

How to rebuild a molecular only‏
by dingding830106＠sina.com 09 Jul '11

09 Jul '11

Dear All, There are 8 moleculars in an asymmetry unit, but only one molecular should be rebulit, what can I do ? Thanks a lot! Wei Feng

2 1

Re: [phenixbb] Elbow and smiles to make CO coordinates
by Yarrow Madrona 08 Jul '11

08 Jul '11

Thanks! Didn't realize it was so easy. > For getting carbon monoxide coordinates: In coot, go to File > Get > Monomer... > Type in "CMO" and click OK. > > A carbon monoxide molecule will po up in the middle of the window > (provided everything is set up correctly). > > If you are trying to fit an electron density, place that density in the > middle of the window before you retrieve the ligand. > > Hope that helps. > > Cheers! > MM > > On Jul 8, 2011, at 1:41 PM, Yarrow Madrona wrote: > >> >> Hi, >> >> I am trying to use Elbow and smiles to generate carbon monoxide >> coordinates. I thought it would be easy enough. I made two waters in >> coot >> and changed the the oxygens to C and O and set the distance around 1.13 >> angstrom. Saved these coordinates as a PDB. I then used this in Elbow >> but >> it gives me a methyl group. Not sure what is going on. Any help would >> be >> appreciated. Thank you! >> >> -Yarrow >> >> -- >> Yarrow Madrona >> >> Graduate Student >> Molecular Biology and Biochemistry Dept. >> University of California, Irvine >> Natural Sciences I, Rm 2403 >> Irvine, CA 92697 >> >> >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://phenix-online.org/mailman/listinfo/phenixbb > > > _______________________________________________ > phenixbb mailing list > phenixbb(a)phenix-online.org > http://phenix-online.org/mailman/listinfo/phenixbb > > -- Yarrow Madrona Graduate Student Molecular Biology and Biochemistry Dept. University of California, Irvine Natural Sciences I, Rm 2403 Irvine, CA 92697

1 0

Elbow and smiles to make CO coordinates
by Yarrow Madrona 08 Jul '11

08 Jul '11

Hi, I am trying to use Elbow and smiles to generate carbon monoxide coordinates. I thought it would be easy enough. I made two waters in coot and changed the the oxygens to C and O and set the distance around 1.13 angstrom. Saved these coordinates as a PDB. I then used this in Elbow but it gives me a methyl group. Not sure what is going on. Any help would be appreciated. Thank you! -Yarrow -- Yarrow Madrona Graduate Student Molecular Biology and Biochemistry Dept. University of California, Irvine Natural Sciences I, Rm 2403 Irvine, CA 92697

3 2

AutoSol question
by Laurie Betts 08 Jul '11

08 Jul '11

If I want to run Autosol with a set of previously determined heavy atom positions (seleno-Met protein) that I might have edited, and NOT search for any more sites, just use those sites for phasing and density modification and model building, I can't seem to find in any of the parameter buttons a way to do this - every time I try it seems to just re-run hyss and find the same sites again. It seems to find multiple Se sites on top of each other that are closer than the minimum distance I specified in the input. Thanks LAurie Betts

2 1

NCS Question
by Yuri 07 Jul '11

07 Jul '11

Dear PHENIX community, Could someone please point me toward a good paper on NCS? Cheers, -- Yuri Pompeu

3 2

SA omit map
by Tjaard Pijning 06 Jul '11

06 Jul '11

Hi, I would like to generate an SA-omit Fo-Fc map for a ligand bound to protein. Using the GUI I selected the AutoBuild-Create Omit Map module and set the following : - data.mtz - protein.pdb (no ligand, no solvent) = start model - ligand.pdb (just ligand) = omit map atoms - omit map type = simulated annealing - omit region = omit around pdb The resulting map (/OMIT/resolve_composite_map.mtz) shows density for both the protein and the ligand. When I feed this into the CCP4 module FFT to generate an nF1-mF2 map (with n=1 and m=1) I still get density for both the protein and ligand. However I need only the ligand Fo-Fc density. What is the right way to do this ? Thanks, Tjaard Ing. Tjaard Pijning Research Assistant Protein Crystallography Group Nijenborgh 7 9747 AG Groningen The Netherlands (31)(0)50 363 43 85

4 6

Re: [phenixbb] phenix.mr_rosetta: Inconsistent query numbering
by Thomas C. Terwilliger 06 Jul '11

06 Jul '11

Hi Pietro, There is something that phenix.mr_model_preparation cannot read in Japanin.hhr (attached): phenix.mr_model_preparation Japanin.hhr gives the same error. I'll copy Gabor Bunkoczi and see if he can help with this. All the best, Tom T >> Dear all, >> >> my run of phenix.mr_rosetta stops with: >> >> ******************* ERROR ENDING *************** >> Inconsistent query numbering >> >> ******************* ERROR ENDING *************** >> >> I attach script and logfiles. >> >> Any suggestions? >> >> Thanks >> >> Pietro >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://phenix-online.org/mailman/listinfo/phenixbb >>

3 2

Validation question: R-factor (all) and # of unique reflections (all)
by Young-Jin Cho 05 Jul '11

05 Jul '11

Hi everyone, While I am trying to deposit, I am faced with several things that I am not sure: FYI. I am in the <Refinement Statistics> section (Auto Dep Input Tool) in the PDB bank. Most values are came from pdb file that was made from phenix.refine but it asks me "R-factor (all)" that satisfy the limits established by high & low resolution cut off and "number of unique reflections (all)". My simple question is if phenix has those kind of information, if so where I can find them. Another question is if I changed water (deletion or changed occupancy) based on 'precheck', should I run phenix.refine again to get R/Rfree values in correct? I am not sure if these are captured by PDB automatically or not. (seemingly they do like that at this point but not 100% sure yet.) With many thanks~ Young-Jin

4 3

Re: [phenixbb] Validation question: R-factor (all) and # of unique reflections (all)
by Pavel Afonine 05 Jul '11

05 Jul '11

Hi Young-Jin, > I feel like I should calculate R(all). As mtz file is usually big, > how can I extract Fobs and Fmodel from phenix generated mtz file? If you feel you want it, then go ahead, and this command might be of help here: phenix.mtz.dump data.mtz --column-data-format=spreadsheet -c Yet another (easier) alternative is to run a dummy phenix.refine run like this: phenix.refine data.mtz model.pdb xray_data.r_free_flags.ignore_r_free_flags=true main.number_of_mac=1 strategy=none output.prefix=junk then take the reported R-factor (Rwork will be equal to Rfree, so take any of them), and discard the files. All the best! Pavel.

1 0

phenix.mr_rosetta: Inconsistent query numbering
by Pietro Roversi 05 Jul '11

05 Jul '11

Dear all, my run of phenix.mr_rosetta stops with: ******************* ERROR ENDING *************** Inconsistent query numbering ******************* ERROR ENDING *************** I attach script and logfiles. Any suggestions? Thanks Pietro

1 0

Ligand woes
by José Trincão 05 Jul '11

05 Jul '11

Hello all, I have been having very weird problems refining a structure with several ligands. The weirdest problem is that if I try to run ligand fit giving just the name of the ligands, no problem - it finds the ligands (FAD and MTE). I save the pdb with all the ligands in and phenix.refine complains that MTE is not recognised, I have to run ready-set. I do it. Run phenix.refine again, and now it complains with the following: Number of atoms with unknown nonbonded energy type symbols: 116 "HETATM18439 C1' MTE A3003 .*. C " "HETATM18441 C2' MTE A3003 .*. C " "HETATM18443 C3' MTE A3003 .*. C " "HETATM18444 O3' MTE A3003 .*. O " "HETATM18474 N9A FAD A3005 .*. N " "HETATM18475 C8A FAD A3005 .*. C " "HETATM18476 N7A FAD A3005 .*. N " "HETATM18477 C5A FAD A3005 .*. C " "HETATM18478 C6A FAD A3005 .*. C " "HETATM18480 N1A FAD A3005 .*. N " ... (remaining 106 not shown) Any help will be greatly appreciated. Jose José Trincão, PhD CQFB@FCT-UNL 2829-516 Caparica, Portugal "It's very hard to make predictions... especially about the future" - Niels Bohr

3 2

NCS Restraints
by Yuri 01 Jul '11

01 Jul '11

Hi Pavel, >It's rather packing.. Say, if more than one copy packs in the > asymmetric >unit then you have NCS, roughly speaking.. So is it correct to say that in cases where you have 1 molecule per ASU, there will be no NCS only crystallographic symmetry?? > In other words, if for the same enzyme one crystal goes to 3.1 and > another 1.5 A could one have significant NCS but not the other? >I'm not sure I understood this... May some one else did? I was trying to understand if this phenomenon would be noticeable depending on the resolution only... You and Nathaniel pretty much answered this now... -- Yuri Pompeu

2 1

NCS restraints
by Yuri 01 Jul '11

01 Jul '11

Hello, I have some fundamental questions about NCS: 1- I think its somewhat of a consensus that at higher res. it should not be used to avoid missing out "real" detailed discrepancy. Is this true? 2- Should one try to use it in lower res data (3 A and worse)? 3- Is the presence of NCS macromolecule specific, or crystal specific? In other words, if for the same enzyme one crystal goes to 3.1 and another 1.5 A could one have significant NCS but not the other? Thanks or any clarification -- Yuri Pompeu

3 3

phenix.ligand_identification problem
by Andrzej 29 Jun '11

29 Jun '11

Dear All, I'm having problem with running ligand_identification script. All well up to the point where the ligand library is created and then for each ligand following error is displayed: Running job sequence 13, ligand 14, in /share/scratch/an... Sorry, resolve was not able to finish on the first cycle I suspect it has something to do with resolve size but there is no keyword like resolve_size for ligand_identification. I've tested my input files with 1.7.1-743 and dev-738 release - same results. Thanks in advance for the support! Regards, Andrzej -- Dr. Andrzej Lyskowski (andrzej[.]lyskowski[(a)]acib[.]at) ACIB Austrian Centre of Industrial Biotechnology GmbH Petersgasse 14, A-8010 Graz, Austria http://acib.at/ c/o IMB Institute of Molecular Biosciences, University of Graz Humboldtstraße 50,III/194, A-8010 Graz, Austria http://strubi.uni-graz.at/ +43 316 380 1993 tel +43 316 380 9897 fax

2 1