- phenixbb - phenix-online.org

(no subject)
by Rajendra Sharma 15 Nov '17

15 Nov '17

1 0

Max_reasonble_bond_distance error in phene.resl.space.refine
by Rune Thomas Kidmose 15 Nov '17

15 Nov '17

I have been playing around with phenix.real.space.refine as part of a small pipeline I am making for EM structures. Testing various deposited EM structures and their corresponding maps, within my pipeline, I sometimes run in to the below error: RuntimeError: Bond distance > max_reasonable_bond_distance: 79.8663 > 50: Followed by a list of distances between, what I am assuming is, atom numbers and the actual distance. I can get past it simply by changing the max_reasonable_bond_distance for the run to, in this case, 80. But I tend to see, that the final results gets significantly worse when I use this bypass. This is based on very similar models going through the pipeline, where the ones which yields the error, scores much worse when using the above mentioned bypass, compared to the models that showed no max_reasonable_distance_error. So I was wondering: A: What does "max reasonable bond distance" actual mean? B: Any ide what I can do, to avoid running into this bond distance problem? I know the latter is asking a lot, as I have not disclosed any of the details surrounding how the error arose etc. but if such info is needed, I would prefer to do so off list and then return with the solution on list, if possible. Maybe I was hoping for a setting that ignores the outliers, without changing the max value. But that all depends on, what the setting actual is, which I dont know Cheers Ruki

2 4

Twin refinement problem
by Leonid Sazanov 15 Nov '17

15 Nov '17

Dear all, there seems to be a problem with twin refinement of some of our datasets, and I wanted to see if that is something specific to our data or if anybody else has seen it. These crystals have been reproduced many times and are always twinned in P21 space group with twin operator: -h,-k,h+l and twin fraction in the range 0.45 - 0.49. Data always refines fine in Refmac and in many cases in Phenix as well, giving similar twin fraction in both programs. However, with some datasets (at 3.4-3.6 A resolution) in Phenix there is a problem: refinement of twin fraction is unstable and jumps each macro-cycle between about 0.40 (wrong?, giving high R-factors) and 0.47 or so (true?, giving much lower R-factors). This can be clearly seen from log file extract: cycle stage r-work r-free twin_fraction 0 : 0.3002 0.3084 1_adp: 0.2511 0.2701 0.49 2_adp: 0.2724 0.2837 0.42 3_adp: 0.2296 0.2465 0.47 4_adp: 0.2464 0.2587 0.40 5_adp: 0.2230 0.2420 0.47 6_adp: 0.2393 0.2554 0.40 7_adp: 0.2225 0.2423 0.47 8_adp: 0.2416 0.2590 0.40 9_adp: 0.2184 0.2376 0.47 10_adp: 0.2380 0.2556 0.40 end: 0.2393 0.2568 0.39 ------------------------------------------- The same dataset refines with stable twin fraction = 0.48 in refmac. When Phenix accepts higher twin_fraction with its lower R-factors, it does cycles 0 and 1 of twin fraction refinement once, while when it fails to update and falls back to low twin_fraction, it does this twice, as shown below. Updates OK: =========================== Bulk-solvent and scaling ========================== ============================= updating all scales ============================= start: r(all,work,free)=0.2383 0.2381 0.2565 n_refl.: 191608 twin_fraction=0.40 re-set all scales: r(all,work,free)=0.2926 0.2926 0.2950 n_refl.: 191608 twin_fraction=0.40 cycle 0: update twin fraction: r(all,work,free)=0.2910 0.2910 0.2933 n_refl.: 191608 twin_fraction=0.48 update scales twin start: r(all,work,free)=0.2910 0.2910 0.2933 n_refl.: 191608 twin_fraction=0.48 reset f_part, k_(total,mask): r(all,work,free)=0.2910 0.2910 0.2933 n_refl.: 191608 twin_fraction=0.48 bulk-solvent and scaling: r(all,work,free)=0.2249 0.2246 0.2473 n_refl.: 191608 twin_fraction=0.48 cycle 1: update twin fraction: r(all,work,free)=0.2249 0.2246 0.2472 n_refl.: 191608 twin_fraction=0.47 update scales twin start: r(all,work,free)=0.2249 0.2246 0.2472 n_refl.: 191608 twin_fraction=0.47 reset f_part, k_(total,mask): r(all,work,free)=0.2910 0.2911 0.2934 n_refl.: 191608 twin_fraction=0.47 bulk-solvent and scaling: r(all,work,free)=0.2249 0.2246 0.2472 n_refl.: 191608 twin_fraction=0.47 r_work=0.2246 r_free=0.2472 twin_fraction=0.47 twin_law=-H,-K,H+L ============================ Set refinement target ============================ Falls back: =========================== Bulk-solvent and scaling ========================== ============================= updating all scales ============================= start: r(all,work,free)=0.2216 0.2214 0.2408 n_refl.: 191608 twin_fraction=0.47 re-set all scales: r(all,work,free)=0.2855 0.2855 0.2848 n_refl.: 191608 twin_fraction=0.47 cycle 0: update twin fraction: r(all,work,free)=0.2854 0.2855 0.2848 n_refl.: 191608 twin_fraction=0.48 update scales twin start: r(all,work,free)=0.2854 0.2855 0.2848 n_refl.: 191608 twin_fraction=0.48 reset f_part, k_(total,mask): r(all,work,free)=0.2854 0.2855 0.2848 n_refl.: 191608 twin_fraction=0.48 bulk-solvent and scaling: r(all,work,free)=0.2186 0.2184 0.2377 n_refl.: 191608 twin_fraction=0.48 cycle 1: update twin fraction: r(all,work,free)=0.2186 0.2184 0.2377 n_refl.: 191608 twin_fraction=0.48 update scales twin start: r(all,work,free)=0.2186 0.2184 0.2377 n_refl.: 191608 twin_fraction=0.48 reset f_part, k_(total,mask): r(all,work,free)=0.2854 0.2855 0.2848 n_refl.: 191608 twin_fraction=0.48 bulk-solvent and scaling: r(all,work,free)=0.2186 0.2184 0.2377 n_refl.: 191608 twin_fraction=0.48 start: r(all,work,free)=0.2783 0.2782 0.2848 n_refl.: 191509 twin_fraction=0.48 re-set all scales: r(all,work,free)=0.2783 0.2782 0.2848 n_refl.: 191509 twin_fraction=0.48 cycle 0: update twin fraction: r(all,work,free)=0.2783 0.2782 0.2848 n_refl.: 191509 twin_fraction=0.48 update scales twin start: r(all,work,free)=0.2783 0.2782 0.2848 n_refl.: 191509 twin_fraction=0.48 reset f_part, k_(total,mask): r(all,work,free)=0.2783 0.2782 0.2848 n_refl.: 191509 twin_fraction=0.48 bulk-solvent and scaling: r(all,work,free)=0.2330 0.2329 0.2486 n_refl.: 191509 twin_fraction=0.48 cycle 1: update twin fraction: r(all,work,free)=0.2306 0.2305 0.2458 n_refl.: 191509 twin_fraction=0.40 update scales twin start: r(all,work,free)=0.2306 0.2305 0.2458 n_refl.: 191509 twin_fraction=0.40 reset f_part, k_(total,mask): r(all,work,free)=0.2781 0.2780 0.2846 n_refl.: 191509 twin_fraction=0.40 bulk-solvent and scaling: r(all,work,free)=0.2295 0.2293 0.2439 n_refl.: 191509 twin_fraction=0.40 r_work=0.2374 r_free=0.2522 twin_fraction=0.40 twin_law=-H,-K,H+L ============================ Set refinement target ============================ With some other datasets Phenix just gets stuck at apparently in-correctly low twin fraction - each cycle it tries to update to 0.47-0.48 (and R factors are better in this case) but then it always falls back to 0.36-0.40 again. So I wonder if anybody else has seen this with their data? Many thanks for any info! Leonid -- Prof. Leonid Sazanov IST Austria Am Campus 1 A-3400 Klosterneuburg Austria Phone: +43 2243 9000 3026 E-mail: sazanov(a)ist.ac.at

1 0

Postdoctoral research position
by Demetres D. Leonidas 15 Nov '17

15 Nov '17

Postdoctoral research position: molecular biologist/bioinformatician A 2-year Post-Doctoral position is available immediately in the Laboratory of Molecular Biology and Genomics in the Department of Biochemistry and Biotechnology of the University of Thessaly, located in the city of Larissa. The selected candidate will work in the framework of a national Research Infrastructure program on Synthetic Biology. Research in the Laboratory of Molecular Biology and Genomics focuses on various aspects of insect olfaction, reproduction and evolution, aiming at the development of new approaches for control of insects of agricultural or medical importance. We are seeking a highly motivated and creative scientist with a track record of successful research. A good bioinformatics background is strongly desired. The ideal candidate should be able to analyze raw sequencing data, perform genome assembly and annotation and have a good understanding of comparative and functional genomics. At the same time, he/she should be able to design and perform experiments in support of the bioinformatics analyses. Candidates should have a PhD and should have training in one or more of the following disciplines: bioinformatics, molecular biology, biochemistry, cell biology, and genome-wide technologies. Familiarity with insect rearing is desirable but not essential. Curriculum Vitae, a brief letter describing interests and prior experience and names of three references should be sent to: kmathiop(a)bio.uth.gr. Research position: MSc Bioinformatician A 2-year position for an MSc Bioinformatician is available immediately in the Laboratory of Molecular Biology and Genomics in the Department of Biochemistry and Biotechnology of the University of Thessaly, located in the city of Larissa. The selected candidate will work in the framework of a national Research Infrastructure program on Synthetic Biology. Research in the Laboratory of Molecular Biology and Genomics focuses on various aspects of insect olfaction, reproduction and evolution, aiming at the development of new approaches for control of insects of agricultural or medical importance. We are seeking a highly motivated individual with a BSc in Biology or related science and an MSc in Bioinformatics. The successful candidate will be involved in raw sequencing data analysis, genome assembly and annotation and comparative and functional genomics studies. His/her research could also lead to a PhD. Curriculum Vitae, a brief letter describing interests and prior experience and names of three references should be sent to: kmathiop(a)bio.uth.gr. -- ----------------------------------------------- Dr. Demetres D. Leonidas Professor of Biochemistry Department of Biochemistry & Biotechnology University of Thessaly Biopolis 41500 Larissa, Greece --------------------------------------------- Tel. +302410 565278 Tel. +302410 565297 (Lab) Fax. +302410 565290 E-mail: ddleonidas(a)bio.uth.gr http://www.bio.uth.gr ORCID ID orcid.org/0000-0002-3874-2523 ----------------------------------------------- --- Αυτό το e-mail ελέγχθηκε για ιούς από το πρόγραμμα Avast antivirus. https://www.avast.com/antivirus

1 0

bgroup refinement of ZN ions
by Ursula Schulze-Gahmen 13 Nov '17

13 Nov '17

I have 2 Zn ions which are defined as HETATM, in my structure. I am using tls and bgroup refinement for the whole structure. It seems like the Zn ions are treated as one group. They always have the same B-value. I would like to refine the B-value of the ZN ions separately, but refine the rest of the structure with bgroup refinement. Is that possible? How would I do this? Ursula -- Ursula Schulze-Gahmen, Ph.D. Project Scientist UC Berkeley, QB3 360 Stanley Hall #3220 Berkeley, CA 94720-3220 (510) 643 9491

1 0

mmtbx.prepare_pdb_deposition sequence file
by Ursula Schulze-Gahmen 09 Nov '17

09 Nov '17

I have problems running 'mmtbx.prepare_pdb_deposition' on the cif file created in phenix. I think the problem may be the sequence file. The coordinates are from a complex of 4 proteins plus an RNA. I created a text file with 5 chains, each chain sequence in one letter code. I am getting the following error code when I try to run the program: Traceback (most recent call last): File "/home/programs/x86_64-linux/phenix/1.12-2829/phenix-1.12-2829/build/../modules/cctbx_project/mmtbx/command_line/prepare_pdb_deposition.py", line 99, in <module> run(sys.argv[1:]) File "/home/programs/x86_64-linux/phenix/1.12-2829/phenix-1.12-2829/build/../modules/cctbx_project/mmtbx/command_line/prepare_pdb_deposition.py", line 69, in run alignment_params=params) File "/home/programs/x86_64-linux/phenix/1.12-2829/phenix-1.12-2829/modules/cctbx_project/iotbx/pdb/hierarchy.py", line 1052, in as_cif_block_with_sequence assert len(chain.residue_groups) + chain.n_missing_start + chain.n_missing_end == len(sequence) There seems to be something wrong with chain N, but I am not sure what. Is it ok to have a longer sequence in the sequence file than in the coordinates? Ursula -- Ursula Schulze-Gahmen, Ph.D. Project Scientist UC Berkeley, QB3 360 Stanley Hall #3220 Berkeley, CA 94720-3220 (510) 643 9491

1 0

chain and residue ID in phenix.erraser
by Ursula Schulze-Gahmen 08 Nov '17

08 Nov '17

I am trying to run phenix.erraser, but I am not able to define the single residue for rebuilding correctly. I am getting error messaheg: Please make sure that you include both the chain ID and residue number, for example 'A15'. I want to rebuild nucleotide G32 in chain N, and defined the residue as N32. That is not correct. What would be the correct syntax for this case? Ursula -- Ursula Schulze-Gahmen, Ph.D. Project Scientist UC Berkeley, QB3 360 Stanley Hall #3220 Berkeley, CA 94720-3220 (510) 643 9491

3 3

Question about writing cif files at the last round of refinement
by Boaz Shaanan 07 Nov '17

07 Nov '17

3 5

kicked map calculations
by Ursula Schulze-Gahmen 07 Nov '17

07 Nov '17

I am trying to calculate kicked maps in the phenix GUI, but I don't get the option, nit in the mtz files, and not in the CCP4 map option. Could that have something to do with this error message I am getting when starting Phenix: Gtk-Message: Failed to load module "canberra-gtk-module" Gtk-Message: Failed to load module "pk-gtk-module" Or are there any settings I need to change? Ursula -- Ursula Schulze-Gahmen, Ph.D. Project Scientist UC Berkeley, QB3 360 Stanley Hall #3220 Berkeley, CA 94720-3220 (510) 643 9491

2 1

data fill in in phenix refine?
by Ursula Schulze-Gahmen 26 Oct '17

26 Oct '17

I have data from a crystal with anisotropic diffraction. After scaling and merging using the webserver staraniso, I would like to refine using phenix refine. Could you tell me if phenixrefine will try to fill in the missing reflections that were removed in the staraniso procedure? If so, is there a a way to prevent this fill in? Following is the warning message for the Staraniso server: We are grateful to Andrew Leslie (MRC-LMB) for identifying a potentially undesirable 'feature' of CCP4/REFMAC (and possibly other refinement software), namely the 'fill-in' of unobserved data with *DF*calc which makes it effectively *incompatible* with the MTZ file output by STARANISO. Our advice is therefore to use *BUSTER <http://www.globalphasing.com/buster>* for refinement using the downloaded MTZ file, unless the user can find a way to disable that feature in the alternative software. In general, any treatment of STARANISO's output data that adds missing reflections out to the resolution limit may result in a very high proportion of unobserved data compared with data that have had a cut-off applied according to an isotropic criterion. Where missing reflections have been added in this way the user should take care to verify that any applications used to process the data subsequently can handle it correctly. Ursula -- Ursula Schulze-Gahmen, Ph.D. Project Scientist UC Berkeley, QB3 360 Stanley Hall #3220 Berkeley, CA 94720-3220 (510) 643 9491

3 3