- phenixbb - phenix-online.org

Autobuild
by Heather L Condurso 06 Sep '13

06 Sep '13

I'm a bit confused about what is going on with my autobuild run. I solved a structure from anomalous data to 2.7A but many loops are missing. I ran phaser with the current model and native data that extends to 2.4A and got a solution and started autobuild from there. I used mostly the default settings and rebuild_in_place=False. After a pretty short time I realized the run was over and at the bottom of the log I see this: Cycle 3 R/Rfree=0.28/0.35 Built=488 Placed=467 AutoBuild_build_cycle AutoBuild Run 162 Fri Sep 6 14:33:39 2013 Deciding if we should continue... Recent R change per cycle: 0.0 Ending these build cycles as the recent change in R per cycle has been 0.0 and the required value of change is -0.005 and the R is somewhat acceptable at 0.276670644792 All done with build cycles in this region AutoBuild_set_up_build AutoBuild Run 162 Fri Sep 6 14:33:40 2013 All omit/non-omit regions completed finished AutoBuild Run 162 Fri Sep 6 14:33:40 2013 Finishing AutoBuild Run 162 Facts written to AutoBuild_run_162_/AutoBuild_Facts.dat AutoBuild Run 162 Summary of model-building for run 162 Fri Sep 6 14:33:40 2013 Files are in the directory: /Users/heather/Desktop/M/AutoBuild_run_162_/ Starting mtz file: /Users/heather/Desktop/M/ImportRawData_run_4_/A-3_xx_PHX.mtz Sequence file: /Users/heather/Desktop/M/m.pir Best solution on cycle: 2 R/Rfree=0.26/0.33 Summary of output files for Solution 1 from rebuild cycle 2 --- Model (PDB file) --- pdb_file: /Users/heather/Desktop/M/AutoBuild_run_162_/overall_best.pdb --- Model (CIF file) --- cif_pdb_file: /Users/heather/Desktop/M/AutoBuild_run_162_/overall_best.cif --- Refinement log file --- log_refine: /Users/heather/Desktop/M/AutoBuild_run_162_/overall_best.log_refine --- Model-building log file --- log: /Users/heather/Desktop/M/AutoBuild_run_162_/overall_best.log --- Model-map correlation log file --- log_eval: /Users/heather/Desktop/M/AutoBuild_run_162_/overall_best.log_eval --- 2FoFc and FoFc map coefficients from refinement 2FOFCWT PH2FOFCWT FOFCWT PHFOFCWT --- refine_map_coeffs: /Users/heather/Desktop/M/AutoBuild_run_162_/overall_best_refine_map_coeffs.mtz --- Data for refinement FP SIGFP PHIM FOMM HLAM HLBM HLCM HLDM FreeR_flag --- refine_data: /Users/heather/Desktop/M/AutoBuild_run_162_/overall_best_refine_data.mtz --- Density-modified map coefficients FWT PHWT --- denmod_map_coeffs: /Users/heather/Desktop//AutoBuild_run_162_/overall_best_denmod_map_coeffs.mtz You might consider making one very good model now with: phenix.autobuild \ data=/Users/heather/Desktop/M/AutoBuild_run_162_/overall_best_refine_data.mtz \ model=/Users/heather/Desktop/M/AutoBuild_run_162_/overall_best.pdb \ rebuild_in_place=True \ seq_file=/Users/heather/Desktop/M/m.pir My full sequence is 325 amino acids so only about 75% is build in. I am hoping autobuild can help find at least some of these missing residues. Why would the change in R be zero? What is meant by "You might consider making one very good model now with:"? I tried running autobuild with these new files and rebuild_in_place=true and get the same type of message. Is this the best that autobuild can do? I also ran the phaser solution through refinement and can significantly lower the Rs from 0.26/0.33 to 0.22/0.30. I appreciate any advice on how best to proceed. I am using 1.8.2-1472 if that matters. Thanks, Heather

2 3

Postdoctoral Researcher Opportunity in the application of free-electron lasers to structural biology.
by Edward Snell 06 Sep '13

06 Sep '13

Postdoctoral Researcher in application of free-electron lasers to structural biology. Open from Sept 2013 until filled. Employer: Hauptman Woodward Medical Research Institute (Position is situated in Germany) Outline of Position: This position is supported by a new NSF grant devoted to the application of free-electron X-ray lasers to structural biology. The position will be based in Hamburg, Germany at DESY CFEL, but is supported by a grant at SUNY Buffalo. There will be time for research collaborations in addition to support duties. The ideal candidate will hold a PhD in experimental Biophysics or similar. Major Responsibilities: 1) Application of X-ray free electron laser to structural biology and development of sample-delivery devices. 2) International travel to annual meetings being held in the U.S. 3) Collaboration with partner institutions and labs in the U.S. 4) Interact via web conferencing with the principal investigator and scientific director of the project Please send resume and cover letter to: Jill Szczesek Director of Operations and Academic Affairs Hauptman Woodward Medical Research Institute 700 Ellicott Street Buffalo, New York 14203 phone: 716-898-8597 email: jszczesek(a)hwi.buffalo.edu<mailto:[email protected]> Edward Snell Ph.D. Assistant Prof. Department of Structural Biology, SUNY Buffalo, Senior Scientist, Hauptman-Woodward Medical Research Institute 700 Ellicott Street, Buffalo, NY 14203-1102 Phone: (716) 898 8631 Fax: (716) 898 8660 Skype: eddie.snell Email: esnell(a)hwi.buffalo.edu Telepathy: 42.2 GHz Heisenberg was probably here!

1 0

PDBTools and TLS
by Mark Saper 03 Sep '13

03 Sep '13

Hi, I couldn't find the answers to these questions in the documentation. 1. If I am refining TLS and individual isotropic ADPs, does phenix.refine derive iso-ADP from the ANISOU records presently in the file, or does it use the iso-ADP from the ATOM record and recombine it with the TLS values? 2. Consequently if I reset the B's for a range of residues on the ATOM cards only, will they be used rather than the ANISOU record to derive the ADP? Is it better to use PDBTools, specify the iso-ADP and also convert to ANISOU? I am doing this for only a selection of the residues. 3. Should phenix.refine input files retain the header lines containing the previously cycle TLS parameters? If so PDB Tools doesn't save them. PDB Editor does. Which is the preferred tool? 4. If one explicitly changes atom ADPs with Coot using Measures>Residue info... the output file reflects the change only in the ATOM lines, but the ANISOU line retains the previous 6 parameters. It seems as if coot should also change the ANISOU to reflect anisou (last 3 parameters=0). Comments or answers? Many thanks in advance! -Mark _______________________________________________ Mark A. Saper, Ph.D. Department of Biological Chemistry, University of Michigan Medical School 3040 Chemistry Building | saper(a)umich.edu | +1 (734) 764-3353

2 1

Error determining reference matches
by Dan McNamara 31 Aug '13

31 Aug '13

Hi all, I am encountering a strange crash when running phenix.refine through the command line with reference model restraints turned on. The references point to two different PDB files, one labeled as chain A and one labeled as chain B. This error occurs whether I have my refinement model in P1 with 192 chains or P2(1) with 96 chains. The chains are A-EZ (192) or A-BH (96). It appears to fail to align the reference models to the sequences in the refinement model. This error does not occur with other smaller target structures under 25 chains and reference models used with the same installation of phenix. I am hoping for any advice on why this might be happening or how I might get around it. ================== Extract refinement strategy and selections ================= Refinement flags and selection counts: individual_sites = True (401472 atoms) torsion_angles = False (0 atoms) rigid_body = True (401472 atoms in 192 groups) individual_adp = False (iso = 0 aniso = 0) group_adp = True (201984 atoms in 25824 groups) tls = False (0 atoms in 0 groups) occupancies = False (0 atoms) group_anomalous = False n_use = 401472 n_use_u_iso = 401472 n_use_u_aniso = 0 n_grad_site = 0 n_grad_u_iso = 0 n_grad_u_aniso = 0 n_grad_occupancy = 0 n_grad_fp = 0 n_grad_fdp = 0 total number of scatterers = 401472 *** Adding Reference Model Restraints *** determining reference matches automatically... Traceback (most recent call last): File "/auto_nfs/joule2/programs/phenix/phenix-installer-dev-1457/phenix-dev-1457/build/intel-linux-2.6-x86_64/../../phenix/phenix/command_line/refine.py", line 11, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/auto_nfs/joule2/programs/phenix/phenix-installer-dev-1457/phenix-dev-1457/phenix/phenix/refinement/command_line.py", line 92, in run master_params=customized_master_params) File "/auto_nfs/joule2/programs/phenix/phenix-installer-dev-1457/phenix-dev-1457/phenix/phenix/refinement/driver.py", line 501, in __init__ log=self.log) File "/auto_nfs/joule2/programs/phenix/phenix-installer-dev-1457/phenix-dev-1457/cctbx_project/mmtbx/torsion_restraints/reference_model.py", line 147, in __init__ self.get_reference_dihedral_proxies() File "/auto_nfs/joule2/programs/phenix/phenix-installer-dev-1457/phenix-dev-1457/cctbx_project/mmtbx/torsion_restraints/reference_model.py", line 460, in get_reference_dihedral_proxies log=self.log) File "/auto_nfs/joule2/programs/phenix/phenix-installer-dev-1457/phenix-dev-1457/cctbx_project/mmtbx/torsion_restraints/reference_model.py", line 323, in process_reference_groups moving_chain = mod_h) File "/auto_nfs/joule2/programs/phenix/phenix-installer-dev-1457/phenix-dev-1457/cctbx_project/mmtbx/torsion_restraints/utils.py", line 426, in _ssm_align ssm_alignment = ccp4io_adaptbx.SSMAlignment.residue_groups(match=ssm) File "/auto_nfs/joule2/programs/phenix/phenix-installer-dev-1457/phenix-dev-1457/ccp4io_adaptbx/__init__.py", line 215, in residue_groups return cls( match = match, indexer = indexer ) File "/auto_nfs/joule2/programs/phenix/phenix-installer-dev-1457/phenix-dev-1457/ccp4io_adaptbx/__init__.py", line 176, in __init__ self.pairs.append( ( get( f, indexer1 ), get( s, indexer2 ) ) ) File "/auto_nfs/joule2/programs/phenix/phenix-installer-dev-1457/phenix-dev-1457/ccp4io_adaptbx/__init__.py", line 173, in get assert identifier in indexer, "Id %s missing" % str( identifier ) AssertionError: Id ('A', 1, ' ') missing Best, Dan

3 2

coot and mtz interpretation
by Heather L Condurso 30 Aug '13

30 Aug '13

I'm having an issue that when I launch coot from the phenix refine gui it loads the pdb and mtz fine, but the FoFc map colors are green and royal blue. When I open coot first and then load the pdb and mtz files, the FoFc map colors are green and red, but the blobs that were green in the phenix loaded maps are now red and what was royal blue is now green. Has anyone else had this issue? I am running a mac with OS X 10.6.8 and coot 0.6.1 and phenix 1.8.2-1309. I have an identical mac running all the same versions and everything loads green and red as expected. Is there a parameter file somewhere that tells coot how to display the phenix maps that may provide a clue as to what is going on? I appreciate any suggestions, Heather

3 4

a question on NAG-ASN cif file
by fn1＠rice.edu 29 Aug '13

29 Aug '13

Hi Nigel, I tried your nag_asn.cif file. I did not include hydrogen in the refinement, so i always get the error message that NAG-ASN could not be applied. So i remark the angle definition that contains hydrogen and the same for the plane definition, as follows, #NAG-ASN 1 C1 2 ND2 2 HD2 119.000 3.000 #NAG-ASN plane1 2 HD21 .020 by doing this, the program could run through. However, when i check the geometry output .geo file, only the added plane information is seen, no chirality of C1 in NAG is seen. Did i do something wrong? Thank you! Fengyun

2 1

a question on NAG-ASN cif file
by fn1＠rice.edu 28 Aug '13

28 Aug '13

Hi everyone, I am refining a glycoprotein model, and I have applied NAG-ASN link to the ASN. But the beta-conformation of NAG could not be retained after refinement for some sites. I read the cif file and find that there are only angle definitions for NAG-ASN. The dihedral and charility information are missing. Is there a way that I could retain the charility of C1 of NAG? Thank you! Fengyun

3 3

ALS Call for General User Proposals - DEADLINE SEPTEMBER 4 , 2013
by Banumathi Sankaran 28 Aug '13

28 Aug '13

> > The deadline for Jan/July 2014 Collaborative > Crystallography proposals will be *Sep 4, 2013. * > > ---------------------------------------------------------------------------------------- > 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: > > Please follow the instructions for proposal submission at: > http://www-als.lbl.gov/index.php/user-information/user-guide/58.html > Scroll down to "*Structural Biology beamlines (includes protein SAXS)*" > and click on "New Proposal." Enter your proposal information. > > > > Regards, > Banumathi Sankaran >

1 0

processor usage controlin DEN
by miyatake 24 Aug '13

24 Aug '13

Dear phenixbb, Hello, I run DEN refinements using GUI (dev-1439). I wonder whether I can control the number of the processor used during the DEN refinement. I know that we can use nproc in command line, but I'd like use it in GUI. Best regards, H. Miyatake

2 4

phenix.real_space_refine
by Heather L Condurso 23 Aug '13

23 Aug '13

I tried to run the phenix.real_space_refine that was discussed in the CCN but I get the following error. import mmtbx.refinement.real_space.driver ImportError: No module named driver I downloaded the latest nightly build dev-1463. Even when asking for the help I get the same error. I did however successfully create some FEMs which look lovely. Heather Condurso UF-Department of Chemistry condurso(a)ufl.edu

3 3

Post-doc position for SAXS at the MAX IV Laboratory
by Marjolein Thunnissen 23 Aug '13

23 Aug '13

The MAX IV laboratory in Lund has opened a position for a Postdoctoral fellow to strengthen the biological SAXS environment at the current SAXS beam line at the MAX-II ring. The position will be for two years and is embedded in the BIOSTRUCT-X EU program. Full information of the position is available at: http://admin.lu.se/o.o.i.s?id=22598&Dnr=554436&Type=E. More information can be obtained from Tomás Plivelic (tomas.plivelic(a)maxlab.lu.se<mailto:[email protected]>). Information about the MAX IV laboratory can be obtained at the homepage www.maxlab.lu.se<http://www.maxlab.lu.se> . best regards Marjolein Thunnissen ________________________________________________________________________ Marjolein Thunnissen Phone +46-(0)76-632 0417 Associate Professor Fax +46-(0)46-22 24692 Dept of Biochemistry and Structural Biology, Lund University http://www.mps.lu.se<http://www.mps.lu.se/> PO-Box 124 S-221 00 Lund, Sweden Scientific coordinator MX (The MAX IV Laboratory): I911 and BioMAX ________________________________________________________________________

1 0

Find peaks and holes documentation
by Morten Grøftehauge 22 Aug '13

22 Aug '13

Hi, I wanted to try phenix.find_peaks_holes with anomalous structure factors from XPREP but I ran into a few problems. And then I tried to fix them by reading the manual and that was really hard too. So first the missing manual. If I google it I get this page http://www.phenix-online.org/documentation/find_peaks_holes.htm which is for the GUI. When I go back to the documentation http://www.phenix-online.org/documentation/phenix_documentation.html I can't find the link for it (which is also why I was using Google to find it in the first place). I also checked the source to see if the link was somehow misnamed. Anyway, running it from the command line with default settings is not hard. phenix.find_peaks_holes model.pdb reflections.hkl but I run into this error message Sorry: Unresolved amplitude/intensity ambiguity: /path/path/path/reflections.hkl SHELX reflection files may contain amplitudes or intensities. Please append =amplitudes or =hklf3 or =intensities or =hklf4 to the file name argument or parameter to resolve the ambiguity. and I don't how to do that. If I add .hklf3 at the end of the file name I get the same error. If I go reflections.hkl=hklf3 I just get a message that I need to supply a reflection file. What is the syntax? Cheers, Morten -- Morten K Grøftehauge, PhD Pohl Group Durham University

2 4

coot add residue issue
by Heather L Condurso 21 Aug '13

21 Aug '13

Dear all, I have been trying to add a residue to my model in coot for hours now. The add terminal residue function adds a residue far from the end of my chain and I've tried moving it manually and it keeps going back and it never makes a bond between the two amino acids. Is there a trick in coot or a better way to do this in phenix? Thanks, Heather Condurso UF-Department of Chemistry

5 5

peptide ligand refinement problems
by Mengbin Chen 20 Aug '13

20 Aug '13

Dear Phenix Users, My colleague has got into trouble when he was refining a peptide ligand with a coumarin ring attached via a peptide bond. The peptide part is composed of natural amino acids, which does not require a cif file. On the other hand, a cif file has been created for the coumarin ring. Unfortuantely, after the refinement run, the peptide and the coumarin ring fell apart, and we could not figure out why this is happening. We are using the Phenix-dev-1370 version. Any help would be appreciated! Thank you in advance, Mengbin -- Mengbin Chen Department of Chemistry University of Pennsylvania

2 4

Re: [phenixbb] HKL2000 behaved weirdly
by Mengbin Chen 20 Aug '13

20 Aug '13

Sorry for being late in replying, but thank you all for your suggestions! The resolution is ~3 angstroms, so it was quite easy to pick more than 100 peaks for indexing. So I tried to index ~5-6 images at the same time as you suggested, the Bravais Lattice window appeared (and I reset the beam center too), but the unit cell dimensions were ~10 angstroms along one axis, which was not appropriate since the protein is around 40 kDa. I tried to manually pick peaks for indexing, still cell dimensions were too small to be properly assigned. I feel the crystal has got severe problems of multiple lattices and high mosaicity at the same time, which made it hard to index. I am still playing around with the dataset to see if I can get any luck with it, so I'll let you know if I can work this out. Mengbin On Fri, Aug 16, 2013 at 12:18 PM, Noinaj, Nicholas (NIH/NIDDK) [F] < noinajn(a)niddk.nih.gov> wrote: > Mengbin, > > > This happens in HKL2000 with tough datasets. The data can sometime appear > to be very nice, possibly some high mosaicity, but still not process > properly. I have seen this many times. As mentioned already, the next > step is to start trying different things to get the data to process. > WIthout knowing more about resolution, etc, here is a few things to try. > > 1- Check and the double check the beam center. Then, check it again! > With troublesome datasets, this is EXTREMELY important. In HKL2000, > everytime you hit Abort Refinement, you have to reset the beamcenter again. > SO monitor the beam center and everytime it gets too far away from the > actual numbers, reset it. > > 2- you want to try to have roughly 100 peaks for indexing at minimum (if > possible), so change the resolution limits to 15 ang - 4 ang, then retry > indexing. IF this does not work, try 15-5ang, then 12-5, then 10-5, etc. > > 3- If #2 doesn't work, then you can alway try indexing using multiple > frames, esp if you don't have a enough peaks on the first frame. To do > this, please see the HKL2000 manual, or let me know, i can send you some > screenshots or give you a call to guide you via phone. here, you can use > 5-10 frames for indexing if needed. > > 4- as a last option, you can always manually select your peaks. i have had > this work nicely in the past, but again, last option. > > 5- you could try other data processing programs too such as MoSFLM, XDS, > Xia2, etc. > > If you still don't get it working, i am more than happy to take a quick > look. Just send me one of your frames to analyse. Of course, everything > will be strictly confidential. Thanks in advance and good luck! > > > > > Cheers, > Nick > > > > ________________________________ > From: Mengbin Chen [mengbinc(a)sas.upenn.edu] > Sent: Thursday, August 15, 2013 3:12 PM > To: PHENIX user mailing list > Subject: [phenixbb] HKL2000 behaved weirdly > > Dear All, > > I used HKL2000 to index my data, but weird things happened: after I hit > peak search, the data were selected to index, and I chose primitive > triclinic to continue. After I hit refine, the circles for peak search > disappeared, and if I hit "Bravais Lattice", the table would not show up, > which meant that I could not proceed to refine in the potentially correct > Bravais Lattice. I don't know if anyone else has encountered this problem > before, so any suggestions would be greatly appreciated! > > Thank you in advance, > Mengbin > > -- > Mengbin Chen > Department of Chemistry > University of Pennsylvania > -- Mengbin Chen Department of Chemistry University of Pennsylvania

1 0

Preferred treatment of residues with no electron density
by Mark Saper 20 Aug '13

20 Aug '13

What is the preferred method for handling missing residues or residues with zero occupancy, in terms of PBD submission? I see three options: Remove them from the model set occupancy=0 and flag in REMARKs 475 and 480 keep occupancy=1 and allow the temperature factors to climb I have a disordered region in my molecule which I am restraining to a portion of a structure in which that region is well resolved. There are several residues that have no density in the current refinement, yet I want to keep them in the model. Suggestions? Mark _______________________________________________ Mark A. Saper, Ph.D. Department of Biological Chemistry, University of Michigan Medical School 3040 Chemistry Building | saper(a)umich.edu | +1 (734) 764-3353

3 2

a question on LSTc refinement
by fn1＠rice.edu 19 Aug '13

19 Aug '13

Hi all, I read some previous post on LSTc [SIA(alpha2-6)GAL(beta1-4)NAG] refinement, but i still cannot figure out how to refine it in PHENIX. In the current version of phenix, i still could not find the ALPHA2-6 link. Does anybody have such experience? Another question is that, sialic acid (SIA) is in the monomer library of phenix. But when i run on the pdb directly, it will give the error information like, "unknown bond type for O1A and O1B". Though this error could be solved by provided the cif file generated by phenix.elbow. Thank you very much! Fengyun

3 2

Ensemble refinement, starting water structure and alternative conformations
by Joseph Noel 19 Aug '13

19 Aug '13

Is it important (or at least useful) to remove any waters before initiating ensemble refinement wherein the waters are located in what is obvious electron density for alternative conformations of side chains? Would the presence of such waters possibly result in ensemble refinement NOT modeling the alternative conformations? ______________________________________________________________________________________ Joseph P. Noel, Ph.D. Arthur and Julie Woodrow Chair Investigator, Howard Hughes Medical Institute Professor, The Jack H. Skirball Center for Chemical Biology and Proteomics The Salk Institute for Biological Studies 10010 North Torrey Pines Road La Jolla, CA 92037 USA Phone: (858) 453-4100 extension 1442 Cell: (858) 349-4700 Fax: (858) 597-0855 E-mail: noel(a)salk.edu Publications & Citations: http://scholar.google.com/citations?user=xiL1lscAAAAJ Homepage Salk: http://www.salk.edu/faculty/noel.html Homepage HHMI: http://hhmi.org/research/investigators/noel.html ______________________________________________________________________________________

2 1

Density modification using partial model phase
by Yi xing 17 Aug '13

17 Aug '13

Does mtz file from phenix.refine report FOM? If not, is there a way to have FOM listed in the mtz file. I am trying to use the mtz from the phenix refined model as the input of the ccp4 density modification (DM) or Dmmulti. For example, mtz file from Phaser molecular replacement has all the columns to run the DM. Many thanks for any help or comments! Best, yi

3 2

HKL2000 behaved weirdly
by Mengbin Chen 15 Aug '13

15 Aug '13

Dear All, I used HKL2000 to index my data, but weird things happened: after I hit peak search, the data were selected to index, and I chose primitive triclinic to continue. After I hit refine, the circles for peak search disappeared, and if I hit "Bravais Lattice", the table would not show up, which meant that I could not proceed to refine in the potentially correct Bravais Lattice. I don't know if anyone else has encountered this problem before, so any suggestions would be greatly appreciated! Thank you in advance, Mengbin -- Mengbin Chen Department of Chemistry University of Pennsylvania

2 1

problem with refinement in AutoBuild
by Michael Thompson 15 Aug '13

15 Aug '13

Hello Everyone, I am having a problem running autobuild, which appears to be happening during the refinement step. I start autobuild using : phenix.autobuild data=data.mtz model=model.mtz seq_file=sequence.dat rebuild_from_fragments=True i_ran_seed=124889 find_ncs=False refine_with_ncs=False nproc=4 The runs are failing during the first cycle with the following error message: AutoBuild_build_cycle AutoBuild Run 2 Thu Aug 15 13:37:34 2013 Build cycle 1 of 16 method:rebuild Initial rebuild cycle: just combining previous models or copying over initial PDB file Using NO_LIG_NO_WAT_CloMP_AutoBuild_input_ed.pdb for rebuild, copied to unrefined.pdb_1 Refining model: unrefined.pdb_1 Refinement strategy: ['individual_sites', 'individual_adp', 'occupancies'] No refinement carried out (may need free R flags) Error message from phenix.refine: Refinement failed... Sorry: ******************************************************************************** Failed to carry out AutoBuild_build_cycle: Initial phenix refinement of AutoBuild_run_2_/TEMP0/NO_LIG_NO_WAT_CloMP_AutoBuild_input_ed.pdb with aniso_refinement_data_PHX.mtz failed...quitting... /home/miket/projects/CloMP/Oxidized/APSoct2012/CloMP-TMAOcryo_B2_3/MRwithChainsAB/RefineCChalf/AutobuildPrep/AutoBuild_run_2_/LAST.LOG may have additional information... ******************************************************************************** I also checked the LAST.LOG file that is referenced in the error message. That appears to be a phenix.refine logfile, however it just cuts off after the "refinement start" header, with no information about the cause of the failure (no error message). My input MTZ file does indeed contain R-free flags. I have used this MTZ for refinement with phenix.refine (outside of autobuild) with no problems. Furthermore, earlier in the autobuild logfile, I see: Copying refinement_PHX.mtz to exptl_fobs_phases_freeR_flags.mtz Columns used: FP=FP SIGFP=SIGFP FreeR_flag=FreeRflag Which seems to me like it is recognizing my R-free flags column in the input MTZ. Does anyone have any idea what might be going wrong here?? Thanks, Mike -- Michael C. Thompson Graduate Student Biochemistry & Molecular Biology Division Department of Chemistry & Biochemistry University of California, Los Angeles miket(a)chem.ucla.edu

2 1

Phenix.refine suddenly stopped to refine B-factors
by krimmer＠staff.uni-marburg.de 12 Aug '13

12 Aug '13

Dear all, Phenix.refine suddenly stopped to refine B-factors, even though according to the settings it should refine them (checked box “individual B-factors”). This problem only occurs for one job. The weird thing is that this was always working and I did not do anything different, at least as far as I know. Also, this just happened to a colleague, totally independent from me. We never had this problem before. Here is a part of my log file: |-----overall-----|---macromolecule----|------solvent-------| stage b_max b_min b_ave b_max b_min b_ave b_max b_min b_ave 0 : 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 1_bss: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 1_ohs: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 1_xyz: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 1_adp: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 1_occ: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 2_bss: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 2_ohs: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 2_xyz: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 2_adp: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 2_occ: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 3_bss: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 3_ohs: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 3_xyz: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 3_adp: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 3_occ: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 4_bss: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 4_ohs: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 4_xyz: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 4_adp: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 4_occ: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 5_bss: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 5_ohs: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 5_xyz: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 5_adp: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 5_occ: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 5_bss: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 5_ohs: 46.46 5.37 12.45 29.68 5.37 10.14 46.46 5.58 25.33 The values are not changing at all. The structure is almost done, but it still should change a bit during the refinement cycles. I noticed the problem because I added another water molecule to the structure, and the B-factor stayed at 30.00 which resulted in a negative density blob after refinement. Thank you for your help! Stefan

3 3

How to correct translational psuedosymmetry
by Appu kumar 12 Aug '13

12 Aug '13

Dear all, I have processed a protein-DNA complex data with imosflm. Phenix.xtriage suggest that the is a off origin peak of hight 41.702 hight at fractional coordiantes of 0.500 0.021 0.500 with distance of 40.019 from origin. I am new to such kind of problem, so would you please suggest me how to solve the psuedosymmetry problem encountered in data process. The present unit cells are 55.41 109.72 66.61 90 98.74 90 in space group P1211 with p_value is 3.116e-04. Your help would be much appreciated Thank you Appu

2 1

is phenix trying to polymerize b-octyl glucoside?
by Edward A. Berry 09 Aug '13

09 Aug '13

I'm refining a model with some molecules of the detergent b-octyl glucoside (BOG). I got BOG.cif from a recent ccp4 library, so not sure if this is a suggestion for phenix or ccp4, but: If two BOG molecules are consecutive in the pdb file, phenix stops in setup because it can't find atom HO1 to delete. Even if I had H's in my model there would be no HO1 because O1 is in glycosidic linkage with the octyl chain. I'm guessing phenix somehow realizes there is a sugar and when two are consecutive it autogenerates a link to make a glycosidic bond between them - not appropriate in the case of BOG. easy work around is to intersperse with non-sugar molecules, but if there is a better way, I'll take it. ed RuntimeError: delete_atom_in_place: unknown atom_id: HO1 apply_mod failure: pdbres="BOG R 18 " segid="R " comp id: BOG mod id: DEL-O1 Version: dev Release tag: 1420 Platform: intel-linux-2.6-x86_64 linux

2 2

Ensemble refinement?
by Joseph Noel 08 Aug '13

08 Aug '13

Is anyone having issues executing ensemble refinement in Phenix? I've tried several times now to run things and it always crashes without returning any maps and just returning a set of large PDB files such as: "TEAS-W273E-FPP-A_ensemble_0.6.pdb, TEAS-W273E-FPP-A_ensemble_0.8.pdb, TEAS-W273E-FPP-A_ensemble_0.9.pdb and TEAS-W273E-FPP-A_ensemble_0.1.pdb". The log file always ends with: " Time building chain proxies: 3.01, per 1000 atoms: 0.35 Number of scatterers: 8638 At special positions: 0 Unit cell: (125.75, 125.75, 120.07, 90, 90, 90) Space group: P 41 21 2 (No. 92) Number of sites at special positions: 0 Number of scattering types: 7 Type Number sf(0) S 19 16.00 P 2 15.00 Mg2+ 3 10.00 O 835 8.00 N 717 7.00 C 2765 6.00 H 4297 1.00 sf(0) = scattering factor at diffraction angle 0. AttributeError" I only entering a single chain with solvent and ions and hydrogens. All of the other parameters from the GUI are kept the same. Any suggestions or comments would be much appreciated. Joe ______________________________________________________________________________________ Joseph P. Noel, Ph.D. Arthur and Julie Woodrow Chair Investigator, Howard Hughes Medical Institute Professor, The Jack H. Skirball Center for Chemical Biology and Proteomics The Salk Institute for Biological Studies 10010 North Torrey Pines Road La Jolla, CA 92037 USA Phone: (858) 453-4100 extension 1442 Cell: (858) 349-4700 Fax: (858) 597-0855 E-mail: noel(a)salk.edu Publications & Citations: http://scholar.google.com/citations?user=xiL1lscAAAAJ Homepage Salk: http://www.salk.edu/faculty/noel.html Homepage HHMI: http://hhmi.org/research/investigators/noel.html ______________________________________________________________________________________

1 0