- phenixbb - phenix-online.org

Real space refinement of complete virus structure
by #TOFAYEL AHMED# 13 Mar '18

13 Mar '18

Dear all, I am trying to refine a complete virus pdb structure (possessing icosahedral symmetry) in a cryo-EM map using phenix.real_space_refine. The run fails complaining "Multi-model PDB (with MODEL-ENDMDL) is not supported". So, I followed a previous post in phenixbb: http://www.phenix-online.org/pipermail/phenixbb/2011-May/017079.html and merged multiple models into one but faced following problems: 1. This merged file is unreadable in UCSF Chimera. 2. The secondary structure definition in the PDB file header is not recognizable anymore in Phenix. In my understanding, both the problems are arising due to PDB format's inability to handle more than 62 chains and 99999 atoms in a file. So, what will be the work-around to perform a real space refinement of the complete virus structure at once? Also, how to maintain secondary structure during such refinement? Any help is highly appreciated. Best regards, Tofayel ............................................ Tofayel Ahmed Graduate Student (PhD) | Lab-02n-01 School of Biological Sciences Nanyang Technological University 60 Nanyang Drive, Singapore 637551 +65 8650.0473 |tofayel001(a)e.ntu.edu.sg<mailto:[email protected]>

2 1

AutoSol Br SAD/MAD
by Ian Slaymaker 12 Mar '18

12 Mar '18

Hello Phenixbb, I've been out of the crystallography game for a few years and was hoping for some suggestions. I'm trying to phase a structure (~1200aa protein:DNA complex) with native data out to 2.0A and SAD/MAD data out to 3.0-2.6. I'm trying to use Br/I soaking while I wait for Sel-Met crystals to grow. I'm collecting 1800 frames at a synchrotron and scaling together multiple MAD wavelength sets (30 degree inverse wedges) collected from the same crystal with phenix.scale_and_merge to achieve >85% anomalous data (determined by xtriage). Anomalous signal of the merged data is around 10 at the 3.5 angstrom shell. I'm keeping an eye on the merging statistics, but I'm not sure what the critical stats to monitor are. AutoSol has given me a a few possible solutions, finding between 3 and 8 sites for MAD and SAD sets(I'm letting AutoSol decide how many to find) with Resolve Rfactors around 0.34, but low map skew (<1 and sometimes negative). I've tried using the initial sites to search for additional Br sites, but with no luck, and the density isn't good enough to build anything into. I'm doing hyss searches with various resolution cutoffs and continue to get the same/similar uninterpretable density with very few heavy atom sites found. My questions: What are the critical stats to look at when merging data for anomalous signal? Are there parameters to explore in AutoSol to identify more heavy atom sites? Should I accept that I need to wait for new crystals to grow or keep charging ahead with this data, encouraged by some promising numbers by AutoSol? Thank you in advance for the suggestions. -Ian

2 1

View Solvent Mask
by ben flath 09 Mar '18

09 Mar '18

Hi All I was watching the YouTube tutorial on Polder Maps and noticed that the bulk solvent mask was visualized in coot. How can I visualize my bulk solvent mask? Regards Ben Flath

4 5

Re: [phenixbb] View Solvent Mask
by Flath, Ben 09 Mar '18

09 Mar '18

It seems I am looking for the the compute_box option in polder maps. Ben Sent from my Android -------- Original message -------- From: "Flath, Ben" <bef206(a)mail.usask.ca> Date: 2018-03-09 11:15 AM (GMT-06:00) To: Pavel Afonine <pafonine(a)lbl.gov>, "Flath, Ben" <bef206(a)mail.usask.ca>, phenixbb(a)phenix-online.org Subject: Re: [phenixbb] View Solvent Mask Now the important question How can I modify the mask so it doesn't enter pockets within my model? Add HOH with occupancy set to zero and ignore_zero_occupancy_atoms=False? Or is running polder maps more effective? Thanks again Ben Sent from my Android -------- Original message -------- From: Pavel Afonine <pafonine(a)lbl.gov> Date: 2018-03-09 10:50 AM (GMT-06:00) To: "Flath, Ben" <bef206(a)mail.usask.ca>, phenixbb(a)phenix-online.org Subject: Re: [phenixbb] View Solvent Mask Hi Ben, I was watching the YouTube tutorial on Polder Maps and noticed that the bulk solvent mask was visualized in coot. How can I visualize my bulk solvent mask? phenix.mask model.pdb will calculate the mask. Pavel

1 0

Questions about the Autosol water placing option
by Kaibo Zhang 07 Mar '18

07 Mar '18

Helllo, I was using AutoSol GUI for MR-SAD of my data and trying to turn off "place waters in refinement". However, the output pdb files still got waters in them. I did check the run log file and confirmed both the "place_waters" and the "flood_with_waters" flags were set to false. But in the later session, it still recorded "adding waters in refinement". Is there any parameters I should adjust to avoid adding waters during the Autosol run. Or is it just turned on in the end in order to calculate R_work/free, and doesn't affect your phasing solution? Look forward to your responses. Thank you in advance! Regards Kaibo *Kaibo* PULSe, Purdue University Hockmeyer Hall, Room 321 240 S. Martin Jischke Drive, West Lafayette, IN 47907 Telephone: 765-337-1955 E-mail: zhang644(a)purdue.edu

2 1

Re: [phenixbb] Overlapping ligand-protein side chain density
by van den Bedem, Henry 02 Mar '18

02 Mar '18

If the phe is a gatekeeper, shouldn’t the ligand be refined at partial occupancy; i.e. occupancies not necessarily have to ‘add up’ as you suggest? Maybe this link is helpful too: www.biorxiv.org/content/early/2018/01/24/253419 H From: <phenixbb-bounces(a)phenix-online.org> on behalf of "Phan, Jason" <jason.phan(a)vanderbilt.edu> Date: Wednesday, February 28, 2018 at 1:30 PM To: "phenixbb(a)phenix-online.org" <phenixbb(a)phenix-online.org> Subject: [phenixbb] Overlapping ligand-protein side chain density All, A piece of a ligand and the side chain of a “gate-keeper” Phe occupy the same space with well-defined features observed for both (resolution is 1.7 A). It looks about 50:50. How do you refine both ligand and protein side chain in this case? A couple of phenixbb suggestions for dealing with ligand-ligand overlapping density have been considered. With the first suggestion, the two entities still clashed and moved apart off of their respective densities. The second suggestion is not applicable in this case since the other molecule is not a ligand but part of the protein but it was tested anyway. Although there was no bumping, the occupancies don’t add up, resulting in a big blob of negative density around the ligand piece. ------------------------------ Only two comments: 1. At that resolution, constrained group occupancy refinement should work reasonably well (provided you can model the 2 entities). Then you also do not have clashes between the molecules, because Occ(A)+Occ(B)=1, meaning when one (A) is there, the other one (B) is not. This works with refmac (external keyword file); if you need more sophisticated occupancy re/constraints SHELXL may offer more opportunities. 2. There is no necessity for the two NCS copies of the binding site to look exactly the same (non-equivalent). Maybe there is a good reason/story (accessibility, contacts etc) for one site to be occupied differently than the other one. Best, BR —————————————— Modeling two molecules that occupy overlapping binding sites in a structure simply involves designating them as alternate conformers, with the same chain and residue number, and an occupancy that sums to 1.0. For example, if you have an AMP and an ADP that occupy the same binding site, you would define them as AAMP B 501 BADP B 501 and initially set the occupancies for the atoms in each conformer to the ratio (50:50, 30:70, etc.) that you observe in the density. Refinement in this manner is straightforward in PHENIX. Diana ----------------------------

5 6

Overlapping ligand-protein side chain density
by Phan, Jason 28 Feb '18

28 Feb '18

All, A piece of a ligand and the side chain of a “gate-keeper” Phe occupy the same space with well-defined features observed for both (resolution is 1.7 A). It looks about 50:50. How do you refine both ligand and protein side chain in this case? A couple of phenixbb suggestions for dealing with ligand-ligand overlapping density have been considered. With the first suggestion, the two entities still clashed and moved apart off of their respective densities. The second suggestion is not applicable in this case since the other molecule is not a ligand but part of the protein but it was tested anyway. Although there was no bumping, the occupancies don’t add up, resulting in a big blob of negative density around the ligand piece. ------------------------------ Only two comments: 1. At that resolution, constrained group occupancy refinement should work reasonably well (provided you can model the 2 entities). Then you also do not have clashes between the molecules, because Occ(A)+Occ(B)=1, meaning when one (A) is there, the other one (B) is not. This works with refmac (external keyword file); if you need more sophisticated occupancy re/constraints SHELXL may offer more opportunities. 2. There is no necessity for the two NCS copies of the binding site to look exactly the same (non-equivalent). Maybe there is a good reason/story (accessibility, contacts etc) for one site to be occupied differently than the other one. Best, BR —————————————— Modeling two molecules that occupy overlapping binding sites in a structure simply involves designating them as alternate conformers, with the same chain and residue number, and an occupancy that sums to 1.0. For example, if you have an AMP and an ADP that occupy the same binding site, you would define them as AAMP B 501 BADP B 501 and initially set the occupancies for the atoms in each conformer to the ratio (50:50, 30:70, etc.) that you observe in the density. Refinement in this manner is straightforward in PHENIX. Diana ----------------------------

2 1

Plots of local real-space CC
by Murpholino Peligro 27 Feb '18

27 Feb '18

Dear all I want to get the values for the blue line (i.e. local real-space CC) in the Molprobity multicriterion plot. How can I do that? I want to plot those values with gnuplot. I do not want the whole graph (red and green lines). Thanks Murphy

3 5

hi please unsubscribe me
by Zhenwei Su 27 Feb '18

27 Feb '18

-- Zhenwei Su 4417 Spruce Street, APT 3A Philadelphia, PA 19104 email: zwsuustc(a)gmail.com cell phone: (608) 698-4156

1 0

hi please unsubscribe me
by Zhenwei Su 27 Feb '18

27 Feb '18

-- Zhenwei Su 4417 Spruce Street, APT 3A Philadelphia, PA 19104 email: zwsuustc(a)gmail.com cell phone: (608) 698-4156

1 0