On Friday 27 August 2010, Frank von Delft wrote:
> The requirement sounds extremely suspect: every atom in the structure
> contributes to every reflection, so refining "only some atoms" makes as
> little mathematical sense as refining against "only a subset of
> reflections".
I agree with you that the requirement sounds dubious.
But the specific argument you make is not quite right.
Two common counter-examples are real-space refinement and rigid-body
placement of a known fragment relative to an existing partial model.
> Especially not Bfactors, which are after all where errors get pushed
> into. So if you refine Bfactors for only some atoms, all you're doing
> is making the refinement program shove all the errors in the *whole*
> structure into the Bfactors of a few atoms. You'd have to explain
> extremely carefully why those atoms deserve that -- and even then I
> doubt anybody would believe you.
Indeed.
Ethan
>
> phx.
>
>
>
> On 26/08/2010 22:19, Hailiang Zhang wrote:
> > Thanks a lot Ethan, I will give it a try.
> >
> > Best Regards, Hailiang
> >
> >> On Thursday 26 August 2010 11:56:39 am Hailiang Zhang wrote:
> >>
> >>> Hi,
> >>>
> >>> I want to refine B factors for several residues only (all the other B
> >>> factors and all coordinates fixed, I know it sounds weird but there is a
> >>> reason to try that).
> >>>
> >> Maybe you could tell us what this reason is?
> >>
> >>
> >>> Is there anyway CCP4 can do this? Thanks for any suggestions!
> >>>
> >> Suggestion 1)
> >>
> >> Calculate structure factors for the entire rest of the model.
> >> Include these as F_partial contributing to the refinement of a model
> >> containing only your residues of interest. In this refinement,
> >> refine only the B terms.
> >>
> >> Caveat: I think you will encounter problems with how to handle
> >> the bulk solvent correction. Perhaps that must be included in
> >> F_partial also, and omitted from the subsequence mini-refinement.
> >>
> >> Suggestion 2)
> >>
> >> - Place your residues of interest in a single TLS group.
> >> - Do not assign any other atoms or residues to a TLS group.
> >> - Refine the entire model using refmac in TLS refinement mode.
> >> Choose 5 or 10 cycles of TLS refinement, but 0 cycles of
> >> coordinate/Biso refinement. Disregard all output other than the
> >> refined TLS description for the B factors in your residues of
> >> interest.
> >> - Use TLSANL to expand the TLS parameters back to individual
> >> Biso if you like.
> >>
> >>
> >>
> >>
> >>> Best Regards, Hailiang
> >>>
> >>>
> >> --
> >> Ethan A Merritt
> >> Biomolecular Structure Center, K-428 Health Sciences Bldg
> >> University of Washington, Seattle 98195-7742
> >>
> >>
> >>
>
|