- phenixbb - phenix-online.org

Dealing with malformed ligand geometries in phenix.refine
by Francis E Reyes 26 Aug '11

26 Aug '11

Hi all I've got 3.0A data and I'm ending up with distorted ligand geometry in phenix.refine. The model is complete and fitting the ligand is the final step (refined with individual_sites and individual_sites_real_space as suggested by an earlier posting). I've tried with a wxc_scale of 0.001 and it still looks wonky (tetrahedral C's not being so, etc). My cif is coming from phenix's internal db. Any suggestions? F --------------------------------------------- Francis E. Reyes M.Sc. 215 UCB University of Colorado at Boulder

3 2

maybe not a phenix issue, but with mr_rosetta:
by Francis E Reyes 26 Aug '11

26 Aug '11

What exactly is the nature of the fragment files used for MR/rebuilding? I'm an RNA guy and may have to generate them myself for certain projects. F --------------------------------------------- Francis E. Reyes M.Sc. 215 UCB University of Colorado at Boulder

2 1

unstable phenix?
by Kenneth A. Satyshur 25 Aug '11

25 Aug '11

We are running phenix (743) on a RHEL5 dual quad, multithread, making a composite omit map pf a very large structure (750 AA) and the system keeps getting locked up. Have anyone else seen this type of problem? Does multiprocessing work well on RHEL or just FC versions? We needed to install the FC8 version (/phenix-1.7.1-743/build/intel-linux-2.6-x86_64) to get it to work. So maybe we are using the wrong version? Shoiuld I recompile it for our system? Why is the gui always open when we are running long jobs (many days)? We need to put this stuff in the background and log out so someone else can run. Should scripts be used instead of the gui? Great program but too unstable. kas -- Kenneth A. Satyshur, M.S.,Ph.D. Associate Scientist University of Wisconsin Madison, Wisconsin 53706 608-215-5207

4 4

Weird Hydorgen refinement
by Yuri 23 Aug '11

23 Aug '11

Hi guys, I DEFINETELY ran refinement after I removed the hydrogens. My Rfree went up from 0.159 to about 0.167. Thanks for the tips. Ill give pdb.tools a try... -- Yuri Pompeu

3 2

Reference model restraints in autobuild?
by Francis E Reyes 23 Aug '11

23 Aug '11

Hi all While it is possible to generate refinement parameters to pass into autobuild, how does one pass the reference structure into the autobuild jobs if one were to refine using geometry restraints from the reference structure? (command line autobuild) If I can't supply the reference pdb, is there a phenix cmd utility that'll spit out all the restraints given a pdb so I can supply them in the refinement parameters file? F --------------------------------------------- Francis E. Reyes M.Sc. 215 UCB University of Colorado at Boulder

1 1

Re: [phenixbb] Weird Hydorgen refinement
by Yuri 22 Aug '11

22 Aug '11

Hi Pavel, The pdb did have anisou records for H's. I believe that is because in some previous runs I had only selected chain 'A' - I left out "and not elemnt H" - which I believe caused the H's to get refined anisotropically. I am trying to get rid of the ANISOU records and Ill run refinement again. ps. I am getting ready to deposit the coordinates and I dont want to keep hydrogens in the model, how can I get rid of them? I ve been using pyMol to remove them, but pymol gets rid of all the info -including symmetry -that phenix adds... Thanks -- Yuri Pompeu

3 2

Weird Hydorgen refinement
by Yuri 22 Aug '11

22 Aug '11

Hello everyone. I was trying to refine my model anisotropically. I selected element NA or (chain A and not element H). I looked at the output pdb shows that hydrogens were being refined anisou as well as the protein. I believe its because technically they are also in chain A, but shouldnt the default be never refine H's anisotropically, especially at 1.5A ?. -- Yuri Pompeu

2 1

autosol error in GUI
by r n 22 Aug '11

22 Aug '11

Hi all I downloaded the Phenix developing version 842, autosol gui running gave me the following error, "The key gui_output_dir is not in var_type_dict....command line wizzard needs conversion in the autosol.py script uing translating inputs or else in excluded inputs....wizard setVars routine needs this variable set". I have not had this error earlier version. any help? Thanks ram

2 1

pdb conversion to 3.2 format
by Luca Pellegrini 22 Aug '11

22 Aug '11

Hi, I am trying to use iotbx.pdb_remediator to convert PDB atom names to 3.2 format (as suggested by Nat Echols in a recent post to Coot mail-list). Unfortunately I get the following error message: lp212-mp2281 58% iotbx.pdb_remediator file_name=refine.pdb output_file=refine_3.2.pdb_file Traceback (most recent call last): File "/usr/local/cctbx/cctbx_sources/iotbx/command_line/pdb_remediator.py", line 54, in <module> run(sys.argv[1:]) File "/usr/local/cctbx/cctbx_sources/iotbx/command_line/pdb_remediator.py", line 51, in run remediator.remediator(params) File "/usr/local/cctbx/cctbx_sources/iotbx/pdb/remediation/remediator.py", line 257, in remediator user_dict) File "/usr/local/cctbx/cctbx_sources/iotbx/pdb/remediation/remediator.py", line 95, in build_hash test=os.path.isfile), "rb") TypeError: coercing to Unicode: need string or buffer, NoneType found It looks like a python problem, but I am afraid my insight ends here. Can anyone tell me how to fix this, please? I am not sure if it matters, but iotbx_pdb_remediator was not in my phenix installation (PHENIX-1.7.1-743 on MacOSX 10.6.8), so I installed cctbx separately. Thanks! Luca Luca Pellegrini Department of Biochemistry University of Cambridge 80 Tennis Court Road Cambridge CB2 1GA - UK Email: lp212(a)cam.ac.uk Tel: 0044-1223-760469 Fax: 0044-1223-766002 Sanger building, room 3.59

2 1

MR problem
by intekhab alam 22 Aug '11

22 Aug '11

Hi there I have a data of heterotrimer protein complex at 2.9A resolution. One protein consists of two domains. I tried phaser as well as molrep which gives a solution with only one domain. Rotation function and translation function were found to be fine for these solutions. I tried to find missing domain of the protein after fixing one of the domain or the other partner proteins using whole part or various truncations of missing domain. I also tried to find and build missing domain using Rosetta with the solutions of Molrep or Phaser as template. But, there were no solutions, or the solutions are clashed with other domain or proteins. Furthermore, R and Rfree is 30 and 40, respectively, and I could not reduce them further. There was almost no electron density map in the empty space so that I could not model manually. Plz guide me ,how can i look for the missing domain in the protein. regards -- INTEKHAB ALAM LABORATORY OF STRUCTURAL BIOINFORMATICS KOREA UNIVERSITY, SEOUL

2 1

matplotlib
by Ed Pozharski 20 Aug '11

20 Aug '11

This is really a rather esoteric issue, but I feel I need to mention this to save someone, depending on the situation, anywhere between five minutes of googling or two hours of hair-pulling. In a nutshell, phenix gui may in some circumstances screw up other programs that use matplotlib. I can't track down which specific phenix gui action causes it, but it's safe to assume that whenever matplotlib libs are used the subsequent calls to logarithmic scale plots (and likely other things) throw this error TypeError: coercing to Unicode: need string or buffer, dict found The solution is to erase the ~/.matplotlib/fontList.cache. The likely cause is that phinix gui calls on bundled matplotlib which is different from one I have installed (not to mention that I am using Lucid (because it's LTS) which has python 2.6 and not the 2.7 that is bundled with phenix). However, it still writes into the same ~/.matplotlib folder, thus I end up with incompatible data. Certainly, the problem will be gone when matplotlib gets bumped up to 1.0.1 in next Ubuntu release. This is yet another example of why the standalone installation approach is ideologically objectionable on modern Linux. But of course, the practical advantage gained by not having to package the software for any possible OS flavor/version users may choose outweighs the lower risks of package incompatibility and the reduced size of the packaged product. -- Oh, suddenly throwing a giraffe into a volcano to make water is crazy? Julian, King of Lemurs

5 5

formyl-methionine
by Yuri 19 Aug '11

19 Aug '11

Hello Everyone, Stupid question. I see density at my N-terminus that really looks like formyl-Met. I want to model it and refine. I notice REEl has FME ligand for the formyl-methionine, but I dont want to treat it as a "ligand" since it is covalently attached to SER 2. What is the best way to modify the MET residue? Cheers -- Yuri Pompeu

3 3

Definition of "Saenger Class"
by jens j birktoft 19 Aug '11

19 Aug '11

Hi All, Could somebody direct me to "saenger class" definitions used in base-pair restraints for DNA/RNA. Cheers Jens *+++++++++++++++++++++++++++++++++++++++++++++++++++++++* * * *Jens J. Birktoft e-mail: birktoft(a)nyu.edu <jens.knold(a)gmail.com> New York University, Dept. Chemistry Phone: 212-749-5057 *

2 1

phenix.phaser multiple pdbs
by G Y 17 Aug '11

17 Aug '11

Dear all, I am writing a bash script for running phenix.phaser automatically on an ensemble of multiple pdbs. These pdbs have been superimposed by other tools before running phaser. The following is what I wrote for the ensemble part in params.eff file: ensemble { model_id = a1 coordinates { pdb = model1.pdb pdb = model2.pdb pdb = model3.pdb pdb = model4.pdb rmsd = None identity = 0.9 } It runs well. When it finished, I want to confirm that it really run MR on all four models or on a ensemble of these models. However in .log file, I found at the refinement part it shown me: ---------- REFINEMENT ---------- There is 1 solution to refine Refining solutions 0% 100% |==| DONE Configuring ensembles --------------------- PDB file # 1: model1.pdb This pdb file contains 1 models The (input) RmsD of this model with respect to the real structure is 0.8 **************************** Does this mean that phenix.phaser only run MR with the first model? Does it due to the way I wrote in my params.eff file? What would be the correct way? Many thanks! best, G

3 5

changing project directory
by Jens Kaiser 16 Aug '11

16 Aug '11

Hi all, I have project data on external hard drives. As I use several different file systems, and mount the drives in random order, I cannot rely on the device being mounted in the same mount point. In order to keep paths consistent, I link the project directories into my home. This used to work fine. I recently created a project, and phenix seemed to have used the real path instead of the link-path (That's not the problem, yet). I tried to manually edit ~/.phenix/project_db.phil, but it still tries to write in a non-existing directory. A naive 'grep -R "disk-1" ~/.phenix' did not reveal any obvious file either. Where can I change the directory associated with a project {per,con}sistently? Thanks for your time, Jens

2 1

nucleic acid secondary restraints breaking in 1.7.1-743
by Francis E Reyes 16 Aug '11

16 Aug '11

Hi all Traceback (most recent call last): File "/usr/local/phenix-1.7.1-743/phenix/phenix/command_line/refine.py", line 11, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/usr/local/phenix-1.7.1-743/phenix/phenix/refinement/command_line.py", line 67, in run overwrite = command_line_interpreter.command_line.options.overwrite) File "/usr/local/phenix-1.7.1-743/phenix/phenix/refinement/driver.py", line 112, in __init__ hbond_params=self.params.hydrogen_bonding) File "/usr/local/phenix-1.7.1-743/cctbx_project/mmtbx/secondary_structure/__init__.py", line 494, in create_hbond_proxies log=log) File "/usr/local/phenix-1.7.1-743/cctbx_project/mmtbx/secondary_structure/__init__.py", line 286, in hydrogen_bond_proxies_from_selections log=log) File "/usr/local/phenix-1.7.1-743/cctbx_project/mmtbx/secondary_structure/base_pairing.py", line 415, in create_hbond_proxies if distance_values[i][2] != '_': IndexError: tuple index out of range When specifying a single base pair when using secondary_structure nucleic acids restraints Any ideas? Thanks! F --------------------------------------------- Francis E. Reyes M.Sc. 215 UCB University of Colorado at Boulder

2 6

modelling of disordered protein
by intekhab alam 16 Aug '11

16 Aug '11

Hi I solved an RNA polymerase structure in its native as well as in complex with chemical ligands (2.5A and C2 space group). The structure has been solved by MR a protein having 60% identity. The ASU has a monomer and dimer. Recently i complexed this RNA polymerase (protein A) with a highly flexible protein B and got the crystal. These crytsal also have C2 space group (resolution 2.3A) having similar cell parameters but has psudotranslation (17.34%). I tried MR using rosetta. HHpred server selected one of my deposited structure in PDB (3NAI 100% identity) as template. Rosetta MR gave 2 solutions one of which is dimer-monomer arrangement while the other one is a trimer. The values of these solution are below. Solution 1: trimer REMARK Log-Likelihood Gain: 10662.548 REMARK RFZ=16.5 TFZ=14.4 PAK=0 LLG=959 TFZ==23.2 RFZ=0.0 TFZ=48.9 PAK=0 LLG=3195 TFZ==55.6 (& TFZ==52.6) LLG+=(3191 & 7008) LLG=7024 TFZ==72.0 LLG=10663 TFZ==90.1 Solution 2: Monomer and dimer in ASU REMARK Log-Likelihood Gain: 2120.965 REMARK RFZ=39.7 TFZ=21.9 PAK=0 LLG=811 TFZ==37.7 LLG=2121 TFZ==43.6 REMARK ENSEMBLE ensemble_1 EULER 204.07 178.91 204.37 FRAC 0.499 Map and model fits well in both the cases with R and R free 26 and 30% respectively. My queries: 1. which one will be the correct solution? Is it possible that the protein B has resulted in the trimeric arrangement? 2. i cannot find any continuous 2Fo-Fc map for proteinB, but there are some continuous density at 0.7 which is very ambiguous. how can i improve the map so that i can model some fragments of my proteinB (protein B is of 73 amino acids and is predicted to be highly flexible and disordered. ((biochemical as well as bioinformatical analysis done )) I already have tried autobuild but it fails to build even a partial fragment of the proteinB. Kindly suggest how to proceed , or how can i prove the presence of protein B in the structure. Regards INTEKHAB ALAM LABORATORY OF STRUCTURAL BIOINFORMATICS KOREA UNIVERSITY, SEOUL

2 1

LigandFit error
by Priit Eek 14 Aug '11

14 Aug '11

Hi, I'm getting this error when trying to run LigandFit from the GUI: The key gui_output_dir is not in var_type_dict…this means the command-line Wizard needs a conversion in the ligandfit.py script under 'Translating inputs' or else in excluded_inputs or else the Wizard SetVars routine needs this variable set Is this a bug? I am using the 64-bit dev-842 version on OS X 10.7 Lion. Best regards, Priit

2 1

Unexpected Occupancy refiniement for ligands...
by Yuri 13 Aug '11

13 Aug '11

Hey, So I am finishing up a model at 1.5 A and I tried to refine my ligands (3) and waters (468) in addition to altloc residues. The residues refine as expected e.g OccA + OccB = 1.0 and all atoms in A have OccA and all the atoms in B have OccB. My ligands on the other are being refined as the following: HETATM 6867 C1 MLA L 3 32.357 -28.612 7.978 0.93 23.77 C C HETATM 6868 O1A MLA L 3 31.209 -28.949 7.502 0.85 24.79 C O HETATM 6869 O1B MLA L 3 32.454 -28.160 9.141 0.70 22.70 C O HETATM 6870 C2 MLA L 3 33.553 -28.727 7.095 1.00 25.46 C C HETATM 6871 C3 MLA L 3 34.774 -28.070 7.627 0.92 28.45 C C HETATM 6872 O3A MLA L 3 34.830 -27.664 8.808 1.00 30.15 C O HETATM 6873 O3B MLA L 3 35.775 -27.948 6.840 1.00 30.23 C O HETATM 6874 HC21 MLA L 3 33.321 -28.281 6.139 1.00 30.56 C H HETATM 6875 HC22 MLA L 3 33.760 -29.776 6.943 1.00 30.56 C H Each atom is getting an individual occupancy assigned, which physically is impossible :) Any light? Could it be my selection syntax? Corrigendum: I think I sent a bug email last night about an error when writing mtz files in refinement. It was an accident. My naming convention was bad, thats why it crashed. Sorry -- Yuri Pompeu

3 4

R/Rfree discrepancy Phenix vs Refmac
by Sabine Schneider 12 Aug '11

12 Aug '11

Hello, I am refining a structure at 2.1A solved by MR. For curiosity I'm testing things like 'shaking' coordinates, simulated annealing, refinement Refmac vs Phenix and so on, to see what influence that has on stats, maps etc. For instance after MR I did a bit of shaking the coordinates with pdbset (noise 0.1), followed by simulated annealing in Phenix. Phenix states after SA: Start R-work = 0.2671, R-free = 0.2992 Final R-work = 0.2312, R-free = 0.2666 When I use the output pdb of phenix directly in Refmac (with same mtz as input for Phenix) Refmac tells me: Initial R factor 0.2392 R free 0.2887 So I am quite puzzled about the discrepancy. Or can someone tell me if I made an error in reasoning somewhere? Thanks a lot for the help! Sabine

6 8

getting the crystal form
by Bradley Hintze 11 Aug '11

11 Aug '11

I was wondering if there was a easy way to grab the crystal form for a given PDB code? -- Bradley J. Hintze Duke University Graduate Student Department of Biochemistry

2 1

phenix.grow_density documentation
by Joern Krausze 11 Aug '11

11 Aug '11

Dear all, could anybody please provide me with a documentation for phenix.grow_density? I did'nt find one on the Phenix homepage or elsewhere in the depths of the internet. Thanks in advance, Joern ****************************************** Address: Joern Krausze Molecular Structural Biology Helmholtz Centre for Infection Research Inhoffenstrasse 7 38124 Braunschweig Germany Email: joern.krausze(a)helmholtz-hzi.de Phone: +49 (0)531 6181 7023 (office) +49 (0)531 6181 7020 (lab) ****************************************** Helmholtz-Zentrum für Infektionsforschung GmbH | Inhoffenstraße 7 | 38124 Braunschweig | www.helmholtz-hzi.de Vorsitzende des Aufsichtsrates: MinDir’in Bärbel Brumme-Bothe, Bundesministerium für Bildung und Forschung Stellvertreter: MinDirig Heiko Gevers, Niedersächsisches Ministerium für Wissenschaft und Kultur Geschäftsführung: Prof. Dr. Dirk Heinz; Ulf Richter, MBA Gesellschaft mit beschränkter Haftung (GmbH) Sitz der Gesellschaft: Braunschweig Handelsregister: Amtsgericht Braunschweig, HRB 477

2 2

Large size display system
by jens j birktoft 11 Aug '11

11 Aug '11

Sorry if this topic is bit off-list We are looking for some advice regarding a presentation system that permit the display in high definition of interactive molecular graphic sessions (incl. PHENIX,COOT, PyMOL. CHIMERA etc ) to our research group (10-12 people). Ideally we would like a 32"-36" LCD monitor coupled to a computer equipped with the proper graphics card for the large display as well as 24" monitor . Any suggestions would be appreciated Cheers Jens *+++++++++++++++++++++++++++++++++++++++++++++++++++++++* * * *Jens J. Birktoft, Department of Chemistry, New York University e-mail: birktoft(a)nyu.edu <jens.knold(a)gmail.com> slow-mail: 350 Central Park West, Suite 9F, New York, NY 10025 Phone: 212-749-5057* *+++++++++++++++++++++++++++++++++++++++++++++++++++++++* * * *Jens J. Birktoft e-mail: jens.knold(a)gmail.com slow-mail: 350 Central Park West, Suite 9F, New York, NY 10025 Phone: 212-749-5057 *

3 2

Re: [phenixbb] rosetta fragment files -- corrected command
by Thomas C. Terwilliger 11 Aug '11

11 Aug '11

Hi again David, Sorry, corrected command below. Can you try the version of the script I put in phenix: phenix.adjust_robetta_resid fragments_in fragments_adjusted <start_resno> Let me know if that doesn't do it! All the best, Tom T >> Hi, >> >> I'm also trying to use rosetta with a protein that has more than 650 aminoacids. I have found the answer on last June to the similar question about the fragment files but I couldn't make the script to change the aminoacids numbers work so far >> http://www.phenix-online.org/pipermail/phenixbb/2011-June/017268.html >> >> However I think the script got out of format and I can't make it work. I have python installed however I'm a total beginner on python. >> >> Is it possible that I could have the scrip file so I don't lose the proper format? >> Also, should I replace something inside (like put the file name somewhere) so it works? >> >> Thanks in advance >> >> David Neves >> Institut de Biologie Structurale >> Grenoble France >> >>

1 0

rosetta fragment files
by David NEVES 11 Aug '11

11 Aug '11

Hi, I'm also trying to use rosetta with a protein that has more than 650 aminoacids. I have found the answer on last June to the similar question about the fragment files but I couldn't make the script to change the aminoacids numbers work so far http://www.phenix-online.org/pipermail/phenixbb/2011-June/017268.html However I think the script got out of format and I can't make it work. I have python installed however I'm a total beginner on python. Is it possible that I could have the scrip file so I don't lose the proper format? Also, should I replace something inside (like put the file name somewhere) so it works? Thanks in advance David Neves Institut de Biologie Structurale Grenoble France

2 1