- phenixbb - phenix-online.org

Phenix.refine RuntimeError
by JXQI 10 Jan '10

10 Jan '10

Dear PHENIXBBs, I tried to refine a N-linked glycan using phenix.refine and used the data_link=NAG-ASN for ASN-NAG connection. But I got a RuntimeError of unknown atom_id: HD22. From the mon_lib_list.cif file of Phenix-1.5-2, the atom_id HD22 is already defined. Any suggestions are appreciated. Thanks J.X. QI ********************************************************* Monomer Library directory: "/home/jxqi/152Phenix/phenix-1.5-2/chem_data/mon_lib" Total number of atoms: 11844 apply_cif_link: data_link: NAG-ASN mod_id_1: DEL-O1 mod_id_2: DEL-HD22 residue_selection_1: chain A and resname ASN and resid 170 residue_selection_2: chain A and resname NAG and resid 1 apply_cif_link: data_link: NAG-ASN mod_id_1: DEL-O1 mod_id_2: DEL-HD22 residue_selection_1: chain C and resname ASN and resid 170 residue_selection_2: chain C and resname NAG and resid 1 apply_cif_link: data_link: NAG-ASN mod_id_1: DEL-O1 mod_id_2: DEL-HD22 residue_selection_1: chain E and resname ASN and resid 170 residue_selection_2: chain E and resname NAG and resid 1 Number of models: 1 Model: "" Number of chains: 7 Chain: "A" Number of atoms: 2559 Number of conformers: 1 Conformer: "" Traceback (most recent call last): File "/home/jxqi/152Phenix/phenix-1.5-2/phenix/phenix/command_line/refine.py", line 11, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/home/jxqi/152Phenix/phenix-1.5-2/phenix/phenix/refinement/command_line.py", line 68, in run log=log) File "/home/jxqi/152Phenix/phenix-1.5-2/phenix/phenix/refinement/command_line.py", line 227, in __init__ self.pdb_file_names) File "/home/jxqi/152Phenix/phenix-1.5-2/cctbx_project/mmtbx/utils.py", line 907, in process_pdb_files stop_if_duplicate_labels = stop_if_duplicate_labels) File "/home/jxqi/152Phenix/phenix-1.5-2/cctbx_project/mmtbx/utils.py", line 940, in _process_pdb_file log = self.log) File "/home/jxqi/152Phenix/phenix-1.5-2/cctbx_project/mmtbx/monomer_library/pdb_interpretation.py", line 2995, in __init__ log=log) File "/home/jxqi/152Phenix/phenix-1.5-2/cctbx_project/mmtbx/monomer_library/pdb_interpretation.py", line 2105, in __init__ log=log) File "/home/jxqi/152Phenix/phenix-1.5-2/cctbx_project/mmtbx/monomer_library/pdb_interpretation.py", line 1575, in __init__ prev_mm=prev_mm) File "/home/jxqi/152Phenix/phenix-1.5-2/cctbx_project/mmtbx/monomer_library/pdb_interpretation.py", line 543, in __init__ mod_mod_id=self.mon_lib_srv.mod_mod_id_dict[apply_data_mod]) File "/home/jxqi/152Phenix/phenix-1.5-2/cctbx_project/mmtbx/monomer_library/pdb_interpretation.py", line 819, in apply_mod mod_mon = self.monomer.apply_mod(mod_mod_id) File "/home/jxqi/152Phenix/phenix-1.5-2/cctbx_project/mmtbx/monomer_library/cif_types.py", line 359, in apply_mod result.delete_atom_in_place(mod_atom.atom_id) File "/home/jxqi/152Phenix/phenix-1.5-2/cctbx_project/mmtbx/monomer_library/cif_types.py", line 264, in delete_atom_in_place raise RuntimeError("delete_atom_in_place: unknown atom_id: %s" % atom_id) RuntimeError: delete_atom_in_place: unknown atom_id: HD22 **********************************************************

2 1

Hello
by Pankaj Chauhan 06 Jan '10

06 Jan '10

Hello, While during NCS based restrained refinement, how can we take two chains as single translation unit as the pdb model is a dimer-dimer. Thanks -- Pankaj Chauhan Senior Research Fellow, PhD Scholar Crystallography Laboratory Division of Protein Science & Engg. IMTECH , Chandigarh, INDIA

1 0

coot - display map
by r n 06 Jan '10

06 Jan '10

Hi I am using Phenix GUI, after refinement , tried to load the coot interface in Phenix, it opens up the coot, but did not load the map or model file. I am using coot0.6 version, Is there any place in phenix enviroment that I could mention the path for coot? Thanks

2 4

refinement
by r n 05 Jan '10

05 Jan '10

Hi I am doing refinement using phenix. I did simulated annealing, Rfree goes up to 10% more than before. (40%) and R stays at 35%. If use previous phenix refinment file for same pdb and data , R and Rfree is 25/27%. I thought even different r-free set after simulated annealing and 7 cycles of refinement in phenix should bring Rfree down, but not ? But same data set and model , in CNS, R/Rfree 26/28%. Any suggestions will be appreciated. I am using phenix-dev249 - linux. thanks ram

2 1

text output for model quality and validation report at phenix
by Ravindra D. Makde 04 Jan '10

04 Jan '10

Hi all, I will greatly appreciate the answers to following query: How can I get the text or excel output of the model quality and validation report of a model after phenix refinement? I am specially interested in real space CC of the various residues. I am using phenix.refine_gui at phenix-1.4-162. I have searched all output files, but couldn't find. Thanks in advance. Ravi Ravindra D. Makde, PhD Postdoctoral fellow, Tan lab, The Department of Biochemistry and Molecular Biology, The Pennsylvania State University, University Park, PA 16802 USA

3 2

optimize phenix
by r n 04 Jan '10

04 Jan '10

hi all I am trying to optimize the model, in which secondary structures are distorted and I did pdbtools , geometry regulatization , after optimization still the helices are distorted. Is it any tools in the phenix that could optimize the secondary structures alone in the pdb? Thanks ram

3 4

some sort of occupancy refinement
by Bryan W. Lepore 01 Jan '10

01 Jan '10

question is after the setup here : i have a ligand with multiple conformers. i gave each conformer a different chain label. refining occupancy for the overlapped parts kicks that part out of the density. i have remove_selection for the parts that are 1.0 and they seem OK. a single conformer stays in the density. i tried starting with and without site or occupancy refinement, no help. i do not want group occupancy refinement b/c some parts do not overlap - i.e. they "should" be occupancy 1.0 -- i tried it anyways with constrained_group =... to check. is the chain label an improper way to do this? -bryan p.s: Happy New Year with Blue Moon Eclipse (~10degW->180degE long. eastern HS)!

4 3

high resolution Anomalous map
by vennila Natesan 31 Dec '09

31 Dec '09

Hello all... How do we create anomalous map using high resolution reflections? For example, if we have 2.2A resolution data (i.e 20-2.2A), i want to know how to create the anomalous map for the same dataset but from 2.4-2.2A resolution range. Vennila.

4 8

phenix.refine change: support for residues on symmetry elements
by Ralf W. Grosse-Kunstleve 28 Dec '09

28 Dec '09

The current development version of phenix.refine (dev-268 or later) supports residues located on symmetry elements, e.g. SO4 on a three-fold axis. phenix.refine now computes structure factors slightly differently compared to older versions. The new rules are: - Atoms with occupancy == 1 are subject to automatic handling of the symmetry. Example: An atom near a three-fold axis (distance to symmetry copy <= 0.5 A) is moved onto the axis and constrained in refinement to stay on the axis. The "special position factor" (1/3 is this case) is taken into account internally. The factor does not appear in the occupancy. - Atoms with occupancy != 1 are assumed to take the special position factor into account in the occupancy. Automatic handling of the symmetry is turned off. Example: A disordered SO4 molecule near a three-fold axis should have occupancies of 0.33 for all atoms. - If all atoms in a residue have occupancies != 1, all nonbonded symmetry interactions between the atoms in the residue are disabled. These three rules are meant - to bring phenix.refine in line with other programs and the PDB. - "to do the right thing" in most cases without manual intervention. Please let me know if you see problems with the new way phenix.refine uses occupancies to determine special position constraints and exclusions for nonbonded symmetry interactions. Does anyone need atoms with partial occupancies but constrained to stay on a symmetry axis during refinement? Ralf

6 7

Re: [phenixbb] Phaser SAD question
by Kitty Pattridge 24 Dec '09

24 Dec '09

"Randy J. Read" <rjr27(a)cam.ac.uk> wrote: >Dear Peter, > >I'm interested in your observation that you get better signal from an >anomalous difference map computed with density-modified phases than you get >with an LLG map using the MR model. Presumably the density modification >just started from the model phases as well? We've never tried that, so >we'll have to look into that. Unfortunately, you can't at the moment start >the LLG completion with a "model" specified through structure factors, >though this is on our to-do list. > >With the MR+SAD option, because it is iterative, there are sites found in >subsequent cycles (where the information from some of the anomalous >scatterers has become available) and you end up with a set of reliable >sites that can be significantly larger than the set obtained from the first >LLG map. Does your anomalous difference Fourier with density modified >phases give you a larger apparently reliable substructure than the >iterative completion? You could use phenix.emma to compare the two >substructures and see how much they have in common. > >You can indeed give Phaser the set of sites that you got from the anomalous >difference map. One way to choose the cutoff would be to see how deep the >deepest hole is (i.e. how many standard deviations the map goes in the >negative direction), and choose a larger number in the positive direction. >Or you could just use the peaks above 6-7 sigma. > >I'd be interested in hearing how you get on with this approach, off-line if >you prefer. > >Good luck and Happy Christmas! > >Randy Read > >On Dec 24 2009, Peter Grey wrote: > >>Dear Phenix Experts, >> >>I use the SAD+MR option in Phenix. The ASU is huge and the number of >>scatterres is in the few hundreds range. I try to minimize the number of >>runs of Phaser SAD+MR since they take long time on my computer and would >>appreciate your insight regarding the following : >> >> - I noticed that I get much higher signal and much larger number of sites >> if I use density modified phases to calculate the anomalous difference >> map. Is it possible to input these phases instead of the model to Phaser >> ? >> >> - If not, do you recommend providing Phaser with the list of sites found >> in this difference anomalous map. Should such a list include only the >> sites you can be sure of (say above 7 sigma) or all possible sites (say >> above 5 sigma) ? >> >> - Should I use the (peak height)/(sigma of map) as the first estimate of >> occupancy ? if so do I have to normalize it (strongest site has occupancy >> of 1) or can I leave the height/sigma values as they are ? >> >> >>Many thanks, >> >>Peter. >> > >_______________________________________________ >phenixbb mailing list >phenixbb(a)phenix-online.org >http://phenix-online.org/mailman/listinfo/phenixbb > >

1 0