Re: [ccp4bb] Apparent twinning in P 1 21 1 (refmac)
Yuri, note, the R-factors in case of twinning are not directly comparable. Pavel On Thu, Sep 29, 2011 at 7:25 PM, Yuri Pompeu wrote: > after 1 round of refmac rigid body and restrained refinement with twin law > (estimated alpha 0.47) > I am looking at 0.25 -0.29 Rwork Rfree and overall FOM 0.72. > I also defined NCS restraints... >
Re: [ccp4bb] Apparent twinning in P 1 21 1 (refmac)
after 1 round of refmac rigid body and restrained refinement with twin law (estimated alpha 0.47) I am looking at 0.25 -0.29 Rwork Rfree and overall FOM 0.72. I also defined NCS restraints...
Re: [ccp4bb] Apparent twinning in P 1 21 1
I am now realizing that! I have been refining with twin law enabled, its been somewhat succesful. I have been using phenix and now getting ready to try refmac... regards On Thu, 29 Sep 2011 15:23:13 +0100, Garib N Murshudov wrote: > Why do you detwin? It would not be normal procedure if you have twinning. Molecular replacement programs usually do not have much problem with twinned data and refinement programs can deal with them more or less properly. > when you detwin then errors are increased (As far as I remember proportional to 1/(1-2alpha), if alpha is 0.46 then errors will increase more than 10 times). Moreover it is very likely that twin and pseudo rotation are close to each other and estimated twin fractions may not be accurate. > regards > Garib > > On 29 Sep 2011, at 15:03, Yuri Pompeu wrote: > >> After I ran DETWIN with the estimated 0.46 alpha, my completeness for the detwinned data is now down to 54%!!! >> Is this normal behavior? (I am guessing yes since the lower symmetry untwinned dat is P1 21 1) > > Garib N Murshudov > Structural Studies Division > MRC Laboratory of Molecular Biology > Hills Road > Cambridge > CB2 0QH UK > Email: ga...@mrc-lmb.cam.ac.uk [1] > Web http://www.mrc-lmb.cam.ac.uk [2] -- Yuri Pompeu Links: -- [1] mailto:jen...@mrc-lmb.cam.ac.uk [2] http://www.mrc-lmb.cam.ac.uk/
Re: [ccp4bb] Apparent twinning in P 1 21 1
I am refining in phenix twin law h,-k,-h-l. Without twinning I was around 0.36 Rfree and 0.25 with twinning. I am noticing however that my Rgap keeps increasing slowly... (little concerned now its at 8% - 0.18-0.26 - to 2.4A) Maps look decent for 2.4A I have a lot of clashes however some are just bad waters though. phenix does not do real space when twinning is enabled. Any ideas here? thank you much On Thu, 29 Sep 2011 10:25:17 -0400, Phil Jeffrey wrote: Yuri, Detwinning relies on having both twin-related reflections present to calculate either/both of the the de-twinned data values. Therefore it magnifies incompleteness depending on where your missing data is with respect to the twin operator. I'd recommend against trying to do this with a twin fraction close to 0.5. From the DETWIN docs: Itrue(h1) = ((1-tf)*iTw(h1) -tf*iTw(h2)) / (1-2tf) i.e. tf = twin fraction, so 1/(1-2tf) becomes a large number and it's multiplying a weighted term of the form: (iTw(h1) - iTw(h2)) which becomes a very small number as the twin fraction approaches 0.5. The latter difference can easily be less than sigma(I), and so the signal/noise of your data plummets. Better to use REFMAC and phenix.refine's abilities to compensate for the twin fraction directly in refinement and leave your data as it is. Phil Jeffrey Princeton On 9/29/11 10:03 AM, Yuri Pompeu wrote: After I ran DETWIN with the estimated 0.46 alpha, my completeness for the detwinned data is now down to 54%!!! Is this normal behavior? (I am guessing yes since the lower symmetry untwinned dat is P1 21 1) -- Yuri Pompeu
Re: [ccp4bb] Apparent twinning in P 1 21 1
Yuri, Detwinning relies on having both twin-related reflections present to calculate either/both of the the de-twinned data values. Therefore it magnifies incompleteness depending on where your missing data is with respect to the twin operator. I'd recommend against trying to do this with a twin fraction close to 0.5. From the DETWIN docs: Itrue(h1) = ((1-tf)*iTw(h1) -tf*iTw(h2)) / (1-2tf) i.e. tf = twin fraction, so 1/(1-2tf) becomes a large number and it's multiplying a weighted term of the form: (iTw(h1) - iTw(h2)) which becomes a very small number as the twin fraction approaches 0.5. The latter difference can easily be less than sigma(I), and so the signal/noise of your data plummets. Better to use REFMAC and phenix.refine's abilities to compensate for the twin fraction directly in refinement and leave your data as it is. Phil Jeffrey Princeton On 9/29/11 10:03 AM, Yuri Pompeu wrote: After I ran DETWIN with the estimated 0.46 alpha, my completeness for the detwinned data is now down to 54%!!! Is this normal behavior? (I am guessing yes since the lower symmetry untwinned dat is P1 21 1)
Re: [ccp4bb] Apparent twinning in P 1 21 1
Why do you detwin? It would not be normal procedure if you have twinning. Molecular replacement programs usually do not have much problem with twinned data and refinement programs can deal with them more or less properly. when you detwin then errors are increased (As far as I remember proportional to 1/(1-2alpha), if alpha is 0.46 then errors will increase more than 10 times). Moreover it is very likely that twin and pseudo rotation are close to each other and estimated twin fractions may not be accurate. regards Garib On 29 Sep 2011, at 15:03, Yuri Pompeu wrote: > After I ran DETWIN with the estimated 0.46 alpha, my completeness for the > detwinned data is now down to 54%!!! > Is this normal behavior? (I am guessing yes since the lower symmetry > untwinned dat is P1 21 1) Garib N Murshudov Structural Studies Division MRC Laboratory of Molecular Biology Hills Road Cambridge CB2 0QH UK Email: ga...@mrc-lmb.cam.ac.uk Web http://www.mrc-lmb.cam.ac.uk
Re: [ccp4bb] Apparent twinning in P 1 21 1
After I ran DETWIN with the estimated 0.46 alpha, my completeness for the detwinned data is now down to 54%!!! Is this normal behavior? (I am guessing yes since the lower symmetry untwinned dat is P1 21 1)
Re: [ccp4bb] Apparent twinning in P 1 21 1
You might like to look at this.. It tries to explain likely twinning possibilities in P21. If you get C and P21, then probably a~=c - then Beta can have any value. C222 axes are then always possible with a* +c* , a*-c*, b* all having angles ~ 90 Without twinning you wont get 222 symmetry though. Pointless helps here. Eleanor http://www.ccp4.ac.uk/dist/html/twinning.html On 09/27/2011 08:23 PM, Linda Schuldt wrote: Dear Yuri, in a monoclinic space group an orthorhombic lattice metric can be simulated when one of the following conditions is fulfilled: i) a = c [e.g. in Wittmann& Rudolph (2007) Acta Cryst. D63, 744-749], ii) the beta angle is close to 90° [e.g. in Larsen et al. (2002) Acta Cryst. D58, 2055-2059 ] or iii) c cos beta is about -a/2 [e.g. in Declercq& Evrard, (2002) Acta Cryst. D57, 1829-1835]. The a and b axes of the orthorhombic cell are identical to the monoclinic a and c axes, respectively. The length of the orthorhombic b-axis can also be calculated by "c(monoclinic) cos(beta-90°) = 1/2b(orthorhomic)". I would assume that you have the case iii with a quite high twin fraction. If I recall correctly, Declercq and Evrard have a nice figure in their paper showing the geometric relationship. If not, let me know and I can sent you a figure. Good luck! Linda Yuri Pompeu schrieb: Hello everyone, I have a 2.3A data set that could be scaled in C 2 2 21 and P 1 21 1 Intensity statistics tests indicate twinning (pseudo-merohedral h,-k,-h-l in P 1 21 1) I find a good MR solution and when I try to refine it with the twin law I get fairly good maps and decent Rs 21-28%. I can see features tha were not in the search model Which leads me to think that this a valid solution. The one thing that bothers me however is the fact that my beta angle in P 1 21 1 is 104 (not close to 90) and that the geometry gets worse after refinement? Any suggestions? cheers *** Dr. Linda Schuldt Department of Molecular Biology University of Aarhus Science Park Gustav Wieds Vej 10c DK-8000 Århus C Denmark
Re: [ccp4bb] Apparent twinning in P 1 21 1
These papers describe something similar to what I see. Acta Cryst. (2001). D57, 1829-1835 Acta Cryst. (2009). D65, 388-392
Re: [ccp4bb] Apparent twinning in P 1 21 1
Dear Yuri, in a monoclinic space group an orthorhombic lattice metric can be simulated when one of the following conditions is fulfilled: i) a = c [e.g. in Wittmann & Rudolph (2007) Acta Cryst. D63, 744-749], ii) the beta angle is close to 90° [e.g. in Larsen et al. (2002) Acta Cryst. D58, 2055-2059 ] or iii) c cos beta is about -a/2 [e.g. in Declercq & Evrard, (2002) Acta Cryst. D57, 1829-1835]. The a and b axes of the orthorhombic cell are identical to the monoclinic a and c axes, respectively. The length of the orthorhombic b-axis can also be calculated by "c(monoclinic) cos(beta-90°) = 1/2b(orthorhomic)". I would assume that you have the case iii with a quite high twin fraction. If I recall correctly, Declercq and Evrard have a nice figure in their paper showing the geometric relationship. If not, let me know and I can sent you a figure. Good luck! Linda Yuri Pompeu schrieb: > Hello everyone, > I have a 2.3A data set that could be scaled in C 2 2 21 and P 1 21 1 > Intensity statistics tests indicate twinning (pseudo-merohedral h,-k,-h-l > in P 1 21 1) > I find a good MR solution and when I try to refine it with the twin law I > get fairly good maps and decent Rs 21-28%. I can see features tha were not > in the search model > Which leads me to think that this a valid solution. The one thing that > bothers me however is the fact that my beta angle in P 1 21 1 is 104 (not > close to 90) and that the geometry gets worse after refinement? > Any suggestions? > cheers > *** Dr. Linda Schuldt Department of Molecular Biology University of Aarhus Science Park Gustav Wieds Vej 10c DK-8000 Århus C Denmark
Re: [ccp4bb] Apparent twinning in P 1 21 1
On Tue, Sep 27, 2011 at 11:16 AM, Yuri Pompeu wrote: > I have a 2.3A data set that could be scaled in C 2 2 21 and P 1 21 1 > Intensity statistics tests indicate twinning (pseudo-merohedral h,-k,-h-l > in P 1 21 1) > I find a good MR solution and when I try to refine it with the twin law I > get fairly good maps and decent Rs 21-28%. I can see features tha were not > in the search model > Which leads me to think that this a valid solution. The one thing that > bothers me however is the fact that my beta angle in P 1 21 1 is 104 (not > close to 90) and that the geometry gets worse after refinement? > I've seen this before - the conventions for the C2221 and P21 unit cells are very different, so even if beta=104 in P21, the equivalent C2221 cell can still have all angles equal to 90. And you can definitely have pseudo-merohedral twinning in these circumstances (I did too - PDB ID 3ori). The problem with geometry is a separate issue - probably the automatic weighting not working properly, or an improper fixed weight. -Nat
[ccp4bb] Apparent twinning in P 1 21 1
Hello everyone, I have a 2.3A data set that could be scaled in C 2 2 21 and P 1 21 1 Intensity statistics tests indicate twinning (pseudo-merohedral h,-k,-h-l in P 1 21 1) I find a good MR solution and when I try to refine it with the twin law I get fairly good maps and decent Rs 21-28%. I can see features tha were not in the search model Which leads me to think that this a valid solution. The one thing that bothers me however is the fact that my beta angle in P 1 21 1 is 104 (not close to 90) and that the geometry gets worse after refinement? Any suggestions? cheers