April 2010 - phenixbb - phenix-online.org

individual_sites_real_space?
by Francis E Reyes 07 Apr '10

07 Apr '10

Phenixer's I just upgraded from 1.4-87 and noticed that there's individual_sites_real_space now. The documentation (http://www.phenix-online.org/documentation/refinement.htm ) just shows that it's a parameter phenix.refine. What is this and when is it helpful (bad resolution/ bad model? or high resolution/ good model? or high resolution/bad model?) Thx 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

LSTc refinement in phenix.refine
by JXQI 07 Apr '10

07 Apr '10

Dear all, I am refining a structure complexed with LSTc (alpha 2,6-linked galactose,i.e. SIA-GAL-NAG-GAL).I know that the BETA1-4 or BETA1-3 in the mon-lib could restrain the bond distance,but I could not find the links between the SIA and GAL......any suggestions are appreciated. Best regards Janxon

2 1

polygon and membrane proteins
by Pascal Egea 06 Apr '10

06 Apr '10

Dear Developers, I wonder if there would be a simple way to get the Polygon program implemented in Phenix to use only membrane protein deposited structures as the reference data set. I understand that the reference set is much smaller than the whole menagerie of available structures but it represents an homogenous class of protein (crystallized in detergent). Do you think that this would be meaningful and useful? Thanks in advance, -- Pascal F. Egea, PhD Assistant Professor UCLA, David Geffen School of Medicine Department of Biological Chemistry 314 Biomedical Sciences Research Building office (310)-983-3515 lab (310)-983-3516 email pegea(a)mednet.ucla.edu

4 4

Re Double stranded DNA restraints?
by Thomas C. Terwilliger 06 Apr '10

06 Apr '10

Hi Sena, I would add this command helices_strands_only=True to force autobuild to just try and build helical segments. This probably won't fix the lack of connectivity however. All the best, Tom T Hi Everyone, Is there a way to force phenix (autobuild) into building a contiguous double stranded DNA. Currently, when I am solving a protein-DNA complex using SAD data: 1) its not recognizing that DNA needs to be double stranded, so it is fitting one of the strands correctly (at least part of it), but the other strand is in the protein density region. 2) It breaks down each of the strand. So it is not fitting as a contiguous 30 mer but instead a 10 mer, two 6 mers and a 8 mer. Thanks Sena

1 0

PPC binaries?
by Marcelo Carlos Sousa 05 Apr '10

05 Apr '10

Is Phenix no longer supported on Mac PPC? ____________________________________________ Marcelo C. Sousa, Ph.D. Department of Chemistry and Biochemistry Cristol Chemistry Building. Room 100 215 UCB University of Colorado at Boulder Boulder, CO 80309 Phone: (303) 735 4341 Fax: (303) 492 5894

2 1

(no subject)
by Prof Jennifer Martin (IMB) 04 Apr '10

04 Apr '10

The protein crystal structure I am refining has a redox modification: cysteine sulfenic acid (S-hydroxy cysteine). How do I go about including parameters for this modification in phenix refine? I have found parameters on HIC-UP (under CSO) but these are for XPLOR/ CNS. Many thanks in advance for suggestions, Jenny Martin University of Queensland

3 2

Torsion angels refinement question
by Peter Grey 03 Apr '10

03 Apr '10

Dear Phenix users, I try to refine a big complex at low resolution. I am a little confused about the syntax of using torsion angels gradient-driven refinement. Is it now the default for individual sites refinement instead of cartesian ? If not, should I add : sites.Torsion_angles=all for example ? Could you shed some light on this, Thank you, -- Peter

2 1

Optimizing the Geometry Weight for a 2.76 A structure
by Sam Stampfer 02 Apr '10

02 Apr '10

Hi, I'm refining a 2.76 A structure and need help optimizing the refinement weight on geometry. When I refine at the default refinement weight (wxc = 0.5), phenix seems to generate a much poorer structure (both geometrically AND by Rfree value) relative to lower wxc values which put a stronger weight on geometry. I am trying to pick the best structure to use for my next round of rebuilding. The R-stats and condensed molprobity stats are listed below (at the end of this email). I have a few questions regarding this structure and those stats: 1. What is the relevance of having a LOW difference between Rfree and Rwork? My structure seems to have a higher than expected Rwork value relative to Rfree (aren't they supposed to differ by 0.05 or so, instead of the 0.03 observed in my structure?) Additionally, as I refine at lower wxc values (more weight on geometry) I get an increased Rwork value with an equal or decreased Rfree (relative to the refinement at a higher wxc value). Does that mean my structure is better or worse? 2. How can I judge the output from my refinement? I have looked at Rwork, Rfree, and the molprobity clashscore and overall score values. I included them below, at the end of this email. How do I tell which the best refinement is? Which one would you suggest? I thought the best was wxc = 0.1 since the R-work and Rfree aren't changed much from the start values but the geometry is far better. 3. How do I optimize the weight on geometry in refinement? My PI said that phenix can do it but not for twinned structures (my structure is twinned). Can the new version of phenix do weight optimization for twinned structures? How? If not, what's the best way for me to do it manually? As always, thanks immensely for all your help! -Sam All have the same values for: Start R-work = 0.2002, R-free = 0.2376 *wxc = 0.01 *Final R-work = 0.2104, R-free = 0.2378 Molprobity Clashscore: 19.46 Molprobity Score: 2.29 *wxc = 0.05 *Final R-work = 0.2049, R-free = 0.2377 Molprobity Clashscore: 20.96 Molprobity Score: 2.38 *wxc = 0.1 *Final R-work = 0.2010, R-free = 0.2388 Molprobity Clashscore: 23.50 Molprobity Score: 2.54 *wxc = 0.2 *Final R-work = 0.1969, R-free = 0.2399 Molprobity Clashscore: 27.15 Molprobity Score: 2.79 *wxc = 0.3 *Final R-work = 0.1945, R-free = 0.2411 Molprobity Clashscore: 30.76 Molprobity Score: 2.97 *wxc = 0.4 *Final R-work = 0.1930, R-free = 0.2420 Molprobity Clashscore: 34.54 Molprobity Score: 3.11 *wxc = default (0.5) *Final R-work = 0.1918, R-free = 0.2429 Molprobity Clashscore: 37.73 Molprobity Score: 3.22

6 7

Re: [phenixbb] Optimizing the Geometry Weight for a 2.76 A structure
by det102＠uoxray.uoregon.edu 02 Apr '10

02 Apr '10

Dear Phil, I don't have access to vol 277 of Met Enz here, which is sad since I wrote a chapter in the thing, so I can't comment on Axel's exact words. I agree that the uncertainty in the ability of the free R estimate based on a subset of reflections to predict the mean of the free R calculated from all possible test sets of that size is proportional to 1/N. My point is that the pdf of the difference of two random variables is not the same as the pdf of either of them. You have to know the coupling between them to decide if the difference between them has greater or less precision than the individual variable. For example, let's assume two measurements with a precision of 1% but the difference in the two measurements is 0.5%. This small difference can be significant if the causes of the uncertainties are systematic. The "errors" would cancel out. I believe the fluctuations in free R estimate due to the small size of the test set are due to the particular indices in the test set and therefore mainly systematic. The intensities of some reflections may have been measured better than others and your test set may happen to be more enriched with those, resulting is a slightly lower free R estimate. Since the test set is always the same during refinement the benefit of this enrichment (or detriment if you were unlucky) will be subtracted out when you compare the two free R estimates during your optimization. I expect the difference between two free R estimates (using the same test set) will be a more precise indicator of the true change in free R than either estimate is in predicting the free R itself. Dale Phil Jeffrey wrote: > Although I am totally unclear why the form of the equation, which is > just counting/Poisson statistics, doesn't also apply to the uncertainty > of calculating the average percentage deviation over a test set of size > N. For example a test set of N=1000 could be regarded as 100 different > test sets of size N=10, and I think it's likely that the distribution of > R-free for these 100 mini test sets would be Poisson in form centered > around the R-free for the superset N=1000. > > So, since the change in |Fo-Fc| for each reflection doesn't simply scale > with wxc in structure refinement nor does it inevitably decrease for > every reflection on structure improvement, the change in R-free for any > given change in structure should be related to N in the form of the > equation given, so does reflect the s.d. of R-free as an estimate for > structure improvement. > > No ? > > If not, then I misinterpreted what Brunger mean on p.394 of his 1997 Met > Enz paper because that's certainly what it read like. > > Phil > > det102(a)uoxray.uoregon.edu wrote: >> >> That uncertainty is not quite the same thing. What you describe >> is the uncertainty in the free R due to the small sampling size of >> the test set. It is the spread of free R's obtained when different >> test sets are chosen for the same project and circumstances, and is >> useful when comparing free R's calculated using different test sets, >> or, heaven forbid, trying to compare free R's from different crystals.

1 0

map grid size
by Sean Gay 02 Apr '10

02 Apr '10

Does anyone know if there is a way to increase the map sampling rate for Phenix fofc and 2fofc maps? I want to make a figure for publication in PyMOL and the map_double command helps a bit, but i'd still like more points on the map grids. Sean Sean C. Gay, PhD Postdoctoral Scholar Skaggs School of Pharmacy and Pharmaceutical Sciences University of California, San Diego 9500 Gilman Dr. #0703 La Jolla, CA 92093 Tel. (858)822-7804 Fax. (858)246-0089 Email. scgay(a)ucsd.edu

2 1

Double stranded DNA restraints?
by Rajagopalan, Senapathy 02 Apr '10

02 Apr '10

Hi Everyone, Is there a way to force phenix (autobuild) into building a contiguous double stranded DNA. Currently, when I am solving a protein-DNA complex using SAD data: 1) its not recognizing that DNA needs to be double stranded, so it is fitting one of the strands correctly (at least part of it), but the other strand is in the protein density region. 2) It breaks down each of the strand. So it is not fitting as a contiguous 30 mer but instead a 10 mer, two 6 mers and a 8 mer. Thanks 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.

1 0

Secondary structure restraints for nucleic acids
by Peter Grey 02 Apr '10

02 Apr '10

Dear developers, Does the new mechanism of secondary structure restraints apply to nucleic acids as well ? Are the average bond distances similar and does DSSP works well with RNA/DNA ? Thank you in advance for clarifying this issue, Peter.

2 2

Installation problem
by Tatsuya Kaminishi 02 Apr '10

02 Apr '10

Dear Phenix BB When I tried to install the latest binary version (1.6.1), the "install" script always stopped and dumped an error message as follows: ----- /software/PHENIX/phenix-installer-1.6.1-357/bin/install-from-binary: line 347: 10964 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 /software/PHENIX/phenix-installer-1.6.1-357/ build-binary/intel-linux-2.6-x86_64/lupus/log/binary.log for details ----- And this "binary.log" file looks like this; ========== . . After many lines... . ./reduce/comment_out_carboxylate_H_v2.3.pl ./reduce/Makefile ./reduce/libtbx_refresh.py ========== This never happened when I installed version 1.6-289. Could anybody out there tell me how I can get around this? Or anything important is missing? The system I'm using is: Linux crux 2.6.16.60-0.27-smp #1 SMP Mon Jul 28 12:55:32 UTC 2008 x86_64 x86_64 x86_64 GNU/Linux SUSE Linux Enterprise Desktop 10 (x86_64) VERSION = 10 PATCHLEVEL = 2 Thank you very much for your help in advance! Best regards, kami Tatsuya KAMINISHI, Ph.D. Fucini Laboratory Institute for Organic Chemistry and Chemical Biology Johann Wolfgang Goethe University Max-von-Laue-Strasse 7, N160/B413 D-60438 Frankfurt am Main, Germany Tel: +49 69 798-29165 Fax: +49 69 798-29268 Web: http://user.uni-frankfurt.de/~joharms/akfucini/

2 1

Restraints on dihedrals in phenix.refine
by Gregers Rom Andersen 02 Apr '10

02 Apr '10

hi phenix-bb in the standard phenix.refine scripts there are options for adding user defined restraints on bonds and angles, but not on dihedrals. Is this currently possible in phenix.refine? One application would be to restrain the phi-psi angles of a low resolution structure to phi-psi angles from a high-res structure of the same protein, an option previously made possible for xplor in the USF program XPLO2D gregers Gregers Rom Andersen Associate professor Department of Molecular Biology University of Aarhus Gustav Wiedsvej 10C DK8000 Aarhus C phone +45 89425024 fax +45 86123178 www.bioxray.au.dk/gra

3 2