- phenixbb - phenix-online.org

Re: [phenixbb] [ccp4bb] Lysozyme/Xenon derivatives & Phenix.
by Randy Read 03 Feb '12

03 Feb '12

Hi, I'm cross-posting this to the Phenix-BB, because you did this calculation in Phenix, but since the site completion for SAD phasing would be done with Phaser, the answer is relevant to both packages. Basically, I think Jim Pflugrath's suggestion is correct. Comparing your Xe coordinates to S and Cl coordinates from one of our test SAD phasing datasets, the majority of your extra sites are in the positions of either intrinsic S anomalous scatterers or bound halides. The Phaser completion step found them, even though you only asked for 2 sites initially, because it continues adding new sites as long as there are significant peaks in the SAD LLG maps. The refined occupancies make sense, assuming you were collecting your data around 1A wavelength. At that wavelength f" is about 4-5e for Xe and about 0.25-0.3 for S or Cl, and multiplying (say) 4.5e by occupancies of 0.03-0.08 gives effective f" values of 0.14-0.36e. Of course, Phaser isn't just trying to fit the imaginary part of the scattering, so the occupancies could end up slightly higher because it's trying to account for atoms with 16-17 electrons using partially-occupied atoms with 54 electrons. Knowing that, once Phaser has found the Xe atoms, it will then find the weak anomalous scattering from S and Cl, you could tell Phaser to look in addition for S atoms. (At around 1A wavelength, there's not enough difference between S and Cl, in either their real or anomalous scattering, to matter, so you don't need to look for Cl in addition.) Phaser will reassign atom types based on refined occupancy, so it would likely get most of these right, which would very slightly improve your phasing model (because you would then have the right ratio of real to imaginary scattering). Another advantage to looking for both atom types is that, when you mix atoms with different ratios of real and imaginary scattering, Friedel's law is broken even for the substructure structure factors, which has two consequences. First, the likelihood score will probably be significantly better for the correct enantiomorph than the incorrect one; second, you're likely to end up with a slightly more complete substructure when working with the correct enantiomorph. Note that AutoSol also chose the correct enantiomorph after phasing, by looking at map quality indicators for the two hands. Regards, Randy Read On 2 Feb 2012, at 16:43, Brennan Bonnet wrote: > Hi all, > > I have a strange result using Phenix's AutoSol to look for xenon sites in lysozyme. > > For a few months now I have been trying to produce xenon derivatives of lysozyme using pressures in the range of 50-350psi and time ranges between 5-60min trying to find the "sweet spot" so that this technique can be applied to other proteins. > > For each dataset, I AutoSol using phenix and have it look for 2, 3, and 4 xenon sites. I haven't had much success though and typical values when asked to look for 2 sites have been Rwork/Rfree = 0.5585/0.5892, CC=0.14, Bayes-CC=6.5 and with xenon sites: > > (a, b, c, alpha, beta, gamma, space group) > CRYST1 78.870 78.870 36.940 90.00 90.00 90.00 P 41 21 2 > > (columns after numbering are a, b, c, occ, B factor) > HETATM 1 XE XE Z 1 48.601 67.475 1.088 0.06 6.58 XE > HETATM 2 XE XE Z 2 54.986 76.636 2.746 0.11 7.49 XE > > which, to me, says that there is no binding (or at least nothing obvious). > > I finally got a good result though with 350psi and 5 minutes pressurization. After AutoSol, FOM=0.436, Bayes-CC=33.57, Model CC=0.83, Rwork/Rfree=0.2843/0.3023, 117/125 residues build and placed. And even though I told it to only look for 2 sites, Phenix went ahead and found 22 sites: > > (a, b, c, alpha, beta, gamma, space group) > CRYST1 79.130 79.130 36.920 90.00 90.00 90.00 P 43 21 2 > > (colums after numbering are a, b, c, occ, B factor) > HETATM 1 XE XE Z 1 -36.117 -56.764 -1.825 0.18 9.43 XE > HETATM 2 XE XE Z 2 -54.805 -54.805 0.000 0.14 6.44 XE > HETATM 3 XE XE Z 3 -8.337 -31.695 -16.766 0.08 8.26 XE > HETATM 4 XE XE Z 4 -5.702 -64.959 -35.134 0.05 5.21 XE > HETATM 5 XE XE Z 5 -9.535 -22.397 -31.839 0.06 14.33 XE > HETATM 6 XE XE Z 6 -7.319 -66.343 -35.123 0.05 6.02 XE > HETATM 7 XE XE Z 7 -8.580 -21.153 -30.274 0.05 7.43 XE > HETATM 8 XE XE Z 8 -1.494 -29.583 -2.339 0.06 7.36 XE > HETATM 9 XE XE Z 9 -0.518 -16.483 -16.776 0.05 6.56 XE > HETATM 10 XE XE Z 10 -0.464 -29.847 -4.386 0.06 6.04 XE > HETATM 11 XE XE Z 11 -1.222 -16.234 -18.952 0.04 6.51 XE > HETATM 12 XE XE Z 12 -6.001 -20.380 -4.532 0.05 7.22 XE > HETATM 13 XE XE Z 13 -5.169 -27.522 -10.128 0.05 3.54 XE > HETATM 14 XE XE Z 14 -11.931 -67.039 -4.699 0.05 7.96 XE > HETATM 15 XE XE Z 15 -6.538 -10.408 -3.761 0.04 10.22 XE > HETATM 16 XE XE Z 16 -0.458 -23.897 -31.221 0.05 6.85 XE > HETATM 17 XE XE Z 17 -3.308 -65.181 -8.754 0.03 8.08 XE > HETATM 18 XE XE Z 18 -11.080 -29.129 -2.144 0.04 5.61 XE > HETATM 19 XE XE Z 19 -5.720 -73.029 -2.035 0.04 11.24 XE > HETATM 20 XE XE Z 20 10.236 60.097 33.822 0.03 5.38 XE > HETATM 21 XE XE Z 21 5.980 68.581 3.532 0.03 7.49 XE > HETATM 22 XE XE Z 22 11.107 13.277 27.334 0.03 7.10 XE > > I'm pretty sure there aren't 22 sites for xenon to bind in lysozyme but it looks to me like the top 2 (maybe top 3) are actual sites which is great. > The strange part is that when I use the same data and have Phenix look for 3 or 4 sites it then only finds 3 or 4 very low occupancy sites and the CC's/R's are again all very poor. > Has anyone had this problem or does anyone know what's going on? This is the first success that I'm having with xenon derivatization and it seems to me that if Phenix can find 22 sites when asked to find 2 it should have no problem when asked to find 3 or 4. > > Thanks in advance, > ~Brennan~ ------ Randy J. Read Department of Haematology, University of Cambridge Cambridge Institute for Medical Research Tel: + 44 1223 336500 Wellcome Trust/MRC Building Fax: + 44 1223 336827 Hills Road E-mail: rjr27(a)cam.ac.uk Cambridge CB2 0XY, U.K. www-structmed.cimr.cam.ac.uk

1 0

Re: [phenixbb] Lysozyme/Xenon derivatives & Phenix.
by Terwilliger, Thomas C 03 Feb '12

03 Feb '12

Hi Brennan, The number of sites that you ask for affects the HYSS heavy-atom search procedure, but not too much else. Once some sites are found by HYSS, then autosol uses phaser to identify additional sites (that's why you got 22 sites, not 2). I am a little surprised that changing the number of sites from 2 to 4 makes the difference between finding correct sites and not finding them at all, but it does seem possible. For your experiment, one thing you could do: once you figure out where the "standard" Xe sites are, just feed them in to autosol with a sites_file=xxx.pdb for all your datasets. Then they are all starting out with pretty much the same information. All the best, Tom T ________________________________________ From: CCP4 bulletin board [CCP4BB(a)JISCMAIL.AC.UK] on behalf of Brennan Bonnet [Brennan.Bonnet(a)LIGHTSOURCE.CA] Sent: Thursday, February 02, 2012 9:43 AM To: CCP4BB(a)JISCMAIL.AC.UK Subject: [ccp4bb] Lysozyme/Xenon derivatives & Phenix. Hi all, I have a strange result using Phenix's AutoSol to look for xenon sites in lysozyme. For a few months now I have been trying to produce xenon derivatives of lysozyme using pressures in the range of 50-350psi and time ranges between 5-60min trying to find the "sweet spot" so that this technique can be applied to other proteins. For each dataset, I AutoSol using phenix and have it look for 2, 3, and 4 xenon sites. I haven't had much success though and typical values when asked to look for 2 sites have been Rwork/Rfree = 0.5585/0.5892, CC=0.14, Bayes-CC=6.5 and with xenon sites: (a, b, c, alpha, beta, gamma, space group) CRYST1 78.870 78.870 36.940 90.00 90.00 90.00 P 41 21 2 (columns after numbering are a, b, c, occ, B factor) HETATM 1 XE XE Z 1 48.601 67.475 1.088 0.06 6.58 XE HETATM 2 XE XE Z 2 54.986 76.636 2.746 0.11 7.49 XE which, to me, says that there is no binding (or at least nothing obvious). I finally got a good result though with 350psi and 5 minutes pressurization. After AutoSol, FOM=0.436, Bayes-CC=33.57, Model CC=0.83, Rwork/Rfree=0.2843/0.3023, 117/125 residues build and placed. And even though I told it to only look for 2 sites, Phenix went ahead and found 22 sites: (a, b, c, alpha, beta, gamma, space group) CRYST1 79.130 79.130 36.920 90.00 90.00 90.00 P 43 21 2 (colums after numbering are a, b, c, occ, B factor) HETATM 1 XE XE Z 1 -36.117 -56.764 -1.825 0.18 9.43 XE HETATM 2 XE XE Z 2 -54.805 -54.805 0.000 0.14 6.44 XE HETATM 3 XE XE Z 3 -8.337 -31.695 -16.766 0.08 8.26 XE HETATM 4 XE XE Z 4 -5.702 -64.959 -35.134 0.05 5.21 XE HETATM 5 XE XE Z 5 -9.535 -22.397 -31.839 0.06 14.33 XE HETATM 6 XE XE Z 6 -7.319 -66.343 -35.123 0.05 6.02 XE HETATM 7 XE XE Z 7 -8.580 -21.153 -30.274 0.05 7.43 XE HETATM 8 XE XE Z 8 -1.494 -29.583 -2.339 0.06 7.36 XE HETATM 9 XE XE Z 9 -0.518 -16.483 -16.776 0.05 6.56 XE HETATM 10 XE XE Z 10 -0.464 -29.847 -4.386 0.06 6.04 XE HETATM 11 XE XE Z 11 -1.222 -16.234 -18.952 0.04 6.51 XE HETATM 12 XE XE Z 12 -6.001 -20.380 -4.532 0.05 7.22 XE HETATM 13 XE XE Z 13 -5.169 -27.522 -10.128 0.05 3.54 XE HETATM 14 XE XE Z 14 -11.931 -67.039 -4.699 0.05 7.96 XE HETATM 15 XE XE Z 15 -6.538 -10.408 -3.761 0.04 10.22 XE HETATM 16 XE XE Z 16 -0.458 -23.897 -31.221 0.05 6.85 XE HETATM 17 XE XE Z 17 -3.308 -65.181 -8.754 0.03 8.08 XE HETATM 18 XE XE Z 18 -11.080 -29.129 -2.144 0.04 5.61 XE HETATM 19 XE XE Z 19 -5.720 -73.029 -2.035 0.04 11.24 XE HETATM 20 XE XE Z 20 10.236 60.097 33.822 0.03 5.38 XE HETATM 21 XE XE Z 21 5.980 68.581 3.532 0.03 7.49 XE HETATM 22 XE XE Z 22 11.107 13.277 27.334 0.03 7.10 XE I'm pretty sure there aren't 22 sites for xenon to bind in lysozyme but it looks to me like the top 2 (maybe top 3) are actual sites which is great. The strange part is that when I use the same data and have Phenix look for 3 or 4 sites it then only finds 3 or 4 very low occupancy sites and the CC's/R's are again all very poor. Has anyone had this problem or does anyone know what's going on? This is the first success that I'm having with xenon derivatization and it seems to me that if Phenix can find 22 sites when asked to find 2 it should have no problem when asked to find 3 or 4. Thanks in advance, ~Brennan~

2 1

Picking Rfree in thin resolution shells using command line
by Simon Kolstoe 01 Feb '12

01 Feb '12

Dear phenixbb, I see from a quick google that it is possible to pick my Rfree's using thin resolution shells (coz I've got 20 fold NCS), however as I am someone who tries to avoid the GUI where at all possible, could someone let me know what the command line way of doing this is? Thanks, Simon --------------------------------------------------------------- Dr Simon Kolstoe Laboratory for Protein Crystallography Wolfson Drug Discovery Unit University College London Rowland Hill Street, London NW3 2PF Tel: 020 7433 2765 http://www.ucl.ac.uk/~rmhasek ---------------------------------------------------------------

9 14

labelit.index python error
by Christian Roth 01 Feb '12

01 Feb '12

Hi, I had a few images, some of them are very bad. So itried to index 1 of the pictures, and indexing went very well, but then a strange error appeared which I cannot relate to anything. What might be the problem here? Best Regards Christian labelit.index_dev-968 /hosts/vulcan1/mar345/moschuetz/064c2_4/064c2_4_001.mar2300 Beam center is not immediately clear; rigorously retesting 2 solutions Beam x 172.5 y 172.9, initial score 405; refined rmsd: 0.1445 <Z>=20.4 Beam x 171.8 y 175.6, initial score 335; refined rmsd: 0.2160 <Z>=14.9 LABELIT Indexing results: Beam center x 172.58mm, y 172.91mm, distance 209.93mm ; 80% mosaicity=0.25 deg. Solution Metric fit rmsd #spots crystal_system unit_cell volume :) 2 0.2528 dg 0.105 190 monoclinic mC 248.48 71.93 154.76 90.00 112.70 90.00 2551913 :) 1 0.0000 dg 0.076 203 triclinic aP 71.93 129.30 154.76 111.69 90.25 106.07 1275944 MOSFLM Integration results: Solution SpaceGroup Beam x y distance Resolution Mosaicity RMS :) 2 C2 172.89 172.41 209.92 2.54 0.250000 0.109 1 P1 172.88 172.36 209.85 2.37 0.250000 0.094 show_stack(1): /userdata-2/christian/phenix- dev-968/labelit/command_line/screen.py(84) top_level_indexing show_stack(2): /userdata-2/christian/phenix-dev-968/build/intel- linux-2.6/../../labelit/command_line/index.py(5) <module> libc backtrace (28 frames, most recent call last): /userdata-2/christian/phenix-dev-968/build/intel-linux-2.6/base/bin/python [0x8058241] /lib/tls/i686/cmov/libc.so.6(__libc_start_main+0xe6) [0x4009fb56] /userdata-2/christian/phenix-dev-968/build/intel- linux-2.6/base/bin/python(main+0x17) [0x80582df] /userdata-2/christian/phenix-dev-968/build/intel- linux-2.6/base/bin/python(Py_Main+0xaba) [0x8058f2a] /userdata-2/christian/phenix-dev-968/build/intel- linux-2.6/base/bin/python(PyRun_SimpleFileExFlags+0xb5) [0x80f61fd] /userdata-2/christian/phenix-dev-968/build/intel- linux-2.6/base/bin/python(PyRun_FileExFlags+0xdc) [0x80f52cc] /userdata-2/christian/phenix-dev-968/build/intel- linux-2.6/base/bin/python(PyEval_EvalCode+0x22) [0x80d659e] /userdata-2/christian/phenix-dev-968/build/intel- linux-2.6/base/bin/python(PyEval_EvalCodeEx+0x603) [0x80d63f3] /userdata-2/christian/phenix-dev-968/build/intel- linux-2.6/base/bin/python(PyEval_EvalFrameEx+0x622e) [0x80d5ad6] /userdata-2/christian/phenix-dev-968/build/intel- linux-2.6/base/bin/python(PyEval_EvalFrameEx+0x3dc7) [0x80d366f] /userdata-2/christian/phenix-dev-968/build/intel- linux-2.6/base/bin/python(PyObject_Call+0x46) [0x806162a] /userdata-2/christian/phenix-dev-968/build/intel-linux-2.6/base/bin/python [0x8147760] /userdata-2/christian/phenix-dev-968/build/intel- linux-2.6/base/bin/python(PyEval_EvalCodeEx+0x4ba) [0x80d62aa] /userdata-2/christian/phenix-dev-968/build/intel-linux-2.6/base/bin/python [0x8145ea2] /userdata-2/christian/phenix-dev-968/build/intel-linux-2.6/base/bin/python [0x8064f09] /userdata-2/christian/phenix-dev-968/build/intel-linux-2.6/base/bin/python [0x8089ab0] /userdata-2/christian/phenix-dev-968/build/intel-linux-2.6/base/bin/python [0x8064f09] /userdata-2/christian/phenix-dev-968/build/intel-linux-2.6/base/bin/python [0x8089ab0] /userdata-2/christian/phenix-dev-968/build/intel-linux-2.6/base/bin/python [0x807b94a] /userdata-2/christian/phenix-dev-968/build/intel-linux-2.6/base/bin/python [0x80a3710] /userdata-2/christian/phenix-dev-968/build/intel-linux-2.6/base/bin/python [0x8089ab0] /userdata-2/christian/phenix-dev-968/build/intel-linux-2.6/base/bin/python [0x80a3677] /userdata-2/christian/phenix-dev-968/build/intel- linux-2.6/lib/libboost_python.so [0x4036b67b] /userdata-2/christian/phenix-dev-968/build/intel- linux-2.6/lib/cbflib_ext.so(boost::python::objects::value_holder<iotbx::detectors::Mar345Adaptor>::~value_holder()+0x4a) [0x42a13f22] /lib/tls/i686/cmov/libc.so.6(fclose+0x1a0) [0x400e4640] /lib/tls/i686/cmov/libc.so.6(fclose+0x41) [0x4018f241] [0x4001d400] /userdata-2/christian/phenix-dev-968/build/intel- linux-2.6/lib/boost_python_meta_ext.so(boost_adaptbx_segmentation_fault_backtrace+0x1f) [0x4032231f] Segmentation fault (Python and libc call stacks above) This crash may be due to a problem in any imported Python module, including modules which are not part of the cctbx project. To disable the traps leading to this message, define these environment variables (e.g. assign the value 1): BOOST_ADAPTBX_FPE_DEFAULT BOOST_ADAPTBX_SIGNALS_DEFAULT This will NOT solve the problem, just mask it, but may allow you to proceed in case it is not critical. -- Christian Roth Institut für Bioanalytische Chemie Biotechnologisch-Biomedizinisches Zentrum Fakultät für Chemie und Mineralogie Universität Leipzig Deutscher Platz 5 04103 Leipzig Telefon: +49 (0)341 97 31316 Fax: +49 (0)341 97 31319

2 1

There are still seats available for RapiData 2012, a course on Data Collection and Structure Solving at the NSLS.
by Robert Sweet 31 Jan '12

31 Jan '12

We are extending the deadline because seats are still available for RapiData 2012, the fourteenth 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 22-27 April 2012: 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 tenth time we will hold a short lecture course on the fundamentals of crystallography for roughly five hours on Sunday 22 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/ Biology Dept Brookhaven Nat'l Lab. Phones: Upton, NY 11973 631 344 3401 (Office) U.S.A. 631 344 2741 (Facsimile) =========================================================================

1 0

Rosetta-MR - How long does it take?
by Tri Ngo 28 Jan '12

28 Jan '12

Hi Phenix users, I tried to utilize the Rosetta-MR to solve an Esterase structure (389 residues, 2 molecules/asu) which has 22% sequence identity to the know structure. I used HHPred to prepare the alignment input and Robetta fragment server to generate the fragment files. I also performed a quick test of mr-rosetta by running the command phenix_regression.wizards.test_command_line_rosetta_quick_tests. I think there is no error on my input files or the system. However, the program never finished. It still shows me the "Placing Model..." status. It's running around 48 hours and still show the same log file (at "Running group of place_model" step). I am not sure if it's normal to wait days or weeks to get the work done when running Rosetta-MR? Thank you for your helpful input. Best wishes, Tri -- * TriNgo * Structural Biology Lab School of Medicine - Sungkyunkwan University Phone: 031-299-6150 Cell: 010-7774-8210

2 1

Problems with opening Coot and Firefox from gui under RHEL 6.2
by Brick, Peter 27 Jan '12

27 Jan '12

Hi, I've installed Phenix 1.7.3-928 on a linux workstation by downloading the binary phenix-installer-1.7.3-928-intel-linux-2.6-x86_64-fc3.tar I'm running Red Hat Enterprise Linux Workstation release 6.2. (Mainly 64 bit libraries but some additional 32 bit libraries that are necessary to run Alwyn Jones' 32 bit version of O.) *Help and Firefox* When I click on the Help icon in the Phenix GUI I get the message: [pbrick@pb1 ~/test]$ phenix phenix_gui.htm Using browser /usr/bin/firefox phenix_gui.htm Using browser /usr/bin/firefox /usr/lib64/firefox-3.6/firefox: symbol lookup error: /usr/lib64/libgnome-2.so.0: undefined symbol: g_dgettext But if firefox is already running it opens a new tab and I receive no error message. Firefox is the version shipped with RHEL 6.2: firefox 3.6.24 (I notice that the mozilla website is offering version 9.0.1 for i686 processors.) *Starting Coot* When I click on the Coot icon in the Phenix gui no window opens but I can see the processes: coot, coot_build_phenix_coot.csh and coot-real running. In fact coot-real consumes 100% of the cpu of one of the cores. What I find really weird is if I look at the processes with the system monitor gui and right click and 'kill' coot-real then instead of killing the process the coot window opens. At the same time coot_build_phenix_coot.csh disappears and the cpu usage drops to zero. (Note: I've also tried the fc12 version. I found that the version of pymol distributed with phenix uses 32 bit libraries and I needed to install libXmu.so.6) Any suggestions as to what I might be doing wrong? Peter

5 6

a story of acetoacetate decarboxylase
by Kevin Jin 27 Jan '12

27 Jan '12

Dear all, I wrote a story for acetoacetate decarboxylase catalysis. I hope you will like it. http://www.jinkai.org/AAD_history.html http://www.jinkai.org/AAD_catalysis.html Regards, Kevin

1 0

Bond to symmetry mate
by Eike Schulz 27 Jan '12

27 Jan '12

Dear phenixbb, I am sure this is a trivial problem encountered before and that there is a solution to it however, I was unfortunately not able to solve it myself or with the help of the online documentation. The problem: In one of my structures a nucleic acid chain extends over several symmetry mates. How can I create and refine a bond (in my case 5'-phosphate to 3'-end ) to that symmetry mate? Many thanks in advance Eike

2 1

Twin_law in GUI Phenix Refinement
by Tri Ngo 26 Jan '12

26 Jan '12

Hi phenix users, I have two questions when using Phenix Refinement. 1. I'd like to run phenix refinement with twin data but I don't know how to define twin_law using GUI Phenix refinement module in Linux. 2. I used phenix refinement in Windows to refine a 2.3A structure. My problem is the Rfree is quite high. - I applied twin refinement + TLS + NCS and SA refinement (wxu_scale = 1) Start: r_work = 0.2149 r_free = 0.2833 bonds = 0.016 angles = 1.988 Final: r_work = 0.2029 r_free = 0.2735 bonds = 0.008 angles = 1.448 - Then I use auto-weight to optimize the Rfree. optimize_xyz_weight = True optimize_adp_weight = True wxc_scale = 0.5 wxu_scale = 1 Start: r_work = 0.2152 r_free = 0.2735 bonds = 0.008 angles = 1.322 Final: r_work = 0.2309 r_free = 0.2781 bonds = 0.003 angles = 0.848 I think the bond rmsd and angle rmsd were too small in this final refinement. Do you think it's normal or I need to try more? Thank you very much for your helpful input. Best wishes, -- * TriNgo * Structural Biology Lab School of Medicine - Sungkyunkwan University Phone: 031-299-6150 Cell: 010-7774-8210

4 4

model_vs_data
by Karolina Michalska 26 Jan '12

26 Jan '12

Hi everybody, I 've run model_vs_data for one of my structures and here is what I got: Information extracted from PDB file header: program_name : REFMAC year : None r_work : 0.15793 r_free : 0.19845 high_resolution : 1.45 low_resolution : 36.92 sigma_cutoff : None matthews_coeff : None solvent_cont : None TLS : False (number of groups: 0) After applying resolution and sigma cutoffs: n_refl_cutoff : 26357 r_work_cutoff : 0.1450 r_free_cutoff : 0.1687 Can someone, please, explain why there is such a big difference in R values? When I do refinement in phenix I get practically the same results that Refmac produces. Thanks, Karolina

3 3

links in phenix
by Christian Roth 26 Jan '12

26 Jan '12

Hi, I have a pdb with a sugar chain as ligand. I used phenix ready set to get the links according to the link remarks in the pdb. However when I look into the edits file. All links are defined as geometry restraints and despite it is alpha 1-4 all bond restraints differ slightly in their length. Within the pdb file the definitions for the links were set to alpha1-4. A Similar thing happend with a peptide link between 2 Aminoacids which have not a continously numbering because of a deletion, but it is a simple peptide bond. I modified the .link.edits file accordingly for the sugars like this ---------- refinement.geometry_restraints.edits { apply_cif_link { data_link = ALPHA1-4 atom_selection_1 = name O4 and chain B and resname GLC and resseq 1 atom_selection_2 = name C1 and chain B and resname GLC and resseq 2 } } ------------ And applied a trans link for the peptide bond I added the modified file as additional parameter file in the gui. However in the final .eff file I found just the old restraints from the previous run. In the final .eff file it is stated that the file with my updated links is loaded ------ gui { base_output_dir = "/homes/roth/structure/250/082_c2_5/refinement/phenix" tmp_dir = "/homes/roth/structure/250/082_c2_5/refinement/phenix/phenix/tmp" send_notification = False notify_email = "christian.roth(a)bbz.uni-leipzig.de" add_hydrogens = True skip_rsr = False phil_file = "/homes/roth/structure/250/082_c2_5/refinement/phenix/refmac20link.link.edits" ----- Somehow the old restraits are already in the .eff file and the new ones are ignored how I could I avoid this. Simply closing the GUI an reopening? A second question is a very non standard peptide link due to a SNN residue in the chain. Based on the link a geometry restrain is created througn phenix.link_edits There is a slight difference in the length between the N- terminal and the C-terminal peptide bond. Which is probably fine and the resoution is good enough (roughly 2 Ang.) that is looks fine, but no angle definitions for phi psi and so on are stated. What is the best way to define this two bonds correctly for the subsequent refinement? A third question is a alpha1-4 link which extend to a symmetry related sugar. How I incorporate the symmetry information into the link? A last question regarding the hydrogens in refinement. For rebuilding it is nice to remove the hydrogens. Normally I use the commandline for this purpose. But It would be nice to do it automatically after refinement as a third pdb output in addition to the two standard outputs phenix_output.pdb and the phenix_output_reduce.pdb. Is there an option in the GUI to do this automatically. I thought I have seen this previously but cannot find it anymore I am sorry for this long Mail and so many question. I hope you could help me to clarify some the things. Best Regards Christian

3 5

DEN Refinement in phenix 1.7.3
by Eike Schulz 25 Jan '12

25 Jan '12

Dear phenixbb, I have noticed that in (since?) the latest release it is possible to activate DEN refinement in phenix.refine. I assume this refers to deformable elastic networks right? What is your experience with DEN refinements ? As I could not find any reference to it in the documentation I would be happy to have feedback on when to use and when to avoid it e.g. resolution limits etc. and what to watch out for. Many thanks in advance Eike

3 3

cif files
by Shya Biswas 24 Jan '12

24 Jan '12

Hi all, It appears that the cif file generated from prodrg (REFMAC format) is not recognised by phenix during refinement. Did anyone have similar issues or i am doing something wrong? somehow the geometries obtained from elbow are not right for my ligand so had to use prodrg to generate cif files. thanks, Shya

4 5

Refinement of two co-factors/ligands at the same position
by Kornelius Zeth 23 Jan '12

23 Jan '12

Dear all, how can I refine two different co-factors/ligands ('XXX and YYY', with differing occupancy and cif-Files) located at the same position within an enzyme? What about the chain IDs, segment IDs etc? Thanks and best wishes Kornelius -- Kornelius Zeth ZMBP, Uni-Tuebingen Auf der Morgenstelle 1 72076 Tübingen Germany

2 1

There's BARELY A WEEK left to apply for RapiData 2012, a course on Data Collection and Structure Solving at the NSLS.
by Robert Sweet 21 Jan '12

21 Jan '12

Half of the seats are still available for RapiData 2012, the fourteenth 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 22-27 April 2012: 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 tenth time we will hold a short lecture course on the fundamentals of crystallography for roughly five hours on Sunday 22 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/ Biology Dept Brookhaven Nat'l Lab. Phones: Upton, NY 11973 631 344 3401 (Office) U.S.A. 631 344 2741 (Facsimile) =========================================================================

1 0

phosphate stereochemistry
by Luca Pellegrini 21 Jan '12

21 Jan '12

Hi, I have noticed that Phenix mangles O-P bond lengths and angles in the phosphate groups of my DNA chain (example in pict attached). This only happens to 3 out of 28 phosphates. The refined protein-DNA structure does not present any other unusual stereochemistry issues and refinement runs normally. Could anybody advice on what is going on and how to fix this, please? I am using Phenix-1.7.3-928. Refinement is with default target weights, hydrogens on and does not include simulated annealing. I could make the stereochemistry targets more stringent, but the default values seem strict enough (bond and angle rmsd is 0.005 and 1.225). Thanks, Luca Luca Pellegrini Department of Biochemistry University of Cambridge 80 Tennis Court Road Cambridge CB2 1GA - UK Email: lp212(a)cam.ac.uk Tel: 0044-1223-760469 Fax: 0044-1223-766002 Sanger building, room 3.59

5 7

phenix.refine is auto selecting an old .pdb file for refinement
by Andrews, Forest Hoyt 21 Jan '12

21 Jan '12

I'm having a bit of an issue with phenix.refine. I was receiving the following error "number of groups of duplicate atom labels: 4789 total number of affected atoms: 9578" However when I looked inside my .pdb file there where only 4789 atoms total....No duplicates. I discovered that phenix.refine is auto selecting a .pdb from a previous run. I found this out by not selecting any .pdb file and found phenix.refine refining my data with an older model I use to refine a previous round. How can I stop phenix.refine from auto selecting a .pdb during refinement? __________________________________ From: phenixbb-bounces(a)phenix-online.org [phenixbb-bounces(a)phenix-online.org] on behalf of Pavel Afonine [pafonine(a)lbl.gov] Sent: Friday, January 20, 2012 12:20 PM To: PHENIX user mailing list Subject: Re: [phenixbb] phosphate stereochemistry Hi Luca, what kind of NCS restraints did you use: torsion or Cartesian? I wonder if torsion NCS would naturally avoid this problem? I hope Jeff comments on this. Pavel On 1/14/12 10:51 AM, Luca Pellegrini wrote: > I think I discovered what was happening to my DNA molecule during > refinement, that caused distortion of the phosphate groups. It was > something along the lines of what Dale suggested, I mention it here so > that it might be useful to others. > > I have two DNA molecules in the ASU and I didn't realised that I had > accidentally swapped OP1 and OP2 in 3 out of 28 phosphates (so that, > if the two DNA molecules were superimposed, OP1 ended on top of OP2). > Because I had NCS restraints on during refinement, I guess phenix was > trying to move OP1 in one DNA molecule towards the NCS-equivalent > position of OP1 in the second molecule, hence the altered geometry. > > I have fixed the atom nomenclature problem and now the geometry of my > DNA is ok. I appreciate that most people tend to avoid such mistakes > ;-) but perhaps it would be good if phenix could flag situations such > as these. > > Luca > > On 6 Jan 2012, at 16:15, Pavel Afonine wrote: > >> Hi Luca, >> >> phenix.refine always generates a *.geo file that lists all the >> geometry restraints (bonds, angles, planarities, chiralities, >> dihedral, non-bonded, ncs(if any)) used in refinement for all atoms. >> For each restraint it list current model value, target (library) >> value, etc. >> >> Can you have a look at *.geo file for O-P bonds in question? May be >> this gives a hint about what's going on? >> >> Pavel >> >> On 1/6/12 4:31 AM, Luca Pellegrini wrote: >>> >>> Hi, >>> >>> I have noticed that Phenix mangles O-P bond lengths and angles in >>> the phosphate groups of my DNA chain (example in pict attached). >>> This only happens to 3 out of 28 phosphates. The refined protein-DNA >>> structure does not present any other unusual stereochemistry issues >>> and refinement runs normally. Could anybody advice on what is going >>> on and how to fix this, please? >>> >>> I am using Phenix-1.7.3-928. Refinement is with default target >>> weights, hydrogens on and does not include simulated annealing. I >>> could make the stereochemistry targets more stringent, but the >>> default values seem strict enough (bond and angle rmsd is 0.005 and >>> 1.225). >>> >>> Thanks, >>> Luca _______________________________________________ phenixbb mailing list phenixbb(a)phenix-online.org http://phenix-online.org/mailman/listinfo/phenixbb

1 0

Real_space_refinement at low resolution question
by fulvio saccoccia 20 Jan '12

20 Jan '12

Dear phenix users, I was wondering if it could be correct to exclude the real_space_refinement step working with a 3.3 A resolution map: I noted that in every real_space_refinemet step phenix told me that it was not useful. How can I exclude the routine from the command line? I am using 1.7.3-928 version.

2 1

Real_space_refinement at low resolution question
by fulvio saccoccia 20 Jan '12

20 Jan '12

Dear phenix users, I was wondering if it could be correct to exclude the real_space_refinement step working with a 3.3 A resolution map: I noted that in every real_space_refinemet step phenix told me that it was not useful. How can I exclude the routine from the command line? I am using 1.7.3-928 version. Fulvio Saccoccia, PhD student Dept. of Biochemical Sciences Sapienza University of Rome

1 0

problem converting PDB file to mmcif for deposition
by Salameh, Mohd A., Ph.D. 20 Jan '12

20 Jan '12

Hi Everybody, I encountered a problem recently in converting my PDB file to mmcif using the pdb_extract tool in rcsb. It simply writes out an empty file!! I wonder if any of you have faced such a problem before and know the solution!! I carefully examined the PDB file header and I couldn't really find a reason why this PDB failed to be converted. I only use PHENIX for molecular replacement and refinement. I used the default TLS refinement so I assume the whole structure will be refined as a group. I hope somebody is familiar with this can help me out. Thank you very much, Mohd

2 1

real space correlation
by Jerome Nwachukwu 20 Jan '12

20 Jan '12

Hi: I am trying to run phenix.real_space_correlation from the command line and received the following error message. Is there a specific syntax i should use? Thank you, -Jerome Syntax error: expected "=", found "Date" (file "3_refine_001.pdb", line 1) Traceback (most recent call last): File "/usr/local/phenix-1.7.1-743/cctbx_project/mmtbx/command_line/real_space_correlation.py", line 9, in <module> command_name = "phenix.real_space_correlation") File "/usr/local/phenix-1.7.1-743/cctbx_project/mmtbx/real_space_correlation.py", line 332, in cmd_run track_unused_definitions = True) File "/usr/local/phenix-1.7.1-743/cctbx_project/libtbx/phil/__init__.py", line 1718, in fetch assert source.name == self.name AttributeError: 'NoneType' object has no attribute 'name'

2 2

noisy sa-omit maps in solvent region.
by Francis E Reyes 20 Jan '12

20 Jan '12

Hi phenixers I'm calculating a composite sa-omit map for my entire molecule at 4 A. While my molecule more or less has density (very good for backbone, not so much for sidechains), the solvent region seems to be particularly noisy. Is this the normal experience for composite sa-omit maps at this resolution? I'm particularly interested in what autobuild does when it 'stitches' the individual maps together, is it a direct stitch of the 2fo-fc's after refinement or is there some kind of real space density modification involved? Or a combination of both? Any suggestions appreciated, F --------------------------------------------- Francis E. Reyes M.Sc. 215 UCB University of Colorado at Boulder

2 1

adding exchangeable H to a perdeuterated structure with ready set?
by Maxime Cuypers 20 Jan '12

20 Jan '12

hello, how can i add exchangeable H to a perdeuterated structure with readyset? (other than manualy editing the file please) i see readyset does the opposite. thanks in advance Maxime

2 1

an error in the operation of SF-Tool
by Dialing Pretty 19 Jan '12

19 Jan '12

Dear All, By the RCSB SF-Tool I have tried to convert my PDB file and the corresponding mtz file into mmCIF for the purpose of RCSB depositation, by the Automatic (default) mode. I got the following message: "Converting structure factor to mmCIF format.... Error! No labels are written for MTZ to mmCIF conversion Structure factor conversion is done. Header of the mtz file can be seen from here Click here to see the converted mmCIF file" Obviously there is an error message, although the conversion is succesful. As for I intend to deposit the structure including the mmCIF file into RCSB, I am afraid I have lost some message in the MTZ file and mmCIF file. Will you please tell me how can I add the labels, which have been lost I think, to the mtz file, so that the labels message will also be included into the mmCIF file if it was necessary? Cheers, Dialing

7 18

tool for dispersive difference fourier ?
by Francis E Reyes 17 Jan '12

17 Jan '12

I'd like to make a dispersive difference fourier using phases from a model. Is there a tool to do so? phenix.maps doesn't seem to have the 'map_type' to do this.. Thanks! F --------------------------------------------- Francis E. Reyes M.Sc. 215 UCB University of Colorado at Boulder

1 0

Rosetta rebuilding...
by Francis E Reyes 17 Jan '12

17 Jan '12

Phenixers, I'm curious about using rosetta to rebuild an existing model (phenix.mr_rosetta it seems) [0] What is being rebuilt in rosetta? The whole model? Or is the model divided into chunks and each chunk sent to rosetta? [1] Can a specific region be chosen to be rebuilt? (In the case where most of the model is correct) [2] How are the rosetta models scored ? Specifically, "(d) refine Rosetta models, average density from top 20%, continue rebuilding each Rosetta model using averaged density". What's the scoring function (mapCC, R-Free, LLG??) How's the top 20% obtained? Thanks! F --------------------------------------------- Francis E. Reyes M.Sc. 215 UCB University of Colorado at Boulder

3 2

Only two weeks left to apply for RapiData 2012, a course on Data Collection and Structure Solving at the NSLS.
by Robert Sweet 14 Jan '12

14 Jan '12

Many seats are still available for RapiData 2012, the fourteenth 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 22-27 April 2012: 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 tenth time we will hold a short lecture course on the fundamentals of crystallography for roughly five hours on Sunday 22 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/ Biology Dept Brookhaven Nat'l Lab. Phones: Upton, NY 11973 631 344 3401 (Office) U.S.A. 631 344 2741 (Facsimile) =========================================================================

1 0

Simulated annealing omit maps
by SEKHAR1＠pitt.edu 13 Jan '12

13 Jan '12

Hello friends, I have a couple of questions regarding simulated annealing composite map calculations. (1) Is there a way of specifying the percentage of residue co-ordinates to be omitted each time as we move through the asymmetric unit? I have found, when using CNS earlier, that depending on the structure, the map quality often varies if you change this from the usual 5% to 4% or even 3%, sometimes getting clearer. (2) Does the 'number of boxes' you can specify in the GUI only affect the fine structure/grid spacing of the map? The map probably does not get better all the time as we keep on increasing the number of boxes; when do we reach the point of diminishing returns and what would be the best way to estimate it? Thank you very much and wish you all a Very Happy 2012. K. Chandrasekhar

2 1

phenix.find_ncs
by Ursula Schulze-Gahmen 12 Jan '12

12 Jan '12

I am working with a complex of 3 different proteins, and there are 3 complexes in the asymmetric unit. During phenix.refine I used automatic NCS and it worked fine. I now wanted to determine the NCS parameters with phenix.find_ncs. The output file gives me 3 NCS groups, 2 groups with 3 chains, and 1 group with 2 chains.I cannot figure out why one chain is not included. I tried to run it with only pdb input, or pdb anf map coefficients, same result. Any suggestions hat else to try? Thanks Ursula -- Ursula Schulze-Gahmen, Ph.D. Assistant Researcher UC Berkeley, QB3 356 Stanley Hall #3220 Berkeley, CA 94720-3220

2 4

brainpool research associate position
by intekhab alam 12 Jan '12

12 Jan '12

Dear members Have a look on the advertisement for the position of research associate in structural biology. Interested candidates please reply to khkim(a)korea.ac.kr. Informal enquiries can also be made. Structural Bioinformatics Lab, Department of Biotechnology & Bioinformatics, Korea University, Seoul KoreaBrainpool Research Associate – Full time, fixed term appointment for up to two years The Structural Bioinformatics Lab is seeking a Brainpool Research Associate to carry out structural analysis (x-ray crystallography) and biophysical studies of complexes related to virus-host interaction. The group is led by Professor Kyung Hyun Kim and has over the years produced good results (PNAS, NAR, JGV). The research programme is funded by a grant from the Korean Federation of Science and technology Societies (www.brainpool.or.kr) to which both the applicant and our team apply together. Further details on the research group can be found at http://sbl.korea.ac.kr/. Applicants should have a PhD in Structural Biology, Biophysics or a related area, with at least 5 years of research experience after Ph.D., or high profile publications in the highest impact journals. Salary range will be from $45,000-$55,000 per annum plus plane ticket, moving expense and medical insurance. Initial salary for possibly the two months will be dependent on the skills and experience of the successful applicant. Closing date for completed applications: 5 February 2012. To apply for this position please send emails to khkim(a)korea.ac.kr with your resume. Korea U is an equal opportunities employer and encourages applications from all candidates irrespective of gender, race, disability, sexual orientation, age, religion and belief or another protected characteristic. -- INTEKHAB ALAM LABORATORY OF STRUCTURAL BIOINFORMATICS KOREA UNIVERSITY, SEOUL

1 0

Re: [phenixbb] phenix.find_ncs
by Thomas C. Terwilliger 12 Jan '12

12 Jan '12

Hi Ursula, In that case...sounds like something I need to fix! Would it be possible to send me (off-list, terwilliger(a)lanl.gov) the entire PDB file and I'll have a look? In the meantime, I guess if you run each group separately you have the information that you are looking for... All the best, Tom T >> I just did waht you suggested. I ran phenix.find_ncs on the problem group >> only. This time it found all 3 chains in the group. >> >> Ursula >> >> On Wed, Jan 11, 2012 at 10:16 AM, Thomas C. Terwilliger < >> terwilliger(a)lanl.gov> wrote: >> >>> Hi Ursula, >>> >>> I am not sure why find_ncs did not get the 3rd chain for the third >>> group. >>> I think that when you run phenix.refine it is (for the moment) actually >>> just calling find_ncs to get the NCS groups, so the groups found there >>> should (if everything is working the way I think it should) be identical >>> to the ones you would find with phenix.find_ncs. >>> >>> Here is something to try: Split your model into just protein #3 (you can >>> use phenix.pdbtools and keep="chain C or chain F or chain I" to do >>> this), >>> and also just protein #1 and #2. Then run find_ncs separately on these. >>> Can find_ncs get the ncs for your 3 copies of protein #3 all by >>> themselves? >>> >>> If not...then it indicates that perhaps one of the chains is actually >>> different (maybe very different conformation, maybe mostly missing, >>> maybe >>> different sequence...?) >>> >>> All the best, >>> Tom T >>> >>> >>> >> I am working with a complex of 3 different proteins, and there are 3 >>> >> complexes in the asymmetric unit. During phenix.refine I used >>> automatic >>> >> NCS >>> >> and it worked fine. I now wanted to determine the NCS parameters with >>> >> phenix.find_ncs. The output file gives me 3 NCS groups, 2 groups with >>> 3 >>> >> chains, and 1 group with 2 chains.I cannot figure out why one chain >>> is >>> not >>> >> included. I tried to run it with only pdb input, or pdb anf map >>> >> coefficients, same result. Any suggestions hat else to try? >>> >> >>> >> Thanks >>> >> >>> >> Ursula >>> >> >>> >> -- >>> >> Ursula Schulze-Gahmen, Ph.D. >>> >> Assistant Researcher >>> >> UC Berkeley, QB3 >>> >> 356 Stanley Hall #3220 >>> >> Berkeley, CA 94720-3220 >>> >> _______________________________________________ >>> >> 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 >>> >> >> >> >> -- >> Ursula Schulze-Gahmen, Ph.D. >> Assistant Researcher >> UC Berkeley, QB3 >> 356 Stanley Hall #3220 >> Berkeley, CA 94720-3220 >>

1 0

Phenix for CUDA?
by Schubert, Carsten [JRDUS] 11 Jan '12

11 Jan '12

Hi, are there any development plans for phenix in the works to take advantage of NVIDIA's CUDA platform? Cheers, Carsten

3 2

Question about unrestrained refinement
by Louis Lazar 11 Jan '12

11 Jan '12

Hello all, I had a question about unrestrained refinement. I have several data sets (right near 1 A resolution) that are nearly fully refined, but to finish I would like to run one round of unrestrained refinement. At this point, I want to focus on unrestrained refinement of just the protein, and not any waters. How would the syntax for refining just the chains of the protein (and not water) look in the command line? I also want to include anisotropic refinement of B-factors in this last refinement run. Should it look something like: strategy=individual_sites+occupanices+individual_adp \ adp.individual.anisotropic="not water and not element H" wc=0 wu=0 Then what should I add after this to include just my protein chains and omit waters from the unrestrained refinement? Also, if I wanted to not refine the waters at all in this round of refinement, what should I add to the command line to make that happen? Thank you in advance as always. Regards, Louis

2 1

zn anomalous
by Pengfei Fang 11 Jan '12

11 Jan '12

Dear All, I collected a dataset at 0.9774 angstrom. the structure is already by MR. Can you please teach me how to tell a specific site has Zn-atom anomalous signal or not using phenix? thank you! Pengfei

2 1

phenix.refine request: output den restraints to a file...
by Francis E Reyes 11 Jan '12

11 Jan '12

Hi phenixbb (Jeff probably), I wanted to try a few different protocols with the den refinement. It occurred to me that phenix.refine won't take den restraints as an external file, instead it seems to be generating them on the fly. Can I have an option to output (and also read in) den restraints? I'd like to generate them once and use the same restraints to test a different refinement params. F --------------------------------------------- Francis E. Reyes M.Sc. 215 UCB University of Colorado at Boulder

2 1

Output for den optimization runs
by Francis E Reyes 10 Jan '12

10 Jan '12

Is there a refine.log for each den (gamma, weight) optimization? I can't seem to find one anywhere and I'm curious to see how it's doing. Thanks! F --------------------------------------------- Francis E. Reyes M.Sc. 215 UCB University of Colorado at Boulder

2 5

Chirality of Ribityl chain of FMN inverted
by Yuri 09 Jan '12

09 Jan '12

Hello Developers, When I built Ligand restraints using eLBOW, my resulting .pdb had an inverted chiral center in the ribityl chain of the flavin mononucleotide particularly C4'. In the database the ligand had an S configuration (correct), but after minimization and restraint generation the output had R configuration at that particular center. I wonder if there is a way to enforce the configuration for each center, although I suspect it could be a little bug, because I noticed my restraints have that center annotated as S eventhough it is R in the actual .pdb Cheers, -- Yuri Pompeu

4 4

molecular replacement of protein/DNA complex
by Wei Shi 08 Jan '12

08 Jan '12

Hi all, I have been trying to solve a protein-DNA complex structure using molecular replacement. I suspect two copies of protein bind one piece of the DNA, and the angle between the two copies of protein is somewhere between 130 to 180 degrees. I could get molecular replacement solution using a search model of protein/DNA complex I made with deletion in part of the protein and with the angle between the two proteins about 145 degrees. Only 1 copy of this two-protein/one-DNA complex is expected in the ASU. Below is the statistics for the solution. RFZ=2.9 TFZ=6.9 PAK=0 LLG=111 TFZ==14.6 LLG=111 And when I open the pdb file and generate symmetry mates, no clashes and I could see contacts between the end of DNA, but space instead of contacts between layers of the protein/DNA complex. I suspect that the angle between the two proteins I use in the search model is not quite right. No solution if I search with DNA and protein separately. I am trying to find a way to move the protein a little bit around the original solution and see whether that can help me make the solution better. I just read that “ROTate AROUnd option of the brute rotation serch” and “NMAPdb” might do this? Does anyone have any suggestions about what I could try? Thank you so much! Best, Wei

2 2

Only three weeks left to apply for RapiData 2012, a course on Data Collection and Structure Solving at the NSLS.
by Robert Sweet 06 Jan '12

06 Jan '12

We are offering RapiData 2012, the fourteenth 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 22-27 April 2012: 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 tenth time we will hold a short lecture course on the fundamentals of crystallography for roughly five hours on Sunday 22 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/ Biology Dept Brookhaven Nat'l Lab. Phones: Upton, NY 11973 631 344 3401 (Office) U.S.A. 631 344 2741 (Facsimile) =========================================================================

1 0

CCI-Apps - Cygwin incompatibility after upgrade
by Schubert, Carsten [JRDUS] 06 Jan '12

06 Jan '12

Hi, After upgrading my cci_apps install on Windows from 1.7-718 to 1.7.3-928 the interoperability with cygwin is gone. Calling phenix.refine just gives a 'Command not found' error. Turns out that the little .exe files for each subprogram were replaced by .bat files. That is fine if one just wants to use the DOS CLI, but cygwin is not able to call programs via .bat files as far as I know. Is there a workaround for this problem besides working in the DOS CLI (yuk)? Thanks Carsten

2 2

phenix.metal_coordination
by Young-Jin Cho 04 Jan '12

04 Jan '12

Hi everyone and wish your happy new year. I have a question regarding phenix.metal_coordination. when I typed my pdb file that has two alternative conformations of different species including metal such as a metal(B) and a water(A). It seems like, in this case, phenix.metal_coordination doesn't read the metal. Otherwise, is there other way I can let make .edits file for this pdb file? Thanks always in advance, Young-Jin

3 2

whether my following data are acceptable
by Dialing Pretty 03 Jan '12

03 Jan '12

Dear All, Following is the data range of my x-ray statistics by resolution bin: FOM between 0.84 to 0.86 Phase error between 17.90 to 24.14 scale factor from 0.88 to 1.07 SigmaA from 0.85 to 0.95 #work from 2688 to 2922 #test from 133 to 157 Do these data are within the acceptable range? Cheers, Dialing

4 10

reference model restraints
by Ursula Schulze-Gahmen 24 Dec '11

24 Dec '11

I am refining a 3.1 A structure and are trying to use model restraints for this. There are 3 molecules in the asymmetric unit of my structure. I input the high resolution structure of one molecule of the reference structure to create restraints during the refinement. From looking at the log file, it seems to me that the reference structure is aligned to only one of my molecules, but not the other 2. How can I make sure the restraints are valid for all 3 molecules. Ursula -- Ursula Schulze-Gahmen, Ph.D. Assistant Researcher UC Berkeley, QB3 356 Stanley Hall #3220 Berkeley, CA 94720-3220

4 3

eLBOW generates Mogul geometries
by Nigel Moriarty 24 Dec '11

24 Dec '11

Folks For a few weeks now eLBOW has been able to generate Mogul geometries. If you have CCDC's mogul installed and accessible at the command-line, eLBOW can request Mogul provide an experimentally accurate geometry. Simply use the --mogul option at the command-line and eLBOW will query Mogul for the geometry to use to generate a restraints file and a PDB file. Stay tuned for more details in the Computational Crystallography Newsletter. Happy Holidays, Merry Christmas and Happy New Year Nigel -- Nigel W. Moriarty, PhD 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 0

partial model in autobuild?
by Yarrow Madrona 22 Dec '11

22 Dec '11

Hello, I would like to use autobuild to start from a partial model. I have done Phaser EP with phases from MR. I can see the fourth molecule but it is not well defined and I am not exactly sure where to place it. Is it possible to run autobuild from the partial model to fill in the fourth model? Thanks. -Yarrow -- Yarrow Madrona Graduate Student Molecular Biology and Biochemistry Dept. University of California, Irvine Natural Sciences I, Rm 2403 Irvine, CA 92697

3 2

Re: [phenixbb] Reflection file woes phenix/ccp4/PDB
by Machius, Mischa Christian 21 Dec '11

21 Dec '11

Robbie, Yes, the thought that the anomalous R free flags could have something to do with it occurred to me as well, so I checked. As far as I can see, the test reflections are evenly distributed between the + and - sets; they were not selected from the + or - set only. For the vast majority of the test reflections, both Bijvoet mates were measured and both were flagged properly. For some, however, only one Bijvoet mate was measured and flagged. It probably wouldn't do any harm to consolidate the test flags, discarding the few with only one measured Bijvoet mate, but I really don't see why it should be necessary to do so. For now, I have converted the reflection file into the CNS format and sent to the PDB, hoping that they will be able to deal with that one properly. Bit surprised though that no one has run into that problem before. Thanks for your suggestion. MM On Dec 16, 2011, at 3:54 AM, Robbie P. Joosten wrote: > Hi Mischa, > > Your mtz file seems to contain 'anomalous' free R-flags (that is what I gather from Pavel's post). This is kind of weird. Do you really want I+ in the work set and I- in the test set? That would bias your R-free. If the two flags match (which I really hope), then you can leave one column out before you (re)deposit your data. This will hopefully solves the problem. > > Cheers, > Robbie Joosten > > Netherlands Cancer Institute > > On 12/15/2011 05:40 PM, Machius, Mischa Christian wrote: >> Y'all, >> >> While depositing a structure, the PDB 'complained' about my structure-factor file, saying they don't contain R-free flags or are otherwise incomplete. The comment was prompted by their Refmac R-factor test failing. >> >> The structure-factor file I submitted came from phenix.refine (xxx_data.mtz). I am using version 1.7.3-928. >> >> Indeed, using this file in CCP4 fails with the message: >> "CCP4MTZfile: Mtz column type mismatch: I-obs(+) K-J CCP4MTZfile: Mtz column type mismatch: SIGI-obs(+) M-Q" >> >> The file contains anomalous data. I have no problems with files that do not contain anomalous data. >> >> Needless to say that Phenix uses the structure-factor file just fine when reading it back into, say, its validation utility. >> >> When I use iotbx.reflection_file_reader, I cannot detect any anomalies. >> >> I am at a loss as to how to resolve this issue. I would like to deposit the data as prepared by Phenix so that the PDB tools can use them. >> >> Anybody ran into this issue before? Any ideas for how to resolve it? >> >> Many thanks in advance! >> >> MM >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://phenix-online.org/mailman/listinfo/phenixbb

4 4

Failure in installing Phenix-1.7.3
by arka chakraborty 21 Dec '11

21 Dec '11

Hi all I am facing a problem while installing phenix 1.7.3 on Redhat Enterprise 5 OS. I have tried with fedora 3 binary as instructed in the download page and also with others but nothing works. I have tried installing from the root. Every time the error log is states that GLIB 2.4 is unavailable which is required by user/local/.. Please help! Thanks in advance, ARKO -- *ARKA CHAKRABORTY* *CAS in Crystallography and Biophysics* *University of Madras* *Chennai,India*

3 2

SF-Tool check
by Dialing Pretty 20 Dec '11

20 Dec '11

Dear All, I have done the SF-Tool check for my protein crystal structure, and the result is as following. My 3 questions are (1) What does the following mean and how to solve the problem "Error! No labels are written for MTZ to mmCIF conversion Structure factor conversion is done?" (2)Where there is no "Calculated phenix R_factor(all data)"? (3) Where there is no "Correlation (Fo-Fc)"? I am looking forward to getting your reply. Dialing Error! No labels are written for MTZ to mmCIF conversion Structure factor conversion is done. Header of the mtz file can be seen from here Click here to see the converted mmCIF file Checking Xray Data by PHENIX (model_vs_data) ... PHENIX (model_vs_data) is done. Checking structure factor quality ... PHENIX.xtriage is done. Summary of Results (XXXX : ) Reported PHENIX Calculated phenix High resolution 2.829 2.83 Low resolution 50.937 50.94 Number of Reflection17446 17446 Completeness 99.98 100.0 R_factor (all data) 0.2134 ? R_factor (work set) 0.2106 0.2093 R_free 0.2668 0.2667 Correlation (Fo-Fc) ? ?

2 1

Prevent MacPyMOL from registering as a valid pdb/pml/etc viewer.
by Francis E Reyes 20 Dec '11

20 Dec '11

Hi phenix developers I just installed 1.7.3-928. Apparently, OS X thinks the internally distributed pymol is my preferred pdb viewer (as opposed to my more recent build in /Applications). This may have something to do with the Info.plist registering document types. A small change in Finder for the user, but it was quite maddening to figure out why I had opened a pdb and pseudoatom was not working. Can future builds keep the internally distributed pymol as walled in as possible? Thanks! F --------------------------------------------- Francis E. Reyes M.Sc. 215 UCB University of Colorado at Boulder

2 3