- phenixbb - phenix-online.org

SA Fo-Fc omit map
by Mengbin Chen 12 Dec '12

12 Dec '12

Dear Phenix Users, I have been playing around with Phenix AutoBuild GUI to generate a fo-fc sa omit map of the ligand in a protein structure, and I have the input files as follows: output.mtz (data type:experimental data) ligand.pdb (data type: omit box atoms & ligands) ligand.cif (data type: restraints of the ligand) protein.pdb (data type: starting model) omit map type: simulated annealing omit region: omit around PDB fie I have set the occupancies of all ligand atoms zero and deleted the ligand coordinates in my starting model. After the run finished and I checked the refined coordinates and map.mtz in COOT, I found that there was no electron density around the ligand while the protein model was in electron density, which was exactly opposite of what I wanted before. Furthermore, I have 3 Mg2+ ions in protein.pdb file, all of which moved away from where they were after the refinement run, which was not wanted either. Could anyone tell me how to do the sa omit map of ligand correctly? Thanks in advance! Mengbin -- Mengbin Chen Department of Chemistry University of Pennsylvania

2 1

Re: [phenixbb] Fwd: parameters for averaged kick maps
by Pavel Afonine 10 Dec '12

10 Dec '12

Hi Jose, kick-map related parameters in ordered_solvent scope of parameters are leftover that I forgot to remove back when I was working on it. So you can see it as a minor bug. And yes, you can't change them. Sorry about delayed answer (I saw the email on phenixbb but did note reply immediately and it sank in the flood of other stuff). Pavel On 12/10/12 11:04 AM, Jose M de Pereda wrote: > Dear Phenix Team, > > Last week I sent an e-mail to the phenixbb regarding the parameters > that control de generation of averaged kick maps. Unfortunately, I > have not received any feedback from the community. > > I wonder if you could give me some help on this subject. > > Probably, the calculation of averaged kick maps is already very well > optimized in phenix.refine and phenix.maps. Yet, I am still curious. > > Thanks a lot in advance > All the best > > Jose > > > -------- Original Message -------- > Subject: [phenixbb] parameters for averaged kick maps > Date: Tue, 04 Dec 2012 10:12:46 +0100 > From: Jose M de Pereda <pereda(a)usal.es> > Reply-To: PHENIX user mailing list <phenixbb(a)phenix-online.org> > To: PHENIX user mailing list <phenixbb(a)phenix-online.org> > > > > Dear all, > > I would like to compare averaged kick maps obtained using different “kick sizes”. > Is it possible to define the kick size and the number of kicks in phenix.refine and/or phenix.maps? > > In phenix.refine, the “ordered_solvent” section has keywords to control the generation of kick maps: > kick_map { > kick_size = 0.5 > number_of_kicks = 100 > } > But, it seems to me that these can not be used during the generation of maps or maps’ coefficients. > > Related with this issue, does phenix.refine always use the same kick > size or does it adjust it ad hoc (e.g. based on ML error estimates)? > > Thanks a lot, > Jose > > > -- > ------------------------------------------------------------ > Jose M de Pereda, PhD > Instituto de Biologia Molecular y Celular del Cancer (IBMCC) > Spanish National Research Council - University of Salamanca > Campus Unamuno s/n > E-37007 Salamanca, Spain > Phone: +34-923-294817 > Fax: +34-923-294795 > http://xtal.cicancer.org/ > ------------------------------------------------------------ > >

1 0

Simulated annealing map
by Appu kumar 10 Dec '12

10 Dec '12

Thank you Terwilliger, Thomas for you quick reply, for guiding me on how to make simulated annealing omit map. Now i am able to make it but i have one question. i want to show the SA omit map of ligand with the map of omitted region around ligand. I tried it, i got the map also, but the omitted map of region around ligand is show in recatangular box which is not as great as of ligand inside the omitted region. Thank you

2 1

Please apply for RapiData 2013, a course on Data Collection and Structure Solving at the NSLS.
by Robert Sweet 10 Dec '12

10 Dec '12

We are offering RapiData 2013, the fifteenth offering of our popular course: Rapid Data Collection and Structure Solving at the NSLS: A Practical Course in Macromolecular X-Ray Diffraction Measurement The course will be held 21-26 April 2013: http://www.bnl.gov/RapiData/. Students could be at any level from advanced undergraduate to full professor. The course should accommodate 48 students total. We encourage all students to bring their own specimens for data collection, and to bring old data for the data-reduction and structure-solving tutorials. Please read the Course Announcement at http://www.bnl.gov/RapiData/. You'll see that many experts in the field will be available for lectures and tutorials. You'll find the application materials on the Course Application tab at this site. For the eleventh time we will hold a short lecture course on the fundamentals of crystallography for roughly five hours on Sunday 21 April. The body of the RapiData course really requires that students have a healthy knowledge of crystallography. For potential students who have some experience but are shaky about fundamentals, this course will help. There will be a small additional fee for the fundamentals course, to pay for Saturday night accomodations and food on Sunday morning and noon. Latin American Scientists: Several scholarships are available, from the International Union of Crystallography, to pay partial travel and subsistence costs for Latin-American students and junior faculty (under 40 yrs). Please apply for the course, and then contact R. Sweet (sweet(a)bnl.gov) if you are interested in applying for a scholarship. In accordance with the standards of the International Union of Crystallography, we observe the basic policy of non-discrimination, affirming the right and freedom of scientists to associate in international scientific activity without regard to such factors as citizenship, religion, creed, political stance, ethnic origin, race, colour, language, age, or gender, in accordance with the Statutes of the International Council for Science. At this course no barriers will exist beyond the application procedure that would prevent the participation of bona fide scientists. Please apply or send your students to our course, Bob Sweet, Sonya Kiss, and Alex Soares Course Announcement at http://www.bnl.gov/RapiData/ ========================================================================= 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/ Photon Sciences and Biology Dept Mail Stop, Bldg 463 Brookhaven Nat'l Lab. Phones: Upton, NY 11973 631 344 3401 (Office) U.S.A. 631 344 2741 (Facsimile) =========================================================================

1 0

Real Space CC Raw Data
by Eddie Pryor 08 Dec '12

08 Dec '12

Hi, All, To satisfy a reviewer, I would like to include a supplemental figure that contains, among other things, Real Space CC for each residue. I was wondering where this raw data might be stored so that I can parse and create a graph with it. The only thing I could find was in the Refine_###_info.txt log file, but that only included values less than 0.7. I know that all of the residues / CC values can be listed through the GUI, but with ~900 residues in my ASU, this is something I do not want to do manually! Thanks! --Eddie Pryor

2 2

How to make simmulated annealing omit map
by Appu kumar 07 Dec '12

07 Dec '12

Dear Users, I am refining a enzyme structure with ligand. I want to make the simulated annealing omit map of ligand including with 4 Angestron radius around the ligand. I am seeking your valuable advice to help me out. I know how to make omit map but i have not come accross making simmulated annealing omit map. I would appreciate your help. Thank you Appu

2 1

Is this a Phaser bug or am I misinterpreting ?
by Phil Jeffrey 05 Dec '12

05 Dec '12

Phaser from CCP4 6.3.0 built via fink running on OSX 10.6.8 in 64-bit mode. This is a snippet from a much larger run, ongoing. In multiple instances the order of the LLG and the Z-score in the translation function are not in sync: (42.16, 40.23, 40.18 vs 5.14, 5.24, 4.82) and looks like the Z-scores are still shown in the peak rank order before rescoring, while the LLGs are after rescoring. SET #5 of 6 TRIAL #60 of 115 ---------------------------- Search Euler = 263.0 82.1 85.7, Ensemble = ensemb3 ANNOTATION: RFZ=3.6 TFZ=3.8 PAK=0 LLG=26 LLG=30 Known MR solutions SOLU SPAC C 1 2 1 SOLU 6DIM ENSE ensemb3 EULER 30.6 65.7 246.9 FRAC 0.90 0.00 0.17 BFAC -5.14 Grid sampling: 0.825029 Angstroms Select peaks over 67.5% of top (i.e. 0.675*(top-mean)+mean) Top 360 translations before clustering will be rescored Calculating Likelihood for TF SET #5 of 6 TRIAL #60 of 115 0% 100% |=========================================================================| DONE Scoring 1 randomly sampled translations Generating Statistics for TF SET #5 of 6 TRIAL #60 of 115 0% 100% |==| DONE Top Peaks With Clustering ------------------------- # Rank of the peak after rescoring search points (#) Rank of the peak before rescoring search points LLG Log-Likelihood Gain Z-Score Number of standard deviations of LLG above the mean FSS Fast Search Score Select all peaks There were 159 peaks # (#) Frac X Frac Y Frac Z LLG Z-score Split #Group raw/top 1 6 0.842 0.923 0.779 +42.16 5.14 0 2 52.58/ 52.58 2 3 0.369 0.865 0.553 +40.23 5.24 37 2 53.18/ 53.18 3 10 0.743 0.937 0.261 +40.18 4.82 48 2 50.69/ 50.69 #SITES = 159: OUTPUT TRUNCATED TO 3 SITES

2 1

Re: [phenixbb] log files for refinement (was: completeness using Intensities vs. structure factor)
by Pavel Afonine 05 Dec '12

05 Dec '12

Hi Yarrow, thanks for sending log files. Completeness reported in log files reflects the content of input MTZ data files: *32.log:* (...) F-obs: (...) Miller array info: (...)/reflectionsFreeR_1.38A.mtz:F_New(+),SIGF_New(+),F_New(-),SIGF_New(-),systematic_absences_eliminated Observation type: xray.amplitude Type of data: double, size=325598 Type of sigmas: double, size=325598 Number of Miller indices: 325598 Anomalous flag: True Unit cell: (64.3967, 68.0966, 103.663, 90, 96.1167, 90) Space group: P 1 21 1 (No. 4) Systematic absences: 0 Centric reflections: 4910 Resolution range: 43.7941 1.3786 Completeness in resolution range: 0.901282 Completeness with d_max=infinity: 0.901259 (...) and *35.log:* (...) I-obs: (...) Miller array info: (...)/reflectionsFandIFreeR_1.38A.mtz:IMEAN_New,SIGIMEAN_New,systematic_absences_eliminated Observation type: xray.intensity Type of data: double, size=181940 Type of sigmas: double, size=181940 Number of Miller indices: 181940 Anomalous flag: False Unit cell: (64.3967, 68.0966, 103.663, 90, 96.1167, 90) Space group: P 1 21 1 (No. 4) Systematic absences: 0 Centric reflections: 4910 Resolution range: 43.7941 1.3786 Completeness in resolution range: 0.992202 Completeness with d_max=infinity: 0.992164 (...) The numbers above are reported before anything else is done, that is they reflect exactly what's in input MTZ files. Somehow you managed to loose ~10% of data when converting Iobs into Fobs. How did you do this conversion? In theory this is possible since French&Wilson Iobs->Fobs conversion may result on some number of reflections lost. Note, when phenix.refine does IMEAN->Fobs conversion (35.log file) this does not result in ~10% loss. So based on looking at log files I can say it's not phenix.refine that manipulated your data or reported it wrong. You lost reflections at previous steps. Pavel On 12/4/12 4:57 PM, Yarrow Madrona wrote: > Hi, > > Here are the log files for refinement you asked for. > > 35.log is refined with Intensities. The MTZ file has a slightly different > name but it is the same MTZ. One had Intensities taken out of it because I > didn't think I needed them any longer. > > 32.log is refined with structure factors. Intensities were removed. > > Thanks for your help! > > -Yarrow > > >

1 0

completeness using Intensities vs. structure factor
by Yarrow Madrona 05 Dec '12

05 Dec '12

Hi, (Phenix v.1.8.1-1168) I recently refined a structure using experimental Intensities and compared it to the same structure refined with calculated structure factors. I know if you provide intensities Phenix.refine will convert to structure factors. I am not sure if it runs truncate or what. Anyway, I get this in the log file which matches my completeness in the scala log file: Completeness in resolution range: 0.992202 Completeness with d_max=infinity: 0.992164 When I run using structure factors (Truncate run in Scala) I get the following: Completeness in resolution range: 0.901282 Completeness with d_max=infinity: 0.901259 Perhaps Scala did not merge anomolous and non-anomolous data? This does not may a whole lot of sense to me. Is some data thrown out? Can someone explain? Thank you. -- Yarrow Madrona Graduate Student Molecular Biology and Biochemistry Dept. University of California, Irvine Natural Sciences I, Rm 2403 Irvine, CA 92697

4 5

parameters for averaged kick maps
by Jose M de Pereda 04 Dec '12

04 Dec '12

Dear all, I would like to compare averaged kick maps obtained using different “kick sizes”. Is it possible to define the kick size and the number of kicks in phenix.refine and/or phenix.maps? In phenix.refine, the “ordered_solvent” section has keywords to control the generation of kick maps: kick_map { kick_size = 0.5 number_of_kicks = 100 } But, it seems to me that these can not be used during the generation of maps or maps’ coefficients. Related with this issue, does phenix.refine always use the same kick size or does it adjust it ad hoc (e.g. based on ML error estimates)? Thanks a lot, Jose -- ------------------------------------------------------------ Jose M de Pereda, PhD Instituto de Biologia Molecular y Celular del Cancer (IBMCC) Spanish National Research Council - University of Salamanca Campus Unamuno s/n E-37007 Salamanca, Spain Phone: +34-923-294817 Fax: +34-923-294795 http://xtal.cicancer.org/ ------------------------------------------------------------

1 0