- phenixbb - phenix-online.org

Re: [phenixbb] phenixbb Digest, Vol 98, Issue 2
by Kay Diederichs 14 Jan '14

14 Jan '14

Am 14.01.14 21:00, schrieb phenixbb-request(a)phenix-online.org: > PS. It's possible I'm still misunderstanding the meaning of FRIEDEL'S_LAW=TRUE. > Under what circumstances does XDS output anomalous vs. non-anomalous data? This keyword affects the way the statistics are calculated in CORRECT.LP: with FRIEDEL'S_LAW=FALSE the members of Bijvoet pairs are treated as separate reflections, which almost doubles the number of unique reflections, affecting the reported completeness and R-values, I/sigmas and so on. The reflection output file XDS_ASCII.HKL has unmerged intensities in all cases (whether FRIEDEL'S_LAW=TRUE or FALSE). Phenix should ignore the FRIEDEL'S_LAW setting in XDS_ASCII.HKL. hope this helps, Kay

1 0

XDS data
by Tim Gruene 14 Jan '14

14 Jan '14

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Dear all, I would like to run Phaser-EP through the phenix GUI. When I first uploaded the XDS_ASCII.HKL file, there was an error message about "No valid observed reflections found in this file". Subsequently I uploaded an unmerged hklf4 file. The GUI tells me I need to let phenix know about it being hklf4 by appending '=hklf4'. However, when I do so, there is an error message the file did not exist and the Data file line goes red. Where should I correctly append '=hklf4'? At the bottom, the message states it would be best to recover the original reflection data such as [...] XDS [...]. This suggests that the original XDS_ASCII.HKL could be used. Is this correct in if so, how do I put it into phenix? Best, Tim - -- - -- Dr Tim Gruene Institut fuer anorganische Chemie Tammannstr. 4 D-37077 Goettingen GPG Key ID = A46BEE1A -----BEGIN PGP SIGNATURE----- Version: GnuPG v1.4.12 (GNU/Linux) Comment: Using GnuPG with Icedove - http://www.enigmail.net/ iD8DBQFS1TDEUxlJ7aRr7hoRArKdAKDpXBrs4wsWMvWw0JNo8F746fY76ACfeFuG geGY6T6se8peE9dEGz8ETow= =AaLk -----END PGP SIGNATURE-----

3 6

Refining coordinates
by Subhani Bandara 10 Jan '14

10 Jan '14

Dear all, I have a protein crystallized with metal chelate complexes and some of the atoms of chelators and the metals sit on the two fold axis. Is there a way that I can refine the coordinates of those atoms only along the two fold axis (let's say X) while restraining the coordinate refinement along other two axes (Y and Z)? Thanks Subhani

2 1

fwd: Announcement: International Summer School of Crystallography - May 2014, Hamburg, Germany
by Nathaniel Echols 10 Jan '14

10 Jan '14

---------- Forwarded message ---------- From: "Oberthuer, Dominik" <dominik.oberthuer(a)cfel.de> To: phenixbb(a)phenix-online.org Date: Fri, 10 Jan 2014 16:14:08 +0100 Subject: Announcement: International Summer School of Crystallography - May 2014, Hamburg, Germany Dear colleagues, we are pleased to announce the first International Summer School of Crystallography at the Center for Free-Electron Laser Science at DESY, Hamburg (D). As part of the International Year of Crystallography 2014, Prof. Carmelo Giacovazzo will give a full week of lectures on the basics of crystallography, as well as the basics of modern phasing methods in protein crystallography. The course will cover basic mathematical understanding of crystallographic point and space groups, diffraction experiments, structure factor calculations, systematic absences, determination of space groups and many additional topics. After the discussions of theory, one afternoon will be spent to learn about modern applications of protein crystallography using advanced light sources, such as 3rd generation synchrotron and XFELs. The course will be given mainly for PhD students working in the field of crystallography, but exceptions for motivated PostDocs and MSc students can be made. The course will mostly cover mathematical background of crystallography, hence we are looking for students, having a (bio)physics background. This will NOT be a workshop covering protein crystallographic programs but rather the theory behind those programs. The course will start Sunday (11.05.2014) around lunch and will end Saturday (17.05.2014) around lunch. The registration fee is 100 EUR per person including coffee breaks and lunch. Accommodations can be made at the DESY Guest House: 70 EUR per person for 7 nights (10.-17.05.2014). Application deadline will be 28th of February 2014 and notification of accepted applications will be sent during March 2014. The applications will be managed via the DESY Indico system: https://indico.desy.de//event/ISSC14 Applicants must provide a motivation letter, as well as a CV, directly to the organizer Cornelius Gati (cornelius.gati(a)cfel.de). Places and accommodation are limited to 30 students in total. Please contact Cornelius Gati regarding further questions ( cornelius.gati(a)cfel.de). We are looking forward to welcome you in Hamburg! -- Dr. Dominik Oberthür CFEL - Center for Free-Electron Laser Science Coherent Imaging Division and University of Hamburg Laboratory for Structural Biology of Infection and Inflammation Notkestrasse 85 22607 Hamburg Germany phone: +49 (0)40 8998 6394 dominik.oberthuer(a)cfel.de obertuer(a)chemie.uni-hamburg.de

1 0

Low F-obs Values
by 张怀东 06 Jan '14

06 Jan '14

Hi, everyone I used phenix-1.8.2-1309 to do MR and Autobuild, then Refinement, the max F-obs value of the mtz file from refinement was about 300. However, I used phenix-1.8.4-1496 to do the same work, the max F-obs value was about 3000. The resulting structure models donot seem to be any different. We tried several datasets of different proteins, and the max F-obs value are low with phenix-1.8.2-1309. So is this a real problem or not? Best regards Huaidong Zhang

2 1

I remind you: Please apply for RapiData 2014, a course on Data Collection and Structure Solving at the NSLS.
by Robert Sweet 06 Jan '14

06 Jan '14

We are offering RapiData 2014, the sixteenth 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 27 April - 2 May 2014: 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. In addition to MX there will be opportunity for learning about correlated MX and spectroscopy, and SAXS. You'll find the application materials on the Course Application tab at this site. For the twelfth time we will hold a short lecture course on the fundamentals of crystallography for roughly five hours on Sunday 27 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 an 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 (no more than 30 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 yourself, 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 Biosciences Dept Office and mail, Bldg 745, a.k.a. LOB-5 Brookhaven Nat'l Lab. Phones: Upton, NY 11973 631 344 3401 (Office) U.S.A. 631 344 2741 (Facsimile) =========================================================================

1 0

Postdoc positions (SGC, Diamond): fragments into crystals
by Frank von Delft 02 Jan '14

02 Jan '14

Hello, I have 3 open postdoc positions, to work on technologies around fragment screening in crystals, which is what my (joint) research team at Diamond (beamline I04-1) and SGC (PX group) is trying to make simple and accessible. Details of the project are on my website <http://www.thesgc.org/scientists/groups/oxford/crystallography>; we need smart people to: * (at SGC) test-drive with start-to-end soaking campaigns <http://www.thesgc.org/careers/oxford/crystallography-postdoc-fragment-scree…> * (at SGC) develop crystal engineering to get soakable systems (molecular biology, biochemistry) <http://www.thesgc.org/careers/oxford/crystallography-postdoc-crystal-engine…> * (at Diamond) explore soaking techniques <http://www.diamond.ac.uk/Home/Jobs/Current/DIA0886_CH.html> The application deadlines for the SGC posts are a function of the UK's charming political climate; if you're interested but from outside the EU, alert me with a direct email plus CV. phx. P.S. if you know smart undergrads who'd be interested in a DPhil in that sort of thing -- applied engineering, technology -- send them to that website too. (DPhils are what you get in Oxford - that makes it special. Or "special".)

1 1

Multiple conformation, two ligands
by YANG HAI 31 Dec '13

31 Dec '13

Dear Phenix Community, I tried to model two ligands one of which has two conformations. In other words, the two scenarios are: 1) ligand X with conformation A (XA, 75% occupancy); 2) ligand X with conformation B (XB, 25% occupancy) plus another ligand Y. However, the problem is ligand Y will have steric clash with ligand XA. If I rename ligand X and Y to the same residue name then I do not know how to define the CIF file. I apologize if this question was asked before and I do not know how Phenix or other refinement can deal with multiple conformation a group of ligands. Thanks. Happy new year! Yang Hai Christianson Lab University of Pennsylvania

3 2

phenix morph- no shift problems
by Ming Sun 29 Dec '13

29 Dec '13

Hi Everyone I'm currently trying running phenix.morph model on a large complex. However, after 5 cycles of running, there were no shift at all. And we set data, model, map_coeffs, reference_restraints, morph_main.. Does anyone happen to have similar issues? Any specific parameter I should change? Thanks again -Ming

2 1

Reduce NQH
by jp d 19 Dec '13

19 Dec '13

hi, we were going through the logs recently and see something odd during Reduce N/Q/H analysis the log file looks like this: ================================== NQH flips ================================== r_work=0.3047 r_free=0.3518 Analyzing N/Q/H residues for possible flip corrections... segID AL01 renamed chain A for Reduce N/Q/H analysis segID AL02 renamed chain B for Reduce N/Q/H analysis ... segID AL29 renamed chain Y for Reduce N/Q/H analysis segID AL30 renamed chain Z for Reduce N/Q/H analysis segID AL32 renamed chain a for Reduce N/Q/H analysis segID AL33 renamed chain a for Reduce N/Q/H analysis segID AL34 renamed chain b for Reduce N/Q/H analysis segID AL35 renamed chain b for Reduce N/Q/H analysis ... segID BL28 renamed chain z for Reduce N/Q/H analysis segID BL31 renamed chain z for Reduce N/Q/H analysis segID BL13 renamed chain AA for Reduce N/Q/H analysis segID BL29 renamed chain AB for Reduce N/Q/H analysis ... we noticed that the lower case chain letters start to duplicate themselves and are not unique to each segid. when it starts the two letter chain, everything seems unique again. forgive me if that is expected, it just seems kind of strange this was seen in mutlple log files and we are using version Version: 1.8.4 Release tag: 1496 if you need our files or more details, let me know thanks jpd

2 2

ensemble refinement rfree values
by Arjun Raman 19 Dec '13

19 Dec '13

Dear Phenix Guys, I have recently been trying my hand at the ensemble refinement option in phenix.refine. I generally observe an increase in rfree after performing the ensemble refinement with a well refined (using anisotropic B) structure, but the r-free increases are fairly small. For instance, 18.53% to 18.9% (1.8 Angstrom structure) or 20.45% to 21.1% (1.9 A)...I understand what this technically means (the ensemble refined model isn't as good) but given the noise in the calculation due to the limited number of reflections for the test set (~1000), would someone consider these rfree increases significant? My apologies if this question is a well-answered one...I can't seem to find much discussion about the range of rfrees to be expected or whether a 0.3-0.5% increase is large enough to toss out the refinement. Normally, I would take the lower rfree but because the ensemble refinement procedure gives potentially much more information, this seemed like a relevant question to ask. Thanks so much for your time, Arjun Raman Ranganathan Lab UT Southwestern ________________________________ UT Southwestern Medical Center The future of medicine, today.

3 2

MR-SAD, Yb, combined phases
by Joseph Noel 19 Dec '13

19 Dec '13

Hi All, I am carrying out a structure refinement using experimental phases combined with an MR solution. The experimental phases are being determined from bound Yb3+ cations and the data is collected at one of the Yb absorption peaks of 1.385 A. I do have a few questions about the heavy atom arrangement and then refinement. At this wavelength, I am also seeing peaks after running MR-SAD that are actually not the Yb cations but the sulfurs on several of the proteins Cys and Met residues. Should I try and change this "heavy atom arrangement" to include Yb and S before determining the final refined phases that I would like to then use as restraints during the structure refinement? If so, its not clear how to do this in Phenix without it finding additional sites that it attributes to Yb but which are actually sulfurs. For MR, this is an exact structure. All I am doing is using Yb which substitute for Mg coordination sites as a source of experimental phases to see how it improves refinement and maps. Second, when I do refine using the MLHL target function, I am finding quite a bit of excess positive density around the Yb sites. How do I try and account for this? I have attempted to take care of it with ADP and occupancy refinement but that doesn't seem to be enough. The resolution is 2.2 A and the anomalous signal is quite good (alas that is why I am also seeing the sulfurs clearly in the anomalous difference map). Do I need to change the scattering table since the structure now includes at least 4 Yb cations per asymmetric unit (one 550 amino acid residue protein per asymmetric unit - 68% solvent)? Is it important to refine the anomalous groups? For the latter, one can't refine with phase restraints (MLHL) and anomalous groups together. I have to choose one or the other - no phase restraints and anomalous groups refined, or experimental phase restraints and no anomalous group refinement. My goal is to use additional experimental measurements - amplitudes and phases to improve refinement, map interpretation and ultimately ensemble refinement. Any help, comments, criticisms, etc will be very appreciated! Happy Holidays! Joe ______________________________________________________________________________________ Joseph P. Noel, Ph.D. Arthur and Julie Woodrow Chair Investigator, Howard Hughes Medical Institute Professor, The Jack H. Skirball Center for Chemical Biology and Proteomics The Salk Institute for Biological Studies 10010 North Torrey Pines Road La Jolla, CA 92037 USA Phone: (858) 453-4100 extension 1442 Cell: (858) 349-4700 Fax: (858) 597-0855 E-mail: noel(a)salk.edu Publications & Citations: http://scholar.google.com/citations?user=xiL1lscAAAAJ Homepage Salk: http://www.salk.edu/faculty/noel.html Homepage HHMI: http://hhmi.org/research/investigators/noel.html ______________________________________________________________________________________

3 3

phenix.rosetta_refine on OS X?
by Ben Eisenbraun 19 Dec '13

19 Dec '13

4 5

Please apply for RapiData 2014, a course on Data Collection and Structure Solving at the NSLS.
by Robert Sweet 15 Dec '13

15 Dec '13

We are offering RapiData 2014, the sixteenth 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 27 April - 2 May 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. In addition to MX there will be opportunity for learning about correlated MX and spectroscopy, and SAXS. You'll find the application materials on the Course Application tab at this site. For the twelfth time we will hold a short lecture course on the fundamentals of crystallography for roughly five hours on Sunday 27 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 an 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 (no more than 30 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 yourself, 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 Biosciences Dept Office and mail, Bldg 745, a.k.a. LOB-5 Brookhaven Nat'l Lab. Phones: Upton, NY 11973 631 344 3401 (Office) U.S.A. 631 344 2741 (Facsimile) =========================================================================

1 0

off topic - linux flavors for crystallography
by Todd Jason Green 12 Dec '13

12 Dec '13

Hello All- I have used a trusty box running fedora core 3 for nearly a decade and I am (finally) upgrading to a new workstation. In lieu of Nat's recent mention of new linux installers of phenix, I thought I would try to get a consensus on what linux flavor(s), I should consider installing on my new box. It will be used with emphasis on typical crystallography software including 3D graphics (3 GB NVIDIA Quadro K4000 card). If anyone would like to make a suggestion for an OS, I thank you in advance. Pointers to a web review can work too. Best- Todd

7 7

cyclic peptide refinement
by wei xia 12 Dec '13

12 Dec '13

Hi, I am trying to refine a structure with a cyclic peptide ligand.I create a 'TRANS' link between the C-ter and N-ter residues of the cyclic peptide during refinement.But there is obvious no bond formed between the two residues in the refined structure.I checked the geo file and found that there are only angle restraints but no bond restraints between the two residuesAnyone can help me out?Thanks in advance! Wei P.S. refinement.pdb_interpretation.apply_cif_link { data_link = TRANS residue_selection_1 = chain D and resname MVA and resid 5005 residue_selection_2 = chain D and resname ARG and resid 5001 }

3 3

phenix.find_ncs_from_density density_search.dump_ncs_density doesn't output mtz
by Francis Reyes 11 Dec '13

11 Dec '13

I'd like to see the region of my maps that obey the NCS operators and I was hoping for separate mtz's for each region given I had set density_search.dump_ncs_density = True. However, I'm not getting any output mtzs. Any ideas? Thanks, F --------------------------------------------- Francis E. Reyes PhD 215 UCB University of Colorado at Boulder

2 5

phenix.refine nproc
by de Waal, Parker 10 Dec '13

10 Dec '13

Hello Everyone, I recently compiled phenix 1.8.4-1496 on a CentOS 6.4 server and am having trouble using phenix.refine in parallel using the nproc flag which always runs as a single thread. Please find my input parameters here: https://gist.github.com/anonymous/13b9d4d39ed4d0c9622a I'm wondering if maybe OpenMP isn't configured properly on our server? Or maybe the correct version is not installed when I compliled phenix. If anyone has any input or experience with a problem such as this I would be extremely grateful. Cheers, Parker

3 3

Spatial filtering
by Machius, Mischa Christian 06 Dec '13

06 Dec '13

Potentially silly question: is spatial filtering used for increasing contrast in electron density maps? Is it used for any other kind of data manipulation in macromolecular crystallography, e.g., diffraction data processing, density modification, etc. I'm a bit fuzzy on that subject, and so I said 'potentially silly', because it could be that such approaches are totally common but have different names (such as local averaging/scaling and so on). Many thanks in advance. Cheers, MM

2 1

Submissions for next CCN newsletter
by Nigel Moriarty 06 Dec '13

06 Dec '13

Folks Calling for articles and short communications of interest to structural biologists. The deadline for submission in the January 2014 issue is 15 December 2013 but if you can send me an email we can make exceptions so you can use the holiday break to work on your drafts. The Computational Crystallography Newsletter (CCN) is an electronic newsletter for structural biologists, and is published online every 6 months. Feature articles, meeting announcements and reports can be submitted to me at any time for consideration. Submission of text by email or word-processing files using the CCN templates is requested. Past newsletters and a Microsoft Word template are available at www.phenix-online.org/newsletter. 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

1 1

DEN refinement
by Scott Thomas Walsh 05 Dec '13

05 Dec '13

Dear Phenix Users, Is it possible to use NCS restraints in DEN refinement like we can do in regular phenix.refine refinement? Thank you for your help. Cheers, Scott ************************************************ Scott T. R. Walsh, PhD Assistant Professor University of Maryland IBBR/CBMG 3127E CARB-2 9600 Gudelsky Drive Rockville, MD 20850 USA phone: (240) 314-6478 fax: (240) 314-6225 email: swalsh12(a)umd.edu<mailto:[email protected]>

2 1

new Linux installers
by Nathaniel Echols 05 Dec '13

05 Dec '13

Hi Linux/Phenix users-- We are now building installers on several "new" platforms, including Ubuntu 10.04 and several versions of CentOS (which isn't much different from Fedora, except that it's supported by the virtual machine software we're using). Since we don't actually run these OSes outside of the VMs, we're going to crowd-source the testing a little more than usual. Our goal is to retire our ancient Fedora 3 systems - these are old and slow and take up too much space, but it's also the most compatible with other distributions. I'm hoping that the CentOS 5 builds will be a suitable replacement, or perhaps CentOS 4 (but we can only build the 32-bit installer on that OS). In addition to Ubuntu, I'm particularly interested in hearing whether any of these builds runs smoothly on SuSE and Debian. The virtual machine software supposedly supports these OSes too but I haven't yet managed to install SuSE on it after several attempts. The latest builds are here: https://www.phenix-online.org/download/phenix/nightly/?version=dev-1558&sho… We're missing the 64-bit CentOS 5 build right now, but I hope to get that working in the next few days. thanks, Nat

1 0

Re: [phenixbb] reduce hangs
by Andreas Förster 05 Dec '13

05 Dec '13

Dear all, I'm not sure how much of my conversation with Jeff Headd and the comments from others that I received made it onto the mailing list. The problem with reduce being stuck in an infinite loop and refinement never completing has now a workaround, if no solution: Mark van Raaij pointed out that running refinement repeatedly while varying the number of cores would eventually avoid the infinite loop. I didn't observe that, but while running refinement with a handful of Xeons got stuck, running the same job with the same executable on the same OS with a bucketful of Opterons instead succeeded. As I said, there's no solution but a number of ways to get unstuck if you find yourself in the same situation like Mark and me. Best. Andreas -------- Original Message -------- Subject: Re: [phenixbb] reduce hangs Date: Wed, 4 Dec 2013 09:46:20 -0500 From: Jeff Headd <jjheadd(a)lbl.gov> To: Andreas Förster <docandreas(a)gmail.com> If the different machine is a slightly different architecture it's possible that it may run, even under the same OS version. My guess is that it's down to slight floating point differences in the coordinates throughout the course of refinement, but it's so hard to pinpoint. If you have access to other machines it is probably worth trying. <snip> Jeff

1 0

Coot not allowing real space refinement when connected to phenix
by Jobichen Chacko 04 Dec '13

04 Dec '13

Dear All, I recently upgraded to latest version of coot (version0.7.2). When I try to do the Real space refine zone option, coot is giving the following complaint (Refinement setup failure. Failed to find restraints for PROLYSILE). But it works fine when I open coot separately from a terminal and read the pdb and mtz file. I downloaded this coot along with the ccp4 package and it also works fine when I open from ccp4 interface. Any suggestions will be much appreciated. Jobi

3 2

hi.
by Hexi Zhang 03 Dec '13

03 Dec '13

http://satis.egitmeni.org/dxibo/axzjmylquusc.html Hexi Zhang

1 0

reduce hangs
by Andreas Förster 03 Dec '13

03 Dec '13

Dear all, I'm refining a structure with 11 molecules in the asymmetric unit. Recently I took out NCS (data are to 2.2A and differences between molecules of interest). Refinement now stalls halfway through with two reduce jobs running at 100% CPU. The current job has been sitting like this for four days. The last entry in the log file are a few XYZ individual lines. Refinement includes XYZ, TLS, individual B, weight optimization. I launched 8 cycles. Phenix hangs after 4_bss. I've observed this twice in a row. Any suggestions to overcome this, short of going back to using NCS? This is with Phenix dev-1538 on RHEL 6.3. Thanks. Andreas -- Andreas Förster Crystallization and Xray Facility Manager Centre for Structural Biology Imperial College London

4 3

AutoSol/Heavy atom sites/outside the model
by Egor Tchesnokov 03 Dec '13

03 Dec '13

Thank you for the feedback. WinCoot 0.7/Draw/Cell Symmetry brings up "Show Symmetry" window, which only allows to alter the display but not to coincide the heavy atom sites with the model. Am I not looking at the right menu in Coot ? Thank you for your time and help. Egor

1 0

Autobuild combine
by Heather L Condurso 02 Dec '13

02 Dec '13

Hi all, I am still having issues with Autobuild giving errors during the combine step. I've downloaded the latest build and hope that will solve the problem, but in the mean time I'm having an issue trying to combine the models. I've been following the script: #!/bin/csh -f #COMBINE_MODELS SCRIPT if (-d PDS || -d AutoBuild_run_1_) then echo "Please run in a directory without PDS or AutoBuild_run_1_" exit 1 endif echo "Setting up combine models with a dummy run. NOTE: multiple_models_group_number must be correct" phenix.autobuild fobs.mtz multiple_models=true seq_file=seq.dat combine_only=true multiple_models_group_number=2 multiple_models_number=1 > dummy_autobuild.log echo "Copying files to AutoBuild_run_1_/MULTIPLE_MODELS" mkdir AutoBuild_run_1_/MULTIPLE_MODELS cp coords1.pdb AutoBuild_run_1_/MULTIPLE_MODELS/initial_model.pdb_1_1 cp coords2.pdb AutoBuild_run_1_/MULTIPLE_MODELS/initial_model.pdb_1_2 cp map_coeffs_1.mtz AutoBuild_run_1_/MULTIPLE_MODELS/initial_model.mtz_1_1 cp map_coeffs_2.mtz AutoBuild_run_1_/MULTIPLE_MODELS/initial_model.mtz_1_2 ls AutoBuild_run_1_/MULTIPLE_MODELS/ echo "Running autobuild to combine files in AutoBuild_run_1_/MULTIPLE_MODELS" phenix.autobuild combine_only=true seq_file=seq.dat carry_on=true run=1 > autobuild_combine.log # END OF COMBINE_MODELS SCRIPT # Here are the two log files. I used the expl_fobs_phases_freeR_flags.mtz as fobs.mtz and changed each overall_best.pdb and overall_best_refine_map_coeffs.mtz as the model and map to combine. Thanks, Heather

2 2

AutoSol/Heavy atom sites/outside the model
by Egor Tchesnokov 02 Dec '13

02 Dec '13

Hello, Level: beginner 1. Data collected at wavelength of 0.9537 2. Ran Xtriage: anomalous data seem to extend to 2.2 A 3. Ran Phaser-MR to get partial model (1 copy in AU) 4. Ran AutoSol using partial model and anomalous data AutoSol finds heavy atom sites (FOM=0.69; Overall score=75.6+/-3.3). Coot displays model/density without heavy atoms Heavy atoms are displayed outside the model (using zoom out). Question: what did I do wrong? Thank you for your time and help. Egor

2 1

Red patches of electron density
by Singh, Bishal 01 Dec '13

01 Dec '13

Hi, I am solving a crystal structure at 1.2 Å with completeness 97.3% in the highest resolution shell. There are few small patches of red electron densities on the main chain of a flexible loop at a contour level of ~3sigma of difference ED map. Refinement parameters were following: XYZ coordinates, individual B-factor, optimize X-ray/stereochemistry weight, optimize X-ray/ADP weight and individual Anisotropic ADP for amino acids excluding hydrogen. Target weights are as default in phenix.refine. I already have observed two alternative conformation for this region in another crystal structure. Therefore I am thinking that these red patches are because of lower occupancy of these residues. Please suggest me what shall I do. I shall be grateful for your words. With regards, Bishal Singh University of Heidelberg, Germany

2 1

Iterative build omit map error
by Heather L Condurso 27 Nov '13

27 Nov '13

SO just as my omit map run should have finished, it gave an error instead. I'm running 1.8.4-1496 on a Mac with OSX 10.6.8 and 8GB of memory. Here is the end of the log file: Splitting work into 32 jobs and running with 1 processors using sh background=False in /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0 Final job will be run with sh with background=True Starting job 1...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_1.log Starting job 2...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_2.log Starting job 3...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_3.log Starting job 4...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_4.log Starting job 5...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_5.log Starting job 6...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_6.log Starting job 7...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_7.log Starting job 8...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_8.log Starting job 9...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_9.log Starting job 10...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_10.log Starting job 11...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_11.log Starting job 12...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_12.log Starting job 13...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_13.log Starting job 14...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_14.log Starting job 15...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_15.log Starting job 16...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_16.log Starting job 17...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_17.log Starting job 18...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_18.log Starting job 19...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_19.log Starting job 20...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_20.log Starting job 21...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_21.log Starting job 22...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_22.log Starting job 23...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_23.log Starting job 24...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_24.log Starting job 25...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_25.log Starting job 26...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_26.log Starting job 27...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_27.log Starting job 28...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_28.log Starting job 29...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_29.log Starting job 30...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_30.log Starting job 31...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_31.log Starting job 32...Log will be: /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_32.log ******************************************************************************** Sorry, a subprocess has failed... END OF LOG FILE /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_32.log : /TEMP0/Build_composite_refined_3.pdb R/Rfree=0.32/0.39 Maps from refinement will be filled (2FOFCWT PH2FOFCWT) Set refine map coeffs file pdb_file_refine_map_coeffs to refine_map_coeffs_3.mtz Rebuilding side chains from Build_composite_refined_3.pdb Refining model: edit_model.pdb approx_equal eps: 1e-06 approx_equal multiplier: 10000000000.0 0.0197876048373 approx_equal ERROR 0.0197856373852 approx_equal ERROR Refinement strategy: ['individual_sites', 'individual_adp', 'occupancies'] Moving refinement file AutoBuild_run_32_/TEMP0/Build_composite_refined_side_3_001_f_model.mtz to AutoBuild_run_32_/TEMP0/refine_fmodel.mtz Moving refinement file AutoBuild_run_32_/TEMP0/Build_composite_refined_side_3_001.mtz to AutoBuild_run_32_/TEMP0/refine_map_coeffs_3.mtz Model: AutoBuild_run_32_/TEMP0/Build_composite_refined_side_3.pdb R/Rfree=0.33/0.38 Maps from refinement will be filled (2FOFCWT PH2FOFCWT) Set refine map coeffs file pdb_file_refine_map_coeffs to refine_map_coeffs_3.mtz ******************************************************************************** Failed to carry out AutoBuild_parallel_omit: ******************************************************************************** Sorry, a subprocess has failed... END OF LOG FILE /Users/stevebruner/Desktop/APS_102313/MdnJB/AutoBuild_run_22_/TEMP0/RUN_FILE_32.log : /TEMP0/Build_composite_refined_3.pdb R/Rfree=0.32/0.39 Maps from refinement will be filled (2FOFCWT PH2FOFCWT) Set refine map coeffs file pdb_file_refine_map_coeffs to refine_map_coeffs_3.mtz Rebuilding side chains from Build_composite_refined_3.pdb Refining model: edit_model.pdb approx_equal eps: 1e-06 approx_equal multiplier: 10000000000.0 0.0197876048373 approx_equal ERROR 0.0197856373852 approx_equal ERROR Refinement strategy: ['individual_sites', 'individual_adp', 'occupancies'] Moving refinement file AutoBuild_run_32_/TEMP0/Build_composite_refined_side_3_001_f_model.mtz to AutoBuild_run_32_/TEMP0/refine_fmodel.mtz Moving refinement file AutoBuild_run_32_/TEMP0/Build_composite_refined_side_3_001.mtz to AutoBuild_run_32_/TEMP0/refine_map_coeffs_3.mtz Model: AutoBuild_run_32_/TEMP0/Build_composite_refined_side_3.pdb R/Rfree=0.33/0.38 Maps from refinement will be filled (2FOFCWT PH2FOFCWT) Set refine map coeffs file pdb_file_refine_map_coeffs to refine_map_coeffs_3.mtz ******************************************************************************** Can this be salvaged? It took 5 days to run. I found a section about this in the documentation, but I can't find the maps for each region overall_best_denmod_map_coeffs.mtz_OMIT_REGION_1 etc. I can only find working_best_denmod_map_coeffs.mtz_OMIT_REGION_1 etc. Will these maps work? Can you tell what the error means so that I can avoid this problem in the future? Thanks, Heather

2 1

zero B-factors
by Edna Patricia Matta, Ms 25 Nov '13

25 Nov '13

Hello, Recently I did submit a structure to the PDB but they found that some atoms had zero value for the B-factors. I am new in this filed and I do not know what to do. There have been reports of problems with some versions of phenix about this. I used Phenix 1.8. I repeted the refinement doing restrain NCS-related B-factors but I do not find where to check for the the B-factor values in the the log file. Somebody can help me? Thanks a lot, Edna

7 15

Re: [phenixbb] [ccp4bb] Ligandfit - problem with ligand_start
by Terwilliger, Thomas C 25 Nov '13

25 Nov '13

Hi Danilo, The most likely problem is that the ligand_start file has to have enough atoms in it to define a fixed group in the ligand. The way that ligandfit works is that it breaks a ligand up into rigid groups connected by rotatable bonds. The ligand_start file has to contain one complete rigid group. So probably if you supply a little more of the ligand then it might run. Just to test this, can you try this test of using a ligand_start file: phenix_regression.wizards.test_resolve resolve test_fit_start And then if that gives an "OK", try this: cd test_fit_start phenix.ligandfit ligand=side.pdb ligand_start=start.pdb resolve_map.mtz quick=true Does that run? If you the look at LigandFit_run_1_/TEMP0/resolve_ligand_fit.log Does that say >ligand_start start.pdb Starting configuration of ligand will be read from: start.pdb If all that works, then try adding more atoms to your ligand_start file. Let me know if that doesn't do it! All the best, Tom T Dear all, I am working on a membrane protein covalently bound to a molecular antanna: it is known that this molecule binds to lysine residue but I do not know how many and which lysine residues it binds. 20 diffraction datasets of this protein-ligand complex have been obtained and now, I would quickly localize the ligand using the Fo-Fc map of each data set and using the information on the covalent bound protein-ligand. Ligandfit tool (PHENIX) seems to be indicated to do this; to use the information on the covalent bound, I am using the ligand_start keyword with a pdb containing a ghost atom (however present in ligand model) perfectly superposed to the lysine atom that should bind the ligand. The command used is: phenix.ligandfit data=prot.mtz model=prot.pdb ligand=lig.pdb ligand_start=lig_start.pdb input_labels="FOFCWT PHFOFCWT" \ refine_ligand=True \ nproc=32 \ cif_def_file_list=lig.cif description: - prot.mtz (data) - prot.pdb (protein without ligand) - lig.pdb (ligand containing ghost atom) - lig_start.pdb (ghost atom superpose to NZ of a lysine) - lig.cif (restrain of lig.pdb) Strangely, no ligand is found at the end of the process even reducing ligand_cc_min to 0.01. I have run the same command by using an other protein where an other ligand has been correctly fitted but, also in this case, no ligand has been detected. Conversely, without the use of ligand_start, ligandfit properly localizes the ligand. I'm doing some mistake in the use of ligand_start? Do you know an other tool to perform a ligand fitting in these conditions? Thanks for your answers. Danilo

1 0

Ligandfit - problem with ligand_start
by Danilo Belviso 25 Nov '13

25 Nov '13

Dear all, I am working on a membrane protein bound to a antenna molecule: I know that the molecule binds to lysine residue but I do not know how many and which lysine residues it binds. 20 diffraction datasets of protein-ligand complex have been obtained and now, I would quickly localize the ligand using the Fo-Fc map of each data set and using the information on the covalent bound protein-ligand. Ligandfit tool (PHENIX) seems to be indicated to do this; to use the information on the covalent bound, I have thought to use the ligand_start keyword with a pdb containing a ghost atom (however present in ligand model) perfectly superposed to the lysine atom that should bind the ligand. The command used is: phenix.ligandfit data=prot.mtz model=prot.pdb ligand=lig.pdb ligand_start=lig_start.pdb input_labels="FOFCWT PHFOFCWT" \ refine_ligand=True \ nproc=32 \ cif_def_file_list=lig.cif description: - prot.mtz (data) - prot.pdb (protein without ligand) - lig.pdb (ligand containing ghost atom) - lig_start.pdb (ghost atom superpose to NZ of a lysine) - lig.cif (restrain of lig.pdb) At the end of the processing, ligand is not found. So, I have tryed with ligand_cc_min=0.01, but also in this case ligandfit does not give result. I have thought that no ligand is present. However, to be sure that all was correct, I have run the same command by using an other protein where an other ligand has been correctly fitted. Also in this case, no result has been obtained. Conversely, without the use of ligand_start, ligandfit properly localizes the ligand. I'm doing some mistake in the use of ligand_start? Thanks for your answers. Danilo

1 0

using reference model with polyA chain.
by Francis Reyes 19 Nov '13

19 Nov '13

I'm refining a structure ( at low res ), so I trimmed off all the side chain atoms from my model while my reference model still has these chains. I'm receiving errors like... ----------Correcting rotamer outliers to match reference model---------- ** evaluating rotamers for working model ** A 502 ILE is missing some sidechain atoms, **skipping** A 504 LYS is missing some sidechain atoms, **skipping** A 505 VAL is missing some sidechain atoms, **skipping** I don't care about the side chain atoms, but I'd still like the reference model restraints to provide secondary structure. Does "**skipping**" mean that the reference model is no longer being considered for the particular residue? Or will I still retain my secondary structure from the reference model? Thanks, F --------------------------------------------- Francis E. Reyes PhD 215 UCB University of Colorado at Boulder

2 1

completeness problem in phenix.refine
by Yarrow Madrona 19 Nov '13

19 Nov '13

Hi, I recently merged two data sets, one having a phase and the other Amplitudes. For some reason my completeness is always 50%. I had anomalous data but the completeness is around 100% after processing and when I do an mtz dump. I'm confused why I get a completeness of 50%. Thank you for any help you can provide. -Yarrow -- Yarrow Madrona Graduate Student Molecular Biology and Biochemistry Dept. University of California, Irvine Natural Sciences I, Rm 2403 Irvine, CA 92697

2 2

Error: "Error running main function (INPUT): LLG completion peak selection method not recognized"
by Egor Tchesnokov 18 Nov '13

18 Nov '13

Hello, I am trying to run Phaser/SAD or Phaser/MR_SAD (GUI using defaults), but I get the following error: "Error running main function (INPUT): LLG completion peak selection method not recognized". Thank you for your help Egor

2 1

MSE out density after phenix.refine
by Pavel Afonine 15 Nov '13

15 Nov '13

Hello, anyone experiencing problems related to MSE residues being pushed out density after phenix.refine run should update to next Phenix nightly build, which should fix this problem. Please let me know otherwise off-list. Pavel

1 0

MR_Rosetta prerefine question
by Heather L Condurso 12 Nov '13

12 Nov '13

Dear all, I am trying to use the rosetta prerefinement fuction of MR_rosetta to complete a homology model with missing stretches of amino acids. Everything runs fine and residues are being added to fill the gaps using the fragment files, but rosetta is also altering parts of the model which should already be well placed. Is there a way to tell Rosetta only to prerefine certain residues? I've been following the command line suggestion: phenix.mr_rosetta \ seq_file=seq.dat \ data=coords1.mtz \ search_models=coords1.pdb \ run_prerefine=True \ number_of_prerefine_models=1000 \ fragment_files = test3.gz \ fragment_files = test9.gz \ rescore_mr.relax=False \ rosetta_models=100 \ ncs_copies=2 \ space_group=p212121 \ use_all_plausible_sg=False \ nproc=200 \ group_run_command=qsub Thanks, Heather Condurso condurso(a)bc.edu

4 5

Re: [phenixbb] Phenix.Rosetta_refine on Ubuntu 12.04 LTS
by qwertgfdsa78 11 Nov '13

11 Nov '13

Hi All, I tried Jan's tutorial and finally I was able to install and ran phenix.rosetta with test set 1yjp successfully. However, when I tried my data, I always got *"terminate called after throwing an instance of 'std::bad_alloc'"* Google suggests maybe memory is not enough. I used VirutalBox and assign 8gb memory to it. My protein has ~3000 residues. I guess this will use more memory but how much RAM I need? Or the error is due to some other problem? Thank for any inputs. Lee > ------------------------------ > > Message: 3 > Date: Tue, 05 Nov 2013 14:26:00 +0100 > From: Jan Gebauer <jan.gebauer(a)uni-koeln.de> > To: PHENIX user mailing list <phenixbb(a)phenix-online.org> > Subject: Re: [phenixbb] Phenix.Rosetta_refine on Ubuntu 12.04 LTS > Message-ID: <5278F1E8.9020602(a)uni-koeln.de> > Content-Type: text/plain; charset="utf-8" > > Hi David, hi all, > > I now managed to at least get a test refinement working. What I mostly > changed was the Ubuntu version. > I have no idea why, but I managed to get it to work on Ubuntu 10.04. > > In principle I installed 10.04 freshly in a virtualbox and from their it > worked rather straight-forward. > > Here is how I proceeded in detail for those of you who are interested: > > 1) Setting up a Virtualbox 4.3.2 r90405 on Ubuntu 12.04.3. > 2) Installing the 64bit version of Ubuntu 10.04 in a virtualbox (4 > processors / 90% speed cap / 4GB Ram /VT-x enabled) > 3) Install VBoxAdditions > 4) Update ubuntu (sudo apt-get update) > 5) Installing g++ and zlib library (sudo apt-get install g++ zlib1g-dev) > 6) Installing binary version of phenix-dev-1525 to > /usr/local/science/phenix-dev-1525 (sudo ./install > --prefix=/usr/local/science) > 7) Unpack rosetta_2013_wk42_bundle to > /usr/local/science/rosetta_2013_wk42_bundle > 8) set rossetta variables by inserting them at the start of > /usr/local/science/pheninx-dev-1525/phenix_env.sh (just in line 3ff): > export PHENIX_ROSETTA_PATH=/usr/local/science/rosetta_2013wk42_bundle > export ROSETTA3_DB=/usr/local/science/rosetta_2013wk42_bundle/main/database > export > PATH=$PATH:/usr/local/science/rosetta_2013wk42_bundle/main/source/bin > 9) assign ownership to active user (sudo chown gebauer:gebauer > /usr/local/science/ -R) > 10) Compile it (phenix.build_phenix_interface nproc=4) > It should run without problems > > Testing: > 1) Create a new directory (mkdir test && cd test) > 2) Download testset (suggested by Nat: phenix.fetch_pdb --mtz 1yjp) > 3) Run rosetta_refine (phenix.rosetta_refine 1yjp.pdb 1yjp.mtz nproc=4) > In my case rFree dropped from 18.21 to 14.90. It's not a low resolution > dataset, but it show that it general works and lowRes would potentially > take to long for a mere testing > > My real data are still being processed so I have no idea how well it > will works... only thing I found out already: > 4.5 GB RAM for 4 processes calculating ~500Aa via rosetta seems not to > be enough. I thereby reduced my "real run" to only one processor... > > Regards, > Jan > > >

3 3

structure factor labels in phenix.refine
by Tanner, John J. 08 Nov '13

08 Nov '13

I have two questions about which structure factor amplitudes are refined against in phenix.refine. 1. Which structure factor amplitude is used in refinement when the labels in the .def file are the following? labels = "F-obs(+),SIGF-obs(+),F-obs(-),SIGF-obs(-)" 2. Is anomalous data used in the refinement when the labels are as follows? labels = "F,SIGF,DANO,SIGDANO,ISYM" I ask these questions because I am wondering which structure factor amplitudes to deposit, and I am trying to rationalize differences between R-factors from refinement and Generate Table 1. Thanks. Jack Tanner John J. Tanner Professor of Biochemistry and Chemistry University of Missouri-Columbia 125 Chemistry Building Columbia, MO 65211 Phone: 573-884-1280 Fax: 573-882-2754 Email: tannerjj(a)missouri.edu<mailto:[email protected]> http://faculty.missouri.edu/~tannerjj/tannergroup/tanner.html

4 4

few questions about resolving new structure through MR
by Zhihong Yu 07 Nov '13

07 Nov '13

Hi, all I'm a rookie in resolving a brand new structure. I have some questions for my current case and look forward to some suggestions. Now I’m working on a protein like this: N-ter(55aa)—domainA(110aa)—linker(30aa)—domainB(150aa)—C-ter(20aa), I got a diffraction data just to 3.5Å, and there is no complete homology structure in pdb bank, but only a homology structure (named as structureX later) for domainB with ~30% sequence identity, so I have some questions as following: 1. Is it possible to find a resolution through MR approach using structureX as a search model? Especially considering that the resolution is only 3.5Å. Currently I just tried once using phaser and refine the structure for 3 cycyles, I can get a R/Rfree of 0.45/0.55, and it looks like most of backbone in the structureX, especially those within helix or sheet, can be well described by 2Fo-Fc density. Is this primary result promising or not? 2. If it’s possible, what’s the general optimal procedure I should follow? Really thanks for any advice and suggestions! Zhihong

2 3

Issue with SeMet residues
by Mark Wilson 07 Nov '13

07 Nov '13

Dear PHENIX users and developers, I am refining a structure containing SeMet residues at atomic (1.05 Å) resolution. The model was initially refined in Refmac5 and behaved as expected (i.e. reasonable R factors and agreement with maps). I've used this final Refmac5 model as a starting point for refinement in PHENIX 1.8.4 with standard settings, as below. phenix.refine data.mtz \ refinement.input.xray_data.labels="F, SIGF" model.pdb \ main.use_experimental_phases=false \ refinement.main.number_of_macro_cycles=5 output.prefix="model_refine"> phenix_refine_1.log Strangely, the conformations of the MSE residues are changed quite significantly after refinement in PHENIX, departing both from electron density and from the prior model. Any help would be appreciated. Best regards, Mark Mark A. Wilson Associate Professor Department of Biochemistry/Redox Biology Center University of Nebraska N118 Beadle Center 1901 Vine Street Lincoln, NE 68588 (402) 472-3626 mwilson13(a)unl.edu

3 3

Issues using reference model with Autobuild
by Dan McNamara 07 Nov '13

07 Nov '13

Hello all, I'm hoping to troubleshoot a problem with autobuild and the use of a reference model. Autobuild has been always failing once it gets to the initial refinement step and attempts to add reference model restraints. I've included my autobuild parameterization. Refinement is where it dies and the contents of LAST.LOG and my error log file are below too. I'm not sure how to get the reference model file location properly interpreted. I even tried manually moving the reference PDB into the AutoBuild_run_X_/ and TEMP0/ folders as they're generated to no avail. Best, Dan phenix.autobuild \ data=/home/.../xyz.mtz \ refinement_file=/home/.../xyz.mtz \ input_refinement_labels="FP SIGFP FreeRflag" \ model=/home/.../xyz.pdb \ seq_file=/home/.../xyz.dat \ s_annealing=False \ rebuild_in_place=True \ n_cycle_rebuild_max=3 \ refine_eff_file=/home/.../xyz.eff \ refine_with_ncs=True \ model_building.reference_model=/home/.../def.pdb \ #I have also tried locations like ../def.pdb and def.pdb that give Errno 21 or Errno 2 below nproc=200 \ general.run_command="/home/.../xyz.sh" AutoBuild_run_1_/LAST.LOG: ================== Extract refinement strategy and selections ================= Refinement flags and selection counts: individual_sites = True (30132 atoms) torsion_angles = False (0 atoms) rigid_body = True (30120 atoms in 24 groups) individual_adp = False (iso = 0 aniso = 0) group_adp = True (30132 atoms in 4092 groups) tls = True (30120 atoms in 147 groups) occupancies = True (12 atoms) group_anomalous = False n_use = 30132 n_use_u_both = 30120 n_use_u_iso = 12 n_use_u_aniso = 0 n_grad_site = 0 n_grad_u_iso = 0 n_grad_u_aniso = 0 n_grad_occupancy = 0 n_grad_fp = 0 n_grad_fdp = 0 total number of scatterers = 30132 *** Adding Reference Model Restraints *** errors.txt: Sorry: ******************************************************************************** Failed to carry out AutoBuild_build_cycle: Initial phenix refinement of AutoBuild_run_1_/TEMP0/NO_LIG_NO_WAT_xyz_ed.pdb with aniso_refinement_data_PHX.mtz failed...quitting... /home/.../AutoBuild_run_1_/LAST.LOG may have additional information... Cannot open file for reading: "/" [Errno 21] Is a directory: '/' ******************************************************************************** Sorry: ******************************************************************************** Failed to carry out AutoBuild_build_cycle: Initial phenix refinement of AutoBuild_run_2_/TEMP0/NO_LIG_NO_WAT_xyz_ed.pdb with aniso_refinement_data_PHX.mtz failed...quitting... /home/.../AutoBuild_run_2_/LAST.LOG may have additional information... Cannot open file for reading: "." [Errno 21] Is a directory: '.' ******************************************************************************** Sorry: ******************************************************************************** Failed to carry out AutoBuild_build_cycle: Initial phenix refinement of AutoBuild_run_3_/TEMP0/NO_LIG_NO_WAT_xyz_ed.pdb with aniso_refinement_data_PHX.mtz failed...quitting... /home/.../AutoBuild_run_3_/LAST.LOG may have additional information... Cannot open file for reading: "d" [Errno 2] No such file or directory: 'd' ********************************************************************************

2 2

Sigmas of F+ and F- equal in _data file
by Tanner, John J. 07 Nov '13

07 Nov '13

I have a question about the _data.mtz file that phenix.refine outputs (version 1.8-1069). It appears that sigmaF+ equals sigmaF- for each reflection in the _data.mtz file, although they were unequal in the input .mtz file. You can compare the input mtz and the output _data.mtz file by viewing these links: http://faculty.missouri.edu/~tannerjj/input.jpeg http://faculty.missouri.edu/~tannerjj/samesigma.jpeg Thanks. Jack Tanner John J. Tanner Professor of Biochemistry and Chemistry University of Missouri-Columbia 125 Chemistry Building Columbia, MO 65211 Phone: 573-884-1280 Fax: 573-882-2754 Email: tannerjj(a)missouri.edu<mailto:[email protected]> http://faculty.missouri.edu/~tannerjj/tannergroup/tanner.html

2 1

real-space refinement
by "Weiergräber, Oliver H." 06 Nov '13

06 Nov '13

Hello, beginning with phenix 1.8.3, the individual_sites_real_space option will not run globally if, e.g., R < 0.25. Is there a way to force the original behaviour? Regards, Oliver ================================================ PD Dr. Oliver H. Weiergräber Institute of Complex Systems ICS-6: Structural Biochemistry Tel.: +49 2461 61-2028 Fax: +49 2461 61-9540 ================================================ ------------------------------------------------------------------------------------------------ ------------------------------------------------------------------------------------------------ Forschungszentrum Juelich GmbH 52425 Juelich Sitz der Gesellschaft: Juelich Eingetragen im Handelsregister des Amtsgerichts Dueren Nr. HR B 3498 Vorsitzender des Aufsichtsrats: MinDir Dr. Karl Eugen Huthmacher Geschaeftsfuehrung: Prof. Dr. Achim Bachem (Vorsitzender), Karsten Beneke (stellv. Vorsitzender), Prof. Dr.-Ing. Harald Bolt, Prof. Dr. Sebastian M. Schmidt ------------------------------------------------------------------------------------------------ ------------------------------------------------------------------------------------------------

2 3

Is the protein-ligand complex structure model right??
by Wei Shi 06 Nov '13

06 Nov '13

Hi all, I am trying to solve a protein-ligand complex structure. The resolution of the data is 2.8Å. The molecular replacement solution is found in space group P21 21 21. The current best model: start final ------------------------------------------------ R-work: 0.2836 0.2484 R-free: 0.3236 0.3066 RMS(angles): 0.60 0.63 RMS(bonds): 0.002 0.002 ==============MolProbity validation ============== Ramachandran outliers: 2.4% (Goal: < 0.2%) Ramachandran favored: 89.2% (Goal: > 98%) Rotamer outliers: 4.7% (Goal: 1%) C-beta outliers: 0 (Goal: 0) Clashscore: 10.89 Overall score: 2.64 Four ligands are placed for every dimer of the protein. And we are not sure whether we put the ligands right, since first it's thought that only two ligands bind every dimer of the protein. I created a SA omit map (composite). The density for the most of the N-terminal DNA binding domain is very bad, and the density for C-terminal ligand binding domain is much better. And, I could see density for all 4 ligands. Does this necessarily mean that the model for the N-terminal DNA binding domain is not right or should I try other types of omit map (simple, refined or iterative building)? Also, I am wondering what's the difference between different types of omit map (simple, refined, iterative building) and whether it's the most reliable method to tell whether the structure model is right. So, I am not sure whether the structure model I have is correct. It's suspected that the molecular solution is in the wrong space group or the data is twinned and only two ligands bind every dimer of the protein......But, I am not sure what I should check or what to do to improve the model. Thank you so much! Best, Wei

1 0

Is the protein-ligand complex structure model right??
by Wei Shi 05 Nov '13

05 Nov '13

Hi all, I am trying to solve a protein-ligand complex structure. The resolution of the data is 2.8Å. The molecular replacement solution is found in space group P21 21 21. The current best model: start final ------------------------------------------------ R-work: 0.2836 0.2484 R-free: 0.3236 0.3066 RMS(angles): 0.60 0.63 RMS(bonds): 0.002 0.002 ==============MolProbity validation ============== Ramachandran outliers: 2.4% (Goal: < 0.2%) Ramachandran favored: 89.2% (Goal: > 98%) Rotamer outliers: 4.7% (Goal: 1%) C-beta outliers: 0 (Goal: 0) Clashscore: 10.89 Overall score: 2.64 Four ligands are placed for every dimer of the protein. And we are not sure whether we put the ligands right, since first it's thought that only two ligands bind every dimer of the protein. I created a SA omit map (composite). The density for the most of the N-terminal DNA binding domain is very bad, and the density for C-terminal ligand binding domain is much better. And, I could see density for all 4 ligands. Does this necessarily mean that the model for the N-terminal DNA binding domain is not right or should I try other types of omit map (simple, refined or iterative building)? Also, I am wondering what's the difference between different types of omit map (simple, refined, iterative building) and whether it's the most reliable method to tell whether the structure model is right. So, I am not sure whether the structure model I have is correct. It's suspected that the molecular solution is in the wrong space group or the data is twinned and only two ligands bind every dimer of the protein......But, I am not sure what I should check or what to do to improve the model. Thank you so much! Best, Wei

1 0

Green density for ligands in places where the ligand is in density
by Wei Shi 05 Nov '13

05 Nov '13

Hi guys, I am working with a protein-ligand complex structure. I generated the ligand pdb and cif with Phenix.eLBOW. Then, I fit the ligand to density (green or continuous blue denisty) in Coot. Then, I did a refinement of the protein-ligand complex model in Phenix. I opened the refinement results in Coot, I don't know why in some places where the ligand is in density, it always shows green density in Coot... And, from the pdb file, the occupancy for all ligands are 0? Do any of you happens to know whether this is unusual and what might go wrong? Thank you so much! Best, Wei

3 4

Re: [phenixbb] Restraints generation for TaBr
by Rebecca Oot 05 Nov '13

05 Nov '13

Hi Nigel, Thanks so much for the response, the information and for the files. Best, Rebecca

1 0