February 2014 - phenixbb - phenix-online.org

FAQ high res refinement
by Eike Schulz 27 Feb '14

27 Feb '14

> > Dear all, > > I am currently refining a rather well resolved protein at 0.97Å. However, even > after many cycles of refinement the R-factors are stuck at around Rwork/Rfree > 14.5/16.5. Considering number of aa, Rmerge (3.5%) and resolution this seems > to be too high. I assume 10 25 macrocycles should be enough for effective > weight refinement in phenix, however no improvement of the R-factors. > > A quick run with default settings in competitor software 'R' yields >2% better > R-factors. > > I would appreciate suggestion how to improve my refinement protocol, values > differing from default are listed below. > > Thank you very much in advance > > Eike > > > > > > > refine { > adp { > individual { > isotropic = element H > anisotropic = not element H > } > group_adp_refinement_mode = one_adp_group_per_residue \ > *two_adp_groups_per_residue group_selection > } > } > main { > apply_overall_isotropic_scale_to_adp = False > nqh_flips = False > ordered_solvent = True > place_ions = True > number_of_macro_cycles = 10 > hydrogen_bonds = True > use_convergence_test = True > target = auto *ml mlhl ml_sad ls > random_seed = 3326496 > wavelength = 0.885611 > nproc = 4 > } > hydrogens { > refine = individual *riding Auto > } > tls { > find_automatically = False > one_residue_one_group = False > } > ordered_solvent { > mode = second_half filter_only every_macro_cycle \ > *every_macro_cycle_after_first > h_bond_min_mac = 1 > h_bond_min_sol = 1 > h_bond_max = 6 > refine_occupancies = True > new_solvent = isotropic *anisotropic > b_iso_min = 0 > } > peak_search { > max_number_of_peaks = 1303 > } > target_weights { > optimize_xyz_weight = True > optimize_adp_weight = True > wxc_scale = 0.2 > wxu_scale = 0.2

2 1

Option for renumber the pdb
by Rojan Shrestha 27 Feb '14

27 Feb '14

Hello, I have found the program phenix.pdbtools. Maybe, it can renumber residues of pdb. I have found the options renumber_residues options but I am unable to see the options for starting residues. I have pdb starts from 1 and want to renumber 1 into 250. Can somebody tell me how to put the option for starting residue? Regards, Rojan

5 5

scaling in phenix.fobs_minus_fobs_map
by Michael Thompson 27 Feb '14

27 Feb '14

Hello, I have a very simple question about phenix.fobs_minus_fobs_map, and it's implementation in the GUI. I just want to make sure I am not skipping an important step in generating my isomorphous difference map. I'm trying to create an Fo-Fo map for two data sets at pretty different resolutions, using phases from a model. I truncated the higher-res data so the high resolution limit matches the lower-res data set. Are these 2 mtz files sufficient for input into phenix.fobs_minus_fobs_map, or do I need to use another utility to scale the two data sets together prior to calculating the map coefficients? I suppose another way to ask my question would be, does phenix.fobs_minus_fobs_map perform scaling of the two Fobs arrays before calculating the new map? Thanks for clarifying, Mike -- Michael C. Thompson Graduate Student Biochemistry & Molecular Biology Division Department of Chemistry & Biochemistry University of California, Los Angeles miket(a)chem.ucla.edu

2 1

problem phenix ensemble
by Firas FADEL 25 Feb '14

25 Feb '14

Hello, I am using phenix ensemble to refine an X-ray structure of a protein (40kDa) at 1.05 A in complex with a ligand. The occupancy of the ligand is 73%. After refining the structure, I found out that the program confused the position of the ligand and put some water molecules. Do you have any suggestion or a command to avoid this problem. If you need details I can send you the mtz and pdb files. Thanks in advance. Firas Fadel Firas FADEL Institut de Génétique et de Biologie Moléculaire et Cellulaire, IGBMC Equipe Alberto PODJARNY, 1 Rue Laurent Fries, 67404 Illkirch-France Email : fadel(a)igbmc.fr<mailto:[email protected]> Tel: + 33 (0) 3 88 65 33 13<tel:%2B%2033%20%280%29%203%2088%2065%2033%2013>

1 0

Announcement: QuantumBio Release of QM-based X-ray Crystallography Plugin v1.0 for Phenix
by Lance Westerhoff 25 Feb '14

25 Feb '14

Hello All- We are pleased to announce the release of our linear scaling, quantum mechanics DivCon Discovery Suite plug-in v1.0 for the Phenix, x-ray crystallography platform. As was detailed in the recent Phenix Newsletter, when this plugin is installed, the conventional refinement functional along with the stereochemical restraints or CIFile for the ligand(s) and for the surrounding active site(s) are replaced with QM data in "real time" during the refinement in order to increase structural accuracy and decrease artificial strain. To see the full press release along with references and links for further information, please take a look at the following website: http://www.quantumbioinc.com/news/show/PhenixAnnouncement To encourage adoption, we have waived the license fee for non-profit, academic use and we have posted our commercial evaluation license as well. These licenses are available at the following link: http://www.quantumbioinc.com/products/software_licensing Certainly, let me know if you have any questions or if you would like to try out the technology on your own structures. -Lance ____________________ Lance M. Westerhoff, Ph.D. President and General Manager QuantumBio Inc. WWW: http://www.quantumbioinc.com Email: lance(a)quantumbioinc.com

1 0

revisiting "Problem with loop fitting "
by Alejandro Buschiazzo 24 Feb '14

24 Feb '14

Dear phenix community I am attempting to use phenix.fit_loops for the first time and getting an error message: "No suitable gap found matching residues" I found a short exchange exactly about this (Jan 2011), but eventually Tom Terwilliger asked for the files directly to see if he could solve it. I can't seem to find a later summary of the solution. (apparently, the most probable reason was some kind of mismatch between PDB numbering and sequence file...which I shouldn't have... I'm probably missing some silly here - I have been fiddling around with the sequence file in many ways...with no success) many thanks for any help! alejandro -- Alejandro Buschiazzo, PhD Research Scientist Unit of Protein Crystallography Institut Pasteur de Montevideo Mataojo 2020 Montevideo 11400 URUGUAY Phone: +598 25220910 int. 120 Fax: +598 25224185 http://www.pasteur.edu.uy/pxf

2 4

Questions about phenix.mr_rosseta
by Yonglin Hu 24 Feb '14

24 Feb '14

Hello List, I am very new to using Phenix/Rosetta to solve MR. I just finished a job by following the procedures described on phenix-online.org using Phenix GUI. The program reported that no solution was found, but when I went through the RUN_FILE_*.log, I found that despite all the log files ended with the following sentences: Stage to relax: refined_rescored_rosetta_solution No solutions to relax No autobuilt solutions to iterate... Finishing up mr_rosetta... one solution had a Rfree of 0.48 with correct number of molecules in ASU. When I refined that solution using phenix, Rfree went down further, proving to be a correct one. So my question is: why phenix.mr_rosetta did not report this solution as a good one? My second question is regarding model rebuild using phenix/rosetta. No matter I used the GUI or command lines copied from phenix-online like the following: phenix.mr_rosetta \ seq_file=~/MA/Mpa.pir \ data=~/MA/MA.mtz \ search_models=~/MA/refine_6.pdb \ already_placed=True \ fragment_files = ~/MA/MA-3mer \ fragment_files = ~/MA/MA-9mer \ rescore_mr.relax=False \ rosetta_models=20 \ ncs_copies=6 \ space_group="C 1 2 1" \ use_all_plausible_sg=False \ nproc=4 the program always fails. The log files have the following messages: Sorry, a subprocess has failed... ...... Sorry: Failed to run REBUILD...see log file in MR_ROSETTA_2/GROUP_OF_ROSETTA_REBUILD_1/RUN_1/REBUILD_IN_SETS_1/RUN_4/WORK_1/rebuild.log They pointed to rebuild.log files, but these files don't contain error message and they ended with the following seemingly innocuous lines: protocols.relax.FastRelax: CMD: accept_to_best -3682.86 0.994359 0.994359 0.44 core.scoring.rms_util: WARNING: In CA_rmsd, residue range 1 to 402 requested but only 401 protein CA atoms found. core.scoring.rms_util: WARNING: In CA_rmsd, residue range 1 to 402 requested but only 401 protein CA atoms found. protocols.relax.FastRelax: CMD: endrepeat -3682.86 0.994359 0.994359 0.44 protocols::checkpoint: Deleting checkpoints of FastRelax protocols::checkpoint: Deleting checkpoints of FastRelax protocols.jd2.JobDistributor: S_MA_R_0005 reported success in 4573 seconds protocols.jd2.JobDistributor: no more batches to process... protocols.jd2.JobDistributor: 5 jobs considered, 5 jobs attempted in 20716 seconds I would appreciate it if you could tell me what could be the causes of the failures. Thank in advance, Yonglin

1 0

stop_for_unknowns=False in new version ?
by Michaela Renner-Schneck 19 Feb '14

19 Feb '14

Hi, I haven’t been running phenix.refine for a couple of weeks and meanwhile we have updated to the latest version. Now I need to re-run an old job but unfortunately the pdb-interpretation option "stop_for_unknowns=False“ which I needed for this seems to be no longer valid. Is there a way to get this option work again ? thx in advance Michaela

2 1

New position
by Gillette, William (NIH/NCI) [C] 19 Feb '14

19 Feb '14

Dear colleagues, We have an opening for a structural biologist here in Frederick, MD. If not interested, please forward this email to those who might be. Thank you. Best regards, Bill <http://jobs.leidos.com/job/Frederick-Lead-Structural-Biology-603181-%28NCI%…> http://jobs.leidos.com/job/Frederick-Lead-Structural-Biology-603181-%28NCI%… ______________________________________________________ Bill Gillette, Ph.D. (Contractor) Sr. Scientist, Protein Expression Laboratory Cancer Research Technology Program Leidos Biomedical Research, Inc. Frederick National Laboratory for Cancer Research (FNL) P.O. Box B, Frederick, MD 21702 Phone: 301-846-1029 William.Gillette(a)nih.gov<gillettew(a)mail.nih.gov> Notice: This communication may contain privileged or other confidential information. If you are not the intended recipient, or believe that you have received this communication in error, please do not print, copy, retransmit, disseminate or otherwise use the information. Please indicate to the sender that you have received this email in error and delete the copy you received.

1 0

phenix.composite-omit-map and PBScluster
by Mark A Saper 19 Feb '14

19 Feb '14

Hi, This is a technical question for cluster experts. I tried to run phenix.composite-omit-map which can divide the calculations into separate jobs, on a very large cluster controlled by the PBS system. I used this command of which the first 2 lines specify to divide the work up into 16 processes, and each process will use 1 of the CPU's in the cluster. phenix.composite_omit_map nproc=16 background=False technology=pbs \ qsub_command="qsub -V -l procs=1,walltime=00:30:00,mem=2000mb,qos=flux -j oe -M saper(a)umich.edu -q flux -m abe " \ omit_type=anneal map_type=2mFo-Dfc output.file_name=composite_omit_map9.mtz \ reflections3_5percent.mtz Refine_122.pdb This is the end of the log file before it aborted: ----------Calculating individual omit maps---------- |-----------------------------------------------------------------------------| | target_work(ml) = 4.65539 r_work = 0.1750 r_free = 0.2066 | |-----------------------------------------------------------------------------| And this is the error log: Traceback (most recent call last): File "/usr/cac/rhel6/lsa/phenix/phenix-dev-1624/build/intel-linux-2.6-x86_64/../../phenix/phenix/command_line/composite_omit_map.py", line 238, in <module> run(sys.argv[1:]) File "/usr/cac/rhel6/lsa/phenix/phenix-dev-1624/build/intel-linux-2.6-x86_64/../../phenix/phenix/command_line/composite_omit_map.py", line 166, in run out=out) File "/usr/cac/rhel6/lsa/phenix/phenix-dev-1624/phenix/phenix/automation/composite_omit_map.py", line 395, in __init__ preserve_exception_message=True) File "/usr/cac/rhel6/lsa/phenix/phenix-dev-1624/cctbx_project/libtbx/easy_mp.py", line 502, in parallel_map display_stderr = not preserve_exception_message, File "/usr/cac/rhel6/lsa/phenix/phenix-dev-1624/cctbx_project/libtbx/queuing_system_utils/processing/__init__.py", line 405, in PBS raise RuntimeError, "PBS does not support synchronous submission" RuntimeError: PBS does not support synchronous submission I know that PBS doesn't support synchronous submission. Perhaps each job (or map region) requires more than one processor? Thanks, Mark _________________________________ Mark A. Saper, Ph.D. Associate Professor of Biological Chemistry, University of Michigan c/o Biophysics, Chem 3040 930 N University Ave Ann Arbor MI 48109-1055 U.S.A. saper(a)umich.edu phone (734) 764-3353 fax (734) 764-3323

2 1