November 2012 - phenixbb - phenix-online.org

Simulated anneal
by Qiang Chen 14 Nov '12

14 Nov '12

Hi all, I have a 3.05Å structure. I refine it to Rwork=25.2 Rfree=28.81. I use Molprobity to check the structure. The clash score is 100th and Molprobity score is 99th. I was suggested to run simulated anneal to help reduce Rwork and Rfree. I tried Cartesian and Torsion and different seed. The best result I got is: torsion anneal, start_temp=2500,final_temp=300, random_seed=4250146. The Rwork:25.8, Rfree=31.2. I checked the new model and found that some places in a domain were not correctly traced. The density for that domain is not good and have several flexible loops. My question is: does simulated anneal always work and reduce Rwork/Rfree? At the final stage of refinement, does it necessary to use simulated anneal again? Thanks! yamei The information in this e-mail is intended only for the person to whom it is addressed. If you believe this e-mail was sent to you in error and the e-mail contains patient information, please contact the Partners Compliance HelpLine at http://www.partners.org/complianceline . If the e-mail was sent to you in error but does not contain patient information, please contact the sender and properly dispose of the e-mail.

3 2

phenix.pdbtools - renumber specific chains
by Antony Oliver 13 Nov '12

13 Nov '12

Apologies if this has been covered before - but can someone point me to a way of renumbering only specific chains (i.e. F, G, H etc). I have tried using phenix.pdbtools with either selection = "chain F or chain G or chain H" or keep "chain A or chain B or chain C or chain D" - but in both cases phenix renumbers every single residue! With thanks, Antony. --- 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 tel (office): +44 (0)1273 678349 tel (lab): +44 (0)1273 677512

3 3

Setting up MR_Rosetta
by Julio Lenin 08 Nov '12

08 Nov '12

Hello guys, I am currently trying to setup MR_Rosetta but I am running into some trouble. When I run the tutorial for Defensin HNP-3 from the GUI I always get an error during the ... process (see image at this dropbox link). I know that my path for phenix and rosetta are correct since I can run the phenix_regression.wizards.test_command_line_rosetta_quick and phenix_regression.wizards.test_command_line_rosetta_quick test_rosetta_rebuild. If I try to trace down the error, I can only get the same info that I got from the GUI error message (vide supra). As additional info, here are somme details of my setup: MacBookPro 4,1 running OS 10.6.8 Core2Duo 2.5 Ghz 6 Gb RAM Rosetta 3.2 compiled with flags bin mode=release Phenix 1.8.1-1168 And my input settings for the GUI run (dropbox) Now, this might just be confusing but there was one command line test that my setup failed: phenix_regression.wizards.test_command_line_rosetta test_mr_rosetta The error is kinda cryptic (see dropbox copy of the error) but it seem to be a discrepancy in the scoring results not a missing program or any other kind of misconfiguration. Of course, I might be wrong. Any advice? Best regards. -- Lenin Domínguez-Ramírez "And for all the current and rather silly emphasis on structural biology,understanding enzymes means understanding catalysis and catalysis is concerned with kinetics, not structure: as Jeremy Knowles aptly remarked, studying the photograph of a racehorse cannot tell you how fast it can run." Cornish-Bowden, 1995 Lenin Domínguez-Ramírez, PhD "What are we to think of an organism that loses practically all of its cellular water at ordinary temperatures and pressures and ceases to metabolize, but upon being rehydrated resumes all of the characteristics it previously exhibited?...we should be reluctant to call such an organism dead (unless we are willing to entertain the possibility of its resurrection) because it resumes an active life upon simple restoration of water."

3 2

Sugar refine in phenix
by Jan van Agthoven 05 Nov '12

05 Nov '12

Hi everyone, This is a basic question. I would like to refine the sugar ends of my protein. The sugar polymers are made of NAG and BMA and whenever I use phenix.ready_set and subsequently phenix.refine it results in a structure of disconnnected sugar monomers. Indeed the cif file generated by phenix.ready_set does'nt give any geometry constraints for connecting the monomers. How do create a cif file that tells phenix.refine how to refine the sugars? Thanks, Jan

2 4

ligand fitting, density and refinement
by Terwilliger, Thomas C 02 Nov '12

02 Nov '12

Hi Huaidong, I'm not quite sure exactly what you did...could you possibly elaborate a little? My confusion comes from two things. First, this phrase: "did the refinement with the output files from ligandfit", is confusing to me because I would not think that you would want to use any output files from ligandfit except for the ligand PDB file. Second, I agree, I can't see why ligandfit would find good density for your ligand if it is not present in a difference map calculated from the original data. All the best, Tom T ________________________________________ From: phenixbb-bounces(a)phenix-online.org [phenixbb-bounces(a)phenix-online.org] on behalf of 张怀东 [zhanghuaidong(a)hust.edu.cn] Sent: Thursday, November 01, 2012 8:11 PM To: PHENIX user mailing list Subject: Re: [phenixbb] Ligandfit_SigmaF_obs Hi,Tom T We followed the instruction after the ligandfit, and did the refinement with the output files from ligandfit. The electron density for the ligand is clear, but in the refinement with the original mtz file, the density is not there. So how to explain the difference? Is this a bug? Thanks! Huaidong > -----原始邮件----- > 发件人: "Terwilliger, Thomas C" <terwilliger(a)lanl.gov> > 发送时间: 2012年11月1日星期四 > 收件人: "PHENIX user mailing list" <phenixbb(a)phenix-online.org> > 抄送: > 主题: Re: [phenixbb] Ligandfit_SigmaF_obs > > Hi Ed, > > Yes, good point. I will be changing the way the wizards read in and store reflection data (sometime soon I hope!) and at that point this should be moot because the wizards won't write out dummy intermediate files. > > All the best, > Tom T > > ________________________________________ > From: phenixbb-bounces(a)phenix-online.org [phenixbb-bounces(a)phenix-online.org] on behalf of Ed Pozharski [epozh001(a)umaryland.edu] > Sent: Thursday, November 01, 2012 9:44 AM > To: PHENIX user mailing list > Subject: Re: [phenixbb] Ligandfit_SigmaF_obs > > On Thu, 2012-11-01 at 14:23 +0000, Terwilliger, Thomas C wrote: > > The sigma values are not useful, and are just set to 1.0 in that file. > > Tom, > > I've seen this issue before, and I am just wondering what could be the > rationale to reset sigmas to 1.0? Seems like one would have to code for > this to happen, and even if the sigma values are not used for map > generation, why isn't the column simply omitted but instead replaced > with obviously incorrect values? > > Ed. > > -- > Oh, suddenly throwing a giraffe into a volcano to make water is crazy? > Julian, King of Lemurs > > _______________________________________________ > 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

Ligandfit_SigmaF_obs
by 张怀东 02 Nov '12

02 Nov '12

Hello, I am a newcomer. I need help. I have done phenix.ligandfit, and the result contained a mtz file named resolve.mtz. While sigma_F_obs in the file are all 1.0. Is there something wrong with the file? Is this a bug? Could I do phenix.refine with the resolve.mtz added R_flag, next? sincerely, Huaidong Zhang

3 4

refining with f' and f" values
by Subhani Bandara 01 Nov '12

01 Nov '12

Hi Everyone, I have a structure with heavy metals which I would like to refine anomalous scattering factors. I could get the f' and f" values from Autosol, when I gave the atom name and wavelength. But when I ran refinement with anomalous scattering giving .mtz file and pdb file from pHENIX refinement I got the error "Invalid combination-F-obs not anomalous, provide anomalous reflections or fix all f" values. Then I used reflection file editor of GUI to create a mtz file that has anomalous data using the original reflection file and the mtz file of refinement. But still it gives an error. I must be doing something wrong, but cannot figure out what that is. Thanks in advance for any help. Subhani

2 4

Rwork Rfree Bonds Angles
by Hexi Zhang 01 Nov '12

01 Nov '12

Dear All, During Phenix refinement, it request the Rwork and Rfree graduately decrease (at least no increase). But there is no request that the Bonds (should be RMS?) and the Angles (should also be RMS?) should graduately decrease. Am I right? I am looking forward to getting your reply. Hexi

2 1

"H" atom refining occupancy with previous altconf residue in sequence
by Alexander Scouras 01 Nov '12

01 Nov '12

After a series of residues (1+) with alternate conformations, the H atom of the next residue seems to lock its occupancy to that of the final conformation of the previous residue. Hence I'm getting several H's with 0.1-0.5 occupancy on single conformation residues, where the entire rest of the residue has occupancy 1. This happens after every altconf series, and showed up in two separate structures I've worked on this week. Within a series of 2+ residues with altconfs, the H's refine correctly with their altconf groups. So it isn't just that all subsequent residue's H's are getting overridden. This behavior seems to have started with 1.8.1-1168, or at least previous refinements of the same structure didn't do this. 64 bit build on Retina MBP running OS X 10.8.2. I couldn't google any reference to this behavior, so assuming it's a new bug I can send along source files. Below I've pasted a couple residues from one of my PDBs. I added a ***** to mark the offending H with the wrong occupancy. -Alex Scouras ATOM 668 N CMET A 29 -17.970 6.464 59.720 0.21 10.72 N ANISOU 668 N CMET A 29 1641 1300 1133 -328 -379 -50 N ATOM 669 CA CMET A 29 -17.578 5.189 60.316 0.21 9.05 C ANISOU 669 CA CMET A 29 1545 882 1012 -248 -350 -140 C ATOM 670 C CMET A 29 -17.248 5.317 61.797 0.21 7.83 C ANISOU 670 C CMET A 29 1363 666 947 173 -134 166 C ATOM 671 O CMET A 29 -17.496 4.395 62.578 0.21 7.92 O ANISOU 671 O CMET A 29 1187 751 1071 371 -145 338 O ATOM 672 CB CMET A 29 -16.376 4.609 59.573 0.21 9.09 C ANISOU 672 CB CMET A 29 1561 834 1059 -243 -283 -203 C ATOM 673 CG CMET A 29 -16.681 4.210 58.145 0.21 9.05 C ANISOU 673 CG CMET A 29 1529 756 1154 -28 -331 -235 C ATOM 674 SD CMET A 29 -18.049 3.042 58.009 0.21 13.05 S ANISOU 674 SD CMET A 29 2356 1325 1279 -218 -284 -50 S ATOM 675 CE CMET A 29 -17.505 1.757 59.124 0.21 12.06 C ANISOU 675 CE CMET A 29 1866 1485 1232 -630 -667 420 C ATOM 676 H CMET A 29 -17.333 6.878 59.318 0.21 12.87 H ATOM 677 HA CMET A 29 -18.320 4.557 60.216 0.21 10.86 H ATOM 678 HB2CMET A 29 -15.671 5.275 59.551 0.21 10.91 H ATOM 679 HB3CMET A 29 -16.069 3.819 60.044 0.21 10.91 H ATOM 680 HG2CMET A 29 -16.916 5.004 57.640 0.21 10.86 H ATOM 681 HG3CMET A 29 -15.894 3.794 57.761 0.21 10.86 H ATOM 682 HE1CMET A 29 -18.164 1.059 59.140 0.21 14.48 H ATOM 683 HE2CMET A 29 -16.666 1.409 58.813 0.21 14.48 H ATOM 684 HE3CMET A 29 -17.401 2.130 60.002 0.21 14.48 H ATOM 685 N MET A 30 -16.667 6.445 62.187 1.00 6.77 N ANISOU 685 N MET A 30 1361 564 647 -30 -56 18 N ATOM 686 CA MET A 30 -16.416 6.699 63.606 1.00 7.66 C ANISOU 686 CA MET A 30 1362 946 601 -42 -173 2 C ATOM 687 C MET A 30 -17.759 6.724 64.367 1.00 7.62 C ANISOU 687 C MET A 30 1310 944 641 -9 -24 -22 C ATOM 688 O MET A 30 -17.852 6.212 65.486 1.00 9.05 O ANISOU 688 O MET A 30 1571 1151 718 84 98 231 O ATOM 689 CB MET A 30 -15.620 7.987 63.790 1.00 7.65 C ANISOU 689 CB MET A 30 1299 894 714 26 -120 -166 C ATOM 690 CG MET A 30 -14.229 7.965 63.231 1.00 8.62 C ANISOU 690 CG MET A 30 1464 1072 741 31 -163 -10 C ATOM 691 SD MET A 30 -13.130 6.745 64.006 1.00 9.04 S ANISOU 691 SD MET A 30 1555 961 919 105 -180 -94 S ATOM 692 CE MET A 30 -12.812 7.527 65.563 1.00 9.77 C ANISOU 692 CE MET A 30 1540 1284 889 150 -190 -62 C ATOM 693 H MET A 30 -16.395 7.062 61.653 0.21 8.12 H ***** ATOM 694 HA MET A 30 -15.872 5.965 63.962 1.00 9.19 H ATOM 695 HB2 MET A 30 -16.099 8.709 63.352 1.00 9.18 H ATOM 696 HB3 MET A 30 -15.550 8.173 64.739 1.00 9.18 H ATOM 697 HG2 MET A 30 -14.278 7.758 62.284 1.00 10.35 H ATOM 698 HG3 MET A 30 -13.831 8.840 63.354 1.00 10.35 H ATOM 699 HE1 MET A 30 -12.228 6.968 66.079 1.00 11.73 H ATOM 700 HE2 MET A 30 -12.397 8.378 65.406 1.00 11.73 H ATOM 701 HE3 MET A 30 -13.644 7.648 66.026 1.00 11.73 H

4 10