Dear Andreas and all,
refining multiple models together would sound like the most
straightforward solution in a case like that. But (please correct me
if I am wrong), generally, such trick can only work in P1:
obviously, any symmetry operation (but lattice translation) applied to
the shifted model will generate a wrong structure...
Furthermore, is there a way to do an "ensemble refinement " in phenix
(or refmac)? In order to do that correctly one must be able to switch
off all interactions between individual models, while keeping
interactions within each model. Apparently, it can be done in CNS,
although I could not find an example of an input file - can anyone
share an example ?
and, by the way, what stops refmac or phenix developers to account for
lattice translation defects during refinement?
Peter
On 7 Jan 2009, at 11:35, Andreas Förster wrote:
> Hey Stephen,
>
> how about simply putting three models, separated by 16A, into your
> original unit cell and refining them together. You'd have to guess
> their fractional occupancies from the heights of the Patterson peaks
> (and make them add up to one, obviously).
>
> This approach (in contrast to correcting the intensities) was the
> consensus when I talked to Peter Zwart and Tom Terwilliger about a
> similar problem. As my problem still hasn't been solved (probably
> because of flawed data processing due to overlapping spots), I'd be
> very keen on hearing about your progress.
>
> http://scripts.iucr.org/cgi-bin/paper?S0907444908016648
> describes both approaches to lattice-translocated crystals
> (correcting intensities and refining several models concurrently).
>
> http://scripts.iucr.org/cgi-bin/paper?ba5111
> is the best review on the subject.
>
>
> All the best.
>
>
> Andreas
>
>
> Stephen Hare wrote:
>> Dear All,
>> We are currently working on a structure of apparent P21 symmetry
>> which has been solved by molecular replacement. The data are to
>> 2.7Å but the Rfree will not drop below 30%. The density is clear
>> for the model we have, however there is extra density that suggests
>> a shift of the structure by 16Å in either direction - resulting in
>> three possible overlapping positions for the structure. We assume
>> this is the result of twinning.
>> The unit cell dimensions are 102.7Å, 83.0Å, 115.3Å, 90°, 101.8°,
>> 90°. Examining the data with phenix.xtriage also suggests pseudo
>> translational symmetry with a separation of 16Å. A Patterson peak
>> at 0.097, 0.000, -0.096 is approximately 30% of the origin peak,
>> while a second peak of double the translation at 0.192, 0.000,
>> 0.195 is 7% of the origin peak. The structure contains a dimer of
>> dimers with an NCS 2 fold axis almost perpendicular to the
>> crystallographic 2 fold. This NCS axis almost coincides with the
>> diagonal between the A and C axes. A twin axis along the A C
>> diagonal (l,-k,h) could explain the observed extra density, however
>> this is not possible because A and C are different lengths.
>> As a result of the NCS axis running almost perpendicular to the
>> observed P21 axis, it is possible to merge the reflections in a
>> larger orthorhombic unit cell - dimensions 137.1Å, 83.3Å, 169.8Å
>> although here the Rmerge is higher and it is not possible to get a
>> molecular replacement solution.
>> Is it possible to define the (l,-k,h) twinning operator in our
>> original unit cell? or have we missed the actual unit cell?
>> Or....something else?
>> Steve
>> Stephen Hare PhD
>> Post doctoral research associate
>> Jefferiss Research Laboratories
>> Wright-Fleming Institute
>> Division of Medicine
>> Imperial College London
>> Norfolk Place
>> London W2 1PG
>> UK
>> Phone: +44 (0) 20 7594 3908 Fax: +44 (0) 20 7594 3906
>
> --
> Andreas Förster, Research Associate
> Paul Freemont & Xiaodong Zhang Labs
> Department of Biochemistry, Imperial College London
|