- phenixbb - phenix-online.org

CCP4 6.3.0 released
by ronan.keegan＠stfc.ac.uk 16 Jul '12

16 Jul '12

Dear All, We're very pleased to announce the release of the latest version of the CCP4 Software Suite. Version 6.3.0 (Settle) is now available from the CCP4 download website: http://www.ccp4.ac.uk/download.php The release is available for Linux, Mac OSX and Windows platforms, and delivers a number of new programs, updates and features: 1. New Components and Updates: New: Aimless (Scaling and Merging, replacement for Scala) Prosmart (Generation of external restraints for Refmac) Nautilus (Automatic building of RNA/DNA) Ample (Automatic Ab initio modelling for Molecular Replacement) Dimple (Automatic difference map calculation) Zanuda (Space group checker/validator) CMapCut (Preparation of electron-density search model for MR) Comit (Composite omit maps) Gesamt (Structure aligner not using secondary structure assignments) ViewHKL (Graphical viewer for reflection data) Updated: Mosflm v7.0.9 (Data processing/reduction) iMosflm v1.0.7 (Mosflm GUI) Pointless v1.6 (Laue group determination) Refmac v5.7 (Macromolecular refinement) Phaser v2.5.1 (Molecular Replacement and Experimental Phasing) MrBump (Automatic Molecular Replacement, now includes phaser.sculptor models) Xia2 (Automated data reduction) CCP4mg v2.6.2 (CCP4 Molecular Graphics, much improved on features and stability) CCP4 Graphical User Interface: - new automation module - provides interfaces for new and updated modules - ViewHKL option for viewing reflection (mtz) files More information on updates is available from: http://www.ccp4.ac.uk/html/CHANGESinV6_3.html 2. ARP/wARP co-distribution Release 6.3.0 marks a milestone in collaboration between CCP4 and ARP/wARP teams by releasing ARP/wARP Model Building software v7.3 jointly with CCP4 Software Suite. In order to get ARP/wARP Software, academic users need only to select the appropriate package from download pages or with our new Package Manager, and type in their name and e-mail address where required. Commerical users will be required also to provide their reference number, obtainable separately from EMBL-EM. ARP/wARP 7.3 includes the following major updates and changes (on behalf of Victor Lamzin): - The protocols for the use of NCS in protein chain tracing have been improved and expanded to a wider range of resolutions. - Common ligands can now be automatically identified and fit into a specified electron density. - The algorithm for modelling partial ligands (cocktail screening) has been replaced with a new one that applies a number of sophisticated numerical features - The ARP/wARP graphics front-end, ArpNavigator, automatically assigns secondary structure and has an improved visual presentation of a static model or real time protein or ligand building. - Supported computer platforms are Mac powerpc, Mac Intel and Linux (including 32 and 64-bit versions). - The ARP/wARP installer has been considerably updated on all platforms. - Joint software release and authentication with the CCP4 is now possible. - CCP4 6.3.0 and Refmac 5.7.0029 are the recommended versions to use with ARP/wARP 7.3. 3. 64-bit precompiled packages Starting with version 6.3.0, CCP4 distributes both 32-bit and 64-bit precompiled binaries for Linux. Note that 64-bit codes will not run on 32-bit architectures, and 32-bit codes will run on 64-bit setups only if 32-bit libraries are installed (not a default for most Linuxes). THEREFORE: PLEASE MAKE SURE THAT YOU HAVE SET THE 32-bit/64-bit SWITCH ON TOP OF THE DOWNLOAD PAGES CORRECTLY BEFORE ALL! 4. "Friendly" Coot and CCP4mg distributions Coot is known to be less portable than the rest of CCP4 Program Suite. CCP4mg is more portable and requires special builds only for older Linuxes. Our new Package Manager will try to automatically choose Coot and CCP4mg bundles that are most appropriate for your system. In addition, our new download pages will offer you a range of precompiled Coot and CCP4mg packages, from where you may choose ones that work for you. 5. New distribution mechanism We hope that you will find our new Download Pages more intuitive and convenient than before. We have worked hard to improve our Windows installer, to make it working faster and smoother. On Linux and Mac platforms, we recommend using Package Manager, a small graphical application that allows you to select the desired components, and automatically download and install them on your system. Note that Package Manager features resumable downloads, which may be useful for users with slow or less reliable Internet connection. We would like to thank all of the developers who have contributed into CCP4 6.3.0, and all of those who have helped in testing it. The following publication should be used to cite the use of CCP4: M. D. Winn, C. C. Ballard, K. D. Cowtan, E. J. Dodson, P. Emsley, P. R. Evans, R. M. Keegan, E. B. Krissinel, A. G. W. Leslie, A. McCoy, S. J. McNicholas, G. N. Murshudov, N. S. Pannu, E. A. Potterton, H. R. Powell, R. J. Read, A. Vagin and K. S. Wilson (2011) Overview of the CCP4 suite and current developments, Acta Cryst D67 235-242. CCP4 is supported by the BBSRC UK through grant BB/F0202281, as well as by the number of industrial users. CCP4 6.3.0 is out now, and we hope that you will find it a helpful advance in crystallographic computing. Best wishes, CCP4 Core Team -- Scanned by iCritical.

1 0

Van der waals interaction in phenix
by subhash bihani 16 Jul '12

16 Jul '12

Hi all, In one of my structure I have to put some water molecule in alternate confimation with amino acid residues and in another case with ligand. How do I do it in Phenix. Regards -- Subhash Bihani Scientific Officer Solid State Physics Division Bhabha Atomic Research Centre Trombay, Mumbai-85 Ph. 25594688 (o) alternate email: subhashbihani(a)yahoo.co.in

3 2

default treatment of Intensities and structure factors in phenix.refine
by Yarrow Madrona 16 Jul '12

16 Jul '12

Hi, I recently noticed from reviewing the phenix.refine log file that phenix is scaling the input intensities using the French-Wilson Method. However, I have already run Truncate in CCP4 and have structure factors in my MTZ. So it appears that Phenix.refine ignores structure factors and uses intensities regardless. The reason I am concerned is that Nat recently mentioned that the default for the sigma I cutoff is zero. However, The default in Truncate I thought was -3. If I understand correctly, the treatment in truncate allows weak reflections to be scaled and converted to structure factors. It seems that setting the I sigma cutoff to zero would cause you to loose some useable data. Currently, the sigma_fobs_rejection_criterion and sigma_ibos_rejection_criterion = None. This is the default I guess which would be 0. So my questions are: 1. Will setting: sigma_iobs_rejection_criterion = -3 result in Phenix using a sigma I cutoff of -3? 2. How do I tell phenix to use the structure factors created by Truncate? Thanks for your help. I apologize if I have some misconceptions as I am still relatively new to crystallography. -Yarrow -- Yarrow Madrona Graduate Student Molecular Biology and Biochemistry Dept. University of California, Irvine Natural Sciences I, Rm 2403 Irvine, CA 92697

5 7

average temperature factor
by Ursula Schulze-Gahmen 15 Jul '12

15 Jul '12

Is there an easy way to calculate average temperature factors in phenix or another program after ADP + TLS refinement in phenix? Ursula -- Ursula Schulze-Gahmen, Ph.D. Assistant Researcher UC Berkeley, QB3 356 Stanley Hall #3220 Berkeley, CA 94720-3220

6 8

output of labelit.check_pdb_symmetry
by "Weiergräber, Oliver H." 15 Jul '12

15 Jul '12

Hello, I am using labelit.check_pdb_symmetry in order to evaluate a crystal structure with sg P21 (2 copies/asu) and all angles indistinguishable from 90 degs for potential higher (i.e. orthorhombic) symmetry. Comparison of refinement statistics clearly indicates that the most likely scenario is perfect pseudomerohedral twinning of a true P21 crystal (Rfree=0.21). Accordingly, labelit.check_pdb_symmetry reports a striking discrepancy between Rsymop(Iobs) and Rsymop(Icalc): 8 % vs 31 %, i.e. there is an additional (apparent) 2-fold symmetry in the data which cannot be explained by the model. In an attempt to force a solution for comparison and visualization purposes, I have increased the r_symop_tol parameter to 0.35; the program starts to evaluate possible space groups with fTF, yielding scores (these are correlation coefficients, right?) way above 90 %. Still, in the end it reports there are no higher symmetry solutions. Why is this? Does the program check for the contrast between the best solution and the others? Changing peak_cutoff to a "discriminative" value did not help... Best regards, Oliver ------------------------------------------------------------------------------------------------ ------------------------------------------------------------------------------------------------ Forschungszentrum Juelich GmbH 52425 Juelich Sitz der Gesellschaft: Juelich Eingetragen im Handelsregister des Amtsgerichts Dueren Nr. HR B 3498 Vorsitzender des Aufsichtsrats: MinDir Dr. Karl Eugen Huthmacher Geschaeftsfuehrung: Prof. Dr. Achim Bachem (Vorsitzender), Karsten Beneke (stellv. Vorsitzender), Prof. Dr.-Ing. Harald Bolt, Prof. Dr. Sebastian M. Schmidt ------------------------------------------------------------------------------------------------ ------------------------------------------------------------------------------------------------ Kennen Sie schon unseren neuen Film? http://www.fz-juelich.de/film Kennen Sie schon unsere app? http://www.fz-juelich.de/app

6 9

Re: [phenixbb] phenixbb Digest, Vol 80, Issue 6
by Yuri 13 Jul '12

13 Jul '12

On Fri, 13 Jul 2012 12:13:15 -0700, phenixbb-request(a)phenix-online.org wrote: > Send phenixbb mailing list submissions to > phenixbb(a)phenix-online.org > > To subscribe or unsubscribe via the World Wide Web, visit > http://phenix-online.org/mailman/listinfo/phenixbb > or, via email, send a message with subject or body 'help' to > phenixbb-request(a)phenix-online.org > > You can reach the person managing the list at > phenixbb-owner(a)phenix-online.org > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of phenixbb digest..." > > > Today's Topics: > > 1. Re: segmented rigid body refinement (Nathaniel Echols) > 2. Re: segmented rigid body refinement (Pavel Afonine) > 3. Re: average temperature factor (David Cobessi) > 4. Re: segmented rigid body refinement (Gino Cingolani) > 5. Re: average temperature factor (Pavel Afonine) > 6. Re: output of labelit.check_pdb_symmetry (Nicholas Sauter) > 7. Re: average temperature factor (Nathaniel Echols) > > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Thu, 12 Jul 2012 09:40:41 -0700 > From: Nathaniel Echols <nechols(a)lbl.gov> > To: PHENIX user mailing list <phenixbb(a)phenix-online.org> > Subject: Re: [phenixbb] segmented rigid body refinement > Message-ID: > <CALeAa1NE=F6Quwvsm=ijGf_j-QMxwN1KqStpwiE6NJ7bmyKnsw(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > On Thu, Jul 12, 2012 at 9:33 AM, Gino Cingolani < > Gino.Cingolani(a)jefferson.edu> wrote: > >> The latter....something more analytical to split individual chains >> into >> smaller groups >> (e.g. like to define tls groups). >> > > Airlie McCoy just added some functionality to do this in Phaser using > normal mode analysis. It is still somewhat experimental and in flux > (i.e. > the program name isn't even finalized), but I will email you off-list > with > details. > > -Nat >

1 0

segmented rigid body refinement
by Gino Cingolani 12 Jul '12

12 Jul '12

Hi all, I'm refining a large oligomeric structure at low res (~6A). I'd like to define groups to use for 'segmented rigid body refinement'. Is there any rational way to define 'segments'? My macromolecule is already an oligomer on its own, so, intuitively, each protomer is a segment. Thanks, 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 Website: http://www.cingolanilab.org ****************************************************************************** "Nati non foste per viver come bruti, ma per seguir virtute e canoscenza" ("You were not born to live like brutes, but to follow virtue and knowledge") Dante, The Divine Comedy (Inferno, XXVI, vv. 119-120) ________________________________________ From: phenixbb-bounces(a)phenix-online.org [phenixbb-bounces(a)phenix-online.org] on behalf of Nathaniel Echols [nechols(a)lbl.gov] Sent: Thursday, July 12, 2012 7:44 AM To: PHENIX user mailing list Subject: Re: [phenixbb] RFZ is higher than TFZ On Thu, Jul 12, 2012 at 12:38 AM, 周文昌 <wenchangyu2006(a)gmail.com> wrote: > I had one problem using phaser-MR to solve my structure. Phaser-MR gives me > one solution with TFZ=15.4 and LLG=165 which seems to be promising but > phenix.refine does not give me a low R values (0.4884/0.5424). I went back > and checked the solution from phaser-MR, I found RFZ=18.8, higher than TFZ. > Does anyone think this is a problem to solve my structure? > > The space group is C121, I believe there are 4 copies of my search model in > ASU. What is the resolution, and what is the sequence identity of the search model? I don't know whether the high RFZ means anything, but the first thing I would check is whether the packing of the MR solution makes sense; Phaser will do a good job preventing molecules from overlapping (unless you changed it to use much less conservative criteria), but it won't guarantee that the molecules pack into a nice continuous lattice. If there are big gaps, it may not be a real solution. I think we need to add an entire FAQ page for high R-factors. However, there are an increasing number of tools and/or tricks in Phenix to deal with problems like this - in approximate order of runtime: - run phenix.refine for more cycles - at least 20. See this paper for an example of why this can be important: http://journals.iucr.org/d/issues/2012/07/00/kw5044/index.html - use phenix.morph_model (refer to previously mentioned paper) - use phenix.den_refine (or DEN in CNS) - use MR-Rosetta Also, if you're willing to share your data with us, we're collecting difficult-to-refine structures for benchmarking new methods. (It will of course be kept confidential.) -Nat _______________________________________________ phenixbb mailing list phenixbb(a)phenix-online.org http://phenix-online.org/mailman/listinfo/phenixbb

3 5

v1.8 bug ? memory allocation crash in map calculation (phenix.refine, phenix.maps)
by Phil Jeffrey 12 Jul '12

12 Jul '12

This is in F432 in a=b=c=340A at 2.85A resolution. Mac OSX 10.5.8, the machine has 8Gb memory of which something like 6Gb is hanging around begging to be used. 1,327,104,000 if that's in bytes doesn't seem all that excessive. On the other hand, 340*340*340= 39,419,600 so it's perhaps a little greedy. Limits set by shell: phil% limit cputime unlimited filesize unlimited datasize unlimited stacksize 65532 kbytes coredumpsize 0 kbytes memoryuse unlimited descriptors 256 memorylocked unlimited maxproc 266 phenix.refine does the same thing right at the end, which is kind of a shame since it takes quite a long time to run. phenix.maps does the same thing even if I ask it just to generate F,phase in MTZ format. [snip] Compute maps. Python(632) malloc: *** mmap(size=1327104000) failed (error code=12) *** error: can't allocate region *** set a breakpoint in malloc_error_break to debug show_stack(1): /usr/local/phenix-1.8-1069/cctbx_project/cctbx/miller/__init__.py(3067) double_step_filtration show_stack(2): /usr/local/phenix-1.8-1069/cctbx_project/mmtbx/maps/__init__.py(416) map_coefficients_from_fmodel show_stack(3): /usr/local/phenix-1.8-1069/cctbx_project/mmtbx/maps/__init__.py(543) __init__ show_stack(4): /usr/local/phenix-1.8-1069/build/mac-intel-osx/../../cctbx_project/mmtbx/command_line/maps.py(325) run show_stack(5): /usr/local/phenix-1.8-1069/build/mac-intel-osx/../../cctbx_project/mmtbx/command_line/maps.py(402) <module> Bus error (Python call stack 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.

3 8

RFZ is higher than TFZ
by 周文昌 12 Jul '12

12 Jul '12

Dear all, I had one problem using phaser-MR to solve my structure. Phaser-MR gives me one solution with TFZ=15.4 and LLG=165 which seems to be promising but phenix.refine does not give me a low R values (0.4884/0.5424). I went back and checked the solution from phaser-MR, I found RFZ=18.8, higher than TFZ. Does anyone think this is a problem to solve my structure? The space group is C121, I believe there are 4 copies of my search model in ASU. thanks, Wenchang

2 1

sigma_iobs_rejection_criterion
by Yarrow Madrona 11 Jul '12

11 Jul '12

Hello, I am refining a structure using a.eff file that I generated in Phenix version 1.7 GUI and later ran in the command line. When I try to run the same refinement in Phenix V.1.8 I get the following error: For French and Wilson scaling, sigma_iobs_rejection_criterion must be set between -4.0 and -1.0 I see that it was set to zero previously. My questions are: 1. Is this the same observed criterion that is reported in the RCSB Deposition? 2. To be honest I don't understand what this is. Is it a limit for the error in Iobs, below which the reflection will be discarded? 3. What is a reasonable setting for this input and how does it change the refinement? Thank you (I appreciate your help as I know these may be very basic questions). -Yarrow -- Yarrow Madrona Graduate Student Molecular Biology and Biochemistry Dept. University of California, Irvine Natural Sciences I, Rm 2403 Irvine, CA 92697

3 4

Phenix autosol error 2L MAD
by Graeme Winter 10 Jul '12

10 Jul '12

Hi Folks, 2L MAD data set where SHELXC/D/E works fine I (repeatedly) get this error, through GUI or scripts... resolution: 1.7 recommended_resolution: 1.9 datafile_for_hyss: /Users/graeme/projects/test/12287_build/AutoSol_run_10_/TEMP0/ANISO_subset_PHX_I(+)_WAVE1.sca_ano_1.sca aniso_corrected: True aniso_b_min: 13.2 aniso_range_of_b: 6.8 RUNNING HYSS on ANISO_subset_PHX_I(+)_WAVE1.sca_ano_1.sca from PDS as miller arrays *************ERROR ENDING ******************* ******************************************************************************** Failed to carry out AutoSol_hyss: HYSS failed...Probably no iso or ano differences... ******************************************************************************** *************ERROR ENDING ******************* The data set works fine with a third (radiation damaged) data set. There is also a massive anomalous signal in the data (2L) as it stands. here's the PHIL: autosol { seq_file = 12287.seq sites = 8 atom_type = Se crystal_info { unit_cell = 51.59 51.59 157.69 90 90 90 space_group = "P 43 21 2" solvent_fraction = 0.52 chain_type = PROTEIN resolution = 0 change_sg = True } wavelength { data = subset.mtz lambda = 0.97966 f_prime = -9.9 f_double_prime = 3.9 labels = 'I(+)_WAVE1 SIGI(+)_WAVE1 I(-)_WAVE1 SIGI(-)_WAVE1' } wavelength { data = subset.mtz lambda = 1.00000 f_prime = -3.3 f_double_prime = 0.5 labels = 'I(+)_WAVE2 SIGI(+)_WAVE2 I(-)_WAVE2 SIGI(-)_WAVE2' } } Have I done something wrong? Thanks, Graeme

2 1

special position
by Karolina Michalska 10 Jul '12

10 Jul '12

Hi all, I'm refining a structure with a phosphate ion on a 3-fold axis. I set occupancy manually to 0.33 and after refinement the ion has huge B-factor (close to 100) and there is a lot of positive electron density. I found somewhere that phenix recognizes special positions and there is not need to adjust occupancy manually. so I set it to 1 again and re-refined. This time my phosphate sits in a mixture of negative and positive density with B factor around 90 (for comparison Refmac refines it to ~20). Now I'm confused which procedure to follow. Karolina

2 1

I/O error in writing very large B-factors
by Phil Jeffrey 09 Jul '12

09 Jul '12

This occurs at the end of refinement. Quite why it's trying to *read* from a string with Fortran-esque overflow formatting while in a writing step I don't know, but the PDB file gets truncated at REMARK IF THIS FILE IS FOR PDB DEPOSITION: REMOVE ALL FROM THIS LINE UP. The structure is not idyllic, clearly, but I'd like to see the extent of the problem rather than just the first atom. Got to be pretty easy to trap for and map to 999.99 either on the way in or out. Writing refined structure to PDB file: /Users/phil/Structures/reference_models/start_refine_001.pdb n_use = 21264 n_use_u_iso = 21264 n_use_u_aniso = 0 n_grad_site = 0 n_grad_u_iso = 21264 n_grad_u_aniso = 0 n_grad_occupancy = 0 n_grad_fp = 0 n_grad_fdp = 0 total number of scatterers = 21264 Traceback (most recent call last): File "/usr/local/phenix-1.7.3-928/build/mac-intel-osx/../../phenix/phenix/command_line/refine.py", line 11, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/usr/local/phenix-1.7.3-928/phenix/phenix/refinement/command_line.py", line 105, in run refine_object.write_refined_pdb(f_pdb) File "/usr/local/phenix-1.7.3-928/phenix/phenix/refinement/driver.py", line 1379, in write_refined_pdb ignore_hd = not self.neutron_refinement) File "/usr/local/phenix-1.7.3-928/cctbx_project/mmtbx/model_statistics.py", line 589, in __init__ ignore_hd = ignore_hd) File "/usr/local/phenix-1.7.3-928/cctbx_project/mmtbx/model_statistics.py", line 504, in __init__ molprobity_scores=True) File "/usr/local/phenix-1.7.3-928/cctbx_project/mmtbx/model.py", line 1063, in geometry_statistics molprobity_scores = molprobity_scores) File "/usr/local/phenix-1.7.3-928/cctbx_project/mmtbx/model_statistics.py", line 112, in __init__ self.clashscore_obj.analyze_clashes(hierarchy = pdb_hierarchy) File "/usr/local/phenix-1.7.3-928/cctbx_project/mmtbx/validation/clashscore.py", line 182, in analyze_clashes self.pdb_hierarchy = pdb.hierarchy.input(pdb_string=input_str).hierarchy File "/usr/local/phenix-1.7.3-928/cctbx_project/iotbx/pdb/hierarchy.py", line 1044, in __init__ source_info=source_info, lines=flex.split_lines(pdb_string)) File "/usr/local/phenix-1.7.3-928/cctbx_project/iotbx/pdb/__init__.py", line 582, in input return ext.input(source_info=source_info, lines=lines) ValueError: string, line 5148: ATOM 17 N ILE A 268 -16.878 -22.262 49.415 1.00****** A N ------------------------------------------------------------^ not a floating-point number. Phil Jeffrey Princeton

2 1

Iridium hexamine problems with eLBOW
by Benoît Masquida 07 Jul '12

07 Jul '12

I experience problems generating a correct CIF file for iridium hexamine (IRI in the eLBOW database). The geometry is supposed to be octahedral and eLBOW outputs a prism regardless of the options checked in the GUI. eLBOW actually does the same for the octahedral Mg(OH2)6. Does someone know how to rule this out? Thanks to everybody. Ben -- Benoît Masquida Chargé de Recherche CNRS IBMC Strasbourg UPR 9002 ARN b.masquida(a)ibmc-cnrs.unistra.fr +33 (0)3 88 41 70 45 +33 (0)6 70 23 79 50 +33(0)3 88 60 22 18 (fax)

5 6

Subject: no kidding - we need small crystals!
by jens j birktoft 05 Jul '12

05 Jul '12

Hi everybody, Any suggestions regarding how to generate small!! crystals? Yes, in this case that is what we want While we can generate crystals of size 100-500 microns, this time we something in the range of .5-5 microns. Any ideas and suggestions would be greatly appreciated, Cheers Jens

7 6

ask refinement question including an aligo as the ligand
by crystallogrphy 05 Jul '12

05 Jul '12

Hi, I want to refine my structure including a sugar oligomer as the ligand in Phenix GUI. I first put each subunit into the density and merge them together as the oligomer in Coot. Then I merged the oligomer coordinates to the protein file. Then I run phenix.elbow to get the cif file. Finally, I run phenix.refine. However, the oligomer units were pushed away from each other. I guess I need to modify the cif file to link the oligo units together. Does anyone know how I can I do it? thanks!

2 1

Re: [phenixbb] phenixbb Digest, Vol 80, Issue 1
by Matthew J. Whitley 05 Jul '12

05 Jul '12

Hello, One simple thing that might be the problem: As I'm sure you know, oligosaccharides are formed via intermolecular dehydration reactions. If I understand you correctly, you built the oligosaccharide by loading the component monomers into the density in Coot and then joined them together. Whenever two sugars are joined together, three atoms are given off as a water molecule. Did you actually delete the appropriate oxygen and two hydrogens from your coordinates for each glycosidic linkage in your structure? If you didn't, then it doesn't matter whether you told Coot or Phenix that the monomers or linked, you'll still have atoms trying to occupy the same space and thus the refinement will blow them apart. Not that I have EVER made that mistake myself . . . Similar topics have come up on the Phenix BB before. For other possibly useful information, see: http://phenix-online.org/pipermail/phenixbb/2012-April/018390.html http://www.phenix-online.org/pipermail/phenixbb/2009-October/013951.html Good luck! Matthew --- Matthew J. Whitley, Ph.D. Postdoctoral Research Fellow Angela Gronenborn Lab Department of Structural Biology University of Pittsburgh School of Medicine On Wed, July 4, 2012 7:55 pm, phenixbb-request(a)phenix-online.org wrote: > > > Message: 14 > Date: Wed, 4 Jul 2012 15:28:09 -0400 > From: crystallogrphy <qunwan1(a)gmail.com> > To: PHENIX user mailing list <phenixbb(a)phenix-online.org> > Subject: [phenixbb] ask refinement question including an aligo as the > ligand > Message-ID: > <CANcEHNi95ss+Tst8+8JVTK3=B598HjC_NNg5N7_=6-whdMTXpQ(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Hi, > > I want to refine my structure including a sugar oligomer as the ligand in > Phenix GUI. > > I first put each subunit into the density and merge them together as the > oligomer in Coot. Then I merged the oligomer coordinates to the protein > file. Then I run phenix.elbow to get the cif file. Finally, I run > phenix.refine. However, the oligomer units were pushed away from each > other. I guess I need to modify the cif file to link the oligo units > together. Does anyone know how I can I do it? > > thanks!

1 0

need subscription
by Afshan Begum 03 Jul '12

03 Jul '12

Dear Sir I have subscribe successfully to the phenixbb but i have not receive or involved in the discussion board could you kindly help ? Best Regards AFSHAN

1 0

Question about mr-sad
by Edward A. Berry 02 Jul '12

02 Jul '12

in doing SAD phasing with the option: input_partpdb_file=MR.pdb Is it assumed MR.pdb is already located in the cell, or is it a model to be used in MR assisted by information from the heavy atom sites? Thanks, eab

3 2

Re: [phenixbb] Default anisotropic B-factor refinement in phenix.refine?
by Filip YABUKARSKI 02 Jul '12

02 Jul '12

Hello again, First thank you very much for the response But in this case after the second run of phenix.refine i will lose the TLS contribution and the anisoU lines (?). How am I supposed to continue with TLS then (given it lowers my R factors i would be happy to keep it) and be sure that i do not refine anisotropic ADPs? Or you mean that i can do the 2 consecutive runs (with and without TLS) in the end of the refinement? Sorry if I got something wrong Filip On 2012-07-02 18:58, Pavel Afonine wrote: > Hi Filip, > > phenix.refine manual may need an update.. > > At 3A resolution if you run two consecutive refinements first with > TLS and the next one without TLS then in the second run anisotropic > ADPs will be converted to isotropic automatically. This is done by > phenix.refine to prevent accidental refinement of individual > anisotropic ADPs in such cases. The keyword and threshold are: > switch_to_isotropic_high_res_limit=1.5. > > Pavel > > On 7/2/12 9:47 AM, Filip YABUKARSKI wrote: >> Hello, >> >> I recently started to refine a structure with phenix.refine for a >> first time. It is a 3A dataset with NCS of 8. For the refinement i >> wanted to use TLS groups as well as individual B-factors >> (data/parameter ratio should be fine with ncs of 8) and after several >> cycles of refinement i get the Rfac of 20 and Rfree of 24 which is >> quite good (i also use simulated annealing, Secondary structure >> restraints and NCS restraints). Re-reading the refinement manual I >> realized that by using a PDB file which comes from a phenix.refine >> with TLS in the next refinement i am probably doing an anisotropic >> B-factor refinement for each atom, which in my case is not justified. >> I found a post on ccp4bb from 19 Fev 2009 discussing this issue but >> for me it is still not very clear. By reading the PDB i got from >> phenix.refine (version 1.7.3) it is not very clear (at least for me) >> if the refinement was with TLS groups and individual isotropic >> B-factors or with anisotropic B-factors. Basically my question is: if >> the input PDB has already anisoU records from a previous phenix.refine >> with TLS groups is the new phenix.refine going to use isotropic >> individual B-factors or anisotropic ones? (Even if in my case one can >> assume that i can probably use anisotropic B-factor refinement (24000 >> unique reflections with ncs of 8 vs. 6000 atoms in total) this is >> rarely the case for non-atomic resolution datasets) I would be very >> happy if someone gives an explanation. Many thanks >> >> Cheers >> Filip

2 1

Default anisotropic B-factor refinement in phenix.refine?
by Filip YABUKARSKI 02 Jul '12

02 Jul '12

Hello, I recently started to refine a structure with phenix.refine for a first time. It is a 3A dataset with NCS of 8. For the refinement i wanted to use TLS groups as well as individual B-factors (data/parameter ratio should be fine with ncs of 8) and after several cycles of refinement i get the Rfac of 20 and Rfree of 24 which is quite good (i also use simulated annealing, Secondary structure restraints and NCS restraints). Re-reading the refinement manual I realized that by using a PDB file which comes from a phenix.refine with TLS in the next refinement i am probably doing an anisotropic B-factor refinement for each atom, which in my case is not justified. I found a post on ccp4bb from 19 Fev 2009 discussing this issue but for me it is still not very clear. By reading the PDB i got from phenix.refine (version 1.7.3) it is not very clear (at least for me) if the refinement was with TLS groups and individual isotropic B-factors or with anisotropic B-factors. Basically my question is: if the input PDB has already anisoU records from a previous phenix.refine with TLS groups is the new phenix.refine going to use isotropic individual B-factors or anisotropic ones? (Even if in my case one can assume that i can probably use anisotropic B-factor refinement (24000 unique reflections with ncs of 8 vs. 6000 atoms in total) this is rarely the case for non-atomic resolution datasets) I would be very happy if someone gives an explanation. Many thanks Cheers Filip

2 1

Postdoctoral Fellow & Scientific Officer / Higher Scientific Officer Positions at ICR, London, UK
by Sebastian Guettler 30 Jun '12

30 Jun '12

Dear PHENIX Community, We are currently looking for two motivated researchers to join our new, exciting research programme on the regulation and functions of ADP-ribosyltransferases: one Scientific Officer / Higher Scientific Officer and one Postdoctoral Fellow. Please see the job descriptions for the two positions and instructions on how to apply below. Best wishes, Sebastian Postdoctoral Fellow & Scientific Officer / Higher Scientific Officer Positions at ICR, London, UK In the Divisions of Structural and Cancer Biology of the Institute of Cancer Research (ICR; http://www.icr.ac.uk), our new team led by Dr. Sebastian Guettler has a vacancy for a Scientific Officer or early-career Higher Scientific Officer and a Post-doctoral Training Fellow. Research in the Divisions of Structural and Cancer Biology focuses on the functional, biochemical, and structural characterisation of cellular processes relevant to cancer and cancer therapy. The Division of Structural Biology, where the laboratory will be physically based, has managed facilities for protein crystallography (crystallisation robots, x-ray diffractions systems), cryo-electron microscopy and protein production with expertise in multi-subunit expression systems (bacterial, yeast and insect cell expression, including a 60-litre fermentor). The Division is also well equipped with instrumentation for biophysical analysis (e.g., ITC, fluorescence-based methods, multi-angle light scattering). The Division of Cancer Biology provides a state-of-the-art infrastructure for mammalian cell culture, imaging, mammalian genetics, chemical biology and proteomics. Position 1: Post-doctoral Training Fellow (job reference number 1257974) Closing date of vacancy: 6 July 2012 Starting date of position: mid-October / early November 2012 Division: Structural Biology & Cancer Biology Team: Guettler Group Type of contract: fixed term Length of contract: 3 years in the first instance Salary range: £27,536 to £33,852 p.a. inclusive (depending on previous post-doctoral experience) Work location: London (Chester Beatty Laboratories, Chelsea) The successful applicant will undertake biochemical, crystallographic and functional studies on members of the ADP-ribosyltransferase (or “PARP”) family of proteins and their regulators. Our goal is to gain insight into the control and cellular roles of these enzymes with a specific interest in ADP-ribosyltransferase-dependent signalling processes linked to cancer. Applicants should possess a PhD (or equivalent) in biochemistry or molecular biology and a sound knowledge of protein purification for structural analyses. Experience with recombinant DNA techniques, protein expression and purification as well as protein crystallography is essential. Experience in mammalian cell culture, enzymology or protein production in insect cells (MultiBac) would be of advantage. For further details and instructions on how to apply, please refer to our online vacancy and recruitment site at http://www.icr.ac.uk/jobsearch, quoting job reference number 1257974 for this position. While applications need to be placed via our online system, informal inquiries can be sent to Dr. Sebastian Guettler at sebastian.guettler(a)icr.ac.uk. Position 2: Scientific Officer / Higher Scientific Officer (job reference number 1260528) Closing date of vacancy: 20 July 2012 Starting date of position: mid-October / early November 2012 Division: Structural Biology & Cancer Biology Team: Guettler Group Type of contract: fixed term Length of contract: 3 years in the first instance Salary range: £23,000 to £29,545 p.a. inclusive (depending on skills and experience) Work location: London (Chester Beatty Laboratories, Chelsea) The objective of this post will be to provide support within Dr. Guettler's laboratory for research into the structure and function of proteins. The post holder will work as a member of a multidisciplinary team investigating the functions and regulation of ADP-ribosyltransferases. The work will involve the production and purification of proteins and protein complexes, their characterisation, crystallisation and functional analysis. Applicants should minimally possess a BSc (or equivalent) in molecular biology, biochemistry or cell biology. Experience in molecular cloning, recombinant protein expression in E. coli and large-scale, crystallisation-grade protein purification is essential. Experience with protein production in insect cells (MultiBac) and mammalian cell culture would be of advantage. For further details and instructions on how to apply, please refer to our online vacancy and recruitment site at http://www.icr.ac.uk/jobsearch, quoting job reference number 1260528 for this position. While applications need to be placed via our online system, informal inquiries can be sent to Dr. Sebastian Guettler at sebastian.guettler(a)icr.ac.uk. The Institute of Cancer Research: Royal Cancer Hospital, a charitable Company Limited by Guarantee, Registered in England under Company No. 534147 with its Registered Office at 123 Old Brompton Road, London SW7 3RP. This e-mail message is confidential and for use by the addressee only. If the message is received by anyone other than the addressee, please return the message to the sender by replying to it and then delete the message from your computer and network.

1 0

k_sol and b_sol
by Luca Pellegrini 28 Jun '12

28 Jun '12

Hello, Where do I find the minimised values of k_sol and b_sol in phenix.refine? For some reason, they are not included in REMARK 3 of the pdb file after refinement, and I can't find them in the log file. Thanks, Luca Luca Pellegrini Department of Biochemistry University of Cambridge 80 Tennis Court Road Cambridge CB2 1GA - UK Email: lp212(a)cam.ac.uk Tel: 0044-1223-760469 Fax: 0044-1223-766002 Sanger building, room 3.59

2 1

Incorrect solvent modelling in Phenix.refine?
by David Briggs 28 Jun '12

28 Jun '12

Hi all, I'm refining a couple of structures of the same protein (well, point mutants) and I am consistently seeing large (6sig) negative difference map peaks within some cavities within my protein, present in both the datasets. As far as I understand, this is because the bulk solvent scaling parameters are either incorrect, or being incorrectly applied. I have tried a few things like applying the "Refine solvent mask (slow)" options, all to no avail. Resolution is 2.2Å and 2.5Å. Experimental phasing, Xtriage reports no abnormalities in the data. Is there anything that I am missing? Do I need to worry about this? Regards, Dave ============================ David C. Briggs PhD Father, Structural Biologist and Sceptic ============================ University of Manchester E-mail: david.c.briggs(a)manchester.ac.uk ============================ Webs : http://flavors.me/xtaldave Twitter: @xtaldave Skype: DocDCB ============================

2 1

Error in phenix.autobuild
by Tanner, John J. 27 Jun '12

27 Jun '12

Can anyone shed light on the error below? The situation is that I have a ~90% complete model from refinement (r_work = 0.2415, r_free = 0.3167), and I want phenix.autobuild to build the rest of the structure. I issued this command: phenix.autobuild data=../ccp4/DrPRODH3-6_1_scala2.6Ang.mtz model=Refine_9/DrPRODH_refine_9.pdb seq_file=../seq/DrPRODH-noMet1.fasta rebuild_in_place=False nproc=3 I got this error: [tannerjj@europa AutoBuild_run_2_]$ tail -n20 AutoBuild_run_2_1.log Output space group: P 1 21 1 (No. 4) Columns of data in AutoBuild_run_2_/refinement_PHX.mtz: H K L FreeR_flag F_DrPRODH3-6 SIGF_DrPRODH3-6 Unit cell: (46.806, 95.875, 70.504, 90, 107.733, 90) Space group: P 1 21 1 (No. 4) Using AutoBuild_run_2_/refinement_PHX.mtz for refinement Correcting data file /titan/tanner/DrPRODH/ccp4/DrPRODH3-6_1_scala2.6Ang.mtz for anisotropy and writing to AutoBuild_run_2_/aniso_data_PHX.mtz Getting overall aniso B from datafile... *************ERROR ENDING ******************* **************************************** AutoBuild Input failed iotbx InternalError: /net/longnose/scratch2/phenix/phenix-1.7.3-928/cctbx_project/iotbx/mtz/dataset.cpp(97): IOTBX_ASSERT(strlen(label) < sizeof(column_ptr->label)) failure. **************************************** *************ERROR ENDING *******************

3 4

occupancy group definition
by Christian Roth 26 Jun '12

26 Jun '12

Dear All, I have tried an occupancy refinement for a partially modifed residue. I have split the residue in altloc A and B and the modification (PMSF) is attached on B. I included the PMS moiety as part of the residue and gave them altloc B. I created the cif link. The refinement resulted in a strange mixture of occupancys for PMS around 0.8 and both Serins 0.23 an 0.77. So PMS is not part of the residue. I am not sure about the correct syntax for the constrained group. At the moment I have used constrained_group { selection = chain A and resseq 130 and altloc A selection = chain A and resname PMS and altloc B } However the Serine pat of altloc B refines to a value not compatible to the PMS values or the Altloc A value of the serine. How do I incorporate the altloc B for the serine? Thanks in advance for your help Christian

2 3

geometry weight help
by 周文昌 25 Jun '12

25 Jun '12

Dear all, Last week, I played a little bit with geometry weight in phenix refinement. When I check the log file, what I found is that if I use wxc_scale and fix_wxc in my refinement, actually phenix will ignore the wxc_scale, while if I only use wxc_scale, phenix will follow the value specified in the command line, but the wxc value will change from one macro cycle to another cycle. Below shows one of my examples in which the wxc_scale values are ignored when I use wxc_scale and fix_wxc at the same time. Anyone has similar problems? Thanks for all your considerations. ################## Command line arguments: "../new-unique.mtz" "../3tx8.pdb" "fix_rotamers=true" "strategy=individual_adp+individual_sites" "wxc_scale=0.01" "refinement.input.xray_data.labels=FP,SIGFP" "fix_wxc=10" "--overwrite" refinement.target_weights.wxc_scale = 0.01 refinement.target_weights.fix_wxc = 10 wxc_scale = 0.01 fix_wxc = 10 | XYZ refinement: T = Eexperimental * wxc * wxc_scale + Echem * wc | | wxc = 10.000000 wxc_scale = 1.000 wc = 1.000 | ... ################## Wenchang

3 2

Diffractive Imaging Postdoctoral Scientist Opportunity
by Nicholas Sauter 25 Jun '12

25 Jun '12

I'd like to alert subscribers to an exciting postdoctoral scientist opportunity, using diffractive imaging to examine organelle-scale biological structures at the Linac Coherent Light Source. The job posting is listed at: http://academicemployment.chance.berkeley.edu/DetailsJobSearch.cfm?recordID… The advertised position is strongly oriented toward the development of computational methods. The successful candidate will join the group of Cheryl Kerfeld at UC Berkeley, and collaborate daily with scientists at LBNL. Informal inquiries may be directed to Nick Sauter, nksauter at lbl.gov Peter Zwart, phzwart at lbl.gov -- Nicholas K. Sauter, Ph. D. Computer Staff Scientist, Physical Biosciences Division Lawrence Berkeley National Laboratory 1 Cyclotron Rd., Bldg. 64R0121 Berkeley, CA 94720-8118 (510) 486-5713

1 0

dihedral restraints problem with alternate conformations
by Benoît Masquida 25 Jun '12

25 Jun '12

I am having an issue refining 2 superimposed MES molecules A and B (morpholino ethane sulfonic acid, buffer seen in the density). phenix.refine works fine up to the end where it crashes because of the following error: Duplicate dihedral restraints: tor id: var_10 residue: MES, conformer "A" atoms: "HETATM 3993 O1 AMES T 1 .*. O " "HETATM 3994 C2 AMES T 1 .*. C " "HETATM 3995 C3 AMES T 1 .*. C " "HETATM 3996 N4 AMES T 1 .*. N+1" I can't see where is the conflict since all atoms have distinct labels and I provided the CIF file generated with phenix.elbow. What should I change to circumvent this problem? Thanks Ben

2 1

mr_rosetta problem
by Andreas Forster 24 Jun '12

24 Jun '12

Dear all, especially Tom T and Laurie Betts: I've recently installed phenix-dev-1084 and Rosetta 3.4 on RHEL 6.2. No problem. When I test the rosetta commands with phenix_regression.wizards.test_command_line_rosetta_quick, the first 12 of 39 tests succeed but test_prerefine fails and the process aborts. This is the error in the log file: /csb/soft/Linux64/share/rosetta3.4/rosetta_source/bin/mr_protocols.default.linuxgccrelease: error while loading shared libraries: libdevel.so: cannot open shared object file: No such file or directory libdevel.so exists, in the right place I venture to say: ls -l /csb/soft/Linux64/share/rosetta3.4/rosetta_source/build/src/release/linux/2.6/64/x86/gcc/4.4/libdevel.so -rwxr-xr-x. 1 csbadmin csb 6437 Jun 23 18:44 /csb/soft/Linux64/share/rosetta3.4/rosetta_source/build/src/release/linux/2.6/64/x86/gcc/4.4/libdevel.so Besides the tests failing, the actual mr_rosetta run fails as well, as one would expect, with the same "missing library" error. This problem has been reported before (http://www.phenix-online.org/pipermail/phenixbb/2011-May/017157.html), but I can't find a solution. Is there one? Thanks. Andreas

3 3

Phenix version 1.8 released
by Paul Adams 22 Jun '12

22 Jun '12

The Phenix developers are pleased to announce that version 1.8 of Phenix is now available. Binary installers for Linux, Mac OSX, and Windows 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: phenix.multi_crystal_average, phenix.morph_model, phenix.table_one, phenix.find_peaks_holes - Improved layout for AutoBuild and reflection file editor GUIs - Editor for custom bond/angle/plane restraints in phenix.refine - Project groups: manage defaults for similar projects General: ======== - New very fast bulk-solvent and overall anisotropic scaling procedure - Experimental Windows installer (32-bit only) - Improved support for anomalous difference arrays (F SIGF DANO SIGDANO ISYM) - New commands: phenix.den_refine, phenix.find_peaks_holes - Overall Molprobity score now calculated in phenix.model_vs_data and the validation/refinement GUIs Refinement: =========== - Change in default NCS restraints behavior: - use torsion-angle parameterization - restraints on NCS-related B-factors turned off - Automatic matching of chains to reference model; allows multiple chains to be restrained to a single reference chain - Support for multiple reference model files - New method for filling of missing Fobs using simple density modification Molecular Replacement: ====================== - phaser (version 2.5): - TNCS support enabled by default for both MR and SAD Experimental Phasing and Model Building: ======================================== - AutoSol: - TNCS support for SAD phasing Miscellaneous: ============== - phenix.pdbtools: - ability to set atomic charge field (used for X-ray scattering factors and VDW radii in nonbonded restraints) - phenix.find_peaks_holes (new command): - standalone map peak search program: highlight difference map features - includes analysis of anomalous maps, flag suspicious "water" molecules - phenix.maps: - simplified execution on command line: "phenix.maps model.pdb data.mtz" - generates default map coefficients (2mFo-DFc, mFo-DFc, anomalous if available) and 2mFo-DFc CCP4-format map. - phenix.fetch_pdb: - added --maps option: generates 2mFo-DFc and mFo-DFc map coefficients - phenix.get_cc_mtz_pdb: - added atom_selection keyword: compute CC for a subset of atoms - phenix.real_space_correlation: - simplified execution on command line: "phenix.real_space_correlation model.pdb data.mtz" - calculates CC of model to 2mFo-DFc map - added atom_selection keyword: compute CC for a single selection (instead of per-atom or per-residue statistics) 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 ] --

2 2

Run time error while refining
by Subhani Bandara 22 Jun '12

22 Jun '12

Hi everyone, I got the following error message while running refinement in mac, using PHENIX-1.8-1069. This occured when refined after editing Ramachandran outliers in coot, without applying Ramachandran restraints. Not sure wether it has something to do with this. RuntimeError : Estimated memory allocation for cubicles exceeds max_number_of_bytes: This may be due to unreasonable parameters: cubicle_edge=5 space_span=(3528.41,3528.41,3442.53) n_cubicles=(706,706,689) max_number_of_bytes=2147483648 Any suggestions! Thanks Subhani

2 2

error TLS groups
by Christian Roth 22 Jun '12

22 Jun '12

Dear Phenix developers, If I want use the automatic TLS definition the following error appears. At least the first error line with invalid iteral does also appear int the last verification step of the validation against the supplied sequence file. I have a PMS modified Serine could that cause the problem, because with other pdb's it works. Best Regards Christian RuntimeError : Error in subprocess! Original error: invalid literal for int() with base 10: '' Original traceback: File "/userdata-2/christian/phenix-1.8-1069/cctbx_project/libtbx/thread_utils.py", line 233, in run return_value = self._target(self._args, self._kwargs, self._c) File "/userdata-2/christian/phenix-1.8-1069/cctbx_project/mmtbx/command_line/find_tls_groups.py", line 901, in __call__ xray_structure=self.xray_structure) File "/userdata-2/christian/phenix-1.8-1069/cctbx_project/mmtbx/command_line/find_tls_groups.py", line 729, in find_tls ssm.find_automatically(log=out) File "/userdata-2/christian/phenix-1.8-1069/cctbx_project/mmtbx/secondary_structure/__init__.py", line 334, in find_automatically self.sec_str_from_pdb_file = self.find_sec_str(log=log) File "/userdata-2/christian/phenix-1.8-1069/cctbx_project/mmtbx/secondary_structure/__init__.py", line 381, in find_sec_str allow_none=True) File "/userdata-2/christian/phenix-1.8-1069/cctbx_project/iotbx/pdb/secondary_structure.py", line 425, in process_records secondary_structure = annotation(records=records) File "/userdata-2/christian/phenix-1.8-1069/cctbx_project/iotbx/pdb/secondary_structure.py", line 62, in __init__ self.helices = parse_helix_records(records) File "/userdata-2/christian/phenix-1.8-1069/cctbx_project/iotbx/pdb/secondary_structure.py", line 217, in parse_helix_records helix_class=string.atoi(line[38:40]), File "/userdata-2/christian/phenix-1.8-1069/build/intel- linux-2.6/base/lib/python2.7/string.py", line 403, in atoi return _int(s, base) Traceback (most recent call last): File "/userdata-2/christian/phenix-1.8-1069/phenix/wxGUI2/Programs/Refine/TLS.py", line 140, in OnFindTLS self.find_tls_groups() File "/userdata-2/christian/phenix-1.8-1069/phenix/wxGUI2/Programs/Refine/TLS.py", line 103, in find_tls_groups callback=None) File "/userdata-2/christian/phenix-1.8-1069/cctbx_project/wxtbx/process_control.py", line 348, in run_function_as_process_in_dialog dlg.handle_error() File "/userdata-2/christian/phenix-1.8-1069/cctbx_project/wxtbx/process_control.py", line 302, in handle_error %s""" % (str(exception), traceback)) RuntimeError: Error in subprocess! Original error: invalid literal for int() with base 10: '' Original traceback: File "/userdata-2/christian/phenix-1.8-1069/cctbx_project/libtbx/thread_utils.py", line 233, in run return_value = self._target(self._args, self._kwargs, self._c) File "/userdata-2/christian/phenix-1.8-1069/cctbx_project/mmtbx/command_line/find_tls_groups.py", line 901, in __call__ xray_structure=self.xray_structure) File "/userdata-2/christian/phenix-1.8-1069/cctbx_project/mmtbx/command_line/find_tls_groups.py", line 729, in find_tls ssm.find_automatically(log=out) File "/userdata-2/christian/phenix-1.8-1069/cctbx_project/mmtbx/secondary_structure/__init__.py", line 334, in find_automatically self.sec_str_from_pdb_file = self.find_sec_str(log=log) File "/userdata-2/christian/phenix-1.8-1069/cctbx_project/mmtbx/secondary_structure/__init__.py", line 381, in find_sec_str allow_none=True) File "/userdata-2/christian/phenix-1.8-1069/cctbx_project/iotbx/pdb/secondary_structure.py", line 425, in process_records secondary_structure = annotation(records=records) File "/userdata-2/christian/phenix-1.8-1069/cctbx_project/iotbx/pdb/secondary_structure.py", line 62, in __init__ self.helices = parse_helix_records(records) File "/userdata-2/christian/phenix-1.8-1069/cctbx_project/iotbx/pdb/secondary_structure.py", line 217, in parse_helix_records helix_class=string.atoi(line[38:40]), File "/userdata-2/christian/phenix-1.8-1069/build/intel- linux-2.6/base/lib/python2.7/string.py", line 403, in atoi return _int(s, base) -- Christian Roth Institut für Bioanalytische Chemie Biotechnologisch-Biomedizinisches Zentrum Fakultät für Chemie und Mineralogie Universität Leipzig Deutscher Platz 5 04103 Leipzig Telefon: +49 (0)341 97 31316 Fax: +49 (0)341 97 31319

2 1

phenix.mr_model_preparation
by Alexander Batyuk 22 Jun '12

22 Jun '12

Dear all, I'm trying to run phenix.mr_model_preparation in phenix 1.8-1069 with a .hhr file but it fails with Reading homology search file: 4671015.hhr Sorry: Uninterpretable: 'Q ss_pred cccccchHHHHHHHHHHHHHHHHHhhhchheEEE etc. Could somebody point me to the right direction, please? Thank you and best wishes, Alex -- Alex Batyuk The Plueckthun Lab www.bioc.uzh.ch/plueckthun

2 1

phenix.refine and torsion-angle ncs
by Antony Oliver 21 Jun '12

21 Jun '12

Dear phenix.bb How do I specifically exclude particular amino acids from torsion-angle ncs? >From the GUI, it is not obvious. In a couple of my chains, one or two amino acids are in definite crystal contacts, and are therefore in a different rotamer/conformation compared to the other ncs-related copies. Unfortunately, I can't seem to find the right option / keyword to stop the refinement job moving these side-chains to agree with the rotamers found in the other molecules of the asymmetric unit. Many thanks, Tony. --- Dr Antony W Oliver Senior Research Fellow CR-UK DNA Repair Enzymes Group Genome Damage and Stability Centre Science Park Road University of Sussex Falmer, Brighton, BN1 9RQ email: antony.oliver(a)sussex.ac.uk<mailto:[email protected]> tel (office): +44 (0)1273 678349 tel (lab): +44 (0)1273 677512

3 4

Re: [phenixbb] bad geometry problem
by Qian Sarah 21 Jun '12

21 Jun '12

> > > Hi Pavel Thank you so much! The information is VERY detailed and I learnt a lot! Thanks to phenixbb! > > 2012/6/21 Pavel Afonine <pafonine(a)lbl.gov> > >> Hi Sarah, >> >> thanks for sending me the data. I spent 15 minutes trying a few things >> and here is the statistics I get: >> >> r_work = 0.2454 r_free = 0.2700 >> MOLPROBITY STATISTICS. >> ALL-ATOM CLASHSCORE : 13.60 >> RAMACHANDRAN PLOT: >> OUTLIERS : 0.89 % >> ALLOWED : 5.33 % >> FAVORED : 93.78 % >> ROTAMER OUTLIERS : 8.64 % >> CBETA DEVIATIONS : 0 >> >> which isn't that bad compared to the starting point: >> >> r_work = 0.2469 r_free = 0.2613 >> MOLPROBITY STATISTICS. >> ALL-ATOM CLASHSCORE : 48.87 >> RAMACHANDRAN PLOT: >> OUTLIERS : 2.49 % >> ALLOWED : 5.33 % >> FAVORED : 92.18 % >> ROTAMER OUTLIERS : 29.47 % >> CBETA DEVIATIONS : 10 >> >> Here is what I did (it's just one of possible ways of doing this; other >> scenarios may works as good as this one or better!): >> >> - add H atoms: >> >> phenix.reduce model.pdb > model_H.pdb >> >> - then run a quick geometry regularization: >> >> phenix.pdbtools model_H.pdb --geom >> >> - then strip off H: >> >> phenix.reduce -trim model_H_regularized.pdb > start.pdb >> >> - then run refinement using NCS and Ramachandran plot restraints: >> >> phenix.refine start.pdb def_data.mtz main.ncs=true --overwrite >> output.prefix=def2 main.number_of_mac=5 ramachandran_restraints=true >> rama_potential=emsley strategy=individual_sites+**individual_adp >> >> - then take outcome of previous refinement and run the same refinement >> plus weights optimization: >> >> phenix.refine def2_001.pdb def_data.mtz main.ncs=true --overwrite >> output.prefix=def3 main.number_of_mac=5 ramachandran_restraints=true >> rama_potential=emsley strategy=individual_sites+**individual_adp >> optimize_xyz=true nproc=6 >> >> Also, looking at the original maps I see a lot of places where simple >> fixing in Coot is just obvious, for both - rotamers and Ramachandran >> outliers. So you should do it first of all. >> >> Also, use Ramachandran restraints ONLY to preserve good starting model >> and NOT to fix problems. >> >> All in all, geometry problems you were concerned about are all easily >> fixable. One of the reason why these problems arise is that the map is >> really looking like 3A resolution map, may be even worse; plus there are >> places in the model that you need to fix manually (refinement would not be >> able to fix them). >> >> I will send you the files off-list. >> >> >> Pavel >> >> On 6/19/12 1:17 AM, Qian Sarah wrote: >> >>> Hi, this may not be a new problem, but I tried every stretagy that I >>> know from bb, but the geometry is still very bad. So I wish I can get some >>> help from here. >>> Crystal information: >>> SAD data. Resolution 3.0. There is no anisotropy in the data. >>> R/Rfree : 22.7/26.8 (Rfree-flag is always from the Refine_data.mtz file >>> generated after the first refinement.) >>> RMS(bonds):0.02 >>> RMS(angles): 1.453 >>> Ramachandran outliers: 3.0% >>> Ramachandran favored: 91.2% >>> Rotamer outliers: 18.7% !! >>> Clashscore: 31.84 >>> overall score: 3.47 >>> I tried different strategies, like autobuild and then refine; set >>> wxu/wxc=0.3/0.3; turn on "Ramachandran restrains", "fix bad sidechain >>> romaters"etc, however the paramters are still very bad, especially the >>> rotamer outliers. Then I rebuild them manually in coot,after refinement the >>> parameters are even worse. >>> p.s. phenix version: phenix-1.7.3-dev-1046. >>> Thank you very much! >>> best, >>> Sarah >>> >> >> >

1 0

fixing weights
by Morten Groftehauge 21 Jun '12

21 Jun '12

Hi all, I tried specifying wxc, wxu, wxu_scale and wxc_scale as you can see from the log file: target_weights { wxc_scale = 0.50001 wxu_scale = 0.49999 fix_wxc = 6 fix_wxu = 0.35 } but the refinement proceeds with the default values for _scale: ========================== Target weights: x-ray data ========================= |-----------------------------------------------------------------------------| | XYZ refinement: T = Eexperimental * wxc * wxc_scale + Echem * wc | | wxc = 6.000000 wxc_scale = 1.000 wc = 1.000 | | angle between x-ray and geometry gradient vectors: None (deg) | | | | ADP refinement: T = Eexperimental * wxu * wxu_scale + Eadp * wu | | wxu = 0.350000 wxu_scale = 0.500 wu = 1.000 | | angle between Xray and ADP gradient vectors: None (deg) | |-----------------------------------------------------------------------------| I am running dev-1082, used the GUI for phenix.refine. Cheers, Morten -- Morten K Grøftehauge, PhD Pohl Group Durham University

3 3

R-factor calculation
by Abhinav Kumar 20 Jun '12

20 Jun '12

Hi, Is it possible to specify a different high resolution cutoff (than what's in data file) in phenix.model_vs_data for R-factors calculation? -- Thanks, Abhinav JCSG@SSRL, SLAC (650) 926-2992

2 4

bad geometry problem
by Qian Sarah 20 Jun '12

20 Jun '12

Hi, this may not be a new problem, but I tried every stretagy that I know from bb, but the geometry is still very bad. So I wish I can get some help from here. Crystal information: SAD data. Resolution 3.0. There is no anisotropy in the data. R/Rfree : 22.7/26.8 (Rfree-flag is always from the Refine_data.mtz file generated after the first refinement.) RMS(bonds):0.02 RMS(angles): 1.453 Ramachandran outliers: 3.0% Ramachandran favored: 91.2% Rotamer outliers: 18.7% !! Clashscore: 31.84 overall score: 3.47 I tried different strategies, like autobuild and then refine; set wxu/wxc=0.3/0.3; turn on "Ramachandran restrains", "fix bad sidechain romaters"etc, however the paramters are still very bad, especially the rotamer outliers. Then I rebuild them manually in coot,after refinement the parameters are even worse. p.s. phenix version: phenix-1.7.3-dev-1046. Thank you very much! best, Sarah

3 3

Source code build - overriding compiler choice
by Graeme Winter 18 Jun '12

18 Jun '12

Hi Folks, On my macbook with lion and a source code build, I always set --compiler=gcc for cctbx and this works fine. Just tried a source code build of phenix which failed because it picked up clang, so how do I override this in the install script? I didn't see the option from the install script: Graemes-MacBook-Pro:phenix-installer-dev-1078 graeme$ ./install --help Usage: install [options] Options: --prefix=<path> Determines PHENIX installation path, which will be <path > /phenix-[version] --overwrite Forces overwriting of installed packages --nproc=number Number of processors to use in compilation --base-only Do not compile PHENIX sources (cctbx, phaser, resolve, etc.) --no-gui Do not compile and install GUI components --no-apps Don't create Mac app launcher icons --app-prefix=<path> Path for Mac app launchers (default: /Applications) --noopt Turn off optimization during compilation --debug Turn on debugging during compilation --openmp Try to use OpenMP if supported by compiler --force-compile Force compilation of source --binary Force a binary install if binaries and source are both available --source Force a source install if binaries and source are both available --no-packages Do not install other packages such as PHASER --alias-mtype=<string> Force installation from a different machine type e.g. --alias-mtype=intel-linux-2.4 [for binary installations only] --try-unsupported Try source installation on an unsupported platform --help Print this help message Thanks, Graeme

2 3

MR_Rosetta
by Heather Condurso 16 Jun '12

16 Jun '12

Hello all, I am working on a structure for which a good molecular replacement model doesn't really exist ( highest sequence similarity is 24% for a protein with completely different function and even worse for proteins with a more similar function). As I try and solve the phase with Se-Met derived protein, I thought I would try and use MR_rosetta and see what happens. I have also tried to use pyrosetta to generate better models, but the structure alignments suggest my protein to be primarily composed of beta sheets so the models from ab initio methods aren't great. I am using the GUI and have given my sca file, a sequence file, the 3 and 9 fragment files, and both the global and local alignment hhr files. I thought that phenix would use the hhr alignments to pull models directly from the pdb, but without giving a model pdb, MR_rosetta fails. So my question is how do I input the models? Should I download each of the pdbs and clean them up a bit? How many should I include? I unde! rstand the more I include the more computing time this will take. Can you give any recommendations as to other parameters to change keep the computing time down initially to identify the best model before extensive time is used? I appreciate any help you may be able to provide, Heather Condurso UF-Department of Chemistry condurso(a)ufl.edu

3 6

alternate conformation clashes
by Timothy Springer 14 Jun '12

14 Jun '12

We have a 1.4 Å structure with two pseudo-symmetry related A and B chains that form a super b-sheet. An 8-residue loop has two quite different alternate, backbone conformations. When we refine the occupancy of these loops, they come out to 0.51 and 0.49, consistent with them being mutually incompatible (i.e., forced to be 0.50 and 0.50 in the ideal situation). Indeed, clashes would exist between the A chain A alternate and the B chain A alternate, but when these are A and B alternates, there is no clash. Since these are in separate chains, Phenix has no problem with an A and A clash. There is a second loop, non-contiguous in sequence, that also has clashes with one of these. We refine each of the two clashing loops as the A conformat ion. Since these are separate loops, Phenix does not seem to care and seems to treat the A conformation of one loop as if it is compatible with the B conformation of the other loop. This is great that Phenix is so forgiving, but my concern is that shouldn't by convention we reassign the conformations so that the A conformations are all compatible within the same monomer? Is there some formal convention on this, and will PDB ask us to reassign these? In a closely related 1.7 Å lattice, we have a single molecule per asymmetric unit, and now the relationship is between two symmetry-related molecules. The two conformations of loops are similar as in the 1.4 Å structure. Here, we have to keep the same alternate conformation name, and thus have a clash between A and A loops 1 but not A and B loops 1. Again there is a clash between A conformation of loop 1 and A conformation of loop 2, and we could make this go away by switching A and B in loop 2. I find it to be less confusing if we switch the conformation names in loop 2. What is the consensus view? Tim

2 1

cif_as_mtz assertion error
by Eric Williams 14 Jun '12

14 Jun '12

I'm trying to run phenix.cif_as_mtz on a large number of PDB entries, and several entries give me an assertion error. The example below is from 8icz, but the error always seems to be the same. Traceback (most recent call last): File "/media/funky/phenix-1.8-1069/build/intel-linux-2.6/../../cctbx_project/mmtbx/command_line/cif_as_mtz.py", line 544, in <module> run(sys.argv[1:]) File "/media/funky/phenix-1.8-1069/build/intel-linux-2.6/../../cctbx_project/mmtbx/command_line/cif_as_mtz.py", line 135, in run incompatible_flags_to_work_set=command_line.options.incompatible_flags_to_work_set) File "/media/funky/phenix-1.8-1069/build/intel-linux-2.6/../../cctbx_project/mmtbx/command_line/cif_as_mtz.py", line 197, in process_files r_free_flags = r_free_flags).mtz_object() File "/media/funky/phenix-1.8-1069/cctbx_project/mmtbx/utils.py", line 2001, in __init__ skip_twin_detection = True) File "/media/funky/phenix-1.8-1069/cctbx_project/mmtbx/utils.py", line 2113, in get_r_factor twin_laws = twin_laws) File "/media/funky/phenix-1.8-1069/cctbx_project/mmtbx/utils.py", line 1556, in fmodel_simple assert f_obs.sys_absent_flags().data().count(True)==0 Any ideas? Thanks. :) Eric Williams PhD candidate Intelligent Systems Program University of Pittsburgh

3 8

run time crashing
by CPMAS Chen 14 Jun '12

14 Jun '12

It crashes during rotamer fix at the initial stage. If you like, I can send you the necessary file for you to debugging. > > Message: 7 > Date: Thu, 14 Jun 2012 08:24:44 -0700 > From: Pavel Afonine <pafonine(a)lbl.gov> > To: PHENIX user mailing list <phenixbb(a)phenix-online.org> > Subject: Re: [phenixbb] phenixbb Digest, Vol 79, Issue 4 > Message-ID: <4FDA023C.3070306(a)lbl.gov> > Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" > > Hi, > > how long does it take to reproduce this problem in your case? So far I > have examples that take forever before they crash. So if yours takes a > few minutes I would be very interested to have it for debugging. > > Also, Jeff is going to try something else that may eliminate this problem. > > Thanks, > Pavel > > *************************************************** Charles Chen Research Associate University of Pittsburgh School of Medicine Department of Anesthesiology ******************************************************

2 1

Re: [phenixbb] phenixbb Digest, Vol 79, Issue 4
by CPMAS Chen 14 Jun '12

14 Jun '12

I thought I am the only one have this problem. I have encountered the same problem with the newest version of phenix. My installation is on MacOS X 10.6.8. If more information needed, I am happy to provide. Traceback (most recent call last): File "/usr/local/phenix-1.8-1069/build/mac-intel-osx-x86_64/../../phenix/phenix/command_line/refine.py", line 11, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/usr/local/phenix-1.8-1069/phenix/phenix/refinement/command_line.py", line 146, in run call_back_handler=call_back_handler) File "/usr/local/phenix-1.8-1069/phenix/phenix/refinement/driver.py", line 1510, in run ncs_manager = self.ncs_manager) File "/usr/local/phenix-1.8-1069/phenix/phenix/refinement/strategies.py", line 575, in refinement_machine log=log) File "/usr/local/phenix-1.8-1069/cctbx_project/mmtbx/torsion_restraints/torsion_ncs.py", line 884, in fix_rotamer_outliers monitor_clashscore=False) File "/usr/local/phenix-1.8-1069/cctbx_project/mmtbx/refinement/real_space.py", line 336, in refine geometry_restraints_manager = geo_box) File "/usr/local/phenix-1.8-1069/cctbx_project/mmtbx/refinement/real_space.py", line 250, in __init__ lbfgs_exception_handling_params = lbfgs_exception_handling_params) File "/usr/local/phenix-1.8-1069/cctbx_project/cctbx/maptbx/real_space_refinement_simple.py", line 137, in __init__ exception_handling_params=lbfgs_exception_handling_params) File "/usr/local/phenix-1.8-1069/cctbx_project/scitbx/lbfgs/__init__.py", line 258, in run log) File "/usr/local/phenix-1.8-1069/cctbx_project/scitbx/lbfgs/__init__.py", line 126, in run_c_plus_plus f, g = target_evaluator.compute_functional_and_gradients() File "/usr/local/phenix-1.8-1069/cctbx_project/cctbx/maptbx/real_space_refinement_simple.py", line 193, in compute_functional_and_gradients compute_gradients=True) File "/usr/local/phenix-1.8-1069/cctbx_project/cctbx/geometry_restraints/manager.py", line 716, in energies_sites flags=flags, sites_cart=sites_cart, site_labels=site_labels) File "/usr/local/phenix-1.8-1069/cctbx_project/cctbx/geometry_restraints/manager.py", line 650, in pair_proxies min_cubicle_edge=self.min_cubicle_edge) File "/usr/local/phenix-1.8-1069/cctbx_project/cctbx/geometry_restraints/__init__.py", line 1295, in __init__ shell_asu_tables=shell_asu_tables) RuntimeError: Estimated memory allocation for cubicles exceeds max_number_of_bytes: This may be due to unreasonable parameters: cubicle_edge=5 space_span=(2696.14,2555.98,2650.54) n_cubicles=(540,512,531) max_number_of_bytes=1073741824 > > Message: 1 > Date: Mon, 11 Jun 2012 11:44:07 -0700 > From: Nathaniel Echols <nechols(a)lbl.gov> > To: PHENIX user mailing list <phenixbb(a)phenix-online.org> > Subject: Re: [phenixbb] runtime error in phenix.refine > Message-ID: > <CALeAa1PBCCzKPJTMoDsFbsZMSO=SMB_19ns+mfvC3HozuUcBaQ(a)mail.gmail.com > > > Content-Type: text/plain; charset=ISO-8859-1 > > On Sun, Jun 10, 2012 at 10:04 AM, Eugene Valkov <eugene.valkov(a)gmail.com> > wrote: > > I'm encountering the following during the second macro-cycle in > > phenix.refine in latest 1.8-1069 dev build for Linux 64-bit: > > > > RuntimeError: Excessive number of cubicles: > > ?This may be due to unreasonable parameters: > > ? ?cubicle_edge=5 > > ? ?space_span=(8467.78,8440.85,8372.67) > > ? ?n_cubicles=(1694,1689,1675) > > ? ?max_number_of_bytes=101542821888 > > > > Any thoughts/suggestions much appreciated. > > If I'm interpreting the code correctly, this means that the number of > cubicles will lead to an integer format overflow - which I find hard > to believe if you're using a 64-bit system. Could you please send me > the rest of the error message, along with the output of the command > "boost_adaptbx.show_platform_info"? > > thanks, > Nat > > > -- *************************************************** Charles Chen Research Associate University of Pittsburgh School of Medicine Department of Anesthesiology ******************************************************

3 2

Re: [phenixbb] phenixbb Digest, Vol 79, Issue 4
by Kay Diederichs 14 Jun '12

14 Jun '12

"A 64bit operating system" can only run on hardware that supports 64bits. Here, 64bits refers to addresses in memory - arrays larger than 4GB need 64bit addressing because 32bits do not allow high enough numbers to index such arrays. It has nothing to do with the maximum value of an integer - that is set by the programmer. A programmer can (depending on programming language and compiler) usually decide whether integer numbers should be of size 4 bytes or 8 bytes. It is just by coincidence that 8 bytes are 64 bits. Here it seems to me that a 8byte integer would be needed to hold the number of cubicles. best, Kay > Date: Mon, 11 Jun 2012 11:44:07 -0700 > From: Nathaniel Echols<nechols(a)lbl.gov> > To: PHENIX user mailing list<phenixbb(a)phenix-online.org> > Subject: Re: [phenixbb] runtime error in phenix.refine > Message-ID: > <CALeAa1PBCCzKPJTMoDsFbsZMSO=SMB_19ns+mfvC3HozuUcBaQ(a)mail.gmail.com> > Content-Type: text/plain; charset=ISO-8859-1 > > On Sun, Jun 10, 2012 at 10:04 AM, Eugene Valkov<eugene.valkov(a)gmail.com> wrote: >> I'm encountering the following during the second macro-cycle in >> phenix.refine in latest 1.8-1069 dev build for Linux 64-bit: >> >> RuntimeError: Excessive number of cubicles: >> ?This may be due to unreasonable parameters: >> ? ?cubicle_edge=5 >> ? ?space_span=(8467.78,8440.85,8372.67) >> ? ?n_cubicles=(1694,1689,1675) >> ? ?max_number_of_bytes=101542821888 >> >> Any thoughts/suggestions much appreciated. > > If I'm interpreting the code correctly, this means that the number of > cubicles will lead to an integer format overflow - which I find hard > to believe if you're using a 64-bit system. Could you please send me > the rest of the error message, along with the output of the command > "boost_adaptbx.show_platform_info"? > > thanks, > Nat

2 1

MR error when fix one component and search for the other one
by fn1＠rice.edu 13 Jun '12

13 Jun '12

Hi everyone, I am trying to fix one component from the previous MR search and try to search for the other one, euler and frac are taken from the pdb file of the previous run (MR.1.pdb). Here are my input file and error message. Thanks in advance! Fengyun automr { build=false ensemble { ensembleID = "mol1" copies_to_find = 1 coords = in.pdb RMS = 0.85 } ensemble { ensembleID = "mol2" copies_to_find = 0 coords = in.pdb RMS = 0.85 } component { seq_file="mut.seq" component_type = *protein nucleic_acid component_copies = 4 } fixed_ensembles { fixed_ensembleID_list="mol2" fixed_euler_list = 206.66 146.75 96.21 fixed_frac_list = 0.688 1.011 0.005 } } i get the error message like this, ENSEMBLE: mol1 , 1 PDB file(s) ENSEMBLE: mol2 , 1 PDB file(s) Resolution set to 3.5 Specifying composition through protein sequence file /shared.scratch/ni/mut-p222-all/bb-3.5/AutoMR_run_5_/r2/mut.seq 4 Adding Ensemble mol2 as fixed solution with: Euler angles: [206.66, 146.75, 96.21] Coords (orthogonal A): [0.688, 1.011, 0.005] *************ERROR ENDING ******************* **************************************** AutoMR Input failed Python argument types in SOLU.addSOLU_6DIM_ENSE(InputMR_AUTO, str, list, bool, list, float, bool, bool, bool) did not match C++ signature: addSOLU_6DIM_ENSE(phaser::SOLU {lvalue}, std::string, scitbx::vec3<double>, bool, scitbx::vec3<double>, double, bool, bool, bool, bool, scitbx::vec3<double>, double) **************************************** *************ERROR ENDING *******************

2 1

runtime error in phenix.refine
by Eugene Valkov 11 Jun '12

11 Jun '12

I'm encountering the following during the second macro-cycle in phenix.refine in latest 1.8-1069 dev build for Linux 64-bit: RuntimeError: Excessive number of cubicles: This may be due to unreasonable parameters: cubicle_edge=5 space_span=(8467.78,8440.85,8372.67) n_cubicles=(1694,1689,1675) max_number_of_bytes=101542821888 Any thoughts/suggestions much appreciated. Thanks, Eugene -- Dr Eugene Valkov Division of Structural Studies MRC Laboratory of Molecular Biology Hills Road, Cambridge, CB2 0QH Email: evalkov(a)mrc-lmb.cam.ac.uk Tel: +44 (0)1223 402261

3 3

How to run normal mode analysis under command mode?
by fn1＠rice.edu 11 Jun '12

11 Jun '12

Hi everyone, I'd like to perturb my model with normal mode analysis. Does phenix have a way to do this? Thanks in advance! Fengyun

6 11

Gordon conference - last chance!
by Anastassis Perrakis 11 Jun '12

11 Jun '12

Dear all, The registration deadline for the Gordon conference on diffraction methods in Structural Biology (15-20 July), is at the end of the week (June 17th). This is you last chance to register and get some of the places left!!! See the program and application information at: http://www.grc.org/programs.aspx?year=2012&program=diffrac Ana and Tassos

1 0