I am surprised that this is such a mess. Small molecule crystallographers frequently encounter special positions and have been able to handle them correctly for many years. A well known small molecule refinement program that I don't need to advertise here has been able to generate the special position constraints for the atom coordinates and anisotropic displacement parameters fully automatically for all special positions in all space groups for the last 20 years, no special user action is required.
George On 01/02/2012 12:40 PM, herman.schreu...@sanofi.com wrote: > Dear Ed, > > I thought the clumsy way of handling special positions you describe was > abandoned after the Xplor days, that even rounding errors could cause > special atoms to fly apart, but you may be right. The Mn ion might not > really be special, but might have two minima, one closer to one protein > molecule and one closer to the symmetry related protein molecule. > Supposedly, refmac switches off vd Waals repulsions if the combined > occupancies of the the atoms involved does not exceed 1.0. This should > take care of the problem but apparently does not. Forcing the ion to > have to same position as its symmetry mate may not always be the best > representation of the reality in a crystal. > > The main reason I switched to Buster a few years ago is that in refmac > residual repulsions remain, even if the summed occupancies of atoms in > alternate positions do not exceed 1.0. > > Best, > Herman > > -----Original Message----- > From: CCP4 bulletin board [mailto:CCP4BB@JISCMAIL.AC.UK] On Behalf Of Ed > Pozharski > Sent: Monday, January 02, 2012 7:45 AM > To: CCP4BB@JISCMAIL.AC.UK > Subject: Re: [ccp4bb] Refmac and metal on a two-fold? > > On Sun, 2012-01-01 at 12:14 -0600, Dima Klenchin wrote: >> With Garib's help, I have forced the atom into its position by using >> external distance restrains of zero length against the same >> symmetry-related atom. The cause is unclear because the same program >> handles special positions in another structure just fine. > > Thanks - good to know. Perhaps a good practice is to take care of > special positions manually every time just in case. > >>From what I understand, the way refinement software generally treats > special positions (other than automatically resetting occupancy to 1/z) > is to exclude vdw repulsion to symmetry-related atom which in this case > would always push the atom out of position. Exclusion is usually > applied to the atom when it is within the cutoff distance (0.1A?) from > the corresponding symmetry element. If for whatever reason the atom > shifts, the vdw will push it out. Why would such atoms shift varies - I > suspect in some cases it is due to hydrogen bonding restraints, in > others it could be due to noisy density. > > > > -- > Oh, suddenly throwing a giraffe into a volcano to make water is crazy? > Julian, King of Lemurs > -- Prof. George M. Sheldrick FRS Dept. Structural Chemistry, University of Goettingen, Tammannstr. 4, D37077 Goettingen, Germany Tel. +49-551-39-3021 or -3068 Fax. +49-551-39-22582
