May 2010 - phenixbb - phenix-online.org

Weight optimization
by Jianghai Zhu 29 May '10

29 May '10

Hi, I am using the optimize_wxc=true and optimize_wxu=true in my last refinement. The R and Rfree went down but the RMSDs were way too large. I wonder if phenix.refine still targets the geometry when it optimizes the R and Rfree. Here are the results from almost two identical runs except the weight optimization. 1) optimize_wxc = False optimize_wxu = False wxc_scale=0.08 wxu_scale=1.66 results: R=0.2294 R-free=0.2709 angl=1.087 bond=0.006 2) optimize_wxc = True optimize_wxu = True wxc_scale=0.08 wxu_scale=1.66 results: R=0.2141 R-free=0.2628 angl=2.894 bond=0.033 I am using the dev-412 version. Thanks. -- Jianghai

4 8

offtopic: Problem with Pymol showing maps
by Tatyana Sysoeva 29 May '10

29 May '10

Hi all, if somebody is still not on the memorial day weekend, I have a short question off-topic for Pymol-users I have several maps from Phenix and want to draw some pictures in Pymol. I am opening the files, renamed to ccp4 or map, or not, and all I am getting is just a box, which is not even a unit cell boundary. If you know some nice simple program to do nice pictures of the electron densities- share your knowledge please! Sincerely Tanya

3 3

Super: Rotation and translation
by J. Fleming 29 May '10

29 May '10

Hello, I've been using 'super' to align the c-alphas within one domain of two different models of the same protein while leaving the second domain to move freely. I can do this using: super (A and resid 1-100 and name ca), (B and resid 1-100 and name ca) What I would like to do is determine the rotation and translation of model B's second domain as it moves away from A's second domain. I can output a matrix after running super using: print cmd.get_object_matrix("A") How do I go from the matrix to an angle and distance? Thanks for any thoughts, Jon

2 1

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