Hi all, I attempted to improve the density around a helix and short linker region using an average kick map but it didn't make things look much better. I then tried creating a composite iterative-build omit map from phenix.autobuild, but in this case the electron density around the helix/linker in question largely disappeared! Is there another approach I can use that might improve my maps, or should I start questioning my data? The dataset has >95% completeness to 2.7 A so I am very curious as to why I am missing the backbone now, where before I was only missing the side chain densities. One thing I am wondering is if the problem could have something to do with symmetry; the part of the map I am looking at may be close to a two-fold crystallographic symmetry axis and I have been told that these regions tend to be very noisy and hard to model--comments? Kyle Kyle T. Dolan Department of Biochemistry and Molecular Biology The University of Chicago [email protected]
Hi Kyle, It's always worth questioning your data and your model, so yes I'd examine both carefully. I am guessing that this is molecular replacement solution. (It isn't as common to have a loop that disappears during model-building with experimental phasing because normally you would not build it in the first place unless the density was there. My suggestions would be similar in either case though.) So (still guessing) your starting MR model had some parts that are correct, some that are not, and some that are close. You carry our MR and place your model. You calculate a sigmaA-weighted map, probably after refinement. This map shows great density in many places, medium density in others, little density in some. You rebuild the parts that are clear to you, but some are not so clear, including the helix and loop in question. You refine some more, and now you want to improve the density in this unclear region. Then as you say below, you calculate an average kick map, and it doesn't look so good in this unclear region. You calculate an iterative build omit map and the density in this region is nearly gone. A comment or two first: 1. Yes, it is possible that the presence of a symmetry element very near the problematic area can make the density noisier. Normally this happens only right at the symmetry element however, not even a few A away. 2. The averaged kick maps and iterative build omit maps are ways to remove model bias. In general (with some exceptions) they are not great ways to "improve" maps. Generally your best maps for building are going to be your current 2Fo-Fc map (e.g., from phenix.refine) or a density modified map using model information (e.g., overall_best_denmod_map_coeffs.mtz from autobuild). The omit and kick maps are to help you determine whether you have overinterpreted your density. So if this is all about right, here is my suggestion: try to see if there is any evidence whatsoever for density in this region. A. Go back to the MR solution (very beginning of this structure solution, before any refinement except rigid-body). As this model has not been refined, you can delete any portion, refine the rest, and there will be no model bias in the map. B. So: delete the region in question completely from this unrefined starting model. C. Then refine the model, calculate a map, and look at the map in the region of your helix and linker. If there is density (even not great density) that clearly matches the structure of the helix and the linker, then you know for sure that this helix and linker belong there. If not...perhaps they are disordered or somewhere else. D. If no density yet, you can , iterative build the model, do whatever you want with this model to make it better except do not put the helix and linker back in...and then look at the map in the region of the helix and linker. Still you have not biased the map, so if you can find them they are really there. If you cannot see them...now it is time to call them "disordered" until you get further evidence. I hope that helps! All the best, Tom T On Aug 9, 2010, at 3:08 PM, Kyle Dolan wrote:
Hi all,
I attempted to improve the density around a helix and short linker region using an average kick map but it didn't make things look much better. I then tried creating a composite iterative-build omit map from phenix.autobuild, but in this case the electron density around the helix/linker in question largely disappeared! Is there another approach I can use that might improve my maps, or should I start questioning my data? The dataset has >95% completeness to 2.7 A so I am very curious as to why I am missing the backbone now, where before I was only missing the side chain densities. One thing I am wondering is if the problem could have something to do with symmetry; the part of the map I am looking at may be close to a two- fold crystallographic symmetry axis and I have been told that these regions tend to be very noisy and hard to model--comments?
Kyle
Kyle T. Dolan Department of Biochemistry and Molecular Biology The University of Chicago [email protected] _______________________________________________ phenixbb mailing list [email protected] http://phenix-online.org/mailman/listinfo/phenixbb
Thomas C. Terwilliger Mail Stop M888 Los Alamos National Laboratory Los Alamos, NM 87545 Tel: 505-667-0072 email: [email protected] Fax: 505-665-3024 SOLVE web site: http://solve.lanl.gov PHENIX web site: http:www.phenix-online.org ISFI Integrated Center for Structure and Function Innovation web site: http://techcenter.mbi.ucla.edu TB Structural Genomics Consortium web site: http://www.doe-mbi.ucla.edu/TB CBSS Center for Bio-Security Science web site: http://www.lanl.gov/cbss
Hi Tom, On 8/10/10 7:36 AM, Tom Terwilliger wrote:
1. Yes, it is possible that the presence of a symmetry element very near the problematic area can make the density noisier. Normally this happens only right at the symmetry element however, not even a few A away.
could you please explain why this happens? Thanks! Pavel.
Hi Pavel, I haven't thought about it a lot, but here is my understanding of the effect. A symmetry axis (e.g., 2-fold) basically is the same as having 2 identical cells on top of each other, rotated 180 degrees. Imagine that cell 1 has random density. Then cell 2 has the same random density, rotated 180 degrees. Now add them together. At most points, the density from cell 1 and cell 2 are different, and the RMS of their sum is about sqrt(2) times the RMS of either one. On the other hand, at the symmetry axis, the density from cell 1 and cell 2 are identical, so that they always add constructively and the RMS is 2 times the RMS of either one. As you move away from the axis, the effect is smaller, and it decays away over the resolution of the map. Note that if you have a 3-fold, 4-fold or 6-fold, then the effect is even bigger. All the best, Tom T On Aug 10, 2010, at 9:07 AM, Pavel Afonine wrote:
Hi Tom,
On 8/10/10 7:36 AM, Tom Terwilliger wrote:
1. Yes, it is possible that the presence of a symmetry element very near the problematic area can make the density noisier. Normally this happens only right at the symmetry element however, not even a few A away.
could you please explain why this happens?
Thanks! Pavel.
_______________________________________________ phenixbb mailing list [email protected] http://phenix-online.org/mailman/listinfo/phenixbb
Thomas C. Terwilliger Mail Stop M888 Los Alamos National Laboratory Los Alamos, NM 87545 Tel: 505-667-0072 email: [email protected] Fax: 505-665-3024 SOLVE web site: http://solve.lanl.gov PHENIX web site: http:www.phenix-online.org ISFI Integrated Center for Structure and Function Innovation web site: http://techcenter.mbi.ucla.edu TB Structural Genomics Consortium web site: http://www.doe-mbi.ucla.edu/TB CBSS Center for Bio-Security Science web site: http://www.lanl.gov/cbss
Hi Tom, thanks! Interesting... I wonder if we need to correct Fo-Fc maps for this effect since a few such peaks (especially large ones in case of high order symmetry) may make a weak ligand density even less visible... Pavel. On 8/10/10 8:32 AM, Tom Terwilliger wrote:
Hi Pavel,
I haven't thought about it a lot, but here is my understanding of the effect. A symmetry axis (e.g., 2-fold) basically is the same as having 2 identical cells on top of each other, rotated 180 degrees. Imagine that cell 1 has random density. Then cell 2 has the same random density, rotated 180 degrees. Now add them together. At most points, the density from cell 1 and cell 2 are different, and the RMS of their sum is about sqrt(2) times the RMS of either one. On the other hand, at the symmetry axis, the density from cell 1 and cell 2 are identical, so that they always add constructively and the RMS is 2 times the RMS of either one. As you move away from the axis, the effect is smaller, and it decays away over the resolution of the map. Note that if you have a 3-fold, 4-fold or 6-fold, then the effect is even bigger.
All the best, Tom T
On Aug 10, 2010, at 9:07 AM, Pavel Afonine wrote:
Hi Tom,
On 8/10/10 7:36 AM, Tom Terwilliger wrote:
1. Yes, it is possible that the presence of a symmetry element very near the problematic area can make the density noisier. Normally this happens only right at the symmetry element however, not even a few A away.
could you please explain why this happens?
Thanks! Pavel.
Hi, the refinement of my structure brings me to a new point where I would appreciate any input. - I solved the structure by molecular replacement. - The refinement looked good, I was almost ready to deposit but there always was a fairly big blob of density that i couldn't explain, loosely coordinated to the protein and far from the active site. To try to identify what it could be, i calculated an anomalous difference map (data collected at 1A) and I see a peak right in the middle of this density. I have BaCl2 in my crystallization solution so it is the likely guilty ion! (and with further look, it looks like a loose metal site) So now, I want to refine with Barium in the density and I was wondering what was the best way to do it, and since there is anomalous signal, why not use it? * I updated to the new version of phenix and it reads I+ and I- (while the previous version didn't) but under the refinement settings, the target function ML-SAD is still for developers only. Is it a good idea to use the anomalous Is for refinement, or: * should I re-run Phaser-EP to have phases to use in phenix.refine, or: * should I not bother, the signal is weak anyway (1 barium/ ~400 residues) and continue with my non-anomalous scaled dataset and use the anomalous signal only for a figure to explain why I put the metal here. (resolution of 3.4A, quite anisotropic) I can also try everything but I'd like to understand what would be best. Thanks very much for your help vincent -- Vincent Chaptal Dept. of Physiology at UCLA http://www.physiology.ucla.edu/Labs/Abramson/index.html http://www.physiology.ucla.edu/Labs/Abramson/index.html/ Phone: 1-310-206-1399 IMPORTANT WARNING: This email (and any attachments) is only intended for the use of the person or entity to which it is addressed, and may contain information that is privileged and confidential. You, the recipient, are obligated to maintain it in a safe, secure and confidential manner. Unauthorized redisclosure or failure to maintain confidentiality may subject you to federal and state penalties. If you are not the intended recipient, please immediately notify us by return email, and delete this message from your computer.
Hi Vincent, - I would add that ion and refine its occupancy and anisotropic ADPs... Oh, I see, the resolution is 3.4A, well, then I would try first isotropic ones. - if the data set you are using contains Fobs(+) and Fobs(-), (and not merged Fobs_merged = 0.5*(Fobs(+) + Fobs(-))), I would refine f' and f'' of the anomalous scatterer. After having done this, check if the residual density is clear, and keep in mind that unbalanced syntheses may have Fourier truncation ripples around heavy atom, see relevant discussion here: http://www.mail-archive.com/[email protected]/msg14835.html and some pictures here: http://www.phenix-online.org/presentations/latest/pavel_maps.pdf In general, getting rid of an unmodeled density is good in a sense that it has a potential to improve the density everywhere else and you may better model other some other parts of your structure, so you get a better model overall. Let me know if you need any help with trying the above suggestions. Good luck! Pavel. On 8/11/10 5:56 PM, vincent Chaptal wrote:
Hi,
the refinement of my structure brings me to a new point where I would appreciate any input. - I solved the structure by molecular replacement. - The refinement looked good, I was almost ready to deposit but there always was a fairly big blob of density that i couldn't explain, loosely coordinated to the protein and far from the active site.
To try to identify what it could be, i calculated an anomalous difference map (data collected at 1A) and I see a peak right in the middle of this density. I have BaCl2 in my crystallization solution so it is the likely guilty ion! (and with further look, it looks like a loose metal site)
So now, I want to refine with Barium in the density and I was wondering what was the best way to do it, and since there is anomalous signal, why not use it? * I updated to the new version of phenix and it reads I+ and I- (while the previous version didn't) but under the refinement settings, the target function ML-SAD is still for developers only. Is it a good idea to use the anomalous Is for refinement, or: * should I re-run Phaser-EP to have phases to use in phenix.refine, or: * should I not bother, the signal is weak anyway (1 barium/ ~400 residues) and continue with my non-anomalous scaled dataset and use the anomalous signal only for a figure to explain why I put the metal here.
(resolution of 3.4A, quite anisotropic)
I can also try everything but I'd like to understand what would be best. Thanks very much for your help vincent
participants (4)
-
Kyle Dolan
-
Pavel Afonine
-
Tom Terwilliger
-
vincent Chaptal