- phenixbb - phenix-online.org

Autobuild error
by Nora Cronin 21 Oct '13

21 Oct '13

Whilst running Autobuild SA omit map using 40 segments one OMIT region fail and the program quits. Error as follow: Can you help? Summary of output files for Solution 1 from rebuild cycle 2 --- Model (PDB file) --- pdb_file: /local/ncronin/phenix/APC4/Autobuild_20131018/AutoBuild_run_2_/TEMP0/AutoBuild_run_5_/overall_best.pdb --- Data for refinement FP SIGFP PHIM FOMM HLAM HLBM HLCM HLDM FreeR_flag --- refine_data: /local/ncronin/phenix/APC4/Autobuild_20131018/AutoBuild_run_2_/TEMP0/AutoBuild_run_5_/overall_best_refine_data.mtz AutoBuild_build_cycle AutoBuild Run 5 Fri Oct 18 13:30:32 2013 Ending these rebuild cycles as the limit is reached All done with build cycles in this region AutoBuild_set_up_build AutoBuild Run 5 Fri Oct 18 13:30:32 2013 Traceback (most recent call last): File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/autosol/AutoBaseExtend.py", line 989, in DoNextMethod self.CarryOutBest() # to be obtained after it is finished File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/autosol/AutoBaseExtend.py", line 2146, in CarryOutBest getattr(self,str(self.application_method))() # call this fn File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/wizards/AutoBuild.py", line 1832, in AutoBuild_set_up_build self.AutoBuild_finish_current_region() File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/wizards/AutoBuild.py", line 1909, in AutoBuild_finish_current_region self.AutoBuild_update_ncs() File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/wizards/AutoBuild.py", line 5567, in AutoBuild_update_ncs ],quiet=True,out=f) File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/command_line/find_ncs.py", line 324, in __init__ self.read_ncs_as_file(params,out=out) File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/command_line/find_ncs.py", line 511, in read_ncs_as_file quiet=self.quiet,log=null_out()) File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/command_line/simple_ncs_from_pdb.py", line 546, in __init__ starting_residue_numbers) File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/command_line/simple_ncs_from_pdb.py", line 1326, in find_invariant_groups minimum_size=min_size) File "/common/app/phenix/phenix-1.8.4-1496/cctbx_project/mmtbx/invariant_domain.py", line 35, in __init__ minimum_size=minimum_size) File "/common/app/phenix/phenix-1.8.4-1496/cctbx_project/mmtbx/invariant_domain.py", line 70, in process size, index = self.find_largest( matches, used ) File "/common/app/phenix/phenix-1.8.4-1496/cctbx_project/mmtbx/invariant_domain.py", line 168, in find_largest max_size = flex.max( sizes ) RuntimeError: unidentifiable C++ exception Traceback (most recent call last): File "/common/app/phenix/phenix-1.8.4-1496/build/mac-intel-osx-x86_64/../../phenix/phenix/command_line/autobuild.py", line 230, in <module> run_autobuild=run_autobuild(args=sys.argv[1:]) File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/autosol/wizard_command_line.py", line 564, in __init__ self.RunWizard() File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/autosol/wizard_command_line.py", line 1000, in RunWizard self.runWizard.run() File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/autosol/wizard_command_line.py", line 137, in run self.Wizard.Run(info='script') File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/autosol/AutoBaseExtend.py", line 838, in Run self.DoNextMethod() File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/autosol/AutoBaseExtend.py", line 1003, in DoNextMethod raise AssertionError, line AssertionError: Traceback (most recent call last): File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/autosol/AutoBaseExtend.py", line 989, in DoNextMethod self.CarryOutBest() # to be obtained after it is finished File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/autosol/AutoBaseExtend.py", line 2146, in CarryOutBest getattr(self,str(self.application_method))() # call this fn File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/wizards/AutoBuild.py", line 1832, in AutoBuild_set_up_build self.AutoBuild_finish_current_region() File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/wizards/AutoBuild.py", line 1909, in AutoBuild_finish_current_region self.AutoBuild_update_ncs() File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/wizards/AutoBuild.py", line 5567, in AutoBuild_update_ncs ],quiet=True,out=f) File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/command_line/find_ncs.py", line 324, in __init__ self.read_ncs_as_file(params,out=out) File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/command_line/find_ncs.py", line 511, in read_ncs_as_file quiet=self.quiet,log=null_out()) File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/command_line/simple_ncs_from_pdb.py", line 546, in __init__ starting_residue_numbers) File "/common/app/phenix/phenix-1.8.4-1496/phenix/phenix/command_line/simple_ncs_from_pdb.py", line 1326, in find_invariant_groups minimum_size=min_size) File "/common/app/phenix/phenix-1.8.4-1496/cctbx_project/mmtbx/invariant_domain.py", line 35, in __init__ minimum_size=minimum_size) File "/common/app/phenix/phenix-1.8.4-1496/cctbx_project/mmtbx/invariant_domain.py", line 70, in process size, index = self.find_largest( matches, used ) File "/common/app/phenix/phenix-1.8.4-1496/cctbx_project/mmtbx/invariant_domain.py", line 168, in find_largest max_size = flex.max( sizes ) RuntimeError: unidentifiable C++ exception Best regards Nora Dr N B Cronin X-Ray Laboratory Manager Tel. (+44) 020 7153 5445 Section of Structural Biology Fax. (+44) 020 7153 5456 Chester Beatty Laboratories mobile: (+44) 07787554059 The Institute of Cancer Research 237 Fulham Road London SW3 6JB 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.

3 2

anisotropic refinement based on B values
by Yarrow Madrona 20 Oct '13

20 Oct '13

I'm sure that this has probably already been covered but I am having some problem finding out how to tell phenix to make a decision whether to refine isotropically or anisotropically based on B factor values. Is this something you can do with a python script? Thanks. -Yarrow -- Yarrow Madrona Graduate Student Molecular Biology and Biochemistry Dept. University of California, Irvine Natural Sciences I, Rm 2403 Irvine, CA 92697

3 2

Re: [phenixbb] increasing NCS restraints
by Mark van Raaij 17 Oct '13

17 Oct '13

Oops, 1 vote for making phenixbb default reply to sender only.... Mark van Raaij <mjvanraaij(a)cnb.csic.es> wrote: >Dont worry about it, we are alli diots...even without commuting 4 hrs a day > >Yarrow Madrona <amadrona(a)uci.edu> wrote: > >>Hi Dale, >> >>That is right! And that is what I concluded yesterday yet I somehow forgot >>all that overnight. Starting a new job and commuting 4 hours a day has >>taken a toll on my brain. I am an idiot. >> >>-Yarrow >> >>> >>> I'm puzzled. Didn't you say that the translation of your ncs >>> two-fold was not equal to half the "a" cell edge? If it is not >>> a translation of 1/2 (in fractional coordinates) then it cannot >>> be a crystallographic symmetry element and your space group cannot >>> be P212121. Unless this condition is met, the similarity of your >>> molecules is unimportant. >>> >>> Dale Tronrud >>> >>> On 10/17/2013 12:48 PM, Yarrow Madrona wrote: >>>> I am doing restrained NCS refinement and would like to know how to set >>>> the >>>> target weight of the cartesian restraint in the command line? I also >>>> found >>>> a sigma value in the GUI. I read that I can decrease it in order to >>>> increase the NCS restraints. However, I have seen this under Torsional >>>> restraints. Does this also apply to cartesian restraints? >>>> >>>> I am doing this so that I can look at the differences in position and >>>> B-factors of two molecules in the ASU that I am trying to determine are >>>> related by a perfect or imperfect 2(1) screw axis. (I posted a question >>>> about this on the CCP4 bullitan board.) If the values between molecules >>>> are exactly the same then presumably the space group is actually P212121 >>>> with 1molecule in the asu instead of P21 with two molecules. I already >>>> know that the P21 dimer superimposes almost exactly on the P212121 >>>> monomer >>>> and it's symmetry mate. I just don't know how "perfectly" they align and >>>> I >>>> am hoping this approach will tell me. I hope this makes some sense. >>>> >>>> -Yarrow >>>> >>> _______________________________________________ >>> phenixbb mailing list >>> phenixbb(a)phenix-online.org >>> http://phenix-online.org/mailman/listinfo/phenixbb >>> >>> >> >> >>-- >>Yarrow Madrona >> >>Graduate Student >>Molecular Biology and Biochemistry Dept. >>University of California, Irvine >>Natural Sciences I, Rm 2403 >>Irvine, CA 92697 >> >> >> >>_______________________________________________ >>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

Re: [phenixbb] increasing NCS restraints
by Mark van Raaij 17 Oct '13

17 Oct '13

Dont worry about it, we are alli diots...even without commuting 4 hrs a day Yarrow Madrona <amadrona(a)uci.edu> wrote: >Hi Dale, > >That is right! And that is what I concluded yesterday yet I somehow forgot >all that overnight. Starting a new job and commuting 4 hours a day has >taken a toll on my brain. I am an idiot. > >-Yarrow > >> >> I'm puzzled. Didn't you say that the translation of your ncs >> two-fold was not equal to half the "a" cell edge? If it is not >> a translation of 1/2 (in fractional coordinates) then it cannot >> be a crystallographic symmetry element and your space group cannot >> be P212121. Unless this condition is met, the similarity of your >> molecules is unimportant. >> >> Dale Tronrud >> >> On 10/17/2013 12:48 PM, Yarrow Madrona wrote: >>> I am doing restrained NCS refinement and would like to know how to set >>> the >>> target weight of the cartesian restraint in the command line? I also >>> found >>> a sigma value in the GUI. I read that I can decrease it in order to >>> increase the NCS restraints. However, I have seen this under Torsional >>> restraints. Does this also apply to cartesian restraints? >>> >>> I am doing this so that I can look at the differences in position and >>> B-factors of two molecules in the ASU that I am trying to determine are >>> related by a perfect or imperfect 2(1) screw axis. (I posted a question >>> about this on the CCP4 bullitan board.) If the values between molecules >>> are exactly the same then presumably the space group is actually P212121 >>> with 1molecule in the asu instead of P21 with two molecules. I already >>> know that the P21 dimer superimposes almost exactly on the P212121 >>> monomer >>> and it's symmetry mate. I just don't know how "perfectly" they align and >>> I >>> am hoping this approach will tell me. I hope this makes some sense. >>> >>> -Yarrow >>> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://phenix-online.org/mailman/listinfo/phenixbb >> >> > > >-- >Yarrow Madrona > >Graduate Student >Molecular Biology and Biochemistry Dept. >University of California, Irvine >Natural Sciences I, Rm 2403 >Irvine, CA 92697 > > > >_______________________________________________ >phenixbb mailing list >phenixbb(a)phenix-online.org >http://phenix-online.org/mailman/listinfo/phenixbb

1 0

increasing NCS restraints
by Yarrow Madrona 17 Oct '13

17 Oct '13

I am doing restrained NCS refinement and would like to know how to set the target weight of the cartesian restraint in the command line? I also found a sigma value in the GUI. I read that I can decrease it in order to increase the NCS restraints. However, I have seen this under Torsional restraints. Does this also apply to cartesian restraints? I am doing this so that I can look at the differences in position and B-factors of two molecules in the ASU that I am trying to determine are related by a perfect or imperfect 2(1) screw axis. (I posted a question about this on the CCP4 bullitan board.) If the values between molecules are exactly the same then presumably the space group is actually P212121 with 1molecule in the asu instead of P21 with two molecules. I already know that the P21 dimer superimposes almost exactly on the P212121 monomer and it's symmetry mate. I just don't know how "perfectly" they align and I am hoping this approach will tell me. I hope this makes some sense. -Yarrow -- Yarrow Madrona Graduate Student Molecular Biology and Biochemistry Dept. University of California, Irvine Natural Sciences I, Rm 2403 Irvine, CA 92697

2 2

Re: [phenixbb] Need help with atom selection for anisotropic refinement
by Nathaniel Echols 16 Oct '13

16 Oct '13

I think the problem is with the boolean logic ("and"/"or"). First, the isotropic: optional (element H or element D or resname HOH or chain E or chain F or > chain A or chain B) > This will be processed okay but it's going to select all atoms in chains A/B/E/F in addition to all H/D and waters, which probably isn't what you really want since the anisotropic selection includes chains A and B. I'm not sure what will actually happen here, but it's better not to rely on the atom selections being processed in a particular order (and better not to rely on the code automatically handling the problem of overlapping selections). chain A and chain B and not (element H or element D or resname HOH or chain > E or chain F or chain A resseq 5:11 or chain A resseq 66:72 or chain A > resseq 117:138 or chain A resseq 152:184 or chain A resseq 196:210 or chain > A resseq 249:251 or chain B resseq 3:6 or chain B resseq 66:72 or chain B > resseq 166:210 or chain B resseq 249:251) > The problem here is that you have specified that the atoms match "chain A and chain B", which is impossible. You probably want this: (chain A or chain B) and not (....) -Nat

1 0

phenix.rosetta_refine
by Guenter Fritz 14 Oct '13

14 Oct '13

Dear Phenix developers, I am working on a 3.6 A structure and want to try the latest joint development of phenix and rosetta. Is there are kind of "how_to" available online? Some do's and don'ts? Best, Guenter

3 6

Phenix on Ubuntu12.3
by Swastik Phulera 12 Oct '13

12 Oct '13

Dear All, I am setting up a computer with ubuntu(12.3) installed on it. running "uname -rm" gives me 3.8.0-31-generic x86_64. I am trying to install phenix on this machine, as recommended I downloaded the 64 bit kernel 2.3 (64bit fedora 3 ) command line installer. I then ran ./install and source the required file by: source /usr/local/phenix-1.8.4-1496/phenix_env.sh but some how phenix refuses to start, following is the error phenix Traceback (most recent call last): File "/usr/local/phenix-1.8.4-1496/build/intel-linux-2.6-x86_64/../../phenix/wxGUI2/command_line/main.py", line 22, in <module> run(sys.argv[1:]) File "/usr/local/phenix-1.8.4-1496/build/intel-linux-2.6-x86_64/../../phenix/wxGUI2/command_line/main.py", line 16, in run from wxGUI2 import App File "/usr/local/phenix-1.8.4-1496/phenix/wxGUI2/App.py", line 15, in <module> import wx File "/usr/local/phenix-1.8.4-1496/build/intel-linux-2.6-x86_64/base/lib/python2.7/site-packages/wx-2.8-gtk2-ansi/wx/__init__.py", line 45, in <module> from wx._core import * File "/usr/local/phenix-1.8.4-1496/build/intel-linux-2.6-x86_64/base/lib/python2.7/site-packages/wx-2.8-gtk2-ansi/wx/_core.py", line 4, in <module> import _core_ ImportError: libjpeg.so.62: cannot open shared object file: No such file or directory Any suggestions would be welcome. regards Swastik

2 1

Re: [phenixbb] Problem in reducing R-factor
by Nathaniel Echols 11 Oct '13

11 Oct '13

On Fri, Oct 11, 2013 at 6:19 PM, Ryan Spencer <rspencer(a)uci.edu> wrote: > For high resolution structures, below 1.5 Å, is there a set of > rules > for best practices for refinement? > Not really rules that I'm aware of... just general guidelines, and not well-codified. I've wanted for a long time for us to write some sort of extensive how-to across all resolution ranges, but like everything else, this has to compete with writing code and "real" papers. As Pavel mentioned, these resolution limits are somewhat subjective depending on completeness, how aggressively you process the data, etc. But, a few basic suggestions: - you will probably not want to use real-space refinement after the first couple of rounds - start parameterizing conservatively, then get more aggressive as you near the end of refinement (if nothing else, it's faster this way) - don't use any restraints other than the default geometry restraints plus whatever custom bonds are required by the chemistry (e.g. prosthetic groups) - use hydrogens by the end of refinement, but riding only until you're at subatomic resolution - the automatic weighting is probably going to work great; weight optimization might help for the final round, but probably not necessary - below 1.5Å, try making non-water/non-hydrogen anisotropic (and compare to purely isotropic); below 1.2Å, you can start thinking about making all non-H/D anisotropic. (Also note that heavier atoms can become anisotropic even earlier, so in your case, refining individual anisotropic B-factors for iodine and TLS for everything else would be a reasonable alternative.) - you should probably be able to get rid of all geometry outliers, including clashes - but you may see real (correct) Ramachandran or rotamer outliers at this point - you're inevitably going to see more blobs in the difference map that you can't explain than usual. don't sweat it. - you should definitely see a fair number of alternate conformations - I don't know what is expected statistically but I'd say at least 10% of protein residues will have them - pay attention to waters with close contacts to oxygen - these may be sodium (or other) ions I've attached a few slides that give *very approximate* suggestions - basically just reframing what Pavel and others have said over the years. Although some of these are pretty clear-cut (e.g. you absolutely do not want to use a reference model), you will always need to experiment somewhat to find out what is best for *your* model and data. I think based on your R-factors and statistics that you're headed in the right direction. I'd try to get the clashscore as close to zero as possible, though. I've been toying with different parameters and found that the > weighting didn't make a huge difference, which is probably expected, with a > 1.3 Å dataset. Do stereochemical/ADP weights have any effect? > Do you mean weight optimization? I wouldn't expect it to make much difference. > For the particular data set I'm working with, I'm refining directly > against anomalous data, more specifically against 12 Iodines. Refining the > f' and f" has worked well, generating a beautiful density map. > Good idea. Also consider looking at the anomalous residual map (map_type=anom_residual), which may show weaker anomalous scatterers like chloride or sulfur. (It will be flattened around the iodines if you refine their f'' values.) Refining against the model with hydrogens increased the Rfree/Rwork > values by about 2 points. From what I've read previously I thought refining > hydrogens in a non-"riding" model was only valid below ~1.0 Å. > I'm not sure what you mean here. Did you use the riding model and it increased R-free? Or only when you let them refine individually? The default is now to switch to individual only at 0.9Å or less. > by a few points. For water, I assume that at this resolution there is a > combination of anisotropic (big banana looking density) and isotropic > waters > and that one cannot simply rely entirely on the water_update option. > Yes, at some point you will need to complete the waters yourself. Those bananas may not be simply anisotropic waters, but split sites. You may see other small ligands too. The twin law has worked great, usually dropping Rfree/Rwork values > by 5 points or so. > Just remember to be careful when interpreting the effect of twin law application - I think in this case it's probably quite reasonable, but at lower resolution and with a worse model it can be extremely deceptive. I apologize if any of this is common knowledge and shows my relative > inexperience when it comes to best practices for crystal refinement. > No need to apologize - it is all oral tradition anyways. It is much better to ask questions now than to do something blindly and regret it later. -Nat

1 0

Problem in reducing R-factor
by Singh, Bishal 11 Oct '13

11 Oct '13

Hi all, I am solving a crystal structure at 1.3Å. Completeness of the data-set and Figure of merit at this resolution are 99.5% and 0.62 respectively. Structure has two molecules per asymmetric unit. I used TLS parameters and NCS restraints. WXC_scale and WXU_scale are 0.50 and 1.0 respectively. Rwork and Rfree are 18.09 and 19.62 respectively. I tried with reducing WXC_scale and WXU scale to 0.25 and 0.30 respectively but it did not help. Please suggest me if I need to change some parameters. I shall be grateful for your help. Thanking you, Bishal Singh Graduate student University of Heidelberg, Germany

5 5