- phenixbb - phenix-online.org

phenix-1.7.2-x refine error
by Christoph Parthier 25 Nov '11

25 Nov '11

Dear Phenix developers, Using the current phenix version, refinements of a 'certain structure' always stop with the error message quoted below, while the refinement of other 'similar structures' work flawlessly. In the error case the three scheduled macrocycles of refinement seem to finish (the refined PDB is there) but at there seems to be problem in finishing the job - the output MTZ file has zero byte length and the last line in refine.log is: ============================== Exporting results ============================== This occurs with several 1.7.2-x versions of phenix, I'm not sure about 1.7.1-x. Interestingly, the 'problem structure' can be refined without problems using the same input files in phenix-1.7-650. Any idea whether this is due to an installation problem or a bug? Thanks, Christoph Here's the error message: MemoryError : Traceback: File "/cryst/phenix-1.7.2-869/cctbx_project/libtbx/thread_utils.py", line 229, in run return_value = self._target(self._args, self._kwargs, self._c) File "/cryst/phenix-1.7.2-869/cctbx_project/libtbx/runtime_utils.py", line 72, in __call__ result = self.target() File "/cryst/phenix-1.7.2-869/phenix/phenix/refinement/runtime.py", line 67, in __call__ call_back_handler=callback) File "/cryst/phenix-1.7.2-869/phenix/phenix/refinement/command_line.py", line 115, in run map_manger.write_files() File "/cryst/phenix-1.7.2-869/phenix/phenix/refinement/driver.py", line 1535, in write_files map_coeff_dataset=self.write_mtz_file(), File "/cryst/phenix-1.7.2-869/phenix/phenix/refinement/driver.py", line 1524, in write_mtz_file params = self.map_params.map_coefficients) File "/cryst/phenix-1.7.2-869/cctbx_project/mmtbx/maps/__init__.py", line 502, in __init__ coeffs = map_coefficients_from_fmodel(fmodel = fmodel, params = mcp) File "/cryst/phenix-1.7.2-869/cctbx_project/mmtbx/maps/__init__.py", line 353, in map_coefficients_from_fmodel fill_mode = "dfmodel") File "/cryst/phenix-1.7.2-869/cctbx_project/mmtbx/f_model.py", line 1887, in electron_density_map fill_mode = fill_mode) File "/cryst/phenix-1.7.2-869/cctbx_project/mmtbx/map_tools.py", line 114, in __init__ update_scaling = update_scaling) File "/cryst/phenix-1.7.2-869/cctbx_project/mmtbx/f_model.py", line 1663, in fill_missing_f_obs fmodel=self, fill_mode=fill_mode, update_scaling=update_scaling) File "/cryst/phenix-1.7.2-869/cctbx_project/mmtbx/missing_reflections_handler.py", line 131, in fill_missing_f_obs fmodel = select_by_map_cc(fmodel = fmodel.deep_copy())) File "/cryst/phenix-1.7.2-869/cctbx_project/mmtbx/missing_reflections_handler.py", line 85, in select_by_map_cc fourier_coefficients = coeffs) File "/cryst/phenix-1.7.2-869/cctbx_project/cctbx/miller/__init__.py", line 3688, in __init__ conjugate_flag=True) Platform info: __FILE__ = /net/longnose/scratch2/phenix/phenix-1.7.2-869/cctbx_project/boost_adaptbx/meta_ext.cpp __DATE__ = Sep 23 2011 __TIME__ = 00:13:14 __i386__ __linux __GNUC__ = 3 __GNUC_MINOR__ = 4 __GNUC_PATCHLEVEL__ = 2 boost::python::cxxabi_cxa_demangle_is_broken(): false __GXX_WEAK__ = 1 __VERSION__ = 3.4.2 20041017 (Red Hat 3.4.2-6.fc3) FE_INEXACT = 32 FE_DIVBYZERO = 4 FE_UNDERFLOW = 16 FE_OVERFLOW = 8 FE_INVALID = 1 FE_ALL_EXCEPT = 61 BOOST_LIB_VERSION = 1_48 PY_VERSION = 2.7.2 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) = 12 sizeof(boost::int32_t) = 4 sizeof(boost::uint32_t) = 4 sizeof(wchar_t) = 4 __PTRDIFF_TYPE__ Py_USING_UNICODE sizeof(PY_UNICODE_TYPE) = 2 gnu libc version: 2.6.1 BOOST_ADAPTBX_TYPE_ID_SIZE_T_EQ_UNSIGNED os.name: posix sys.platform: linux2 sys.byteorder: little platform.platform(): Linux-2.6.22.13-GOLIAT-0.3-bigsmp-i686-with-SuSE-10.3-i586 platform.architecture(): ('32bit', 'ELF') floating_point_exceptions.division_by_zero_trapped: False floating_point_exceptions.invalid_trapped: False floating_point_exceptions.overflow_trapped: False number of processors: 8 Memory total: 4,238,028,800 Memory free: 3,715,174,400 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 = Ã

3 3

回复：phenixbb Digest, Vol 72, Issue 10
by dingding830106＠sina.com 25 Nov '11

25 Nov '11

Dear Thomas: I try running phenix.autosol instead of phenix.solve, and it solve this problem. Thanks a lot! Dingwei ----- 原始邮件 ----- 发件人：phenixbb-request(a)phenix-online.org 收件人：phenixbb(a)phenix-online.org 主题：phenixbb Digest, Vol 72, Issue 10 日期：2011年11月23日 01点08分 Message: 14 Date: Tue, 22 Nov 2011 10:08:46 -0700 (MST) From: "Thomas C. Terwilliger" <terwilliger(a)lanl.gov> To: "PHENIX user mailing list" <phenixbb(a)phenix-online.org> Subject: Re: [phenixbb] solve/resolve problem Message-ID: <52791.128.165.72.7.1321981726.squirrel(a)webmail.lanl.gov> Content-Type: text/plain;charset=iso-8859-1 I'm sorry for the problem! I don't recognize this error (the message solve is giving you is just saying that it has finished, unfortunately). Can you possibly (off-list) send me your input script and the output log file and I'll have a look? Or even better, try running phenix.autosol on the data and perhaps that will solve the problem! All the best, Tom T >> >> >> Dear all >> when I use SOLVE to deal with the SAD data, it will stop with a error >> message: >> >> ???STOP at >> /net/longnose/scratch2/phenix/phenix-1.7.2-869/solve_resolve/solve/solve_2.cpp(3133)??? >> Can you tell me what happen? >> >> wding >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://phenix-online.org/mailman/listinfo/phenixbb

1 0

measuring angle
by Hena Dutta 23 Nov '11

23 Nov '11

Hi, Could you please tell me how I can measure angle between two helices. Can COOT or PYMOL measure this? I have another question. In PDBSET, I can use CHAIN command in shell script to write a pdb file with a desired chain ID. What is the command in shell script to write a pdb file with a desired SEGID? As for example:- pdbset xyzin start.pdb xyzout end.pdb <<eof >& pdbset.log ROTATE [INVERT] [MATRIX] values CHAIN A end eof This will change the chain ID to A after rotation. But, I also like to give a segid in the end coordinate. Thank you in advance and have a nice THANKS GIVING WEEKEND Hena

6 8

reference model restrain error
by Stephan Barden 23 Nov '11

23 Nov '11

Dear all, I'm facing a problem with a fresh (yesterday) installation of phenix 1.7.2-869 refining a structure with 6 molecules per ASU at 3.25 A and using a different reference structure with 1 mol/ASU at higher resolution. I wrote an additional parameter file to assign the reference model to each of the six molecules, which worked well with an older version of phenix. After the upgrade neither taking the additional parameter file, nor auto alignment worked for me. The error comming up is always: RuntimeError : Duplicate sequence number and insertion code. Traceback: File "/sw/phenix-1.7.2-869/cctbx_project/libtbx/thread_utils.py", line 229, in run return_value = self._target(self._args, self._kwargs, self._c) File "/sw/phenix-1.7.2-869/cctbx_project/libtbx/runtime_utils.py", line 72, in __call__ result = self.target() File "/sw/phenix-1.7.2-869/phenix/phenix/refinement/runtime.py", line 67, in __call__ call_back_handler=callback) File "/sw/phenix-1.7.2-869/phenix/phenix/refinement/command_line.py", line 75, in run overwrite = command_line_interpreter.command_line.options.overwrite) File "/sw/phenix-1.7.2-869/phenix/phenix/refinement/driver.py", line 209, in __init__ log=self.log) File "/sw/phenix-1.7.2-869/cctbx_project/mmtbx/torsion_restraints/reference_model.py", line 145, in __init__ self.get_reference_dihedral_proxies() File "/sw/phenix-1.7.2-869/cctbx_project/mmtbx/torsion_restraints/reference_model.py", line 611, in get_reference_dihedral_proxies log=self.log) File "/sw/phenix-1.7.2-869/cctbx_project/mmtbx/torsion_restraints/reference_model.py", line 398, in process_reference_groups raise e Has anyone observed a similar problem? Are there any suggestions? Thanks, StephanB.

2 1

solve/resolve problem
by dingding830106＠sina.com 22 Nov '11

22 Nov '11

Dear all when I use SOLVE to deal with the SAD data, it will stop with a error message: “STOP at /net/longnose/scratch2/phenix/phenix-1.7.2-869/solve_resolve/solve/solve_2.cpp(3133)” Can you tell me what happen? wding

2 1

How to reduce clashscore value
by CelinaSocek 22 Nov '11

22 Nov '11

Hi all, I'm now using phenix to refine a small protein about 300 residues long. The resolution of my data is about 2.9 angstrom. Currently I got a refined model with good R-factors about 23/26, RMSD very reasonable and no Ramanchandran outliers. However, I found that the clashscore is extremely high: 63!! And rotamer outliers are also a lot: 8.5%. I wonder, what kind of stretagy I can use to reduce this? Thank you very much! best, Celina

6 16

test_hhr
by Alexander Batyuk 21 Nov '11

21 Nov '11

Dear all, The phenix_regression.wizards.test_command_line_rosetta_quick_tests command hangs at the test_hhr step in phenix-dev-910 on Mac OS X 10.6.8. What could be the reason for that? What kind of log files should I post? The firewall allowed to fetch files from the net. Best wishes, Alex -- Alex Batyuk The Plueckthun Lab www.bioc.uzh.ch/plueckthun

2 2

Re: [phenixbb] How to reduce clashscore value
by CelinaSocek 20 Nov '11

20 Nov '11

Thank you Nat! I will try to fix as many clashes as I can in coot, and hope that I can manage to reduce it to around 20. ;) Also thank other suggestions and discussions! best, Celina > > Which version are you using? We made some changes over the summer > that improve the geometry, especially the clashscore - these are > present in version 1.7.2. If you're still using an earlier version, > this might largely fix the problem. However, if there is > conformational strain in the model, it is difficult to fix severe > clashes without further rebuilding. You should inspect the clashes in > Coot (the Phenix GUI has shortcuts for this, but Coot alone also has > this capability) and fix as many as possible yourself. It is > recommended that you use explicit hydrogen atoms when doing this, > since the reason for clashes is often not obvious if only heavy atoms > are used. (You don't necessarily need to refine with heavy atoms > though - they're just an aid to visualizing clashes.) > > The rotamer problem may be related; there isn't much you can do with > these except fix them manually, but relieving some of the clashes may > help. The dihedral angle restraints simply aren't that strong, so in > areas of weak density it can be difficult to completely prevent > outliers. > > -Nat > > > ------------------------------ > >

1 0

multiple ligand conformations
by Leonard Thomas 20 Nov '11

20 Nov '11

Hello All, One of the groups here has run into a problem in the refinement of their structure or more specifically a ligand bound to the structure. Through a series of other experiments they know that a photo reaction occurs to the compound when exposed to high flux as found at your everyday synchrotron source. The exposure changes the bonding such that there are two conformations with a bit in between obviously. They collected 360 degrees worth of data and if the data is split you can model the two conformations depending which half of the data you use. They are trying to refine with all of the data and modeling of the ligand into the two conformations. In refmac they can model this with no problem but in Phenix the ligand adopts a chemically wrong conformation, sort of a halfway point between the two chemically relevant conformations. They have used Phenix to do similar refinements on other structures but it appears either the last release or the latest release is causing problems. They have started out with everything from 50% occupancy to shifting it towards on one conformation or the other and the ligand always ends up in the chemically wrong position. Data is in the 1.4-1.5 ang. range. Ideas? Thanks in advance. Len Leonard Thomas Ph.D. Macromolecular Crystallography Laboratory Manager University of Oklahoma Department of Chemistry and Biochemistry Stephenson Life Sciences Research Center 101 Stephenson Parkway Norman, OK 73019-5251 lmthomas(a)ou.edu http://barlywine.chem.ou.edu Office: (405)325-1126 Lab: (405)325-7571

2 1

multiple ligand conformations
by Yuri 20 Nov '11

20 Nov '11

Hi, When you say "wrong postion", what does it mean? are the 2 ligands you are trying to refine being "pushed away" from each other? If that is the case first thing I would check is the ALTLOC comlumn of their pdb file entries. Make sure one is ALTLOC A and the other ALTLOC B. I have refined with ligands in four different positions before using phenix before and got satisfactory results... Hope this helps, -- Yuri Pompeu

1 0