- phenixbb - phenix-online.org

Phenix cut out density labels choice
by Mikalai Lapkouski 21 Sep '12

21 Sep '12

Hello, I have a following issue with running Phenix.cutout density either on command line or in GUI. I have multiple columns in the mtz.file and I want to choose map coefficients FOFCWT_d1-8,PHFOFCWT_d1-8 for the difference map to be cut around pdb selection. I run a command: phenix.cut_out_density model.pdb my.mtz cutout_fix_position=True cutout_model_radius=10 cutout_type=model \ map_coeff_labels=FOFCWT_d1-8,PHFOFCWT_d1-8 Problem: The output cutout.mtz has FWT and PHWT (which I did not have at all in the input) and not my difference map coefficients, and obviously the map is not my difference map. In GUI is the same problem, I choose columns I need but it outputs FWT and PHWT again. Please, let me know what I am doing wrong. Thank you very much! Mikalai

3 2

Combine map labels
by Mikalai Lapkouski 21 Sep '12

21 Sep '12

Hello, I have a following question. I have generated several Fo-Fc sa omit maps (FOFCWT, PHFOFCWT) for different portions of a DNA molecule. Now that I have each omit map as a separate file to display the whole DNA omit map I have to load them in coot one by one. Is there a way of how to combine all this maps into one file and have a single name for label columns so that the whole map will open at once? I managed to put them all in one file but column labels are separate: FOFCWT_1, PHFOFCWT_1 FOFCWT_2, PHFOFCWT_2...... Can they all go under FOFCWT, PHFOFCWT? Thank you very much for any suggestion Mikalai

5 5

Specifying union of two chains in pdb_atom_selection
by Francis E Reyes 20 Sep '12

20 Sep '12

Hi all, I issue the following command to pdb_atom_selection but I get 0 atoms selected at the end? $ phenix.pdb_atom_selection refine_to_submit.pdb "chain X and chain B" refine_to_submit.pdb Monomer Library directory: "/usr/local/phenix-1.8-1069/chem_data/mon_lib" Total number of atoms: 3817 Number of models: 1 Model: "" Number of chains: 6 Chain: "B" Number of atoms: 1800 Number of conformers: 1 Conformer: "" Number of residues, atoms: 84, 1800 Classifications: {'RNA': 84} Modifications used: {'rna3p_pyr': 26, 'rna2p_pur': 5, 'rna3p_pur': 44, 'rna2p_pyr': 9} Link IDs: {'rna3p': 70, 'rna2p': 13} Chain: "X" Number of atoms: 1800 Number of conformers: 1 Conformer: "" Number of residues, atoms: 84, 1800 Classifications: {'RNA': 84} Modifications used: {'rna3p_pyr': 26, 'rna2p_pur': 5, 'rna3p_pur': 44, 'rna2p_pyr': 9} Link IDs: {'rna3p': 70, 'rna2p': 13} chain X and chain B 0 atoms selected Thanks! F --------------------------------------------- Francis E. Reyes PhD 215 UCB University of Colorado at Boulder

2 1

Scaling for Isomorphous Difference Maps
by Michael Thompson 20 Sep '12

20 Sep '12

Hello All, I'm currently using the GUI to calculate some Fo-Fo isomorphous difference maps. From looking at the Phenix documentation, it hasn't become apparent to me if the tool that performs this calculation automatically scales the two sets of amplitudes to one another, or if the input mtz files need to be scaled together elsewhere before running the map utility. Any info would be appreciated! Thanks, Mike -- Michael C. Thompson Graduate Student Biochemistry & Molecular Biology Division Department of Chemistry & Biochemistry University of California, Los Angeles miket(a)chem.ucla.edu

3 2

Two versions of TFZ, molecular replacement
by 赵岩 19 Sep '12

19 Sep '12

Hi everyone, I am trying to solve a structure though molecular replacement and I have obtained a suspected solution. The header of the solution as follows, REMARK TITLE [no title set] REMARK Log-Likelihood Gain: 238.884 REMARK RFZ=10.9 TFZ=13.9 PAK=2 LLG=235 TFZ==5.0 LLG=239 TFZ==5.0 REMARK ENSEMBLE 1 EULER 15.53 0.10 344.54 FRAC -0.499 -1.000 -0.499 At first sight, the solution has no problem(TFZ>8), but the second and third TFZ are very low. For a correct solution, the later TFZs always higher than the first TFZ. So I don't konw how to understand the later TFZs and what extent I can refer to them. Thanks a lot! --

5 4

advice for speeding up cif_as_mtz
by Eric Williams 19 Sep '12

19 Sep '12

I'm currently running cif_as_mtz and model_vs_data on the whole PDB as part of my dissertation. Time has become a critical concern, and cif_as_mtz (dev-1141) is running very, very slowly (1-4 minutes per entry). Here's how I'm calling it: phenix.cif_as_mtz <sf_cif_fn> --ignore-bad-sigmas --merge --remove-systematic-absences --map-to-asu --symmetry=<pdb_model_fn> --output-file-name=<sf_mtz_fn> I'm hoping that perhaps someone could help me determine which command line argument(s) could be removed to provide maximum speed-up. Thanks. :) Eric

4 8

B-iso vs. B-aniso
by Yuri 18 Sep '12

18 Sep '12

Dear community, The protein model I am refining has 400 amino acids (3320 atoms). Some real quick calculations tell me that to properly refine it anisotropically, I would need 119,520 observations. Given my unit-cell dimension and space-group it is equivalent to about a 1.24 A complete data set. However, I have had a couple of cases where anisotropic B-factor refinement significantly improved R-work and R-free, while maintaining a reasonable R-gap, for lower resolution models (1.4-1.5 A, around 70,000 reflections). What is the proper way of modelling the B-factors? Any thoughts/opinions/call-outs from the community are welcome. Cheers, -- Yuri Pompeu

4 3

Monitoring MR_ROSETTA progress on a cluster
by Dmitry Rodionov 17 Sep '12

17 Sep '12

Hi Tom, How are you doing? Sorry to bother you again. I got my hands on a login to local cluster (MOAB) and started MR_rosetta there as well. Problem is, I have hard time figuring out what it is doing at this moment and if everything is under control. All I see in the queue are the original submission and a process it spawned less then a minute later. While a day ago there were 30 processes (that's how much I can have). All mr_rosetta.log files end with a list of methods to be run. MR_ROSETTA_1/GROUP_OF_PLACE_MODEL_1/LIST_OF_JOBS_RUNNING says: /sb/home/drodionov/MR_ROSETTA_1/GROUP_OF_PLACE_MODEL_1/RUN_1 /sb/home/drodionov/MR_ROSETTA_1/GROUP_OF_PLACE_MODEL_1/RUN_FILE_1.log MISSING FINAL COUNT: 0 But the file happily exists : GROUP_OF_PLACE_MODEL_1/RUN_FILE_1.log .... Starting job 1000...Log will be: /sb/home/drodionov/MR_ROSETTA_1/GROUP_OF_PLACE_MODEL_1/RUN_1/REBUILD_IN_SETS_1/RUN_FILE_1000.log And that job successfully finished a while ago. Is everything ok with the run? Am I just being impatient? Thanks! Best regards, Dmitry

1 0

Custom restraints editor
by Kristof Van Hecke 17 Sep '12

17 Sep '12

Dear, When I try to define additional 'planarity' restraints via the 'Custom restraints editor', I get the following error (see below): (I'm using the Phenix nightly build dev-1148 on a Mac Book Pro, OS X 10.5.8, but had the same problem with the official release version 1.8) I also tried to add these additional restraints to my restraints .cif file, but then rmsd values increase to unreasonably high values. (I just want to keep 2 chlorine atoms on a phenyl ring better in the same plane) Any ideas please..? Thank you very much Kristof Van Hecke --------------------- untimeError : Syntax error: missing closing quote (python_object, line 1) Traceback (most recent call last): File "/Applications/phenix-dev-1148/phenix/wxGUI2/Programs/Refine/ Dialogs.py", line 296, in OnUpdate new_phil = self.GetPhilObject() File "/Applications/phenix-dev-1148/cctbx_project/wxtbx/ custom_restraints.py", line 646, in GetPhilObject new_phil = self._index.master_phil.format(python_object=self._params) File "/Applications/phenix-dev-1148/cctbx_project/libtbx/phil/ __init__.py", line 1714, in format result.append(object.format(sub_python_object)) File "/Applications/phenix-dev-1148/cctbx_project/libtbx/phil/ __init__.py", line 1714, in format result.append(object.format(sub_python_object)) File "/Applications/phenix-dev-1148/cctbx_project/libtbx/phil/ __init__.py", line 1714, in format result.append(object.format(sub_python_object)) File "/Applications/phenix-dev-1148/cctbx_project/libtbx/phil/ __init__.py", line 1727, in format result.append(object.format(sub_python_object_i)) File "/Applications/phenix-dev-1148/cctbx_project/libtbx/phil/ __init__.py", line 1727, in format result.append(object.format(sub_python_object_i)) File "/Applications/phenix-dev-1148/cctbx_project/libtbx/phil/ __init__.py", line 1094, in format words = type_as_words(python_object=python_object, master=self) File "/Applications/phenix-dev-1148/cctbx_project/libtbx/phil/ __init__.py", line 151, in as_words source_info="python_object") File "/Applications/phenix-dev-1148/cctbx_project/libtbx/phil/ __init__.py", line 60, in tokenize_value_literal tokenizer.settings(contiguous_word_characters="")])) File "/Applications/phenix-dev-1148/cctbx_project/libtbx/phil/ tokenizer.py", line 249, in next where_str(O.source_info, char_iter.line_number))) RuntimeError: Syntax error: missing closing quote (python_object, line 1) Platform info: __FILE__ = /net/patchnose/scratch1/phenix/phenix-dev-1148/ cctbx_project/boost_adaptbx/meta_ext.cpp __DATE__ = Sep 3 2012 __TIME__ = 06:05:56 __i386__ __APPLE_CC__ = 5367 __GNUC__ = 4 __GNUC_MINOR__ = 0 __GNUC_PATCHLEVEL__ = 1 boost::python::cxxabi_cxa_demangle_is_broken(): false __GXX_WEAK__ = 1 __VERSION__ = 4.0.1 (Apple Computer, Inc. build 5367) FE_INEXACT = 32 FE_DIVBYZERO = 4 FE_UNDERFLOW = 16 FE_OVERFLOW = 8 FE_INVALID = 1 FE_ALL_EXCEPT = 61 __SSE2__ = 1 BOOST_LIB_VERSION = 1_51 PY_VERSION = 2.7.3 PYTHON_API_VERSION = 1013 sizeof(bool) = 1 sizeof(short) = 2 sizeof(int) = 4 sizeof(long) = 4 sizeof(std::size_t) = 4 sizeof(void*) = 4 sizeof(long long) = 8 sizeof(float) = 4 sizeof(double) = 8 sizeof(long double) = 16 sizeof(boost::int32_t) = 4 sizeof(boost::uint32_t) = 4 sizeof(wchar_t) = 4 __PTRDIFF_TYPE__ Py_USING_UNICODE sizeof(PY_UNICODE_TYPE) = 2 BOOST_ADAPTBX_TYPE_ID_SIZE_T_EQ_UNSIGNED_LONG boost/TAG: 79991 os.name: posix sys.platform: darwin sys.byteorder: little platform.platform(): Darwin-9.8.0-i386-32bit platform.architecture(): ('32bit', '') floating_point_exceptions.division_by_zero_trapped: False floating_point_exceptions.invalid_trapped: False floating_point_exceptions.overflow_trapped: False number of processors: 2 Memory total: 4,294,967,296 Memory free: unknown import thread: OK "hello" = ['h', 'e', 'l', 'l', 'o'] u"hello" = ['h', '\x00', 'e', '\x00', 'l', '\x00', 'l', '\x00', 'o', '\x00'] u"\u00C5" = ['\xc5', '\x00'] as utf-8 = ['\xc3', '\x85'] LATIN CAPITAL LETTER A WITH RING ABOVE = √Ö

2 2

Sa_Omit_map_DNA_low_resolution
by Mikalai Lapkouski 15 Sep '12

15 Sep '12

Hello, I want to generate a simulated annealing omit map around DNA molecule in my protein-DNA complex for publication. My resolution is very low - 5A. 1. I deleted the whole DNA (947 atoms of DNA against 7543 for protein part) from pdb and ran phenix.refine with simulated_annealing=True. But resulted Fo-Fc in the omitted region is very weak. Is it because of the low resolution and because big part of the model was omitted? 2. Next, I wanted to generate a composit_sa_omit map as in this case DNA would be deleted in small portions and the resulted map would give me a map I can use for publication: phenix.autobuild data=inp.mtz model=whole.pdb (contains DNA) maps_only=True (to skip rebuilding) composite_omit_type=sa_omit omit_selection=DNA.pdb (contains only DNA coordinates). But the job runs too long and still it seems like it is doing rebuilding. Is this a correct command and correct action to create a sa_omit map for the DNA in my complex? Thank you a lot for the help, Nick

3 3