- phenixbb - phenix-online.org

phenix.real_space_refine crashes for very large pdb files
by #AHMED TOFAYEL# 01 Dec '16

01 Dec '16

Dear developers and users, My phenix.real_space_refine run crashes with the error mentioned below when I have used following command to run it on a complete ribosome (containing coordinates for small subunit and large subunit in a single pdb file): phenix.real_space_refine model.pdb map.mrc resolution=4.5 reference_coordinate_restraints.enabled=true reference_coordinate_restraints.selection="chain A and resseq 2:700" reference_coordinate_restraints.sigma=0.05 secondary_structure.enabled=True res.eff > out.out >From the out.out file I can see that, it ran for 5 macro_cycles and started group b-factor refinement which reads (last few lines): |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3809 r_free = 0.3809 target = 0.153923 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 32)----------------| | r_work = 0.3687 r_free = 0.3687 target = 0.144348 restraints weight = 1067.479| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 2; iterations = 37)----------------| | r_work = 0.3724 r_free = 0.3724 target = 0.147257 restraints weight = 482.176| |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 3; iterations = 32)----------------| | r_work = 0.3747 r_free = 0.3747 target = 0.149066 restraints weight = 262.900| |-----------------------------------------------------------------------------| The run crashes at this point.... However, phenix.real_space_refine runs fine when I am refining a protein chain or a complex (for eg: just the large or the small subunit of a ribosome). Does it have to do with my pdb file being so big after merging of both subunits where I have also crossed 99999 atoms and atom serial numbers start from A0000? A related question: What is the latest consensus of using chain IDs and atom serial numbers to use in phenix.real_space_refine for handling large files like ribosome pdbs ? I have found some information on: http://cci.lbl.gov/hybrid_36/ The error in terminal reads: Traceback (most recent call last): File "/usr/local/phenix-1.11.1-2575/build/../modules/phenix/phenix/command_line/real_space_refine.py", line 782, in <module> o = run(sys.argv[1:]) File "/usr/local/phenix-1.11.1-2575/build/../modules/phenix/phenix/command_line/real_space_refine.py", line 518, in __init__ self.caller(func = self.run_macro_cycles, prefix="Macro-cycles") File "/usr/local/phenix-1.11.1-2575/build/../modules/phenix/phenix/command_line/real_space_refine.py", line 547, in caller status = func(prefix=prefix, **kwargs) File "/usr/local/phenix-1.11.1-2575/build/../modules/phenix/phenix/command_line/real_space_refine.py", line 540, in run_macro_cycles log = self.log) File "/usr/local/phenix-1.11.1-2575/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 66, in __init__ self.caller(self.refine_adp, "ADP refinement") File "/usr/local/phenix-1.11.1-2575/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 71, in caller status = func(prefix=prefix) File "/usr/local/phenix-1.11.1-2575/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 317, in refine_adp adp_refinery.refine() File "/usr/local/phenix-1.11.1-2575/modules/cctbx_project/mmtbx/refinement/real_space/adp.py", line 78, in refine b_isos_refined = self.refine_box_with_selected(selection=sel) File "/usr/local/phenix-1.11.1-2575/modules/cctbx_project/mmtbx/refinement/real_space/adp.py", line 56, in refine_box_with_selected fmodel.update_all_scales(update_f_part1=False, apply_back_trace=True) File "/usr/local/phenix-1.11.1-2575/modules/cctbx_project/mmtbx/f_model/f_model.py", line 1452, in update_all_scales self.apply_scale_k1_to_f_obs() File "/usr/local/phenix-1.11.1-2575/modules/cctbx_project/mmtbx/f_model/f_model.py", line 1500, in apply_scale_k1_to_f_obs r_start = self.r_work() File "/usr/local/phenix-1.11.1-2575/modules/cctbx_project/mmtbx/f_model/f_model.py", line 2004, in r_work selection = selection) File "/usr/local/phenix-1.11.1-2575/modules/cctbx_project/mmtbx/f_model/f_model.py", line 1972, in _r_factor f_model = self.f_model_scaled_with_k1_w().data() File "/usr/local/phenix-1.11.1-2575/modules/cctbx_project/mmtbx/f_model/f_model.py", line 1730, in f_model_scaled_with_k1_w data = self.scale_k1_w()*self.f_model_work().data()) File "/usr/local/phenix-1.11.1-2575/modules/cctbx_project/mmtbx/f_model/f_model.py", line 2048, in scale_k1_w selection=selection) File "/usr/local/phenix-1.11.1-2575/modules/cctbx_project/mmtbx/f_model/f_model.py", line 182, in _scale_helper raise RuntimeError("Zero denominator in scale calculation.") RuntimeError: Zero denominator in scale calculation. Best regards, Tofayel

2 2

difference map based on the cryo-EM map
by Smith Lee 01 Dec '16

01 Dec '16

Dear All, Supposed by phenix.real_space_refine I got a PDB of protein-ligand, will you please let me know how to get the subtraction map (difference map?) which only show the density of the ligand, so that I can see whether the ligand was well resolved in the map, or whether the built ligand fits well with the ligand density? Smith

2 1

Incrementing high resolution cutoff in phenix.mtz2map or phenix.cut_out_density
by Andrew Morin 30 Nov '16

30 Nov '16

Hello, I am trying to find a convenient way to generate maps using phenix.mtz2map at increasingly greater values of d_min. For example: I have a .mtz file of a structure solved at 2.4A. I run phenix.mtz2map with d_min=none and obtain a map using all of the resolution data. I would like to also create 10 more maps with d_min = 3.4, 4.4, 5.4, 6.4,...12.4A. An mtz of 1.7A would increment at 2.7, 3.7...11.7A, etc. Might anyone know of a handy way to automate this over a large number of mtz files? ALTERNATIVELY - I am performing a phenix.cut_out_density step prior to this. So if it might be easier to automate a incrementing high resolution cutoff in phenix.cut_out_density, that would also work. Thanks for your help. -Andrew

5 5

Refine only selective chains in a complex using phenix.real_space_refine
by #AHMED TOFAYEL# 30 Nov '16

30 Nov '16

Dear all, Is it possible using phenix.real_space_refine to locally refine a protein or RNA chain (against a cryo-EM map) which is part of a big complex. During the run, I would like to keep most of my complex constant in space except neighboring one or two protein/RNA molecules which I would like to see getting adjusted to remove clashes. I am using the program through command-line. Best regards, Tofayel NTU Singapore

2 2

Can't Run Rosetta_refine
by NISHANT SINGH 29 Nov '16

29 Nov '16

Hi Everyone, I am trying to run Phenix's rosetta refine. I seem to have downloaded and compiled the rosetta package fine, as I got the following statement after compiling was done: scons: done building targets. but I am failing at running: rosetta.build_phenix_interface. It says that command not found. I am attaching my env page with this email. bmbs-imac:rosetta_bin_mac_2016.32.58837_bundle Kumar$ env ROSETTA3_DB=/Applications/rosetta_bin_mac_2016.32.58837_bundle/PHENIX_ROSETTA_PATH//main/database TERM_PROGRAM=Apple_Terminal SHELL=/bin/bash TERM=xterm-256color TMPDIR=/var/folders/cb/5v492std3d31jr42j7h3jm9h0000gp/T/ Apple_PubSub_Socket_Render=/private/tmp/com.apple.launchd.KEC4L2Cjos/Render TERM_PROGRAM_VERSION=387 OLDPWD=/Applications TERM_SESSION_ID=A82EA51D-0920-4F3B-BE03-92AE0C40A44B USER=Kumar SSH_AUTH_SOCK=/private/tmp/com.apple.launchd.tMihawGOC7/Listeners __CF_USER_TEXT_ENCODING=0x1F6:0x0:0x0 ROSETTA_BIN=/Applications/rosetta_bin_mac_2016.32.58837_bundle/PHENIX_ROSETTA_PATH//main/source/bin PATH=/Library/Frameworks/Python.framework/Versions/3.5/bin:/Library/Frameworks/Python.framework/Versions/3.5/bin:/Users/Kumar/miniconda2/bin:/usr/local/bin:/usr/bin:/bin:/usr/sbin:/sbin:/opt/X11/bin PWD=/Applications/rosetta_bin_mac_2016.32.58837_bundle LANG=en_US.UTF-8 XPC_FLAGS=0x0 XPC_SERVICE_NAME=0 SHLVL=1 HOME=/Users/Kumar LOGNAME=Kumar DISPLAY=/private/tmp/com.apple.launchd.BhETYwyInl/org.macosforge.xquartz:0 PHENIX_ROSETTA_PATH=/Applications/rosetta_bin_mac_2016.32.58837_bundle/PHENIX_ROSETTA_PATH/ _=/usr/bin/env bmbs-imac:rosetta_bin_mac_2016.32.58837_bundle Kumar$ env ROSETTA3_DB=/Applications/rosetta_bin_mac_2016.32.58837_bundle/PHENIX_ROSETTA_PATH//main/database TERM_PROGRAM=Apple_Terminal SHELL=/bin/bash TERM=xterm-256color TMPDIR=/var/folders/cb/5v492std3d31jr42j7h3jm9h0000gp/T/ Apple_PubSub_Socket_Render=/private/tmp/com.apple.launchd.KEC4L2Cjos/Render TERM_PROGRAM_VERSION=387 OLDPWD=/Applications TERM_SESSION_ID=A82EA51D-0920-4F3B-BE03-92AE0C40A44B USER=Kumar SSH_AUTH_SOCK=/private/tmp/com.apple.launchd.tMihawGOC7/Listeners __CF_USER_TEXT_ENCODING=0x1F6:0x0:0x0 ROSETTA_BIN=/Applications/rosetta_bin_mac_2016.32.58837_bundle/PHENIX_ROSETTA_PATH//main/source/bin PATH=/Library/Frameworks/Python.framework/Versions/3.5/bin:/Library/Frameworks/Python.framework/Versions/3.5/bin:/Users/Kumar/miniconda2/bin:/usr/local/bin:/usr/bin:/bin:/usr/sbin:/sbin:/opt/X11/bin PWD=/Applications/rosetta_bin_mac_2016.32.58837_bundle LANG=en_US.UTF-8 XPC_FLAGS=0x0 XPC_SERVICE_NAME=0 SHLVL=1 HOME=/Users/Kumar LOGNAME=Kumar DISPLAY=/private/tmp/com.apple.launchd.BhETYwyInl/org.macosforge.xquartz:0 PHENIX_ROSETTA_PATH=/Applications/rosetta_bin_mac_2016.32.58837_bundle/PHENIX_ROSETTA_PATH/ _=/usr/bin/env bmbs-imac:rosetta_bin_mac_2016.32.58837_bundle Kumar$

2 1

kickmaps
by Veronica Fiorentino 25 Nov '16

25 Nov '16

Dear all, I am using the latest Phenix; Calculate Maps does not seem to have an option to for kickmaps as given in the https://www.phenix-online.org/documentation/reference/create_maps.html I have made a screenshot: http://imgur.com/a/upgqy Am I missing something obvious? Alternately, If anyone can share a script I would prefer it! many thanks Veronica

3 4

ensemble refinement crash.
by Rajendra Sharma 24 Nov '16

24 Nov '16

Dear all, I am running ensemble refinement using phenix-1.10. The calculations are running good until the Ensemble reducer where the jobs get killed in some cases automatically without giving output pdb. I am using default setting to run the jobs. pdb_stored_per_block is 50 totaling 500 towards the end. Is the crash happening because of memory? Will reducing the number of pdb per block help? Secondly, how can I restart the killed runs? thanks in advance. Raji *Rajendra SharmaResearch Scholar, PhD.* *Centro Nacional de Biotecnología <http://www.cnb.csic.es/>,* *Universidad Autónoma de Madrid,* *Madrid, Spain.Mobile No. +31-656688390*

2 2

Realspace refinement
by Reza Khayat 23 Nov '16

23 Nov '16

Hi, Does realspace refinement consider temperature factors during refinement? Thanks. Best wishes, Reza Reza Khayat, PhD Assistant Professor Department of Chemistry City College of New York 85 Saint Nicholas Terrace, CDI 2.318 New York, NY 10031 http://www.khayatlab.org/ 212-650-6070

6 15

Postdoctoral positions at the NIH in structural and mechanistic studies on integral membrane proteins
by Anirban Banerjee 19 Nov '16

19 Nov '16

Position Description: The research group of Dr. Anirban Banerjee at the National Institutes of Health (NIH) is seeking candidates for postdoctoral fellows . The broad interests of the lab are in membrane protein structure and function. We combine structural techniques such as macromolecular crystallography with biophysical and biochemical methods such as reconstitution based biochemistry, electrophysiology together with synthetic chemistry to investigate the structural bases of the mechanisms of a number of membrane proteins. The lab has state-of-the-art equipment and facilities. We have dedicated periodic access to the 22-ID beamline at the APS, for crystallographers at the NIH, as well as access to 23-ID and 24-ID beam lines at the APS. Our laboratory is part of the Cell Biology and Neurobiology Branch in the Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD). We also have an affiliation to the National Institute of Neurological Disorders and Stroke (NINDS). Qualifications: The candidates should hold a Ph. D. degree, be author of at least one first author publication and have a strong background in biochemistry or structural biology. Experience with basic molecular biology, protein expression, purification and biochemical characterization are required. Prior experience in membrane protein biochemistry or insect/mammalian expression systems will be added advantages. It is absolutely essential the candidate is motivated and prepared to put in the best efforts for fruition of the project goals. We are a small, friendly lab and ability to work as part of a cohesive group is important. To Apply: Interested candidates can send an e-mail to anirban.banerjee(a)nih.gov with a CV and a summary of previous research experience and future interests.

1 0

Computational Crystallography Newsletter (CCN) submissions
by Nigel Moriarty 18 Nov '16

18 Nov '16

Folks Calling for articles and short communications of interest to structural biologists. The initial deadline for submission in the January 2017 issue is 15 December 2016. However, if you have a brilliant idea during your holiday break please don't hesitate to forward it. The Computational Crystallography Newsletter (CCN) is an electronic newsletter for structural biologists, and is published online every 6 months. Feature articles, meeting announcements and reports can be submitted to me at any time for consideration. Submission of text by email or word-processing files using the CCN templates is requested. Past newsletters and a Microsoft Word template are available at www.phenix-online.org/newsletter. Cheers Nigel --- Nigel W. Moriarty Building 33R0349, Molecular Biophysics and Integrated Bioimaging Lawrence Berkeley National Laboratory Berkeley, CA 94720-8235 Phone : 510-486-5709 Email : NWMoriarty(a)LBL.gov Fax : 510-486-5909 Web : CCI.LBL.gov

1 0