- phenixbb - phenix-online.org

Re: [phenixbb] low resolution data/parameter ratio
by Pavel Afonine 14 Aug '09

14 Aug '09

Hi Gerwald, > Two more quick questions: > how would you identify the best strategy? I would look at Rfree and Rfree-Rwork, and maps of course. To help in deciding what are the "good / reasonable" values, use POLYGON (Acta Cryst. D65, 297-300 (2009)) (phenix.polygon), local real-space correlation (phenix.real_space_correlation). I think Nat bundled them all into his magic tool called phenix.validation with will do this and more (it is GUI based). > I would aim for the lowest number of parameters with reasonably low > Rfree? Yes, this sounds right in general. > Does phenix report the total number of parameters refined anywhere in > the log? It's trivial and not at the same time (at least the way I understand it, of course I might be very wrong). It's trivial to compute the number of refinable parameters given the strategy you choose, but how you will count the restraints which are kind of quasi-observables, and what about weight that if chosen too strong so it can turn the restraints into quasi-constraints? Pavel.

2 2

low resolution data/parameter ratio
by gerwald jogl 14 Aug '09

14 Aug '09

Hi All, I am struggling with low resolution data (3.6 A) and I am concerned about overfitting. I have a ratio of about 3 reflections per atom. Possible refinement scenarios would be a) coordinates+tls, b) coordinates, group_adp, tls or c) coordinates, individual_adp, tls. The last scenario seems to be pushing the limits. I wonder if there are studies, thoughts or opinions out there on which refinement scenario to use. Thanks for any comment on this, Gerwald

2 1

adding sutrates in phenix.refine
by Rune wederkinck Andersen 14 Aug '09

14 Aug '09

Hi everyone I'm running a refinement in phenix.refine. I'm now trying to ad two subtrates AMP and RP5. The phenix.refine seams to reconise AMP but not the RP5 molecules. It seams as if this molecules is well known for the program (phenix-1.4-3/ext_ref_files/chemical_components/r/data_RP5.cif). In my PDB file the two molecules are described like this: ... HETATM 19 C6 AMP A 301 35.852 19.986 36.427 1.00 22.32 A HETATM 20 N6 AMP A 301 35.448 19.736 37.693 1.00 24.44 A HETATM 21 N7 AMP A 301 36.400 17.540 35.616 1.00 23.01 A HETATM 22 C8 AMP A 301 36.853 17.097 34.416 1.00 21.98 A HETATM 23 N9 AMP A 301 37.034 18.126 33.560 1.00 19.84 A HETATM 1 P' RP5 A 302 30.291 24.910 31.815 0.25 10.66 A HETATM 2 C5' RP5 A 302 32.348 25.416 30.301 0.25 14.76 A HETATM 3 O5' RP5 A 302 31.481 24.448 30.868 0.25 12.67 A HETATM 4 C4' RP5 A 302 33.260 24.753 29.307 0.25 15.58 A HETATM 5 C3' RP5 A 302 33.017 25.227 27.872 0.25 15.17 A HETATM 6 O3' RP5 A 302 33.951 26.264 27.522 0.25 14.15 A HETATM 7 C2' RP5 A 302 33.282 23.935 27.083 0.25 16.11 A HETATM 8 O2' RP5 A 302 34.571 23.963 26.478 0.25 13.68 A HETATM 9 C1' RP5 A 302 33.351 22.848 28.160 0.25 16.92 A HETATM 10 O1' RP5 A 302 32.775 23.415 29.296 0.25 16.62 A HETATM 11 O1 RP5 A 302 32.631 21.639 27.867 0.25 20.12 A HETATM 12 O1X RP5 A 302 29.112 25.373 30.809 0.25 10.53 A HETATM 13 O2X RP5 A 302 29.912 23.844 32.781 0.25 10.59 A HETATM 14 O3X RP5 A 302 30.793 26.258 32.567 0.25 9.63 A ATOM 1 N MET B 1 63.305 2.511 45.667 1.00 64.80 B ATOM 2 CA MET B 1 62.528 3.463 46.463 1.00 68.29 B ATOM 3 C MET B 1 62.555 4.863 45.882 1.00 65.73 B ATOM 4 O MET B 1 63.607 5.325 45.446 1.00 64.47 B ... I have tried with phenix.elbow as well but didn't work... Please help... Rune _________________________________________________________________ More than messages–check out the rest of the Windows Live™. http://www.microsoft.com/windows/windowslive/

5 8

Running Autobuild after Autosol
by Katya Heldwein 14 Aug '09

14 Aug '09

Hi, how do I specify running Autobuild after Autosol, in the command line mode? I see the "phenix.autobuild after_autosol" option in the manual, but how will Autobuild know which Autosol run to use? Katya -- Ekaterina Heldwein, Ph.D. Assistant Professor Department of Molecular Biology and Microbiology Tufts University 136 Harrison Ave Boston, MA 02111 phone: 617-636-0858 fax: 617-636-0337 e-mail: katya.heldwein(a)tufts.edu

3 3

CC beamtime application deadline this saturday
by Banumathi Sankaran 14 Aug '09

14 Aug '09

Dear Users, The deadline for October/November/December 2009 Collaborative Crystallography proposals will be August 15, 2009. Thumbnail Description: Through the Collaborative Crystallography Program (CC) at the Advanced Light Source (ALS), scientists can send protein crystals to Berkeley Center for Structural Biology (BCSB) staff researchers for data collection and analysis. The CC Program can provide a number of benefits to researchers: * Obtain high quality data and analysis through collaborating with expert beamline researchers; * Rapid turn around on projects; and * Reduced travel costs. To apply, please submit a proposal through the ALS General User proposal review process for beamtime allocation. Proposals are reviewed and ranked by the Proposal Study Panel, and beamtime is allocated accordingly. BCSB staff schedule the CC projects on Beamlines 5.0.1 and 5.0.2 to fit into the available resources. Only non-proprietary projects will be accepted. As a condition of participation, BCSB staff researchers who participate in data collection and/or analysis must be appropriately acknowledged - typically being included as authors on publications and in PDB depositions. Please consult the website for additional information at: http://bcsb.als.lbl.gov/wiki/index.php/Collaborative_Crystallography ============= How To Apply: To learn more, go to: http://www-als.lbl.gov/als/quickguide/independinvest.html To submit a proposal, go to: http://alsusweb.lbl.gov/. Scroll down to "Structural Biology beamlines (includes protein SAXS)" and click on "New Proposal." Enter your proposal information, with attention to the following details: * For question 3/First choice, select "5.0.1 (Monochromatic);" for question 3/Second choice, select "5.0.2 (MAD)." * Check box (yes) in response to question "Do you want collaborative crystallography (beamline 5.0.1 or 5.0.2 only ?) * In question 9, please describe a specific research project with a clear end point. In order to request CC time for October/November/December 2009 allocation period, proposals must be submitted by August 15, 2009. Please note that the ALS will be shut down in September, 2009. The deadline for CC proposals for the time period January/February 2010 will be November 15, 2009. Regards, Banumathi Sankaran

1 0

Re: [phenixbb] phenix.refine special positions
by Ralf W. Grosse-Kunstleve 14 Aug '09

14 Aug '09

Hi Simon, To add to the previous answers, note that phenix.refine currently doesn't work well if you want to allow symmetry copies of your ligands to overlap. If you move to individual coordinate refinement you have to fix such ligands; i.e. it won't automatically repair bonds distorted by rigid-body refinement. Removing this limitation is on my to-do list. Ralf

1 0

phenix.refine special positions
by Simon Kolstoe 14 Aug '09

14 Aug '09

Dear Pavel et al. Is there any way of defining a residue for rigid body refinement if a couple of the atoms are sitting on a symmetry axis? i.e. can I define the rigid body to be chain L but not the atoms that are giving the problems? I keep on getting the error: ================== Extract refinement strategy and selections ================= Some atoms selected for rigid body refinement are in special positions (listed below). Please remove these atoms from the rigid body selection and run refinement again. "ATOM 4 C4 A12 L 2 .*. C " Sorry: Atoms at special positions are within rigid groups. Thanks, Simon

3 2

Shake
by esmael 13 Aug '09

13 Aug '09

Hi, Would be possible to only add "Shake" to the side-chain atoms during the refinement with Phenix? Thanks, Esmael

2 1

phenix.xtriage says that entire res range contains anomalous signal but chart says differently?
by Francis E Reyes 13 Aug '09

13 Aug '09

Hi all I ran xtriage on my dataset to 4.5.here's the relevant section pertaining to anomalous signal Analyses of anomalous differences Table of measurability as a function of resolution The measurability is defined as the fraction of Bijvoet related intensity differences for which |delta_I|/sigma_delta_I > 3.0 min[I(+)/sigma_I(+), I(-)/sigma_I(-)] > 3.0 holds. The measurability provides an intuitive feeling of the quality of the data, as it is related to the number of reliable Bijvoet differences. When the data are processed properly and the standard deviations have been estimated accurately, values larger than 0.05 are encouraging. Note that this analyses relies on the correctness of the estimated standard deviations of the intensities. unused: - 90.6596 [ 0/1 ] bin 1: 90.6596 - 9.2614 [1371/1382] 0.0375 bin 2: 9.2614 - 7.3520 [1369/1376] 0.0033 bin 3: 7.3520 - 6.4230 [1380/1388] 0.0000 bin 4: 6.4230 - 5.8358 [1353/1359] 0.0016 bin 5: 5.8358 - 5.4176 [1367/1371] 0.0000 bin 6: 5.4176 - 5.0982 [1361/1369] 0.0000 bin 7: 5.0982 - 4.8429 [1378/1382] 0.0000 bin 8: 4.8429 - 4.6321 [1351/1359] 0.0000 bin 9: 4.6321 - 4.4538 [1383/1389] 0.0000 bin 10: 4.4538 - 4.3001 [1390/1393] 0.0000 unused: 4.3001 - [ 0/0 ] The full resolution range seems to contain a useful ammount of anomalous signal. Depending on your specific substructure, you could use all the data available for the location of the heavy atoms, or cut the resolution to speed up the search. If values larger than 0.05 are encouraging why is it telling me that the entire resolution contains a useful amount of anomalous signal? Did I process this incorrectly? Thanks FR --------------------------------------------- Francis Reyes M.Sc. 215 UCB University of Colorado at Boulder gpg --keyserver pgp.mit.edu --recv-keys 67BA8D5D 8AE2 F2F4 90F7 9640 28BC 686F 78FD 6669 67BA 8D5D

6 11

General User Proposal Deadline for PX/SAXS at the ALS: August 15, 2009
by Peter Zwart 12 Aug '09

12 Aug '09

Dear All, August 15, 2009 is the deadline for the October - December 2009 Rapid Access Proposal cycle. (!!!!Please note that the ALS will be shut down in September, 2009.!!!!) Remote data collection is available on all beamlines at the Berkeley Center for Structural Biology (BCSB). The NX remote desktop technology provides the remote user with the full complement of sample visualization, sample manipulation, beamline control, data acquisition and data analysis tools exactly as they would see them if they were stationed at the beamline. This enhanced remote operation capability is coupled with 24hr onsite support by BCSB staff who are able to assist immediately with loading additional samples for remote users or troubleshoot any issues that might arise. Beamlines 8.2.1 and 8.2.2: Both beamlines are now fully remote capable, which means that users are able control sample mounting, data collection and analysis from their home labs. This has dramatically increased the efficiency of data collection and significantly reduced the cost of protein crystallography experiments, since users no longer have to travel to the synchrotron. To facilitate studies on small crystals, a microdiffractometer was installed in the beamline 8.2.1 endstation. The new equipment allows precise sample positioning to within 2 microns, excellent sample viewing of very small crystals, and an off-axis crystal positioning stage. As part of the beamline optics upgrades, a new toroidal focusing mirror was purchased and installed in beamline 8.2.1; new parabolic mirrors are currently being polished. When completed, the new optics will increase the brightness of beamline 8.2.1 by a factor of 4 to 5 and decrease the spot size to 30 microns, greatly facilitating structure solution on large macromolecular complexes. Beamlines 5.0.1, 5.0.2, 5.0.3: On Sector 5, improvements to the existing Berkeley Automounter sample handling system have dramatically improved both total capacity and reliability. Upgrades to the 5.0.1 and 5.0.3 monochromators have been completed. The operational energy of both beamlines has been increased, taking them above the Se K-edge and thus introducing a significant new Se SAD capability for experiments on Se derivatized proteins. The new monochromators also provide significantly better beam focus and performance under the increased power loads of ALS topoff operation. These improvements have been coupled with the introduction of enhanced positional feedback, diagnostic tools, and auto-alignment systems for the beamline components to improve the overall stability and quality of the beam delivered to users. Another mode of access to the BCSB beamlines is through the Collaborative Crystallography (CC) Program. Users apply for beamtime via the general user program, and collaborate with an expert crystallographer who will conduct the experiments and data reduction on behalf of the researchers. Depending on the users, structure solution, model building and refinement can be carried out as well. Please visit http://bcsb.lbl.gov/ for more details about the Center and its beamlines. To find out more, click on: http://www.als.lbl.gov/als/quickguide/independinvest.html Scroll down to "Structural Biology/Small Angle X-Ray Scattering (SAXS)." We invite you to submit a proposal at: http://alsusweb.lbl.gov/4DCGI/WEB_GetForm/PXProposalEntry.shtml/Initialize If you'd like to apply for January-February 2010 beamtime at the Advanced Light Source, please submit a General User proposal by November 15, 2009. If you have any questions or would like to request open beamtime, please e-mail BCSBBeamtime(a)lbl.gov. Please note that executed user agreements must be received by LBNL prior to beamtime. Proprietary fees, if applicable, must be received by LBNL at least five working days prior to scheduled beamtime. -- ----------------------------------------------------------------- P.H. Zwart Beamline Scientist Berkeley Center for Structural Biology Lawrence Berkeley National Laboratories 1 Cyclotron Road, Berkeley, CA-94703, USA Cell: 510 289 9246 BCSB: http://bcsb.als.lbl.gov PHENIX: http://www.phenix-online.org CCTBX: http://cctbx.sf.net -----------------------------------------------------------------

1 0

xtriage
by Francis E Reyes 12 Aug '09

12 Aug '09

Hi all I love the fact that xtriage can measure anomalous signal, but can it also do a quick check for isomorphous differences (with and without anomalous signal) similar to scalepack? Is there a phenix utility that can do this? Thanks FR --------------------------------------------- Francis Reyes M.Sc. 215 UCB University of Colorado at Boulder gpg --keyserver pgp.mit.edu --recv-keys 67BA8D5D 8AE2 F2F4 90F7 9640 28BC 686F 78FD 6669 67BA 8D5D

3 4

duplicated label
by Jianghai Zhu 12 Aug '09

12 Aug '09

Hi, I am using phenix.refine-125 to do a MLHL refinement. I got an error at the end of the refinement while it was trying to exporting results. "Sorry: Duplicate mtz_label_amplitudes:FOFCWT". The output map_coeffs.mtz file is corrupted and the next def file is blank. The refined pdb file seems to be fine. Anyone knows how to fix this problem? Thanks. -- Jianghai

2 1

where are the optimized wxc_scale and wxu_scale when optimize_wxc_scale and optimize_wxu_scale = True?
by Francis E Reyes 07 Aug '09

07 Aug '09

Hi all I went through the time to optimize_wxc and optimize_wxu to help pull out some density for ligands (i.e. all sequence assigned residues are built). However, it would be nice if the resulting .def was updated with the optimized wxc and wxu. Where are those values? I didn't think for wxc it was scale=XX in the following lines.. r_work= 0.2167 r_free= 0.2525 bonds= 0.003 angles= 0.70 scale= 7.00 r_work= 0.2097 r_free= 0.2502 bonds= 0.004 angles= 0.96 scale= 8.00 ... etc Refining with the scale= value that phenix ultimately chose didn't work form me (resulted in awful geometry and high r/rfree). Thanks FR --------------------------------------------- Francis Reyes M.Sc. 215 UCB University of Colorado at Boulder gpg --keyserver pgp.mit.edu --recv-keys 67BA8D5D 8AE2 F2F4 90F7 9640 28BC 686F 78FD 6669 67BA 8D5D

2 1

Re: [phenixbb] Let's talk about maps...
by Dr. Mark Mayer 07 Aug '09

07 Aug '09

Hi Pavel, It would be really nice if the Electron density maps section in the online docs for phenix.refine http://www.phenix-online.org/documentation/refinement.htm#anch70 was updated to decribe all the different map options people have been asking about. Its much harder to find this info by digging through the Q&A on the bb. Mark --

4 4

Re: [phenixbb] Autobuild error (FIXED)
by youssef ben ammar 07 Aug '09

07 Aug '09

Thanks for all replies, Autobuild error linked to the long path was solved based on the suggestion of Ben Eisenbraun (thanks Ben). I reinstalled tcsh, but only that did not solve the problem. So, I deleted csh from the path and kept only the newly installed tcsh after that everything is running perfectly. Youssef

1 0

Re: [phenixbb] phenix.refine 1.3 vs. 1.4 map quality ?
by William G. Scott 06 Aug '09

06 Aug '09

Hi Paul et al: The yellow map is filled and the blue map is no-fill from the phenix- generated mtz file: http://sage.ucsc.edu/~wgscott/mystuff/filled_vs_unfilled.pdf So I think this is clearly what is going on. Thanks. Bill

1 0

Segmentation Fault
by Waight, Andrew 06 Aug '09

06 Aug '09

Hello All, I've recently solved a membrane protein structure largely with the help of phenix refine. I also have been using the autobuild programs recently for ligand fitting and whatnot. The program during this time has been running flawlessly on my machine (Ubuntu 9.04 E8400 3.0ghz etc etc) However this weekend Autobuild and AutoMR spontaneously started throwing this error. (boost_adaptbx_segmentation_fault_backtrace+0x1f) [0xb7de2e27] 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. Memtest 86 says my RAM is good. I have reinstalled phenix and am still getting the same error. I can't seem to find any other information on this error (segmentation faults are tricky I know) and am just wondering if anyone else has seen or solved this type of problem, thanks for your help. Andrew Waight Wang Lab Skirball Institute of Biomolecular Medicine NYU school of medicine </PRE> <html> <body> ------------------------------------------------------------<br /> This email message, including any attachments, is for the sole use of the intended recipient(s) and may contain information that is proprietary, confidential, and exempt from disclosure under applicable law. Any unauthorized review, use, disclosure, or distribution is prohibited. If you have received this email in error please notify the sender by return email and delete the original message. Please note, the recipient should check this email and any attachments for the presence of viruses. The organization accepts no liability for any damage caused by any virus transmitted by this email.<br /> ================================= </body> </html> <PRE>

3 3

Re: [phenixbb] Autobuild error
by youssef ben ammar 06 Aug '09

06 Aug '09

Thanks for your quick reply, I tried again to fix the problem as suggested by Thomas but this did not solve the problem. Actually I have two machines running phenix, one with the old version 1.3, and the other with 1.4-3 version. With the old version I have also coot, ccp4, and arp-warp, and all things are running without any problem in this machine. But in the other computer, I don't have arp-warp in it (ccp4 and coot exist also), and everything was running OK until I updated phenix to 1.4-3. But now even if I downgrade it to 1.3, autobuild does not work. If I compare echo $PATH results from both machines, I can notice that the path in the one with "defected" autobuild is too long. I will review again my path and try to fix it, or see if there are other recently-installed programs that interfere with phenix. Thanks Youssef.

2 1

Autobuild error -
by Gino Cingolani 05 Aug '09

05 Aug '09

Dear Youssef and Tom, I had the same problem. I also emailed the Phenix classified and got a similar response from Tom. In reality, Tom suggestion didn't solve my problem. The funny thing is that in my case PHENIX.Autobuild worked fine for months before suddenly giving the 'Word too long' error. Weird!!! Believe it or not, the only thing I had to do to fix the problem was to uninstall the ARP/wARP package (sorry about that Tassos!). I realized that PHENIX.Autobuilt stopped working properly right when I installed ARP/wARP. I guess there may be a conflict between the two programs (could Phenix read some ARP/wARP library in the wrong directory?), or something like that. I'm not sure. I'm not claiming this is a general solution, but it worked in my case. I hope it will help you! Cheers, Gino ****************************************************************************** Gino Cingolani, Ph.D. Associate Professor Thomas Jefferson University Dept. of Biochemistry & Molecular Biology 233 South 10th Street - Room 826 Philadelphia PA 19107 Office (215) 503 4573 Lab (215) 503 4595 Fax (215) 923 2117 E-mail: gino.cingolani(a)jefferson.edu ****************************************************************************** "Nati non foste per viver come bruti, ma per seguir virtute e conoscenza" ("You were not born to live like brutes, but to follow virtue and knowledge") Dante, The Divine Comedy (Inferno, XXVI, vv. 119-120)

3 2

Re: [phenixbb] Autobuild error
by Thomas C. Terwilliger 05 Aug '09

05 Aug '09

Hi Youssef, I'm sorry for the trouble! If you update to a recent "nightly build" such as 1.4-129 then this should go away. If it still doesn't, then have a look at your path: echo $PATH which is probably quite long....and identify a few things (perhaps ccp4 and coot for example) that you don't need while running phenix...and then add the keyword: remove_path_word_list="coot ccp4" All the best, Tom T remove_path_word_list = None .help = "List of words identifying paths to remove from PATH " "These can be used to shorten your PATH. For example..." "cns ccp4 coot would remove all paths containing these " "words except those also containing phenix. Capitalization" "is ignored." ">> Dear all, >> >> Recently I have problem with Autobuild. The process fails with the >> following >> message: >> Failed to carry out AutoBuild_build_cycle: >> child process stderr output: >> command: 'csh /home/phenix/p21/AutoBuild_run_1_/TEMP0/RUN_FILE_1' >> Word too long. >> Even if I set the path only to things related to phenix, the process fails >> with the same error. >> >> Any idea? >> >> Youssef >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://www.phenix-online.org/mailman/listinfo/phenixbb >>

1 0

Autobuild error
by youssef ben ammar 05 Aug '09

05 Aug '09

Dear all, Recently I have problem with Autobuild. The process fails with the following message: Failed to carry out AutoBuild_build_cycle: child process stderr output: command: 'csh /home/phenix/p21/AutoBuild_run_1_/TEMP0/RUN_FILE_1' Word too long. Even if I set the path only to things related to phenix, the process fails with the same error. Any idea? Youssef

1 0

Using phenix.pdbtools to strip multiple conformations out of a PDB file?
by Stephen Graham 04 Aug '09

04 Aug '09

Dear all, Apologies if this is a dumb question, but I couldn't find the answer in the docs or by asking the googlebrain. Is there a simple way in phenix.pdbtools to strip "alternate conformations" of side chains, etc, from an input PDB file? By removing "alternate conformations" I mean I would like to keep the conformation with the highest occupancy (not necessarily denoted 'A' in the altLoc column) unless all conformations have an equal occupancy, in which case I'd arbitrarily keep conformation 'A'. I think it should be possible by a clever atom_selection command, but the correct syntax has thus-far eluded me! Thanks, Stephen -- Dr Stephen Graham Division of Structural Biology Wellcome Trust Centre for Human Genetics Roosevelt Drive Oxford OX3 7BN United Kingdom Phone: +44 1865 287 549

2 1

Is the twin refine conflict with order_solvent refine
by shi jiahai 04 Aug '09

04 Aug '09

Hi, Dear all, I am refine a twin dataset with phenix 1.4-3 and CCI (2008_02_08_1610) on a CenOS 5.2 workstation. When I run refine with twin law or order_solvent alone, the refinement is good. Everything goes smoothly. However if I enable refinement on both twin law and order_solvent, the error comes during the order_solvent refinement as below. Start model: number = 55 b_iso_min = 5.52 (limit = 1.00) b_iso_max = 60.32 (limit = 60.00) b_iso_mean = 24.19 anisotropy_min = 1.00 (limit = 0.10) occupancy_min = 1.00 (limit = 0.10) occupancy_max = 1.00 (limit = 1.20) dist_sol_mol_min = 2.13 (limit = 1.80) dist_sol_mol_max = 4.38 (limit = 6.00) Traceback (most recent call last): File "/Application/CCI/cci_apps_sources/phenix/phenix/command_line/refine.py", line 5, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/Application/CCI/cci_apps_sources/phenix/phenix/refinement/command_line.py", line 80, in run call_back_after_monitor_collect=call_back_after_monitor_collect) File "/Application/CCI/cci_apps_sources/phenix/phenix/refinement/driver.py", line 1162, in run call_back_after_monitor_collect = call_back_after_monitor_collect) File "/Application/CCI/cci_apps_sources/phenix/phenix/refinement/strategies.py", line 350, in refinement_machine log = log) File "/Application/CCI/cci_apps_sources/mmtbx/mmtbx/solvent/ordered_solvent.py", line 132, in __init__ map_cutoff = self.params.primary_map_cutoff).peaks_mapped() File "/Application/CCI/cci_apps_sources/mmtbx/mmtbx/solvent/ordered_solvent.py", line 324, in find_peaks log = self.log) File "/Application/CCI/cci_apps_sources/mmtbx/mmtbx/find_peaks.py", line 65, in __init__ symmetry_flags = maptbx.use_space_group_symmetry) File "/Application/CCI/cci_apps_sources/mmtbx/mmtbx/twinning/twin_f_model.py", line 2238, in electron_density_map w2 = w2).fft_map( File "/Application/CCI/cci_apps_sources/mmtbx/mmtbx/twinning/twin_f_model.py", line 2126, in map_coefficients "m_gradient" AssertionError Can anybody help me with this? Thanks a lot in advance! Best Regards, Shi Jiahai Department of Biological Sciences, National University of Singapore

2 1

Assertion error while adding water
by Bishal Kumar Singh 03 Aug '09

03 Aug '09

Hi, Following error is showing while adding water. Every command-line is as default except "ordered_solvent = True". I will appreciate your suggestions. "0.7 1.17367178905 0.7 1.2044634853 0.7 1.18455476555 0.7 1.21102571346 0.7 1.20890773595 0.7 1.23028004675 0.7 1.2017635385 0.7 1.2151128962 0.7 1.20672093638 0.7 1.17881942233 0.7 1.21258440895 0.7 1.21100112067 0.7 1.22101928156 0.7 1.20955876073 0.7 1.17722437184 0.7 1.19891104532 0.7 1.18433797592 Traceback (most recent call last): File "/wrkhuge/cci_apps/cci_apps_sources/phenix/phenix/command_line/refine.py", line 16, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/wrkhuge/cci_apps/cci_apps_sources/phenix/phenix/refinement/command_line.py", line 89, in run call_back_after_monitor_collect=call_back_after_monitor_collect) File "/wrkhuge/cci_apps/cci_apps_sources/phenix/phenix/refinement/driver.py", line 1121, in run call_back_after_monitor_collect = call_back_after_monitor_collect) File "/wrkhuge/cci_apps/cci_apps_sources/phenix/phenix/refinement/strategies.py", line 394, in refinement_machine log = log) File "/wrkhuge/cci_apps/cci_apps_sources/mmtbx/solvent/ordered_solvent.py", line 159, in __init__ use_kick_map = self.params.use_kick_maps).peaks_mapped() File "/wrkhuge/cci_apps/cci_apps_sources/mmtbx/solvent/ordered_solvent.py", line 334, in find_peaks log = self.log) File "/wrkhuge/cci_apps/cci_apps_sources/mmtbx/find_peaks.py", line 93, in __init__ crystal_gridding_tags = fft_map.tags() File "/wrkhuge/cci_apps/cci_apps_sources/cctbx/maptbx/__init__.py", line 152, in tags return crystal_gridding_tags(self) File "/wrkhuge/cci_apps/cci_apps_sources/cctbx/maptbx/__init__.py", line 158, in __init__ assert gridding.symmetry_flags() is not None" regards, Bishal Singh

3 2

Autosol error
by Francis E Reyes 31 Jul '09

31 Jul '09

All When issuing phenix.autosol native.data=scala_072809B_3_001.mtz deriv.data=scala_072909A_3_001.mtz deriv.atom_type=Ir deriv.sites=5 deriv.inano=noinano native.labels='IMEAN SIGIMEAN' deriv.labels='IMEAN SIGIMEAN' I get the error **************************************** AutoSol Input failed For SIR datasets please enter two datafiles or datasets.(native +deriv) Please treat MIR as a series of SIR datasets ...to specify a particular set of data you can specify one of the following (command-line format shown, for a script take out the '=' and use the keyword 'group_labels_list' instead): native.labels='IMEAN SIGIMEAN' deriv.labels='IMEAN SIGIMEAN' native.labels='I(+) SIGI(+) I(-) SIGI(-)' deriv.labels='I(+) SIGI(+) I(-) SIGI(-)' **************************************** Any ideas? thanks FR --------------------------------------------- Francis Reyes M.Sc. 215 UCB University of Colorado at Boulder gpg --keyserver pgp.mit.edu --recv-keys 67BA8D5D 8AE2 F2F4 90F7 9640 28BC 686F 78FD 6669 67BA 8D5D

2 2

combining models from multi-start simulated annealing
by James Whittle 30 Jul '09

30 Jul '09

Is there a way to use AutoBuild to combine a set of models created by multi-start simulated annealing? I've used a random number generator to seed different SA runs at two temperatures, as suggested on this bulletin board before ... but can combine the best parts of these models automatically? --James

2 1

Integrating Phenix with Torque
by Jeremy Mann 29 Jul '09

29 Jul '09

A new user of mine wants to run Phenix on two of our clusters. We use Torque/Maui as our batch scheduler and would like to integrated their usage of Phenix. They want to run something like: phenix.autobuild data=data.hkl seq_file=seq.dat model=model.pdb multiple_models=True multiple_models_number=20 n_cycle_rebuild_max=5 nbatch=20 nproc=20 The documentation mentions briefly the ability to use batch schedulers but doesn't give examples of submission scripts. -- Jeremy Mann jeremy(a)biochem.uthscsa.edu University of Texas Health Science Center Bioinformatics Core Facility http://www.bioinformatics.uthscsa.edu Phone: (210) 567-2672

2 1

Phaser RMS - What happens when it is too high
by Peter Grey 29 Jul '09

29 Jul '09

Hi all, In Phaser RMS is a crucial parameter. I had several cases where underestimation of RMS resulted in negative LLG. The cure was simply increasing the RMS. But I wonder, are there negative effects to a too high RMS. What are the dangers of giving Phaser a too high RMS ? Are there any indications in the output for that. This is also important with EM maps as models since in that case estimation of RMS is difficult. Many thanks in advance for your advice and insight, Peter.

2 1

phenix.refine issue
by Dennis Wolan 27 Jul '09

27 Jul '09

I am having an issue with running a simple phenix.refine that I have never encountered before. I get a list of issues and the program quits with: Traceback (most recent call last): File "/usr/local/phenix-1.4-3/phenix/phenix/command_line/ refine.py", line 9, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/usr/local/phenix-1.4-3/phenix/phenix/refinement/ command_line.py", line 68, in run log=log) File "/usr/local/phenix-1.4-3/phenix/phenix/refinement/ command_line.py", line 231, in __init__ reflection_file_server = self.reflection_file_server() File "/usr/local/phenix-1.4-3/phenix/phenix/refinement/ command_line.py", line 566, in reflection_file_server err=self.log) File "/usr/local/phenix-1.4-3/iotbx/reflection_file_utils.py", line 393, in __init__ force_symmetry=self.force_symmetry)) File "/usr/local/phenix-1.4-3/iotbx/reflection_file_reader.py", line 197, in as_miller_arrays base_array_info=base_array_info) File "/usr/local/phenix-1.4-3/iotbx/mtz/__init__.py", line 434, in as_miller_arrays and flex.min(column_group.sigmas()) > 0): RuntimeError: min() argument is an empty array I have no idea if the issue is with my pdb or mtz. I'm sure it's something very simple that I just am not identifying. Thanks in advance, Dennis

2 1

phaser MR
by sbiswas2＠ncsu.edu 27 Jul '09

27 Jul '09

Hi all, I was trying to solve a structure by molecular replacement. I scaled the data in P4222 space group (resolution 2.7A) with two molecules in the assymmetric unit (molecule A and B) I ran phaser with my model and got a Zscore of 5.1. When I look at the map that I got from phaser I could easily see good electron density for both molecules, However upon inspection of the electron density map there were considerable interaction or clashes with molecule B and a symmetry atom. Molecule A had no clashes however with the symmetry atoms. I was wondering if anyone knows how to resolve this. Could it be a problem of space group. The statistics are good for space group P4222 and the I/sigI was good till 2.7A. Any advice is appreciated, Shya

3 6

question about SAD+MR
by crystallogrphy 24 Jul '09

24 Jul '09

Dear all, I have a Se-Met SAD data which diffracted to 2.8A (peak wavelength). The space group is P1. The map from SAD autosol and autobuild still gives crappy map, I can only see one continous alpha-helix density with all other density fragmented and difficult to identify. I suppose it is because of high mosaicity (1.5), low data redundency, and some radiation damage. Now I want to use mocular replace (sequence identify ~0.26) plus SAD data to give better map. 1. The result from MR is shown below and the density is poor as the sole SAD map: REMARK Log-Likelihood Gain: 121.374 REMARK RFZ=3.4 TFZ=100.0 PAK=0 LLG=29 RFZ=3.7 TFZ=5.4 PAK=0 LLG=69 RFZ=6.1 TFZ=4.0 PAK=0 LLG=97 RFZ=5.6 TFZ=5.0 PAK=0 LLG=118 LLG=121 2. I then tried to combine MR model and SAD data set phenix.autosol seq_file=zntb-sd.seq sites=14 atom_type=Se data=output-peak.sca sg="P1" cell="44.355 56.687 58.451 61.061 82.013 82.740" resolution=3.0 input_partpdb_file=AutoMR_run_1_/MR.1.pdb The program failed to run with such error message: Failed to carry out AutoSol_score_ha_solutions: [Errno 2] No such file or directory: 'AutoSol_run_3_/TEMP0/phaser_mtz.hand.mtz' I do not know where to get this phaser.mtz.hand.mtz file or how to generate it. Does anyone know? -- Qun Wan Department of Pharmacology, Case Western Reserve University 10900 Euclid Avenue, Cleveland, OH 44106 email: qunwan1(a)gmail.com lab phone: 216-368-3337

2 1

new gui on osx (10.5)
by Miguel Ortiz Lombardia 24 Jul '09

24 Jul '09

Dear all, We would like to run the new phenix GUIs from a remote computer. Both server and clients are OSX 10.5 (leopard) and they have the same versions of software (X11, python, phenix... all them pretty new versions) When trying to do so, we invariably get this error: This program needs access to the screen. Please run with 'pythonw', not 'python', and only when you are logged in on the main display of your Mac. I have searched the internet for similar problems and I found an old (one year) unresolved question on the same/similar issue: http://www.phenix-online.org/pipermail/phenixbb/2008-July/001129.html In our case, the /usr/local directory is not mounted over the NFS, it is local. As in that case, though, changing the LIBTBX_PYEXE_BASENAME to pythonw (for example in the phenix.autosol_gui wrapper) did not solve the issue. However, it is not clear to me that this variable has to be changed in a single place: it is found in many files. My question is: was this issue solved and I missed that? If not, anyone has a clue on how to solve the problem? Just to make it clear, the GUIs work perfectly locally in the server computer, that is, if not exported. Thanks a lot! Miguel -- Miguel Ortiz Lombardía Architecture et Fonction des Macromolécules Biologiques (UMR6098) CNRS, Universités d'Aix-Marseille I & II Case 932 163 Avenue de Luminy 13288 Marseille cedex 9 France Tel : +33(0) 491 82 55 93 Fax: +33(0) 491 26 67 20 e-mail: miguel.ortiz-lombardia(a)afmb.univ-mrs.fr Web: http://www.pangea.org/mol/spip.php?rubrique2

2 2

sigma cutoff
by dziuba＠amu.edu.pl 24 Jul '09

24 Jul '09

Hi, According to the Phenix manual, phenix.refine uses no sigma cutoff by default. Since I haven't changed this manually, why do I see in my pdb file MIN(FOBS/SIGMA_FOBS) 1.340? Thanks for help, Karolina

2 1

ENSEmble parameters when the model is an EM map
by Peter Grey 24 Jul '09

24 Jul '09

Dear benevolent experts, I am trying to use Phaser for molecular replacement with EM-map as a model. The EM map, originally in a small P1 cell, was put in large P1 using CCP4's MAPROT and a mask that was defined by a certain cut-off level. Could you please advise me how to derive the EXTEnt and CENTre parameters needed for ENSEmble building ? Does the simple way of giving the size of the original P1 as the EXTEnt and its centre as CENTre sound reasonable ? I would be grateful for any advice and comments, Peter.

2 3

anisotropic ADP factor refinement and optimize_wxc
by Joyce, Gordon M. (NIH/NIAID) [F] 22 Jul '09

22 Jul '09

Dear All, I have just recently started using phenix.refine and I hope I can get some advice. I am trying to refine my structure with anisotropic ADP factors while also optimizing the weighting parameters.When I give the following command, anisotropic refinement is carried out but optimization is not carried out as far as I can see (ie the rms bond and angles are not reduced, which is one of my concerns at the moment). phenix.refine protein1.pdb 072109checks.mtz ordered_solvent=true strategy=individual_adp adp.individual.anisotropic="not element H" optimize_wxc=true --overwrite pdb output values REMARK Start: r_work = 0.1952 r_free = 0.2388 bonds = 0.015 angles = 1.646 REMARK Final: r_work = 0.1413 r_free = 0.2448 bonds = 0.015 angles = 1.646 If I give the following command, optimization occurs and rms bond and angle values decrease. phenix.refine protein1.pdb 072109checks.mtz optimize_wxc=true --overwrite REMARK Start: r_work = 0.2048 r_free = 0.2388 bonds = 0.015 angles = 1.646 REMARK Final: r_work = 0.1917 r_free = 0.2352 bonds = 0.009 angles = 1.301 Is it possible to refine both anisotropic B factors and optimize ADP factors in the same run. Any advice would be great. Thanks, Gordon

2 2

Phenix.refine coordinate error estimate discrepancies
by Phil Jeffrey 21 Jul '09

21 Jul '09

There seems to be a large discrepancy between the likelihood-based coordinate error estimate depending on whether the intensities or structure factor magnitudes are used, which doesn't make any sense to me. If I do: % phenix.refine rebuild.pdb data_withI.mtz phenix.refine picks up Imean and SigImean from the MTZ and refines the structure xray_data { file_name = "data_withI.mtz" labels = "IMEAN,SIGIMEAN" r_free_flags { file_name = "data_withI.mtz" label = "FreeR_flag" test_flag_value = 0 } then I get: | maximum likelihood estimate for coordinate error: 1.39 A But if I create an MTZ file without the intensities by using CAD on the original file, and use the same coordinate file: % phenix.refine rebuild.pdb data_withoutI.mtz phenix.refine picks up the structure factors instead: xray_data { file_name = "data_withoutI.mtz" labels = "F,SIGF" r_free_flags { file_name = "data_withoutI.mtz" label = "FreeR_flag" test_flag_value = 0 } } then this seems to be the right ballpark for a 1.75 Angstrom structure.... | maximum likelihood estimate for coordinate error: 0.28 A Similar R-factors, but there is a factor of 10 difference in the scale reported here: I's: r_work= 0.2266 r_free= 0.2570 ksol= 0.41 Bsol= 50.00 scale= 0.097 F's: r_work= 0.2221 r_free= 0.2557 ksol= 0.41 Bsol= 50.00 scale= 0.965 Version: 1.4 Release tag: 3 cctbx tag: 2008_12_07_1353 Platform: mac-intel-osx osx-10.5.6 User: phil Cheers, Phil Jeffrey Princeton

2 1

Assertion Error when using kicked maps in solvent update
by Schubert, Carsten [PRDUS] 21 Jul '09

21 Jul '09

Hi In 1.4.106 I get the following Assertion Error when I use kicked maps (refinement.ordered_solvent.use_kick_maps=True) in the solvent update section: Traceback (most recent call last): File "/raidnet/xdata1/programs/phenix-1.4-106/phenix/phenix/command_line/refi ne.py", line 19, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/raidnet/xdata1/programs/phenix-1.4-106/phenix/phenix/refinement/comman d_line.py", line 89, in run call_back_handler=call_back_handler) File "/raidnet/xdata1/programs/phenix-1.4-106/phenix/phenix/refinement/driver .py", line 1125, in run call_back_handler = call_back_handler) File "/raidnet/xdata1/programs/phenix-1.4-106/phenix/phenix/refinement/strate gies.py", line 358, in refinement_machine log = log) File "/raidnet/xdata1/programs/phenix-1.4-106/cctbx_project/mmtbx/solvent/ord ered_solvent.py", line 205, in __init__ use_kick_map = self.params.use_kick_maps).peaks_mapped() File "/raidnet/xdata1/programs/phenix-1.4-106/cctbx_project/mmtbx/solvent/ord ered_solvent.py", line 394, in find_peaks log = self.log) File "/raidnet/xdata1/programs/phenix-1.4-106/cctbx_project/mmtbx/find_peaks. py", line 98, in __init__ crystal_gridding_tags = fft_map.tags() File "/raidnet/xdata1/programs/phenix-1.4-106/cctbx_project/cctbx/maptbx/__in it__.py", line 152, in tags return crystal_gridding_tags(self) File "/raidnet/xdata1/programs/phenix-1.4-106/cctbx_project/cctbx/maptbx/__in it__.py", line 158, in __init__ assert gridding.symmetry_flags() is not None AssertionError Running from the same .def file w/o kicked maps works fine. Let me know if you need more info to debug this. Cheers, Carsten

2 1

Duration of a phenix.elbow run
by Dirk 21 Jul '09

21 Jul '09

Hi, I'm trying to create restraints for a carbohydrate ligand using phenix.elbow. If I run it without the --opt command everything works fine, but I'm not completely happy with the geometry of my ligand after refinement. I tried to get an optimized cif by using --opt command, but the program just runs up to the following output: MoleculeClass : H:38 C:23 N: 1 O:19 (PDB format) 81 atoms 83 bonds 151 angles 238 dihedrals -------------------- Time Energy delta E RMS Max Step size and then nothing happens for several hours. Does the optimization needs such a long time and if I let it run overnight (or perhaps longer) I will get an output or is something going wrong? The input is a pdb (I tested several) and I tried both the stable phenix release (1.3) and the current version (1.4). Any help would be very appreciated Thanks Dirk

7 9

Unable to produce kicked maps
by Schubert, Carsten [PRDUS] 17 Jul '09

17 Jul '09

Hi , I am trying to make some kicked maps in 1.4.106. At the end of the run all I am getting is Sorry: Duplicate mtz_label_amplitude FOFCWT Any workaround, except not to use kicked maps? Here the diff_params output for easer debugging. Cheers Carsten <<diff.def>>

2 2

autosol
by sbiswas2＠ncsu.edu 16 Jul '09

16 Jul '09

Hi, I was trying to run autosol using gui however it crashed with the following message: Floating-point error (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 Is this again a version problem? I am running version 1.4-6 Any help will be appreciated, Shya

2 1

rigid body refinement in phenix.refine
by Pavel Afonine 16 Jul '09

16 Jul '09

To PHENIX users who asked about this: I recall several people (unfortunately, I don't remember who exactly) asked me about the details of rigid-body refinement implementation in phenix.refine (MZ protocol) and I promised to let them know when the corresponding paper is published. So here it is: J. Appl. Cryst. (2009). 42, 607-615 Automatic multiple-zone rigid-body refinement with a large convergence radius P. V. Afonine, R. W. Grosse-Kunstleve, A. Urzhumtsev and P. D. Adams http://journals.iucr.org/j/issues/2009/04/00/dd5048/index.html Pavel.

1 0

error message
by zachcp 15 Jul '09

15 Jul '09

A Python error was detected. This is almost certainly a bug; please send an email containing this message to phenixbb(a)phenix-online.org so we can fix it. Thanks! SorryException : No array of observed xray data found. Traceback: Traceback (most recent call last): File "/usr/local/phenix-1.4-3/phenix/wxGUI2/Xtriage.py", line 158, in OnExcept raise event.data SorryException: No array of observed xray data found.

3 2

Average B-factor and ML error
by Leo Sazanov 15 Jul '09

15 Jul '09

Hi, phenix.refine outputs "WILSON B" value into PDB header after refinement, which is not the same as average B-factor of all the atoms in the structure, and not exactly the same as Wilson B-factor from SCALA scaling. How exactly is Wilson B-factor for phenix header calculated and would it be possible to output there also the average B-factor of all the atoms in the structure? Also, deposition into PDB requires only Luzzati coordinate errors, while phenix produces only ML based ones. Presumably ML estimate is better than Luzzati?, but it looks like often people just ignore PDB input for errors. It seems that it would it be useful in future to have ML option in the PDB submission process? Or should Phenix produce Luzzati errors as well? (I know that one can calculate av. B-fac in other programs or ask PDB to add some data to the deposition, but if the idea is to streamline the process...) Thanks, -- Dr. Leonid A. Sazanov Research group leader Medical Research Council Mitochondrial Biology Unit Wellcome Trust / MRC Building Hills Road Cambridge CB2 0XY WEB: www.mrc-mbu.cam.ac.uk Tel: +44-1223-252910 Fax: +44-1223-252915 -- Dr. Leonid A. Sazanov Research group leader Medical Research Council Mitochondrial Biology Unit Wellcome Trust / MRC Building Hills Road Cambridge CB2 0XY WEB: www.mrc-mbu.cam.ac.uk Tel: +44-1223-252910 Fax: +44-1223-252915

2 3

Anyone have a phenix.resolve script for after Phaser ?
by Phil Jeffrey 15 Jul '09

15 Jul '09

Although my usual m.o. is SHELX -> SHARP I'm exploring HYSS -> PHASER -> RESOLVE to see if the maps are any more interpretable. I think I've got good scripts for the first two steps but I'm not sure about my resolve script which stands at: phenix.resolve << EOD hklin myp3_auto.mtz LABIN FP=F SIGFP=SIGF PHIB=PHIB FOM=FOM HLA=HLA HLB=HLB HLC=HLC HLD=HLD hklout resolve.mtz solvent_content 0.45 seq_file myp3.pir EOD This runs, but I could not tell from delving into the documentation if this was ideal or not. Thanks Phil Jeffrey (a Luddite and not using the GUI) Princeton

2 2

Re: [phenixbb] refine occupancy of partially overlapping ligands
by Chris Ulens 14 Jul '09

14 Jul '09

Hi, I'm trying to refine the occupancy of 2 ligands that occupy the same binding pocket in 2 *partially* overlapping orientations. Let's assume I would define the ligands in 2 separate "chains" termed A and B. Could the experts point me to the proper syntax to define the correct group restraints to refine the occupancy? Thank you. Chris +++++++++++++++++++++++++++++++ Chris Ulens, Ph.D. Lab of Structural Neurobiology Division of Pharmacology Campus Gasthuisberg, ON1 Herestraat 49, PB 601 B-3000 Leuven Belgium Tel +32 16 345812 Fax +32 16 345699 Telepathy 2.69 GHz +++++++++++++++++++++++++++++++

2 1

can't find phaser_mtz.hand.mtz?
by Pei-Chun Lin 14 Jul '09

14 Jul '09

Hi, I am the beginner of crystallography and have a 3.5A map generated by MAD method using phenix GUI. I already traced and built 2/3 of my protein. I wish to use this partial model combined with SAD to generate a new map and see if the new map can help me to deal with the rest of the structure. Does it sound reasonable to you? However, after running AutoSol. I got this error message: Failed to carry out AutoSol_score_ha_solutions: [Errno 2] No such file or directory: 'AutoSol_run_1_/TEMP0/phaser_mtz.hand.mtz' What should I do now? Can't phenix automatically calculate it? PC _________________________________________________________________ Windows Live™: Keep your life in sync. Check it out! http://windowslive.com/explore?ocid=TXT_TAGLM_WL_t1_allup_explore_012009

2 1

tightening geometry
by Patrick Loll 10 Jul '09

10 Jul '09

Hi, I'm a little baffled. We have a 3 Å structure, and while the R values seem OK, I'd like to tighten up the rmsd values for the geometry. However, altering wxc_scale seems to have no effect. Specifics: In a clean directory, I run the following command phenix.refine X.mtz Y.pdb Z1.cif Z2.cif Z3.cif wxc_scale=0.??? strategy=group_adp (the cif files are libraries for different ligands) When I vary wxc_scale from 0.5 to 0.005, the geometry doesn't change at all, which makes no sense to me: For wxc_scale=0.5 r_work = 0.2672 r_free = 0.2900 bonds = 0.028 angles = 3.151 For wxc_scale=0.05 r_work = 0.2672 r_free = 0.2900 bonds = 0.028 angles = 3.151 For wxc_scale=0.005 r_work = 0.2672 r_free = 0.2900 bonds = 0.028 angles = 3.151 ...you get the idea. These were run with phenix-1.4-3. With phenix-1.4-56, we get subtly different r values, but the bonds and angles numbers are the same as shown above (to 3 decimal places!), and the numbers don't change upon changing wxc_scale. Am I missing something here? Is there a better way to tighten up the geometry? Thanks, Pat --------------------------------------------------------------------------------------- Patrick J. Loll, Ph. D. Professor of Biochemistry & Molecular Biology Director, Biochemistry Graduate Program Drexel University College of Medicine Room 10-102 New College Building 245 N. 15th St., Mailstop 497 Philadelphia, PA 19102-1192 USA (215) 762-7706 pat.loll(a)drexelmed.edu

3 2

Re: [phenixbb] source of fobs for refinement after AutoSol
by Tom Terwilliger 09 Jul '09

09 Jul '09

Hi James, Yes, that is correct. However if you want...you can just go on to AutoBuild and build a new model from scratch and generate a set of free R flags there, supplying your data as data=mydata.sca. I would recommend this anyway because AutoBuild makes a much better model than AutoSol, and including the model from AutoSol doesn't necessarily help (and can make it worse if the AutoSol model is really incomplete, such as from secondary- structure-only model-building). I use the AutoSol model-building just to see if things are going well and then I throw that model away. This is also the default when you use "after_autosol=True" in AutoBuild and when you run the AutoBuild GUI after running AutoSol. It is convenient to use the data and freeR in in exptl_fobs_phases_freeR_flags_*.mtz from AutoSol but you don't have to. The key of course is: once you refine a model and want to use that model further...you must use the set of free R flags used in that original refinement from then on. If you throw away that model and start from scratch you can also make a new set of free R flags. The hexdigest REMARK in your PDB file from phenix.refine helps you make sure that this is done correctly. All the best, Tom T On Jul 9, 2009, at 8:35 AM, James Whittle wrote: > Dear Tom: > > Ok. Thanks for the quick response! > > I should run then run AutoSol again to get this file, and it will > have the HL coefficients in it? > > --James > > > On Jul 9, 2009, at 10:13 AM, Tom Terwilliger wrote: > >> Hi James, >> >> As you are expecting, in AutoSol the file >> exptl_fobs_phases_freeR_flags_*.mtz normally contains the >> experimental Fobs and free R flags for refinement, along with >> phases and HL coefficients from the experimental phasing. >> >> However if an anisotropy correction is applied to the data, then by >> default no refinement is done in AutoSol and no >> exptl_fobs_phases_freeR_flags_*.mtz file is created. This is to >> ensure that refinement is not carried out against anisotropy- >> corrected data (you want to refine against the original data, and >> have phenix.refine apply an anisotropy correction as part of >> refinement). >> >> If you supply >> >> input_refinement_file=my_data.sca >> >> then my_data.sca will be used for refinement and an >> exptl_fobs_phases_freeR_flags_*.mtz will be created. Note that >> my_data.sca can be identical to your input data file if you want. >> I am guessing that in your case an anisotropy correction was >> applied to the data in this way. >> >> I hope that helps! >> -Tom T >> ps. I am thinking of changing the default to make >> input_refinement_file automatically defined if an anisotropy >> correction is applied so that you don't have to do this... >> >> On Jul 9, 2009, at 7:55 AM, whittle wrote: >> >>> Dear all, >>> >>> I have a question regarding the output of AutoSol. >>> >>> I've run AutoSol, using the data directly from scalepack as a .sca >>> file, >>> and would like to now begin refining the structure I've built. >>> >>> There seems to be no exptl_fobs_phases_freeR_flags_*.mtz file in the >>> AutoSol directory. Isn't this the right file to use? Since it's not >>> there, is it proper to use the phaser_2.mtz file for phenix.refine? >>> >>> --James >>> >>> _______________________________________________ >>> phenixbb mailing list >>> phenixbb(a)phenix-online.org >>> http://www.phenix-online.org/mailman/listinfo/phenixbb >> >> >> Thomas C. Terwilliger >> Mail Stop M888 >> Los Alamos National Laboratory >> Los Alamos, NM 87545 >> >> Tel: 505-667-0072 email: terwilliger(a)LANL.gov >> Fax: 505-665-3024 SOLVE web site: http://solve.lanl.gov >> PHENIX web site: http:www.phenix-online.org >> ISFI Integrated Center for Structure and Function Innovation web >> site: http://techcenter.mbi.ucla.edu >> TB Structural Genomics Consortium web site: http://www.doe-mbi.ucla.edu/TB >> CBSS Center for Bio-Security Science web site: http://www.lanl.gov/cbss >> >> >> >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://www.phenix-online.org/mailman/listinfo/phenixbb > Thomas C. Terwilliger Mail Stop M888 Los Alamos National Laboratory Los Alamos, NM 87545 Tel: 505-667-0072 email: terwilliger(a)LANL.gov Fax: 505-665-3024 SOLVE web site: http://solve.lanl.gov PHENIX web site: http:www.phenix-online.org ISFI Integrated Center for Structure and Function Innovation web site: http://techcenter.mbi.ucla.edu TB Structural Genomics Consortium web site: http://www.doe-mbi.ucla.edu/TB CBSS Center for Bio-Security Science web site: http://www.lanl.gov/cbss

2 3

source of fobs for refinement after AutoSol
by whittle 09 Jul '09

09 Jul '09

Dear all, I have a question regarding the output of AutoSol. I've run AutoSol, using the data directly from scalepack as a .sca file, and would like to now begin refining the structure I've built. There seems to be no exptl_fobs_phases_freeR_flags_*.mtz file in the AutoSol directory. Isn't this the right file to use? Since it's not there, is it proper to use the phaser_2.mtz file for phenix.refine? --James

2 1

2010 Gordon Conference on Diffraction Methods
by Andrew Leslie 02 Jul '09

02 Jul '09

The Gordon Research Conference on Diffraction Methods in Structural Biology will be held at Bates College in Lewiston, Maine, from July 18-23, 2010. As in the past, the program will include the latest developments in methodology covering all aspects of macromolecular crystallography, from crystallisation to model refinement and validation. We are currently preparing a preliminary program for the meeting, and would welcome suggestions for particular topics (with or without speakers). Please Email your suggestions to Andrew Leslie (andrew(a)mrc-lmb.cam.ac.uk) or Ana Gonzalez (ana(a)slac.stanford.edu). Andrew and Ana

1 0