Well - there s nothing there to indicate twinning, although as you say
it is hard to be sure with a pseudo translation...
The contrast doesnt look brilliant, but then it often doesnt..
What about phaser? If it gives a contrast between one spacegroup and the
others I always think that is a good sign..
It will take a while though..
Eleanor
Phil Evans program othercell suggests you could reindex to get two axes
equal.. And that is a prequisite for twinning .
C 1 2/m 1 192.3 192.3 117.2 90.0 90.0 92.0 1.98 [k-l,k+l,h]
Same cell
[-k+l,-k-l,h]
Eleanor
Kay Diederichs wrote:
> Eleanor Dodson schrieb:
>> You dont mention any twinning tests?
>
> sorry, I forgot to mention that the twinning tests do not show twinning.
>
> Rather, the actual curves in the "Cumulative distribution of H" lie on
> the "not-twinned-at-all" (i.e opposite) side of the alpha=0 curve (see
> plot below). But I'm pretty sure that this is due to the
> pseudo-translation (almost centering) which results in a high
> proportion of very weak reflections - contrary to what you get with
> twinning.
>
> That's what I get from sfcheck, when run in P212121:
> Perfect twinning test <I^2>/<I>^2 : 3.1699
> Partial Twinning test:
> -h,+l,+k
> Polar angles: 135.00 -89.99 179.99
> Alpha(twin fraction),Npair,Ior,Tol :-0.109 162588 2 0.030
>
> --- Partial Twinning Test : H = !I(h1)-I(h2)!/(I(h1)+I(h2)) ---
>
> Alpha(twinning fraction) = 1/2 - <H>
>
> Reflection related to hkl :
> -h,+l,+k
>
> Cumulative distribution of H
>
> 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
> *---+---+---+---+---+---+---+---+---+---+
> O+++.+ . . . . . . . . !
> O ++.+ . + . . . . . . . !
> !O +++ + . +. . . . . . !
> 0.10!-O-++-+--+---------+-------------------!
> ! O .++ + .+ . . .+ . . . !
> ! O. + + + . +. . . . + . . !
> ! O +.+ + . + . . . . +. !
> 0.20!---O---+-+--+------+-------------------+ alpha=0.4
> ! .O .+ +. +. . +. . . . +
> ! . O . + + .+ . .+ . . . +
> ! . O . +. + . + . . + . . +
> 0.30!------O----+--+----+---------+---------+
> ! . O .+ + . + . . .+ . +
> ! . O . + .+ . +. . . +. +
> ! . .O . +. +. .+ . . . + +
> 0.40!---------O-----+---+------+------------+ alpha=0.3
> ! . . O. .+ .+ . + . . +
> ! . . O. . + . +. . + . . +
> ! . . O . +. + . + . +
> 0.50!------------O------+----+-------+------+
> ! . . .O . .+ . + . . +. +
> ! . . . O . . + . + . .+ +
> ! . . . O. . +. .+ . . + +
> 0.60!---------------O-------+-----+---------+ alpha=0.2
> ! . . . .O . .+ . +. . +
> ! . . . . O . . + . .+ . +
> ! . . . . O. . +. . + . +
> 0.70!------------------O--------+------+----+
> ! . . . . O . .+ . + +
> ! . . . . .O . . + . .+ +
> ! . . . . . O. . +. . ++
> 0.80!-----------------------O-------+-------+ alpha=0.1
> ! . . . . . .O . .+ . +
> ! . . . . . . O . . + . +
> ! . . . . . . O. . +. +
> 0.90!----------------------------O------+---+
> ! . . . . . . . O. .+ +
> ! . . . . . . . .O . + +
> ! . . . . . . . . O. ++
> H *---+---+---+---+---+---+---+---+---+---* alpha=0
> 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
>
> Of course this plot looks _very_ different if I run it in P1 because
> then sfcheck uses -h,k,-l as twinning operator - it then looks like
> perfect twinning.
>
>
>> The L test, now part of the newest ctruncate is pretty good at
>> detecting twinning even with the NCS translation.
>
> this is the L test from ctruncate 1.0.0 : 30/10/08 run in P212121 :
>
> $TABLE: L test for twinning:
> $GRAPHS: cumulative distribution function for |L|:0|1x0|1:1,2,3,4:
> $$ |L| Observed Expected_untwinned Expected_twinned $$
> $$
> 0.000000 0.000000 0.000000 0.000000
> 0.050000 0.048142 0.050000 0.074938
> 0.100000 0.093721 0.100000 0.149500
> 0.150000 0.139205 0.150000 0.223312
> 0.200000 0.184799 0.200000 0.296000
> 0.250000 0.230867 0.250000 0.367188
> 0.300000 0.277173 0.300000 0.436500
> 0.350000 0.323650 0.350000 0.503563
> 0.400000 0.370747 0.400000 0.568000
> 0.450000 0.418233 0.450000 0.629437
> 0.500000 0.466569 0.500000 0.687500
> 0.550000 0.515615 0.550000 0.741812
> 0.600000 0.565706 0.600000 0.792000
> 0.650000 0.616860 0.650000 0.837688
> 0.700000 0.669232 0.700000 0.878500
> 0.750000 0.723200 0.750000 0.914062
> 0.800000 0.779144 0.800000 0.944000
> 0.850000 0.836994 0.850000 0.967938
> 0.900000 0.896933 0.900000 0.985500
> 0.950000 0.957486 0.950000 0.996313
> 1.000000 1.000000 1.000000 1.000000
> $$
>
>
>> And SFCHECK does a good job too.
>> If these are inconclusive I would not assume twinning.
>>
>> Usually you can get solutions for MR with twinned data, but I havent
>> much experience of the signal quality..
>> Can you solve it in P1 then sort out the spacegroup later?
>
> I tried; this is the full story:
> I ran molrep (version 10.2.12 from CCP4 6.1.0) in P1, with NMON=8. It
> uses the pseudo-translation vector and thus places 4 times 2
> molecules. In the "fast" mode (standard RF and TF without rigid body
> refinement) the "contrast" is 1.93/1.77/1.87/14.72 for the
> 1st/2nd/3rd/4th pair of molecules. However the result is different if
> I run it in "medium" (contrast=2.82/4.56/1.55/3.10) or "slow"
> (2.87/11.47/2.01/2.59) mode, and molrep only writes out 4 molecules
> instead of 8, in these two modes.
>
> I therefore suspected a bug and upgraded to molrep 10.2.27 from
> Alexei's webpage. But this version does not find the 8 molecules any
> more:
> Corr_for_fixed_model: 0.166
> WARNING: program can not improve current model
> result is "molrep.crd" with 4 monomers
> so this is quite confusing.
>
> I calculated structure factors from the P1 arrangement that I got in
> "fast" mode, from the 10.2.12 version. But it clearly is
> non-orthorhombic, and should have tetartohedral twinning to account
> for the pseudo-orthorhombic data.
>
> This is where I stopped and thought about asking on CCP4BB.
>
> best,
>
> Kay
>
>
>
>> Eleanor
>>
>>
>>
>> Kay Diederichs wrote:
>>> Dear all,
>>>
>>> we have crystals that nicely diffract to 1.7 A (sharp spots), with
>>> the following characteristics and findings:
>>> a) the data appear as P212121, with axes 117.2 133.6 138.3 (if
>>> reduced in P1, the largest deviation of any angle from 90° is 0.2°);
>>> the odd screw-axis reflections are mostly indistinguishable from
>>> noise; the data do not scale significantly better in P2/P21 (any
>>> setting) or P1.
>>> b) there is a good model available, with coords known from a complex
>>> of this protein with another one; two molecules of this model would
>>> give 64% solvent in P212121 which appears reasonable for a membrane
>>> protein
>>> c) the structure cannot be solved with this model in P212121, nor
>>> can it be solved in P222, P2122, P2212, P2221, P21212, P22121, P21221
>>> d) we conclude that the true space group must be P2 or P21 (with one
>>> of the three possible settings), with almost-perfect twinning. Or it
>>> is P1 with tetartohedral twinning. There are thus still six + one
>>> possibilities.
>>> e) MOLREP tells us
>>> --- Check Patterson for pseudo-translation ---
>>> PST_limit : 0.125 of origin peak
>>> INFO: pseudo-translation was detected.
>>> Origin Patterson peak: P,P/sig : 57.748 257.690
>>> 1 Patterson. peak : p,P/sig : 28.773 128.395
>>> 2 Patterson peak : P,P/sig : 16.551 73.856
>>> 3 Patterson peak : P,P/sig : 8.502 37.936
>>> Peak 1: trans.vector /ort/ : 0.011 55.688 69.399
>>> trans.vector /frac/: 0.000 0.416 0.500
>>> Peak 2: trans.vector /ort/ : 58.554 66.863 0.000
>>> trans.vector /frac/: 0.500 0.500 0.000
>>> Peak 3: trans.vector /ort/ : 58.565 11.385 69.399
>>> trans.vector /frac/: 0.500 0.085 0.500
>>> INFO: translation vector of peak 1 will be used.
>>>
>>> Two molecules (for the orthorhombic spacegroups) may produce only
>>> one pseudo-translation vector. As there is more than one strong
>>> pseudo-translation vector, I conclude that we have at least 3
>>> molecules in the ASU (consistent with monoclinic).
>>>
>>> f) we've calculated all seven molecular replacement searches of d)
>>> in MOLREP. The contrast is very high in all cases. However, Refmac
>>> rigid-body refinement of the solutions, with "Twin refinement"
>>> activated, gives about 51% R and the same for Rfree (give or take 1
>>> %), in all cases.
>>>
>>> I'm wondering: how reliable is a rotation search in the presence of
>>> perfect twinning? Is there any molecular replacement program that
>>> can take a given twinning operator into account in the rotation and
>>> translation search?
>>>
>>> Any other hints what to try?
>>>
>>> best,
>>>
>>> Kay
>>
>
>
|