- phenixbb - phenix-online.org

phenix.superpose_pdbs and nucleic acids.. does it work?
by Francis E Reyes 14 Dec '11

14 Dec '11

Hi phenixer's. I need to do a superposition of two pdbs. I guess phenix.superpose_pdbs is my jiffy but it doesn't seem to work for nucleic acids (tried selection = 'name P' or selection = 'all'). LSQKAB seems to have the same limitation. They are the exact same PDBs (1 to 1 atom/residue mapping) but just shifted and rotated relative to the other. I need a superposition or the transformation matrix for all atoms. Traceback (most recent call last): File "/usr/local/phenix-1.7.1-743/phenix/phenix/command_line/superpose_pdbs.py", line 450, in <module> run(args=sys.argv[1:]) File "/usr/local/phenix-1.7.1-743/phenix/phenix/command_line/superpose_pdbs.py", line 434, in run m = manager(params.extract(), log) File "/usr/local/phenix-1.7.1-743/phenix/phenix/command_line/superpose_pdbs.py", line 171, in __init__ selection_moving = self.selection_moving) File "/usr/local/phenix-1.7.1-743/phenix/phenix/command_line/superpose_pdbs.py", line 306, in _alignment "bugs(a)phenix-online.org for help.") RuntimeError: Failed to identify corresponding atoms from sequence alignment. This is probably a bug - contact bugs(a)phenix-online.org for help. I may just strip it down to phosphate's and rename the residues to ALA but a quick phenix jiffy would be great. Thanks! F --------------------------------------------- Francis E. Reyes M.Sc. 215 UCB University of Colorado at Boulder

2 1

on residue bond deviation
by Dialing Pretty 14 Dec '11

14 Dec '11

Dear All, After a key step refinement byPhenix refine, I find there are 2 residues (total about 450 residues) with bad angles as analyzed by MolProity (deviation ratio 0.46%, goal is 0.1% ), both of the 2 residues are Pro, and both of the deviation angles are C-N-CA. Will you please give the comment on what do the 2 deviations mentioned above mean for my whole protein structure, or should I continue to refine my protein structure? I am looking forward to getting your reply. Dialing

2 1

on the input file for the Phase and build
by Dialing Pretty 14 Dec '11

14 Dec '11

Dear All, For the Phenix Phase and build, the input mtz file is the scala mtz file, and the ouput mtz file will be the mtz file with phase, is this correct? I am looking forard to getting a reply from you. Cheers, Dialing

2 1

strange MR problem
by intekhab alam 12 Dec '11

12 Dec '11

Hi All I have a 3.0A dataset of a protein-protein complex. I used one of the protein structure solved previously as a template and used phaser in CCP4 as well as in phenix. I used the sca file in phenix which gave a solution with 6 monomers in ASU which are packed as a hexamer. In ccp4 phaser i used the converted .mtz file which gave 6 molecules in ASU but with a different packing. In both the cases map fits well in the model and there is some extra density that may correspond to second protein partner of the complex molecule. when i tried to use the MR solution generated in phenix with the mtz file from ccp4 the model shows a lot of clashes. what kind of approach should i take so as to resolve this ambiguity. Regards -- INTEKHAB ALAM LABORATORY OF STRUCTURAL BIOINFORMATICS KOREA UNIVERSITY, SEOUL

2 4

Phenix version 1.7.3 released
by Paul Adams 12 Dec '11

12 Dec '11

The Phenix developers are pleased to announce that version 1.7.3 of Phenix is now available. Binary installers for Linux, and Mac OSX platforms are available at the download site: http://phenix-online.org/download/ Some of the new features in this version are: Graphical interface: ==================== - New GUIs: MR-Rosetta - viewer for diffraction images (also available as phenix.image_viewer) - editable heavy-atom sites in HySS and AutoSol General: ======== - New commands: phenix.ncs_average, phenix.model_vs_sequence, phenix.image_viewer, phenix.find_alt_orig_sym_mate, phaser.MRage - support for CNS v1.3 reflection files Refinement: =========== - Reduced memory requirement for map related calculations - Improved handling of parameter files in the GUI - RNA geometry target now includes pucker & base type-specific Chi 1 parameters - TLS group definitions from PDB file header will be automatically used if no groups are provided by the user - Bug fix in clashscore calculation for large structures with duplicate chain IDs plus SEGIDs - Torsion NCS restraints uses SEGIDs if present in input PDB file Molecular Replacement: ====================== - phaser (version 2.4): - support for pseudo-translational non-crystallographic symmetry (tNCS) in molecular replacement - detects presence of tNCS from inspection of native Patterson - characterizes tNCS parameters (translation, small rotational difference between copies) - accounts for effect of tNCS in molecular replacement calculations - note: support currently limited to one tNCS operator relating two molecules (or sets of molecules) - FAST search method turned on by default, replacing FULL search method - MR possible with a model containing a single atom - B-factor refinement turned on by default in molecular replacement - Bug fix for occasional and irreproducible segmentation faults for 32 bit binary running on 64 bit Linux - Packing function accepts cases where an ensemble with internal point group symmetry is on a special position, and pdb file output deletes atoms overlapped due to the symmetry of special position - phaser.MRage (version 0.1.0) - automated MR: - New command-line switches for common actions, e.g. setting the verbosity, getting PHIL parameters, etc - Define component by sequence, but calculate molecular weight from available models if sequence is omitted - Calculate number of copies to find from Matthews coefficient - Ability to perform homology search (currently only BLAST is supported, either local installation or NCBI service) - alignments automatically generated for "template" if target sequence is known - various options for queueing systems, e.g. qslot - option to automatically write out solutions - fully parallel space group exploration - phenix.find_alt_orig_sym_mate (new command): - For different molecular replacement solutions from the same dataset, finds the copy of moving_pdb closest to reference_pdb with respect to all symmetry operations and alternative origin shifts permitted by the spacegroup in moving_pdb Experimental Phasing and Model Building: ======================================== - AutoSol: - Fixed problems with clean_up=True deleting overall_best.pdb - Allow data file names of any length - Use NCS in scoring by default if ncs copies > 2 - Correct anisotropy and sharpen by default - AutoBuild: - Automatic anisotropy correction and sharpening for all maps - phenix.multi_crystal_average: - Automatic anisotropy correction and sharpening - Automatic filling of reflections to highest resolution of any dataset - Masking heavy-atom sites from averaging - Inputs grouped by crystal Miscellaneous: ============== - phenix.ready_set/phenix.reduce: - Improved handling of alternate conformations - phenix.cif_as_mtz - Support for files containing multiple datasets - Extraction of Hendrickson-Lattman coefficients and anomalous arrays - New map_to_asu and remove_systematic_absences options - phenix.maps: - Added atom selection parameter for simple omit maps (omit.selection) - phenix.ligand_identification: - Implemented methods to generate custom ligand library based on parent protein's SCOP or CATH terms, Pfam accession numbers, GO accession numbers, or InterPro ID - Improved ligand geometry after real-space refinement For a full list of changes see: http://www.phenix-online.org/documentation/CHANGES Please note that this publication should be used to cite use of Phenix: PHENIX: a comprehensive Python-based system for macromolecular structure solution. P. D. Adams, P. V. Afonine, G. Bunkóczi, V. B. Chen, I. W. Davis, N. Echols, J. J. Headd, L.-W. Hung, G. J. Kapral, R. W. Grosse-Kunstleve, A. J. McCoy, N. W. Moriarty, R. Oeffner, R. J. Read, D. C. Richardson, J. S. Richardson, T. C. Terwilliger and P. H. Zwart. Acta Cryst. D66, 213-221 (2010). Full documentation is available here: http://www.phenix-online.org/documentation/ There is a Phenix bulletin board: http://www.phenix-online.org/mailman/listinfo/phenixbb/ Please consult the installer README file or online documentation for installation instructions. Direct questions and problem reports to the bulletin board or: help(a)phenix-online.org and bugs(a)phenix-online.org Commercial users interested in obtaining access to Phenix should visit the Phenix website for information about the Phenix Industrial Consortium. The development of Phenix is principally funded by the National Institute of General Medical Sciences (NIH) under grant P01-GM063210. We also acknowledge the generous support of the members of the Phenix Industrial Consortium. -- Paul Adams Deputy Division Director, Physical Biosciences Division, Lawrence Berkeley Lab Division Deputy for Biosciences, Advanced Light Source, Lawrence Berkeley Lab Adjunct Professor, Department of Bioengineering, U.C. Berkeley Vice President for Technology, the Joint BioEnergy Institute Laboratory Research Manager, ENIGMA Science Focus Area Building 64, Room 248 Tel: 1-510-486-4225, Fax: 1-510-486-5909 http://cci.lbl.gov/paul Lawrence Berkeley Laboratory 1 Cyclotron Road BLDG 64R0121 Berkeley, CA 94720, USA. Executive Assistant: Louise Benvenue [ LBenvenue(a)lbl.gov ][ 1-510-495-2506 ] --

1 0

ncs_average output map
by Ursula Schulze-Gahmen 12 Dec '11

12 Dec '11

I am puzzled about my ncs-averaged output map. I determined the ncs relation using phenix.find_ncs with pdb and mtz input. The program finds 2 NCS groups with 3 molecules each and the rmsds are low. I then read in the ncs_spec file and an mtz file with map coefficients into phenix.ncs_average. I used the FOFCWT, PHFOFCWT coefficients. The output map seems cleaner than the unaveraged map and shows some of the same features, but the density around the 3 molecules is not identical. 2 of them look very similar, but the density around one of them is somewhat different. How is this possible if the map is averaged? Thanks for any suggestions or explanations. Ursula -- Ursula Schulze-Gahmen, Ph.D. Assistant Researcher UC Berkeley, QB3 356 Stanley Hall #3220 Berkeley, CA 94720-3220

2 3

AutoBuild Question
by Page, Rebecca 12 Dec '11

12 Dec '11

Dear phenix community: I ran AutoBuild.dat with a seq.dat file that has 3 separate sequences. Chains A and B were built, but chain C was not, even though there is density for the third chain. It occurred to me that only two of the three chains may have been read in. I seem to recall that seq.dat can be a maximum of two sequences, but am not sure. So, is there a limit to the number of sequences that can be present in seq.dat and if so, what is the limit. Thanks in advance, Rebecca

4 4

phenix.refine_failed
by Xue-Yuan Pei 12 Dec '11

12 Dec '11

Dear Phenix Experts: My phenix.refine has developped the following error (see below). I tried to check my PDB format, there is not B-factor error. so what could be wrong that stoped my phenix.refine without PDB and .def files in the output? Many thanks for your advice. Xue yuan Pei, Ph.D Department of Biochemistry University of Cambridge Cambridge CB2 1GA U.K > Error messages: > Traceback (most recent call last): > File > "/usr/local/xtal/phenix-1.7.2-869/phenix/phenix/command_line/refine.py", > line 11, in <module> > command_line.run(command_name="phenix.refine", args=sys.argv[1:]) > File > "/usr/local/xtal/phenix-1.7.2-869/phenix/phenix/refinement/command_line.py", > line 105, in run > refine_object.write_refined_pdb(f_pdb) > File > "/usr/local/xtal/phenix-1.7.2-869/phenix/phenix/refinement/driver.py", > line 1284, in write_refined_pdb > ignore_hd = not self.neutron_refinement) > File > "/usr/local/xtal/phenix-1.7.2-869/cctbx_project/mmtbx/model_statistics.py", > line 534, in __init__ > ignore_hd = ignore_hd) > File > "/usr/local/xtal/phenix-1.7.2-869/cctbx_project/mmtbx/model_statistics.py", > line 449, in __init__ > molprobity_scores=True) > File "/usr/local/xtal/phenix-1.7.2-869/cctbx_project/mmtbx/model.py", > line 1063, in geometry_statistics > molprobity_scores = molprobity_scores) > File > "/usr/local/xtal/phenix-1.7.2-869/cctbx_project/mmtbx/model_statistics.py", > line 112, in __init__ > self.clashscore_obj.analyze_clashes(hierarchy = pdb_hierarchy) > File > "/usr/local/xtal/phenix-1.7.2-869/cctbx_project/mmtbx/validation/clashscore.py", > line 167, in analyze_clashes > stdin_lines=tmp_r.as_pdb_string()) > File > "/usr/local/xtal/phenix-1.7.2-869/cctbx_project/iotbx/pdb/hierarchy.py", > line 527, in as_pdb_string > siguij=siguij) > RuntimeError: atom B-factor does not fit into F6.2 format: > "ATOM 1 N GLY A 1 " > B-factor: 111147.29 >

3 3

phenix.refine_failed
by Christian Roth 12 Dec '11

12 Dec '11

Hi, indeed your B-Factor shown in the error message for Gly A 1 is unrealistic. Is maybe in your pdb-file a space missing? I think a 8 digit B-Factor is not possible in the pdb format. Best Regards Christian

1 0

Phaser not writting out ligands
by Yarrow Madrona 09 Dec '11

09 Dec '11

Hello, I am having trouble with Phaser not writing out my ligands (HEME and small molecule) to it's output PDB even though they are in the input PDB. I am running through the GUI using Phenix version 1.7.3-921. Phenix refine recognizes HEME from the CCP4 library I think so I am not sure why I at least do not get a HEME in the output PDB. Thanks for your help. -Yarrow -- Yarrow Madrona Graduate Student Molecular Biology and Biochemistry Dept. University of California, Irvine Natural Sciences I, Rm 2403 Irvine, CA 92697 -- Yarrow Madrona Graduate Student Molecular Biology and Biochemistry Dept. University of California, Irvine Natural Sciences I, Rm 2403 Irvine, CA 92697

2 1

Please apply for RapiData 2012, a course on Data Collection and Structure Solving at the NSLS.
by Robert Sweet 09 Dec '11

09 Dec '11

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

map output from ncs_average
by Ursula Schulze-Gahmen 09 Dec '11

09 Dec '11

I used ncs_average to create an averaged fo-fc map. From the manual I thought that the output map is a CCP4 map. I have however difficulties to read it into Pymol as ccp4 map. It is not recognized. Any suggestions? Ursula -- Ursula Schulze-Gahmen, Ph.D. Assistant Researcher UC Berkeley, QB3 356 Stanley Hall #3220 Berkeley, CA 94720-3220

3 4

print_sequence
by Bradley Hintze 09 Dec '11

09 Dec '11

Hi, I am analyzing many PDBs vai a script. Thus, I am feeding print_sequence the full path of the given PDB. This has the unfortunate effect of having a filename more than 80 char. As such the '>filename' line is split into two lines making parsing less than ideal. Is there a way I can avoid having the full path on the '>' line and still give print_sequence the full path? Or can that feature be added? Thanks, -- Bradley

2 1

excluding atoms from refinement
by Crystal VanderZanden 08 Dec '11

08 Dec '11

Does anybody know if there a general code that will allow us to exclude something from refinement? I would like to specify that there is a particular part of my structure that I don't want to refine the position. I've thought about deleting those atoms from the pdb file, but then that phasing information is excluded from the refinement. Any ideas? Thanks, Crystal

3 2

Re: [phenixbb] Autosol problems
by Tom Terwilliger 08 Dec '11

08 Dec '11

Hi Stacy, I'm sorry for the trouble! Normally autosol should not have problems with resolve, so this indicates that either (1) the data are unusual or (2) there is some problem in autosol. Can you possibly send along (terwilliger(a)lanl.gov) the log file from autosol and the LAST.log file if there is one? Also what version of phenix is this? (type "phenix.version") Thanks! -Tom T On Dec 8, 2011, at 8:55 AM, stacy William wrote: > Dear All... > I am new user for Phenix. I am right now working with SAD data. I had run Ausol in GUI interface with my SAD anamolous data. but when i am using the output file for my autobuild program it is giving problems with Resolve... n terminating with an error message suggesting to check resolve file. Do i need to change some parameters in Autosol to produce right output files which in turn can work as input files for Autobuild program > Thanks > Stacy > > 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

1 0

automr wizard: searching for all plausible space groups option
by Vitali Stanevich 08 Dec '11

08 Dec '11

Hi, I'm trying to solve structure by MR with phenix 1.6.1-357 version. I have 100% identity search model and can find solution without problem in ccp4i in P212121 sg. But, I can't do it in phenix.automr command line because program looks only for P222 space group (pointed out by Randy Read). I use .sca from hkl2000 and running program from parameters file. I tried: ... all_plausible_sg_list = True use_all_plausible_sg = True ... all_plausible_sg_list = P222 P2221 P21212 P212121 use_all_plausible_sg = True ... space_group = "P 21 21 21" ... They lead to report that no solution found or program being stalled (last variant). What should I do to make phaser search for all space groups in P222 point group? Thanks, Vitali

3 2

rosetta mr
by Shya Biswas 06 Dec '11

06 Dec '11

Hi all, I was wondering how to generate .ali file that I need to run rosetta MR. I have most of the files for the run except this one. the documentation did not help a lot. thanks, Shya

3 3

pseudo-merohedral twinning
by Zhou, Tongqing (NIH/VRC) [E] 06 Dec '11

06 Dec '11

Hi Pavel, I just collected a 2.85 A dataset with space group P212121, MR went well, however, the refinement was not fine: (eg: Rigid body r=48%, it won't go down anymore). I then ran xtriage, it turned out the my data was pseudo-merohedral twinned. The twin fraction is ~0.02. Twin_law=k,h,-l. It there a way to de-twin the data and then do the refinement or just use the twin_law="k,h,-l" at the command line? Is there a way to specify twin_fraction in the paramenter file? Also re-processed data in P41212 or P43212, the ~25% Rsym indicated the data should still be in P222. Thanks, Tongqing Tongqing Zhou, Ph.D. Staff Scientist Structural Biology Section Vaccine Research Center, NIAID/NIH Building 40, Room 4607B 40 Convent Drive, MSC3027 Bethesda, MD 20892 (301) 594-8710 (Tel) (301) 793-0794 (Cell) (301) 480-2658 (Fax) ****************************************************************************** The information in this e-mail and any of its attachments is confidential and may contain sensitive information. It should not be used by anyone who is not the original intended recipient. If you have received this e-mail in error please inform the sender and delete it from your mailbox or any other storage devices. National Institute of Allergy and Infectious Diseases shall not accept liability for any statements made that are sender's own and not expressly made on behalf of the NIAID by one of its representatives. ******************************************************************

4 3

pdbtools
by Bradley Hintze 06 Dec '11

06 Dec '11

Developers, I am doing some automated analysis on many pdbs. My analysis is contained to a single chain in a given pdb. Aa such, to cut down on computing time I am generating many pdbs cut to the desired chain. When using pdbtools I get the following Sorry: Fatal problems interpreting PDB file: Number of atoms with unknown nonbonded energy type symbols: 46 Please edit the PDB file to resolve the problems and/or supply a CIF file with matching restraint definitions, along with apply_cif_modification and apply_cif_link parameter definitions if necessary. Also note that phenix.ready_set and phenix.elbow are available for creating restraint definitions (CIF files). Is there a way around this without creating a CIF file or another tool to cut pdbs? Thanks, Bradley -- Bradley J. Hintze Duke University Graduate Student Department of Biochemistry

2 1

Rfree flag transfer
by Zhang yu 06 Dec '11

06 Dec '11

Hi, I am using phenix_1.7.1-743 GUI. I would like to transfer Rfree tag from one dataset to the other one, but had the following error, "The uni cell for the Miller array /***/*** (186.326, 103.98, 297.48, 90, 98.465, 90)" is significantly different than the output unit cell (183.486, 103.5, 295.416, 90, 98.828, 90)" Two datasets are from same kind of crystal, but with slight different cells. I guess that is why "Reflection file editor" refused to do it. Could anyone tell me how to avoid the error? Thanks Yu

2 3

autoMR wizard
by Vitali Stanevich 06 Dec '11

06 Dec '11

Hi, I'm trying to solve structure by molecular replacement. We already have structure for this protein and new is exactly the same thing with a ligand soaked. I can solve it without a problem with ccp4 gui (tfz=47.7), but not with phenix.automr, which is disappointing. Can somebody point me out what's wrong with my data input? I did it with this command: phenix.automr output.sca complex_search.pdb identity=90 mass=93500 copies=1 or with .eff file (attached) Thanks, Vitali

2 1

H related bond deviaion
by Dialing Pretty 06 Dec '11

06 Dec '11

Dear All, I did a Phenix refine, and before the refine all the H were added. MolProbity analysis demonstrated that almost everything are acceptable after the refinement. However when I did the geometry analysis by the COOT validation, I find all the residues having a significant deviation bond involving H. with z score higher than 7. I am looking forward to getting an explanation from you on it,including how to avoid it. Cheers, Dialing

3 2

bond deviation deected by Coot
by Dialing Pretty 05 Dec '11

05 Dec '11

Dear All, After Pheneix refine, the quality of the structure improved significantly as analyzed by Molprobity. However the bond deviation increased significantly after the refinement. Will you please tell me which refine function can be used to reduce the bond deviation detected by the Coot geometry analysis? Cheers, Dialing

3 3

rotamer outlier reduction by coot
by Dialing Pretty 05 Dec '11

05 Dec '11

Dear All, I am leaning rotamer outlier reduction by Coot. After select a rotamer and acception, I find some hydrogens in the outlier rotamer will be isolated from the rotamer, and the rotamer will have new hydrogens at its new position. Is any influence by further refine by phenix? After new rotamer acception, I have did the rotamer auto-fit, I find there is a trend on that several oxygens can pack together. How to process this problem? After I save the new PDB file by coot after I selected the new rotamers to replace the outlier rotamers, I did a Phenix refine. There is an error message at the start says "Number of atoms with unknown nonbonded energy typensymbol: 2228". What does it mean and how to solve the problem? Is any simple method by which I can do the rotamer outlier reduction by Coot? I am looking forward to getting your reply. Dialing

3 2

documents explaining parameter settings in phenix-1.7.2
by Ruan Wenjie 04 Dec '11

04 Dec '11

Hi, I want to know, about the new Phenix-1.7.2 version, whether I can find some detailed documents explaining how to set the parameters..since there are so many options and I'm not clear their meaning..Thanks~ best, wenjie

6 9