Hi... I am trying to sharpen a map (x-ray) with the new phenix.auto_sharpen -> phenix.map_sharpening program. Got a script that basically runs phenix.xtriage to get the fraction of solvent content and some things and then I do ... ``` phenix.map_sharpening \ sharpening_method=b-factor \ mtz_in="$MTZ_IN" \ map_coeff_labels="2FOFCWT,PH2FOFCWT" \ map_model.model="$PDB_IN" \ crystal_info.resolution="$RESOL" \ b_factor_sharpening.b_factor_sharpening_crystal_info.solvent_content="$SOLVENT" \ b_factor_sharpening.b_factor_sharpening_crystal_info.is_crystal=True \ > sharpen.log 2>&1 ``` The problem is that I cannot open the mtz file in coot (mtz dump is attached) ... BUT ... I can open the sharpened map (ccp4 format) ... BUT ... the map is not covering the protein completely.. How do I do that? (1.png 2.png just showing the limits of the map) In gui mode I do get the map over the protein...BUT I cannot select the "FOFCWT,PHFOFCWT" labels (3.png ... That's the second part of the plan...) I am using phenix-dev 5989 Thanks a lot for your help
Dear PhenixBB community, I am currently working on solving the structure of protein X (110 amino acids). The protein forms a dimer in solution, and dimerization is functionally relevant. Based on initial analysis, the asymmetric unit (ASU) was estimated to contain one copy [unit cell: 82.699, 82.699, 36.435, 90 , 90, 120]. The data extend to 1.2 Å resolution, but the final R factors are around 0.23/0.25, which seem higher than expected at this resolution, and we have not been able to improve them further. When I examined the structure using symmetry mates in PyMOL, I observed a clear dimer that appears biologically meaningful. This makes me suspect that the ASU estimation might be incorrect—especially for a small protein like this, the ASU could potentially contain a dimer rather than a monomer. However, every time I use a dimer model as the search model for molecular replacement, Phaser consistently returns a monomer solution. I would greatly appreciate any suggestions on how to successfully obtain the dimeric model that matches the crystallographic symmetry. Thank you in advance for your help! Best regards, Liu
Hi,
When a homo-multimeric molecule has point group symmetry, it is very common for part of that symmetry to coincide with the space group symmetry. For instance, if a 222 (D2) tetramer crystallised in space group P222, a tetramer could sit at the intersection of 2 2-folds (1 monomer in the asymmetric unit), on a 2-fold (1 dimer or 2 monomers in the a.u.) or in a general position (4 monomers or one tetramer in the a.u.). In principle you could even have one of each, giving you a total of 7 monomers in the a.u., but all would be forming proper 222 tetramers in the crystal. So in your case the monomer model does match the crystallographic symmetry — you just have to make a dimer for visualisation by adding an appropriate symmetry copy.
Some time ago, Airlie McCoy added a feature to Phaser where, if the search model has point-group symmetry and it lands on a symmetry operator that agrees with that symmetry, it automatically removes the redundant copy that is generated by the space-group symmetry. That’s what has happened in your job, and if you look at the part of the log file dealing with the packing test, you should see some notification that this has happened.
Best wishes,
Randy Read
From: [email protected]
participants (3)
-
Francisco -
princeadam@163.com -
Randy John Read