- phenixbb - phenix-online.org

Re: [phenixbb] Autobuild submitted jobs
by Mark A Saper 27 May '10

27 May '10

Tom -- It turns out that each of the omit region's jobs submits another job, but that latter job is only the "test run_command" job. Why test the run_command if no subsequent jobs are submitted? Our solution, thanks to an extremely savvy cluster manager, was to install the qsub (Torque?) package on each of the computing nodes (117). Also we ran the top-level process on the head (scheduler) node, so as not to tie up a computing node while it waits for its sibling jobs (the omit regions) to finish the calculation. Overall wall-time with a separate node allocated for each omit region was about 50 minutes. Interestingly, although each node has 8 processors, the OMIT region jobs seems to only use one processor! Interestingly when phenix.refine jobs are run on my Macbook Pro, the operating system is savvy enough to spread it over 2 processors. -Mark ______________________________ Mark A. Saper, Ph.D. Associate Professor of Biological Chemistry University of Michigan Biophysics, 3040 Chemistry Building 930 N University Ave Ann Arbor MI 48109-1055 U.S.A. saper(a)umich.edu phone (734) 764-3353 fax (734) 764-3323 http://www.biochem.med.umich.edu/?q=saper http://www.strucbio.umich.edu/

3 2

How to turn off mask optimization
by Schubert, Carsten [PRDUS] 26 May '10

26 May '10

Hi, Is there a way to turn off the mask optimization step in 393? I tried 'main.optimize_mask = False' but that does not work. The program spends a substantial amount of time in repeated mask optimizations and I am trying to shave off a couple minutes for each run... Cheers, Carsten

2 3

Autobuild submitted jobs
by Mark A Saper 26 May '10

26 May '10

Dear Tom and colleagues, We're running Autobuild composite omit on a large cluster, each node has 8 processors. Autobuild submits jobs with the "run_command". And it appears that each of those jobs submit other jobs. Are the latter jobs necessary? That is, are there situtations where they wouldn't be submitted as separate jobs but run as subprocesses? Currently nodes in our cluster don't have capability of qsubing new jobs, but we are in the process of changing this. Thanks, Mark _________________________________ Mark A. Saper, Ph.D. Associate Professor of Biological Chemistry University of Michigan Biophysics, 3040 Chemistry Building 930 N University Ave Ann Arbor MI 48109-1055 U.S.A. saper(a)umich.edu phone (734) 764-3353 fax (734) 764-3323 http://www.biochem.med.umich.edu/?q=saper http://www.strucbio.umich.edu/

2 1

atom selection syntax for ALT atoms?
by Scott Classen 25 May '10

25 May '10

This is the selection for the atom I'm interested in: atom_selection_2 = "name OE1 and chain A and resname GLU and resseq 160" I want to try a refinement with an alternate conformation at GLU160 . Can someone please help me with the proper syntax. Thanks, Scott ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Scott Classen, Ph.D. SIBYLS Beamline 12.3.1 http://bl1231.als.lbl.gov Advanced Light Source Lawrence Berkeley National Laboratory 1 Cyclotron Rd MS6R2100 Berkeley, CA 94720 C) 510.206.4418 O) 510.495.2697 Beamline) 510.495.2134 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

2 1

PDB Deposition Question
by Alexander Hirschi 25 May '10

25 May '10

Hi all - I refined a structure using Phenix 1.4-3 and now I need to submit it to the PDB. When I try to run sfcheck on the .mtz that Phenix generated, sfcheck can't read the file. I also tried to convert the .mtz to .cns using phenix.reflection_file_converter. This time, sfcheck could "not find IOBS or FOBS in your file". Any advice would be greatly appreciated. -Alex

4 4

Installation problem, please help
by Tiberman Danny 24 May '10

24 May '10

Hello, while trying to install I get the following Error message: "======================================================================== Installing from binary bundle ========================================================================== bundle file: bundle-1.6.1-357-intel-linux-2.6-x86_64.tar.gz PHENIX installation source directory set to: /root/phenix-installer-1.6.1-357 PHENIX installation target directory <PHENIX_LOC > set to: /usr/local/phenix-1.6.1-357 PHENIX installation build directory set to: <PHENIX_LOC > /build/intel-linux-2.6-x86_64 PHENIX temporary build directory set to: build-binary/intel-linux-2.6-x86_64/monst/tmp PHENIX installation log directory set to: build-binary/intel-linux-2.6-x86_64/monst/log ==== warning: "/usr/local/phenix-1.6.1-357" already exists ==== ==== warning: cannot determine installation time-stamp, proceeding ==== installing binary components removing existing files in build-binary/intel-linux-2.6-x86_64/monst/tmp/binary finding installed versions finding binary versions binary components up-to-date the binary package will not be installed setting environment configuring PHENIX components /root/phenix-installer-1.6.1-357/bin/install-from-binary: line 347: 19650 Floating point exception$PHENIX_BUILD/base/bin/python $PHENIX_LOC/cctbx_project/libtbx/configure.py --current_working_directory="$PHENIX_BUILD" $PHENIX_COMPONENTS >>$PHENIX_LOG/binary.log 2>&1 Error configuring: see /root/phenix-installer-1.6.1-357/build-binary/intel-linux-2.6-x86_64/monst/log/binary.log for details" Thanks

2 1

Re: [phenixbb] Restraining altloc related atoms during refinement
by John Pascal 20 May '10

20 May '10

>I have my doubts about supporting NCS restraints of alternative >conformations; the two seem to be mutually exclusive. But maybe >I don't understand your situtation. Here are some more details on the refinement situation: it is possible for a piece of duplex DNA to be positioned in two conformations in the crystal, and it appears that both conformations are occupied with roughly the same occupancy. The sequence is not completely palindromic, so at some positions a different base is present depending on the random orientation of the DNA. However, the phosphate and ribose atoms occupy the same positions regardless of the orientation, and I want to enforce this. >I can see that NCS restraints of backbone atoms only are useful, >even if the residue names are different. >The trouble is that atoms in two chains need to be matched >unambiguously while allowing for chain breaks. >I could try to simply remove the requirement that the residue names >have to match and go by residue sequence number + insertion code >only. Would that work in your case? Yes, that sounds like it would work. >A potential downside is that people have to be careful >excluding non-matching sidechain atoms. Perhaps it could be that there is an override option after an error message. Thanks for your help. -John -- John Pascal, PhD ph 215.503.4596 Assistant Professor fx 215.923.2117 Thomas Jefferson University Biochemistry & Molecular Biology 233 South 10th Street, BLSB 804 Philadelphia, Pennsylvania 19107 ----- Forwarded message from John.Pascal(a)mail.jci.tju.edu ----- Date: Thu, 20 May 2010 11:54:33 -0400 From: John Pascal <John.Pascal(a)mail.jci.tju.edu> Reply-To: John Pascal <John.Pascal(a)mail.jci.tju.edu> Subject: Restraining altloc related atoms during refinement To: phenixbb(a)phenix-online.org Hello All, I would like to restrain the positions of ribose/phosphate atoms of altloc 'A' to the positions of those in altloc 'B'. Problem is, altloc A and altloc B are different residue types, so phenix.refine exits stating that there are no common atoms when I enter the following NCS-type restraint: selection: chain I and altloc A and resid 1 and (phosphate or ribose) restrain: chain I and altloc B and resid 1 and (phosphate or ribose) I am using the phenix GUI. Is there a better place to enter this type of restraint (rather than NCS) that does not care if a P atom is from a Cd or from a Gd residue type? thanks. Best, John -- John Pascal, PhD ph 215.503.4596 Assistant Professor fx 215.923.2117 Thomas Jefferson University Biochemistry & Molecular Biology 233 South 10th Street, BLSB 804 Philadelphia, Pennsylvania 19107 ---------------------------------------------------------------- This message was sent using IMP, the Internet Messaging Program. ----- End forwarded message ----- ---------------------------------------------------------------- This message was sent using IMP, the Internet Messaging Program.

2 1

Restraining altloc related atoms during refinement
by John Pascal 20 May '10

20 May '10

Hello All, I would like to restrain the positions of ribose/phosphate atoms of altloc 'A' to the positions of those in altloc 'B'. Problem is, altloc A and altloc B are different residue types, so phenix.refine exits stating that there are no common atoms when I enter the following NCS-type restraint: selection: chain I and altloc A and resid 1 and (phosphate or ribose) restrain: chain I and altloc B and resid 1 and (phosphate or ribose) I am using the phenix GUI. Is there a better place to enter this type of restraint (rather than NCS) that does not care if a P atom is from a Cd or from a Gd residue type? thanks. Best, John -- John Pascal, PhD ph 215.503.4596 Assistant Professor fx 215.923.2117 Thomas Jefferson University Biochemistry & Molecular Biology 233 South 10th Street, BLSB 804 Philadelphia, Pennsylvania 19107 ---------------------------------------------------------------- This message was sent using IMP, the Internet Messaging Program.

3 2

Re: [phenixbb] all reflections are equal, but some more than others
by Frank von Delft 19 May '10

19 May '10

Are you sure? It's the sigma estimation he talked about, and that does assume it's the same foe both strong and weak, doesn't it? ----- Reply message ----- From: "Phil Evans" <pre(a)mrc-lmb.cam.ac.uk> Date: Wed, May 19, 2010 21:22 Subject: [phenixbb] all reflections are equal, but some more than others To: "PHENIX user mailing list" <phenixbb(a)phenix-online.org> If by "scaling" you mean scaling unmerged observations together, I don't think it is right that it assumes a unimodal intensity distribution. In refinement I believe the underlying distributions are assumed to be unimodal, but in data reduction scaling there is no such assumption. That scaling just assumes that the experiment (mainly the diffraction geometry) can be parameterised in such a way that reflections close together in reciprocal space and in measurement time have similar and parameterisable scales. This is independent of the intensity distribution. Phil On 19 May 2010, at 08:44, esko.oksanen(a)helsinki.fi wrote: > Frank, > > Having had a case of such hectic pseudosymmetry, it turns out that the problem originates already from the scaling that assumes a unimodal intensity distribution. The weak reflections then have their sigmas heavily overestimated and the strong ones overestimated. I had to scale the weak and strong reflections separately and ridig-body refine against the weak ones... Of course this type of problem has little to do with issue of treating negative intensities. > > Esko > > Quoting "Frank von Delft" <frank.vondelft(a)sgc.ox.ac.uk>: > >> The philosophical arguments are fine, but is Pavel not justified in asking for real cases where it matters? It's not like he hasn't spent time thinking about it, and it's not like the whole phenix community isn't sitting in his ear about their own favourite missing fix. >> >> I suppose what comes to mind is a case with hectic pseudotranslation causing half the reflections to be systematically almost but not quite zero. But then again, I understand that even if you don't toss the weak ones out, current algorithms don't deal with this well anyway, so it needs special treatment (for now: refine in smaller cell, then rigid-body refine in super-cell). >> >> phx. >> >> >> >> >> >> On 18/05/2010 18:33, Ed Pozharski wrote: >>> On Mon, 2010-05-17 at 16:20 -0700, Peter Zwart wrote: >>> >>>> You make it sound like it is a bad thing. The effect of restraint >>>> weights (ADP, geometry) has most likely a much bigger impact on the >>>> final structure then a small fraction of smallish intensities (*) >>>> >>> Peter, >>> >>> I think discarding data has to be justified. Two points: >>> >>> 1. The fraction of negative intensities is not necessarily small. It >>> depends on resolution cutoff (and you make an excellent point about >>> PTS), but looking at scalepack log-files I can tell that in my hands the >>> fraction is often 10% or more (DISCLAIMER: I belong to I/sigma=1 >>> resolution cutoff cult). >>> >>> 2. Just because these reflections are weak does not mean that they are >>> insignificant. Their contribution to the maps may be small (now I fear >>> another round of "fill-in missing Fobs with Fc for map calculation" >>> discussion), but keeping Fc close to zero for these reflections during >>> refinement seems to be just as important as to keep Fc close to whatever >>> values the strong reflections have. >>> >>> The weak reflections are not fundamentally worse than strong(er) >>> reflections in the same resolution range. They are measured with >>> roughly the same precision. Moreover, the practice of setting negative >>> intensities to zero and then ignoring them in refinement discards those >>> that are barely negative and leaves in those that were (quite randomly) >>> barely positive. >>> >>> You are absolutely right that other factors will have impact on the >>> model. But that does not mean that discarding weak data is justified. >>> Crystallographic refinement is a Rube Goldberg machine, and all the >>> components should be as good as we can make them. Perhaps there could >>> be something better than French&Wilson, but discarding negative >>> intensity reflections is hardly the solution. >>> >>> Cheers, >>> >>> Ed. >>> >>> PS. Personally I have no stake in this, since I always use truncate. >>> >>> >>> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://phenix-online.org/mailman/listinfo/phenixbb >> >> > > > _______________________________________________ > phenixbb mailing list > phenixbb(a)phenix-online.org > http://phenix-online.org/mailman/listinfo/phenixbb _______________________________________________ phenixbb mailing list phenixbb(a)phenix-online.org http://phenix-online.org/mailman/listinfo/phenixbb

1 0

View DM map from AutoSol gui?
by Francis E Reyes 19 May '10

19 May '10

I'm not a GUI kind of person but today I'm messing around with it. After AutoSol (with build=False... errr AutoBuild unchecked), I have an opportunity to 'View sites and experimental map' in coot or PyMOL. Wouldn't it be more useful if the gui loaded the DM maps as well? F --------------------------------------------- Francis Reyes M.Sc. 215 UCB University of Colorado at Boulder gpg --keyserver pgp.mit.edu --recv-keys 67BA8D5D 8AE2 F2F4 90F7 9640 28BC 686F 78FD 6669 67BA 8D5D

2 3

changing TLS groups mid refinement
by Scott Classen 19 May '10

19 May '10

Hello Phenixers, Are there any considerations to changing TLS groups mid refinement? Can I rest assured that after a few rounds of refinement with the new TLS definitions everything will sort itself out or should I do something drastic, like remove all the ANISOU records, reset B factors to some constant value, and start over? Thanks, Scott

14 34

phenix distorts ligand
by Mustafa Koksal 18 May '10

18 May '10

Hi, I have a problem in phenix refinement at 2.25 A resolution, data is 80% complete. I have protein, water, and ions in the .PDB file, and a ligand molecule that is distorted during refinement. I generated the pdb (with H atoms) with ChemDraw3D, edited the atom names etc. I used phenix.elbow to generate the .CIF file. Here is the first problem: Ligand has trans double bonds, which are written out by phenix.elbow as single rotatable bonds although I had correct connect records and H atoms in the .PDB file. I ran phenix.refine, it distorted the ligand, making double bonds not planar and nearly cis. Here is what I have tried: I edited the .CIF file to change these torsion agles to be CONSTANT and 180 or 0 deg. Here is the second problem: phenix.refine still distorted the ligand. Here is what I have tried: I used optimize_wxc/wxu in phenix.refine to optimize the weight, to see if I overweighed the data. Here is the third problem: phenix.refine still distorted the ligand. Is there a way to solve this problem? Is there a way to tighten the constraints for the ligand only? I still do not understand why it even cannot keep the double bonds planar (and trans). I do not want to move back to good old CNS... Thanks to all, Mustafa -- Mustafa Koksal Postdoctoral Fellow University of Pennsylvania Department of Chemistry Roy and Diana Vagelos Laboratories Office: 2090 IAST Office Phone: (215) 898-2227 Mobile Phone: (267) 319-6443 Mail Address: 231 S. 34th St, Philadelphia 19104, USA -- "If you grant me a pen, I shall destroy swords and spears" Mawlana Jalal-ad-Din Rumi --

2 1

geometry_restraints between sym mates?
by Scott Classen 17 May '10

17 May '10

Hi all, Is this possible? # on working model bond { atom_selection_1 = "name SM and chain M and resname SM and resseq 1" atom_selection_2 = "name OE1 and chain A and resname GLU and resseq 313" distance_ideal = 2.5 sigma = 0.2 } # on SYM mol bond { atom_selection_1 = "name SM and chain M and resname SM and resseq 1" atom_selection_2 = "name OE1 and chain A and resname GLU and resseq 285" distance_ideal = 2.5 sigma = 0.2 } I'll find out in a few minutes when I run it through phenix.refine, but I thought to ask anyways. Thanks, Scott

3 4

Looking at omit maps in coot
by xaravich ivan 17 May '10

17 May '10

Hi, Again this is pretty novice and I would apologize in advance. I ran Simulated annealing omit map in PHENIX and it gave me a bunch of files. I started with a partially built model ( I would say most of it was built) and it gave me an edited pdb that I can open in coot but I the composite omit mtz file is not opening in coot. In fact I cannot open any mtz file from the phenix output /OMIT directory and cannot see any map. I am pretty clueless. Thanks, Ivan

2 3

An easier explanation of the regions in omit map
by xaravich ivan 17 May '10

17 May '10

Hi Nat, So I have quite a few regions build by the SA_omit map of phenix. But How do I go about looking at these maps. What am I looking for here. Should I look at the regions where the occupancy was set zero in all these outputs and find if I could build some thing from one of the maps? I guess I am asking, How do I use this information to move forward in refinement? There are 24 regions map and .pdb output in the /OMIT directory. Thanks, Ivan On Fri, May 14, 2010 at 5:52 PM, Nathaniel Echols <nechols(a)lbl.gov> wrote: > On Fri, May 14, 2010 at 5:40 PM, xaravich ivan <xaravich.ivan(a)gmail.com>wrote: > >> Again this is pretty novice and I would apologize in advance. >> I ran Simulated annealing omit map in PHENIX and it gave me a bunch of >> files. I started with a partially built model ( I would say most of it was >> built) and it gave me an edited pdb that I can open in coot but I the >> composite omit mtz file is not opening in coot. In fact I cannot open any >> mtz file from the phenix output /OMIT directory and cannot see any map. >> I am pretty clueless. >> > > Some clarification: > > 1. Which version of Phenix, and which platform? > 2. GUI or command-line? > 3. If you're really making a composite omit map (as opposed to a simple > omit map, e.g. for a ligand), ignore the PDB file output by AutoBuild and > stick to the original. > > One common problem is that on Mac, if you try to open the MTZ file in Coot > (using the "Open MTZ", not auto-open), it will open a window to choose > column labels in the file, but it opens *behind* the main window, which you > need to move to the side temporarily. > > -Nat > > _______________________________________________ > phenixbb mailing list > phenixbb(a)phenix-online.org > http://phenix-online.org/mailman/listinfo/phenixbb > >

2 1

Building DNA
by Rajagopalan, Senapathy 15 May '10

15 May '10

Dear Phenix developers, I am solving a protein-dna complex and Autobuild is able to build both the protein as well as DNA chain (one at a time). However, when I try to do the refinement, I get an error saying "Phenix could not find the geometry restraints in the standard monomer library for one or more residues in the pdb file... DG: 6 copies, DA: 3 copies, DT: 43 copies". It seems like Autobuid uses the nomenclature DA, DC, DG and DT for the DNA bases while Phenix.Refine either doesn't like the nomenclature or the restraint files uses a different one. It still however lets me do the refinement...although I don't know if it just ignores the DNA part. Can someone enlighten me on this. Also, the nomenclature Coot seems to like is Ad, Cd, Gd, Td for the DNA and doesn't like the above phenix style. In order to manipulate in coot, I had to rename all the bases...is it possible to have a format that is compatible both in Phenix and Coot or am I missing something obvious here. Thanks and Have a great weekend! Sena Methodist. Leading Medicine. Recognized by U.S.News & World Report as an "Honor Roll" hospital Named to FORTUNE magazine's "100 Best Companies to Work For" list five years in a row Designated as a Magnet hospital for excellence in nursing Visit us at methodisthealth.com Follow us at twitter.com/MethodistHosp and www.facebook.com/methodisthospital ***CONFIDENTIALITY NOTICE*** This e-mail is the property of The Methodist Hospital and/or its relevant affiliates and may contain restricted and privileged material for the sole use of the intended recipient(s). Any review, use, distribution or disclosure by others is strictly prohibited. If you are not the intended recipient (or authorized to receive for the recipient), please contact the sender and delete all copies of the message. Thank you.

4 3

Improved refinement method for low resolution data
by Florian Nachon 14 May '10

14 May '10

Hi all, Axel Brunger just published a paper in Nature (http://www.nature.com/nature/journal/v464/n7292/abs/nature08892.html) about a new method to improve models despite low-resolution data. It is implemented in CNS 1.3. It might be a bit early, but is there any plan to implement that method in phenix? Florian.

10 21

where do phenix.metal_coordination ideal bond lengths come from?
by Scott Classen 14 May '10

14 May '10

Hello, I have a structure with several water/carboxylate-coordinated metal ions and I am using phenix.metal_coordination to generate the restraints parameter file for phenix.refine. Can someone please tell me where phenix.metal_coordination is getting the ideal bond length/angle values from. Thanks, Scott Classen

2 1

Gordon Conference on Diffraction Methods
by A Leslie 14 May '10

14 May '10

Just a reminder that there is still time to register for the Gordon Conference on Diffraction Methods in Structural Biology, that will be held at Bates College in Lewiston, Maine from 18-23rd July. The program (see below) covers all aspects of macromolecular crystallography, from crystallization to structure solution and refinement, and includes sessions on complementary techniques and the potential of free electron lasers. Many of the leading methods developers will be present and the format of the meeting, with free afternoons, provides ample opportunities for informal discussions. Further details, including a full program and details on the application procedure can be obtained from: http://www.grc.org/programs.aspx?year=2010&program=diffrac We hope to see you there ! Andrew Leslie (Chair) & Ana Gonzalez (Vice Chair) Outline Program: Sunday pm: Macromolecular Structures, Pushing the boundaries Discussion leader: Michael Rossmann 1. ABC Transporters: Doug Rees 2. The Spliceosome : Kyoshi Nagai Monday am: Crystallisation and the diffraction experiment Discussion leader: Zbyszek Dauter 1. Nanolitre crystallisation: Seth Harris 2. Twinning and disorder: Todd Yeates 3. Simulating the diffraction experiment: James Holton 4. Microcrystallography: Gwyndaf Evans Monday pm: Radiation damage and data collection strategies Discussion leader: Vivian Stojanoff 1. Radiation damage, experimental: Elspeth Garman 2. Microbeams: Yanhui Zou 3. Testing data collection strategies with the Pilatus: Marcus Mueller Tues am: Synchrotron developments and automated data processing Discussion leader: Thomas Schneider 1. Synchrotrons, latest developments: Sean Mc Sweeney 2. Synchrotrons, latest developments: Janet Smith 3. Data processing with Xia2/CHEF: Graeme Winter 4. Data processing with AUTOPROC: Clemens Vonrhein Tues pm: Complementary techniques Discussion leader: Ana Gonzales 1. Electron tomography: Sriram Subramaniam 2. SAXS: Hiro Tsuruta 3. Spectroscopy: Carrie Wilmot Weds am: Structure solution and refinement Discussion leader: Paul Adams 1. SAD Phasing: Randy Read 2. Molecular replacement with BALBES: Garib Murshudov 3. PHENIX: Tom Terwilliger 4. TLS: Ethan Merritt Weds pm: Challenging problems/Membrane Proteins Discussion leader: TBC 1. Using cubic lipidic phase: Martin Caffrey 2. Membrane Proteins: Stephen White 3. The Bacterial Ribosome: Venki Ramakrishnan Thurs am: Selected Posters and Map improvement and Model building Discussion leader: Tassos Perrakis 1. COOT: Paul Emsley 2. Modelling nucleic acids: Johan Hattne Thursday pm: Future Methods, FELs, Diffraction imaging Discussion leader: Ed Lattman 1. Anton Barty 2. Serial Crystallography: John Spence and Petra Fromme

1 0

question about free energy calculation in protein complex
by crystallogrphy 14 May '10

14 May '10

Hi, I am trying to calculate the free energy of association for structure of a protein complexes based on the salt bridges, H-bonds and van der waals interaction? Does any one know how to do it? Thanks!

3 2

phenix.model_vs_data
by Young-Jin Cho 13 May '10

13 May '10

Dear Phenix users, This might be trivial questions though... 1. When you did phenix.model_vs_data, I see the discrepancy between R/Rfree. 'After applying resolution and sigma cutoffs' seems more accurate one though, could someone explain why it has such-even small though- discrepancy? 2. Also, when you do phenix.polygon, I used to seeing many branches from old version. With my current version(phenix-1.6.1-343), I normally see only 5 branches. Is there any reasons for that? Thanks in advance as always. Young-Jin

4 6

Refinement with partial model - question
by Peter Grey 13 May '10

13 May '10

Dear Phenix Users, I have a partial model at low resolution. Is there a way let Phenix.refine know the model is not complete and that when a mask is calculated for bulk solvent correction it should take this fact into consideration ? Alternatively, can I use a solvent mask from density modification program in Phenix.refine somehow ? Alternatively, is there a program that will build dummy atoms (zero occupancy) automatically in a map in all peaks above a certain sigma (i.e. atoms just for mask calculation). Thank you in advance, Peter.

14 32

Can simulated annealing remove model bias after MR?
by Francis E Reyes 13 May '10

13 May '10

Just curious.... Thanks F --------------------------------------------- Francis Reyes M.Sc. 215 UCB University of Colorado at Boulder gpg --keyserver pgp.mit.edu --recv-keys 67BA8D5D 8AE2 F2F4 90F7 9640 28BC 686F 78FD 6669 67BA 8D5D

2 1

B factors of zero after refinement
by Kyle Dolan 13 May '10

13 May '10

Hello all, I have been performing standard refinement of atomic coordinates and B factors (ADP) in phenix.refine using NCS and bulk solvent correction without anisotropic scaling. The program has been returning a handful of atoms with B factors of zero after refinement, and I'm not sure why. I tried using adp.set_b_iso to reset all B factors before refinement but some still come back down to zero. Has anyone else noticed something similar, or has a suggestion on what to do? Thanks, Kyle Kyle T. Dolan Department of Biochemistry and Molecular Biology The University of Chicago ktd(a)uchicago.edu

2 1

Phenix Ligand Identification error
by Rajan 13 May '10

13 May '10

Dear All Greetings I was searching for a ligand present in the active site of my protein using Ligand Search in "PHENIX" but i stuck in between with an error : A Python error was detected. This may be a problem with the program settings, an error in your data, or a bug; click "OK" to send a bug report to the PHENIX developers. AttributeError : ligand_identification instance has no attribute 'connections' Traceback: File "/home/programs/linux/phenix- 1.6.1-357/cctbx_project/libtbx/thread_utils.py", line 137, in run return_value = self._target(self._args, self._kwargs, self._c) File "/home/programs/linux/phenix-1.6.1-357/cctbx_project/libtbx/runtime_utils.py", line 55, in __call__ result = self.target() File "/home/programs/linux/phenix-1.6.1-357/phenix/phenix/command_line/ligand_identification.py", line 1769, in __call__ identify_ligands = ligand_identification(args=list(self.args), quiet=False) File "/home/programs/linux/phenix-1.6.1-357/phenix/phenix/command_line/ligand_identification.py", line 308, in __init__ ligand_identification.Run(self, params) File "/home/programs/linux/phenix-1.6.1-357/phenix/phenix/command_line/ligand_identification.py", line 1454, in Run return self.connections['STOP '] I am not able to identify the problem? Please let me know if some one can identify the problem, i don't know whether it is a bug or problem in my files. Thanks for your suggestions in advance. Rajan -- ************************************************ Current Address: Rajan Vyas Research Scholar Deptt. of Biotechnology Panjab University Chandigarh, India 160014 Mob. +919417374197 Fax: +91-172-2625254 ************************************************

5 4