- phenixbb - phenix-online.org

Computational Crystallography Newsletter (CCN) submissions
by Nigel Moriarty 21 Nov '17

21 Nov '17

Folks Calling for articles and short communications of interest to structural biologists. The initial deadline for submission for the January 2018 issue is 20 December 2017. The Computational Crystallography Newsletter (CCN) is an electronic newsletter for structural biologists, and is published online every 6 months. Feature articles, meeting announcements and reports can be submitted to me at any time for consideration. Submission of text by email or word-processing files using the CCN templates is requested. Past newsletters and a Microsoft Word template are available at www.phenix-online.org/newsletter. Cheers Nigel --- Nigel W. Moriarty Building 33R0349, Molecular Biophysics and Integrated Bioimaging Lawrence Berkeley National Laboratory Berkeley, CA 94720-8235 Phone : 510-486-5709 Email : NWMoriarty(a)LBL.gov Fax : 510-486-5909 Web : CCI.LBL.gov

1 0

Bug in 2.10.3 with scrolling
by Mark Saper 20 Nov '17

20 Nov '17

If one scrolls an alignment (not WRAP mode), the names of the sequences scroll but not the sequences with it. Using the scroll bar works correctly. Running 2.10.3 on MacOS 10.12.6 with Java version: 1.8.0_151. Mark _______________________________________________ Mark A. Saper, Ph.D. Associate Professor of Biological Chemistry, U-M Medical School Room 3220B, MSRB III saper(a)umich.edu +1 (734) 764-3353

1 0

Ligand library error in Phenix.Refine
by Gyanendra Kumar 16 Nov '17

16 Nov '17

Hi there, My Phenix.Refine run crashes because of non-availability of dictionary file for "ACP", a non-hydrolyzable ATP. Could this ligand be added to the next nightly build? Phenix.Refine suggested using "ReadySet", but that also failed. I sent the error reports to bugs(a)phenix-online.org, all emails bounced back. Refmac runs fine. Thanks! -Gyan -- Gyanendra Kumar, PhD Associate Scientist, St. Jude Children's Research Hospital, Department of Structural Biology, 262, Danny Thomas Place, MS-311 Memphis, TN 38105 Cell: 631-875-9189 https://www.linkedin.com/in/gyanendrakumar https://scholar.google.com/gyanendrakumar <https://scholar.google.com/citations?user=lnil0IQAAAAJ&hl=en> https://www.researchgate.net/profile/Gyanendra_Kumar3 http://stjude.academia.edu/GyanendraKumar -------------------------------------------------------

3 2

(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