anomalous difference map calculation
Hi, I am in the last refinement steps of a MR model and I want to calculate an anomalous difference map essentially to confirm the presence of a sulfite molecule and to locate vanadium (present in soaking solution). I read that it is necessary to have a mtz file with anomalous data (i.e. F+,F- or I+,I-). However, my data was collected at “normal” wavelenght (0.97) and it was processed with XDS considering Friedls law= true and my mtz file contain the following columns: H K L FP SIGFP. So, can I still calculate a anomalous difference map based on my data? Since I also have a Mo atom in the active site can I try to refine its occupancy by using the option “anomalous groups” in the refinement strategy? Thank you very much, Alexandra
Hi Alexandra, You will have to run the last steps in data processing over and obtain anomalous data in order to calculate that anomalous difference map. Once you have that data you can also refine the occupancy of your Mo atom, however at 0.97 A the anomalous scattering from Mo is not strong (f" is 1.2 electrons). All the best, Tom T ________________________________ From: [email protected] [[email protected]] on behalf of Alexandra Marques [[email protected]] Sent: Wednesday, June 17, 2015 7:04 AM To: [email protected] Subject: [phenixbb] anomalous difference map calculation Hi, I am in the last refinement steps of a MR model and I want to calculate an anomalous difference map essentially to confirm the presence of a sulfite molecule and to locate vanadium (present in soaking solution). I read that it is necessary to have a mtz file with anomalous data (i.e. F+,F- or I+,I-). However, my data was collected at “normal” wavelenght (0.97) and it was processed with XDS considering Friedls law= true and my mtz file contain the following columns: H K L FP SIGFP. So, can I still calculate a anomalous difference map based on my data? Since I also have a Mo atom in the active site can I try to refine its occupancy by using the option “anomalous groups” in the refinement strategy? Thank you very much, Alexandra
Hi Alexandra, I'm sorry, I misinterpreted your question. As Pavel and Tim point out there is no need to do anything special in your case to get the anomalous data. As Pavel mentioned, you can just take your current model and run phenix.refine with that model and your .HKL file and it will produce an anomalous difference map for you. If you have already refined you model with a different data file, then you should provide the free R flags (test set) from the refinement that you have already carried out. You can transfer these flags to your .HKL file with the reflection file utility in the Phenix GUI if you want. All the best, Tom T ________________________________ From: [email protected] [[email protected]] on behalf of Terwilliger, Thomas Charles [[email protected]] Sent: Wednesday, June 17, 2015 7:46 AM To: Alexandra Marques; [email protected] Subject: Re: [phenixbb] anomalous difference map calculation Hi Alexandra, You will have to run the last steps in data processing over and obtain anomalous data in order to calculate that anomalous difference map. Once you have that data you can also refine the occupancy of your Mo atom, however at 0.97 A the anomalous scattering from Mo is not strong (f" is 1.2 electrons). All the best, Tom T ________________________________ From: [email protected] [[email protected]] on behalf of Alexandra Marques [[email protected]] Sent: Wednesday, June 17, 2015 7:04 AM To: [email protected] Subject: [phenixbb] anomalous difference map calculation Hi, I am in the last refinement steps of a MR model and I want to calculate an anomalous difference map essentially to confirm the presence of a sulfite molecule and to locate vanadium (present in soaking solution). I read that it is necessary to have a mtz file with anomalous data (i.e. F+,F- or I+,I-). However, my data was collected at “normal” wavelenght (0.97) and it was processed with XDS considering Friedls law= true and my mtz file contain the following columns: H K L FP SIGFP. So, can I still calculate a anomalous difference map based on my data? Since I also have a Mo atom in the active site can I try to refine its occupancy by using the option “anomalous groups” in the refinement strategy? Thank you very much, Alexandra
Dear Alexandra, I would like to correct the two previous replies: there is no need to reprocess your data. XDS_ASCII.HKL will contain the anomalous information. The keyword 'FRIEDEL'S_LAW' affects the scaling, but not the output format. Actually some people recommend to always process with FRIEDEL'S_LAW=TRUE whether or not you want to extract the anomalous signal, because to get higher multiplicity for scaling (There are caveats if very strong Friedel/Bijvoet differences are present, but these are rare, I believe). You only need to rerun the conversion program and tell it to keep the Bijvoet pairs separate. The fastest way to create the anomalous difference map is to use shelxc for a SAD experiment. shelxc reads the XDS_ASCII.HKL directly and will also treat it correctly. This creates an _fa.hkl file. When you provide this and your pdb-file to the program anode (also available from the SHELX server), it create the anomalous map (ending .pha) that you can read into Coot. When you are familiar with shelx c/d/e and anode, it would take you maybe 1-2 min to create the map. If you are not familiar with the programs, maybe a few minutes more. Take a look at pages 25-27 in http://shelx.uni-ac.gwdg.de/~tg/teaching/anl-ccp4/pdfs/tgruene_chicago2014_s... for a first guide. Best, Tim On 06/17/2015 03:04 PM, Alexandra Marques wrote:
Hi,
I am in the last refinement steps of a MR model and I want to calculate an anomalous difference map essentially to confirm the presence of a sulfite molecule and to locate vanadium (present in soaking solution). I read that it is necessary to have a mtz file with anomalous data (i.e. F+,F- or I+,I-). However, my data was collected at “normal” wavelenght (0.97) and it was processed with XDS considering Friedls law= true and my mtz file contain the following columns: H K L FP SIGFP. So, can I still calculate a anomalous difference map based on my data?
Since I also have a Mo atom in the active site can I try to refine its occupancy by using the option “anomalous groups” in the refinement strategy?
Thank you very much, Alexandra
_______________________________________________ phenixbb mailing list [email protected] http://phenix-online.org/mailman/listinfo/phenixbb Unsubscribe: [email protected]
-- -- Dr Tim Gruene Institut fuer anorganische Chemie Tammannstr. 4 D-37077 Goettingen phone: +49 (0)551 39 22149 GPG Key ID = A46BEE1A
Hi Alexandra, to get anomalous difference map you do not need to do anything special: phenix.refine calculates this map by default as long as your input data are Fobs(+) and Fobs(-) or Iobs(+) and Iobs(-). Refining occupancy of Mo may not be a bad idea given it is likely to be partially occupied. Pavel On 6/17/15 06:04, Alexandra Marques wrote:
Hi,
I am in the last refinement steps of a MR model and I want to calculate an anomalous difference mapessentially to confirm the presence of a sulfite molecule and to locate vanadium (present in soaking solution). I read that it is necessary to have a mtz file with anomalous data (i.e. F+,F- or I+,I-). However, my data was collected at “normal” wavelenght (0.97) and it was processed with XDS considering Friedls law= true and my mtz file contain the following columns: H K L FP SIGFP. So, can I still calculate a anomalous difference map based on my data?
Since I also have a Mo atom in the active site can I try to refine its occupancy by using the option “anomalous groups” in the refinement strategy?
Thank you very much,
Alexandra
_______________________________________________ phenixbb mailing list [email protected] http://phenix-online.org/mailman/listinfo/phenixbb Unsubscribe: [email protected]
participants (4)
-
Alexandra Marques
-
Pavel Afonine
-
Terwilliger, Thomas Charles
-
Tim Gruene